Matlab微分方程,哪一个给出了正确的解?

时间:2014-04-29 18:42:44

标签: matlab ode differential-equations

我正在为不同的角度求解微分方程,我想找到导致特定值的所有角度。我在下面做了一个例子,我有很多解决方案,我想找到产生解决方案的角度说10。

是否有像fsolve这样的方法返回角度?

代码:

xstart = 1;
xend = 5;
table = [];

for angle = 0:30:360
    ustart = [2 sind(angle)]';
    odefun = @(x,u) f(x,u,angle);
    [xout uout]= ode45(odefun,[xstart xend],ustart);
    yend = uout(end,1);
    table(end+1) = yend;

end
disp(table)

我的功能文件:

function uprime = f(x,u,a)
uprime = [u(2)
     1-x*u(1)*sind(a)];

0 个答案:

没有答案