怎样用matlab 解决一个二阶微分方程 边界条件不是初值 如何用matlab求解一个二阶常系数微分方程组
matlab \u4e8c\u9636\u5fae\u5206\u65b9\u7a0b\u6570\u503c\u89e3 \u591a\u4e2a\u8fb9\u754c\u6761\u4ef6\u4e00\u822c\u6765\u8bf4\uff0c\u521d\u503c\u6761\u4ef6\u548c\u65b9\u7a0b\u786e\u5b9a\u4e86\u7684\u8bdd\uff0c\u7ec8\u503c\u6761\u4ef6\u4e5f\u662f\u786e\u5b9a\u4e86\u7684
clc
f=@(z,r)([r(1);-(1/r*r(1)^2*(1+r(1)^2)+1/0.073*(1+r(1)^2)^(3/2))]);
ode45(f,[0 0.003],[0.0017 1.192])
\u7b2c\u4e00\u79cd\u65b9\u6cd5\uff1a\u5229\u7528dsolve\u51fd\u6570\u6c42\u5fae\u5206\u65b9\u7a0b\u7684\u7b26\u53f7\u89e3\uff08\u901a\u89e3\uff09\uff1a\u5bf9\u4e8e\u4e00\u4e9b\u4e0d\u662f\u5f88\u96be\uff0c\u8981\u6c42\u51fa\u901a\u89e3\u7684\u5fae\u5206\u65b9\u7a0b\uff0c\u7528dsolve\u51fd\u6570\u6c42\u89e3\u3002
1\u3001 \u6253\u5f00Matlab\u8f6f\u4ef6-->\u70b9\u51fb\u65b0\u5efa\u811a\u672c\u83dc\u5355\uff0c\u65b0\u5efa\u4e00\u4e2a\u811a\u672c\u6587\u4ef6\u7528\u4e8e\u7f16\u5199\u5fae\u5206\u65b9\u7a0b\u6c42\u89e3\u7a0b\u5e8f\u3002
2\u3001 \u8f93\u5165\u5fae\u5206\u65b9\u7a0b\u6c42\u89e3\u7a0b\u5e8f-->\u70b9\u51fb\u4fdd\u5b58-->\u70b9\u51fb\u8fd0\u884c\u3002
3\u3001\u5728matlab\u7684\u547d\u4ee4\u7a97\u53e3\u5373\u53ef\u770b\u5230\u6c42\u89e3\u7ed3\u679c\uff0c\u662f\u4e00\u4e2a\u5173\u4e8e\u53c2\u6570a\uff0cb\u7684\u8868\u8fbe\u5f0f
\u7b2c\u4e8c\u79cd\u65b9\u6cd5\uff1a\u5229\u7528Matlab\u4e2d\u7684solver\u51fd\u6570(\u5305\u62ecode45\u3001ode23\u3001ode15s\u7b49)\u6765\u6c42\u89e3\u5fae\u5206\u65b9\u7a0b\u7684\u6570\u503c\u89e3\uff0c\u8fd9\u79cd\u65b9\u6cd5\u662f\u6700\u5e38\u7528\u7684\u65b9\u6cd5\uff0c\u5bf9\u4e8edsolve\u51fd\u6570\u96be\u4ee5\u6c42\u89e3\u7684\u65b9\u7a0b\u5c31\u53ef\u4ee5\u5229\u7528\u8fd9\u79cd\u65b9\u6cd5\u6c42\u89e3\u65b9\u7a0b\u7684\u6570\u503c\u89e3\u3002
1\u3001 \u6253\u5f00Matlab-->\u65b0\u5efa\u4e00\u4e2a\u811a\u672c\u6587\u4ef6\u7528\u4e8e\u7f16\u5199\u6c42\u89e3\u7a0b\u5e8f\u3002
2\u3001 \u5728\u811a\u672c\u6587\u4ef6\u4e2d\u8f93\u5165\u6c42\u89e3\u7a0b\u5e8f-->\u70b9\u51fb\u4fdd\u5b58-->\u70b9\u51fb\u8fd0\u884c\u3002
3\u3001 \u8fd9\u91cc\u9700\u8981\u5148\u7f16\u5199\u4e00\u4e2a\u51fd\u6570\u6587\u4ef6\u7528\u4e8e\u8868\u793a\u65b9\u7a0b-->\u70b9\u51fb\u4fdd\u5b58-->\u7f16\u5199\u6c42\u89e3\u7a0b\u5e8f-->\u70b9\u51fb\u4fdd\u5b58-->\u70b9\u51fb\u8fd0\u884c\u3002
4\u3001\u5728figure\u9875\u9762\u5c31\u53ef\u4ee5\u770b\u5230\u6240\u6c42\u89e3\u7684\u5fae\u5206\u65b9\u7a0b\u7684\u56fe\u5f62\u3002
或者用欧拉法,或者说简单点的如下:
>> y=dsolve('D2y=70*t-80*t^2', 'y(0)=0,y(20)=0')
y =
-20/3*t^4+35/3*t^3+146000/3*t
>> y=subs(y,'t','x')
y =
-20/3*(x)^4+35/3*(x)^3+146000/3*(x)
>> y=dsolve('D2y=70*t-80*t^2', 'y(0)=0,y(20)=0')
y =
-20/3*t^4+35/3*t^3+146000/3*t
>> y=subs(y,'t','x')
y =
-20/3*(x)^4+35/3*(x)^3+146000/3*(x)
扩展阅读：matlab好久没用打不开了 ... matlab 国产替代 ... matlab正在暂停卡住 ... 启动matlab黑框闪了一下 ... matlab一运行就崩溃 ... matlab 跳出循环 ... 程序流程图 matlab ... matlab强制停止 ... matlab gui按钮下一步 ...
