| 查看: 326 | 回复: 1 | |||
martinwgc金虫 (小有名气)
|
[求助]
求助：这个matlab的代码有问题吗？ 已有1人参与
|
|
我这里有一个偏微分方程组。希望得到它在t=0,3,6,12,24,48,72,120,168的时候的各个变量的值。为啥我这个运行出来屏幕上啥也没有？请帮我看看代码，谢谢了！ function Model1 clear;clc x0=[0,0,0,100,0.4368,1.1432,0.5389,0.0511,2156.08,946.60,0,0]; [t,x]=ode45(@sac1,[0,3,6,12,24,48,72,120,168],x0); x; function dx=sac1(t,x); global E1max; global K1ad; global E1T; global E2T; global S0; global K2ad; global E2max; global K1r; global K1IG2; global K1IG; global K2r; global K2IG2; global K2IG; global K3r; global K3M; global K3IG; global KG3plus; global KG3minus; global Gcrtri; E1max=0.03257; K1ad=1.238; E1T=1.58; E2T=0.59; S0=100; K2ad=1.865; E2max=0.00102; K1r=18.86; K1IG2=0.0042; K1IG=0.0947; K2r=7.926; K2IG2=147.87; K2IG=0.01509; K3r=228.264; K3M=3.1740; K3IG=0.9617; KG3plus=0.0075; KG3minus=0.1197; Gcrtri=40; dx=zeros(12,1); x(6)=E1max*K1ad*x(5)*x(4)/(1+K1ad*x(5)); x(8)=E2max*K2ad*x(7)*x(4)/(1+K2ad*x(7)); x(9)=K1r*x(6)*x(4)^2/(S0*(1+x(2)/K1IG2+x(1)/K1IG)); x(10)=K2r*(x(6)+x(8))*x(4)^2/(S0*(1+x(2)/K2IG2+x(1)/K2IG)); x(11)=K3r*x(7)*x(2)/(K3M*(1+x(1)/K3IG)+x(2)); dx(4)=-x(9)-x(10); dx(2)=1.056*x(9)-x(11); dx(1)=1.111*x(10)+1.053*x(11)-1.071*x(12); dx(3)=x(12); x(5)=E1T-x(6); x(7)=E2T-x(8); x(12)=(KG3plus*x(1)-KG3minus*x(3))*(1/(1+exp(Gcrtri-x(1)))); |
回复此楼» 猜你喜欢
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有203人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有185人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Matlab 小白问一个简单的 带有变量运算的 问题
已经有3人回复
- Matlab用fsolve求解非线性方程，卡了三四天了跪求帮助！
已经有35人回复
- 求matlab编程代码，求帮助，谢谢！
已经有11人回复
- 请问这段小matlab代码怎么修改？
已经有10人回复
- 请问MATLAB能做出这个函数的图像吗？
已经有8人回复
- Matlab中的fmincon功能求助
已经有16人回复
- matlab解方程组求助攻
已经有4人回复
- matlab一运行就出现的问题
已经有7人回复
- matlab求助？？？
已经有6人回复
- [求助] matlab手写数字识别特征值提取问题
已经有3人回复
- MATLAB-pdepe函数使用问题
已经有5人回复
- 求教matlab线性优化求最小值的问题
已经有7人回复
- 求助：关于matlab中sym syms limit问题。
已经有7人回复
- matlab解决重心法选址问题，运行不了啊，求助各位大神
已经有10人回复
- 请matlab高手看一下这个问题怎样编程才好，附程序代码，谢谢！
已经有9人回复
- 大家爱帮我看看这段MATLAB程序为啥出错了
已经有5人回复
- matlab 拟合反应动力学参数结果很差。大家帮忙看一下
已经有14人回复
- matlab代码 定义函数问题
已经有3人回复
- 同样的MATLAB程序，为何在别人的电脑上运行不了
已经有21人回复
- matlab精度问题
已经有4人回复
- 【求助】求解隐函数
已经有6人回复
【答案】应助回帖
|
把第一个函数里边x后面的分号去掉。以下是显示的结果 x = 1.0e+03 * 0 0 0 0.1000 0.0004 0.0011 0.0005 0.0001 2.1561 0.9466 0 0 0.0178 0.0013 0 0.0828 0.0004 0.0011 0.0005 0.0001 2.1561 0.9466 0 0 0.0255 0.0013 0 0.0759 0.0004 0.0011 0.0005 0.0001 2.1561 0.9466 0 0 0.0354 0.0011 0 0.0671 0.0004 0.0011 0.0005 0.0001 2.1561 0.9466 0 0 0.0470 0.0008 0 0.0570 0.0004 0.0011 0.0005 0.0001 2.1561 0.9466 0 0 0.0591 0.0005 0 0.0464 0.0004 0.0011 0.0005 0.0001 2.1561 0.9466 0 0 0.0658 0.0004 0 0.0405 0.0004 0.0011 0.0005 0.0001 2.1561 0.9466 0 0 0.0736 0.0002 0 0.0336 0.0004 0.0011 0.0005 0.0001 2.1561 0.9466 0 0 0.0782 0.0001 0 0.0295 0.0004 0.0011 0.0005 0.0001 2.1561 0.9466 0 0 |
2楼2015-04-03 15:55:22
投票:
10