Matlab求解方程 y=f(x);f(0)=0; y''=(1+y'^2)^(3/2); 求解区间[0,1] 解得好的再加50分!
1个回答
展开全部
解此方程还需补充一个初始条件。
设此初始条件为f '(0)=-0.15
dy = @(x,y) [y(2);(1 + y(2)^2)^(3/2)];
[X,Y] = ode45(dy,[0 1],[0 -0.15]);
plot(X,Y(:,1),'-o')
设此初始条件为f '(0)=-0.15
dy = @(x,y) [y(2);(1 + y(2)^2)^(3/2)];
[X,Y] = ode45(dy,[0 1],[0 -0.15]);
plot(X,Y(:,1),'-o')
更多追问追答
追问
原来是这样啊,请问这个格式dy = @(x,y) [y(2);(1 + y(1)^2)^(3/2)]; 方括号里面就是等式左右吧,那如果等式左边略复杂比如y(2)*sin(y(2))呢?麻烦你了
追答
如果y''= y'*sin(y')、f '(0) = -0.15 求解区间[0,100]
则
dy = @(x,y) [y(2);y(2)*sin(y(2))];
[X,Y] = ode45(dy,[0 100],[0 -0.15]);
plot(X,Y(:,1),'-o')
推荐律师服务:
若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询