| 查看: 3097 | 回复: 0 | ||
[求助]
matlab求解隐式常微分方程组问题
|
|
我想利用matlab ode15i来求解图1中这个常微分方程组，建立相应的Matlab函数： function fun=c7eqdae(t,y,dy) fun=[dy(1)-sin(y(3)); dy(2)+cos(y(4)); dy(3)-dy(4); 1/(19-32*y(1))-4/3*pi*y(2)*(2*sin(y(4))+1)*(y(1)+2/3*y(2)*cos(y(4))*sin(y(3)))]; 将M矩阵输入给MATLAB工作空间，并在命令窗口中给出如下命令： clear all; y0=[0.3;0.027;pi/2;pi/2];dy0=[1;0;0;0] res=ode15i('c7eqdae',[0.6 0],y0,dy0); plot(res.x,res.y) 代入matlab中，提示我的t步长过小，无法求解。求大神帮忙解答一下，问题在哪？谢谢了  1.png |
回复此楼» 猜你喜欢
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有240人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有86人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
