24小时热门版块排行榜

返回列表

martinwgc

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 183
帖子: 81
在线: 110.2小时
虫号: 1328530
注册: 2011-06-22
性别: GG
专业: 微生物生理与生物化学

[求助] 求助：这个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))));

回复此楼

» 猜你喜欢

求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有240人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复
笼目材料中量子自旋液体基态的证据已经有0人回复
数学教学论硕士可以读数学物理博士吗？已经有0人回复
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生已经有4人回复
澳门大学应用物理及材料工程研究院潘晖教授课题组诚招博士后已经有11人回复
求助NH4V4O10晶体的CIF文件已经有0人回复
英国全奖博士招聘-深度学习与量子物理已经有0人回复
间接带隙半导体有效质量求助已经有0人回复

» 本主题相关价值贴推荐，对您同样有帮助:

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人回复

1楼 2015-03-30 18:13:46

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

超级阿斗

铜虫 (初入文坛)

应助: 3 (幼儿园)
金币: 125.8
帖子: 10
在线: 11小时
虫号: 2907123
注册: 2014-01-01
性别: GG
专业: 激光

【答案】应助回帖

把第一个函数里边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

赞一下(1人)