[未答复] ODE15S函数无法求解的刚性微分方程

[复制链接]
尔疋 发表于 2021-5-2 19:28:26
各位朋友,在下专业是土木,编程太深奥的不懂,现在遇到的问题有关于我的毕业论文(很重要,问了导师,导师也不清楚),自己也想不明白,百度了好多。请求大佬抽出宝贵时间帮助我!(Matlab R2018b)
我整理WORD,直接上图:
1.png
2.png
3.png
程序:
m文件:
[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
/code]
函数文件:
[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
/code]

您需要登录后才可以回帖 登录 | 注册

本版积分规则

相关帖子
热门教程
站长推荐
快速回复 返回顶部 返回列表