【Matlab求解方程y=f(x);f(0)=0;y'-查字典问答网

来自孙肃清的问题

　　【Matlab求解方程y=f(x);f(0)=0;y''=(1+y'^2)^(3/2);求解区间[0,1]解得好的再加50分!】

　　Matlab求解方程y=f(x);f(0)=0;y''=(1+y'^2)^(3/2);求解区间[0,1]解得好的再加50分!

5回答

2019-10-26 20:13

我要回答

请先登录

宋福民

　　解此方程还需补充一个初始条件.

　　设此初始条件为f'(0)=-0.15

　　dy=@(x,y)[y(2);(1+y(2)^2)^(3/2)];

　　[X,Y]=ode45(dy,[01],[0-0.15]);

　　plot(X,Y(:,1),'-o')

2019-10-26 20:15:24

孙肃清

　　原来是这样啊，请问这个格式dy=@(x,y)[y(2);(1+y(1)^2)^(3/2)];方括号里面就是等式左右吧，那如果等式左边略复杂比如y(2)*sin(y(2))呢？麻烦你了

2019-10-26 20:18:09

宋福民

　　如果y''=y'*sin(y')、f'(0)=-0.15求解区间[0,100]则dy=@(x,y)[y(2);y(2)*sin(y(2))];[X,Y]=ode45(dy,[0100],[0-0.15]);plot(X,Y(:,1),'-o')

2019-10-26 20:20:55

孙肃清

　　不好意思，我的意思是原方程如果改成y‘’*sin(y‘’)==(1+y'^2)^(3/2);这样还能像上面那样计算吗？麻烦你了！

2019-10-26 20:23:36

宋福民

　　不能。如sin(y‘’)=(1+y'^2)^(3/2)，则与上面的微分方程不是同一类型。你可查阅隐式微分方程解法ode15i相关的资科。

2019-10-26 20:26:45