各位大佬,在下专业是土木,编程太深奥的不懂,现在遇到的问题有关于我的毕业论文(很重要,问了导师,导师也不清楚),自己也想不明白,百度了好多。请求吧友抽出宝贵时间帮助我!
我整理WORD,直接上图:
除了ODE15S、45,其他的我也尝试了,都不行。我也有尝试修改options,但小白也只能照着网上试试,具体也不知道怎么改、或者能否改。
代码:
function y=eq2(t,n)
lam=0.2;%λ阻尼系数
A=0.5;
y=[0.5*(-lam*n(1)-n(1).^2*A*cos(n(2))+lam*A*sin(n(2))+A*cos(n(2)))/(1+lam^2-4*n(1).^2+3*n(1).^4);...
0.5*(-3*n(1).^4+n(1).^2+3*A*n(1)*sin(n(2))-lam^2+(A*lam*cos(n(2))-A*sin(n(2)))/n(1))/(1+lam^2-4*n(1).^2+3*n(1).^4)];
end
for i=0:0.1:3
tspan=0:0.1:10;
[t,n]=ode15s('eq2',tspan,[j,i]);%初始条件【N,θ】,j=0.1,0.2,……1.2
plot(n(:,2)',n(:,1)');hold on
end
%其中,关于J我并没有设置循环,一般自行赋值检验。
最后,感谢各位吧友赐教哈。
我整理WORD,直接上图:
除了ODE15S、45,其他的我也尝试了,都不行。我也有尝试修改options,但小白也只能照着网上试试,具体也不知道怎么改、或者能否改。
代码:
function y=eq2(t,n)
lam=0.2;%λ阻尼系数
A=0.5;
y=[0.5*(-lam*n(1)-n(1).^2*A*cos(n(2))+lam*A*sin(n(2))+A*cos(n(2)))/(1+lam^2-4*n(1).^2+3*n(1).^4);...
0.5*(-3*n(1).^4+n(1).^2+3*A*n(1)*sin(n(2))-lam^2+(A*lam*cos(n(2))-A*sin(n(2)))/n(1))/(1+lam^2-4*n(1).^2+3*n(1).^4)];
end
for i=0:0.1:3
tspan=0:0.1:10;
[t,n]=ode15s('eq2',tspan,[j,i]);%初始条件【N,θ】,j=0.1,0.2,……1.2
plot(n(:,2)',n(:,1)');hold on
end
%其中,关于J我并没有设置循环,一般自行赋值检验。
最后,感谢各位吧友赐教哈。