版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

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))));

回复此楼

» 猜你喜欢

如何从铁电相到顺电相。已经有1人回复
各位大佬，求一份最新Wien2k版本作为学习用途已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有242人回复
南科大活性流体和软物质课题组诚招2027级博士、硕士研究生和博后已经有42人回复
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向已经有24人回复
PRB投稿，Acknowledgment sent to author状态十天了已经有3人回复
PVA溶解已经有2人回复
基金申请已经有44人回复
CSC与新西兰维多利亚大学PhD奖学金项目已经有0人回复
新西兰Robinson研究所招聘CSC公派访问人员已经有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人)