我正在为不同的角度求解微分方程,我想找到导致特定值的所有角度。我在下面做了一个例子,我有很多解决方案,我想找到产生解决方案的角度说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)];