| 查看: 2363 | 回复: 15 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
matlab解微分方程组
|
||
|
哪位高手指点一下matlab解微分方程组！ 图片是微分方程组。 下面是我根据网上查到的解法写的程序，能运行，但是计算出来的图像完全不对。 主函数： tspan=[0:0.01:2];%求解时间为2秒也可以 tspan=[0 2]表示自动时间步长 x0=[0;0;2.4*10^(-4);20;1.15*10^(-7);0];%初值 [t,x]=ode45(@ff,tspan,x0);%用龙格库塔法求解 figure plot(t,x(:,1)) figure plot(t,x(:,2)) figure plot(t,x(:,3)) figure plot(t,x(:,4)) figure plot(t,x(:,5)) figure plot(t,x(:,6)) M文件 %微分方程 function dx=ff(t,x) dx=zeros(6,1); %x(1)=X1 %x(2)=X2 %x(3)=X3 %x(4)=X4 %x(5)=X5 %x(6)=X6 k1u=1.5*10^(-7) k1b=1.0*10^(-4) k2=1.3*10^(7) k3=7*10^(-1) k4=1.0*10^(10) k5=2.6*10^(9) k6=3*10^(3) Xc=0.5 MMT=0 dx(1)=2*k1u*x(3)+k1b*x(3)^2-k2*x(1)*x(5)+k3*x(2)*x(4)-2*k4*x(1)^2-k5*x(1)*x(2); dx(2)=k1b*x(3)^2+k2*x(1)*x(5)-k3*x(2)*x(4)-k5*x(1)*x(2)-2*k6*x(2); dx(3)=-k1u*x(3)-2*k1b*x(3)^2+k3*x(2)*x(4); dx(4)=-k1u*x(3)-k1b*x(3)^2-k3*x(2)*x(4); dx(5)=-k2*x(1)*x(5)+k6*x(2)^2; dx(6)=(k1u*x(3)+k1b*x(3)^2)*(1-Xc-MMT)*1.1111; end       |
回复此楼» 猜你喜欢
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有204人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
-
笼目材料中量子自旋液体基态的证据
已经有0人回复
-
数学教学论硕士可以读数学物理博士吗？
已经有0人回复
-
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生
已经有4人回复
-
澳门大学 应用物理及材料工程研究院 潘晖教授课题组诚招博士后
已经有11人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 用matlab求解一偏微分方程，求助！
已经有1人回复
- 流体力学一维定常不可压缩流体ns方程，已知压差变化，怎么求解速度u关于x的常微分方程
已经有0人回复
- matlab书籍
已经有103人回复
- 求大神帮忙解题啊！！！
已经有4人回复
- Matlab求解微分方程里面有参数带入怎么求解，求救！！！
已经有0人回复
- 求助CFD软件可以自定义模拟算法
已经有0人回复
- 【紧急求助】 带有微分平方项的ODE求解
已经有11人回复
- 在控制器设计中关于大家如何给定？如果给定是方波怎么给？
已经有2人回复
- 请教高手：PDE工具箱能解两个偏微分方程相互迭代求解吗
已经有2人回复
- 【求助】求解形如y''+a*y'^2+by'+c=0 二阶微分方程
已经有14人回复
- 【资料】一些数学实验,有程序源代码
已经有7人回复
3楼2011-12-02 09:47:19
5楼2011-12-03 11:26:16
7楼2011-12-10 10:43:54
9楼2011-12-10 11:09:40
10楼2011-12-10 11:14:33
12楼2011-12-10 15:04:44
13楼2011-12-10 15:07:03