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

登录注册

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

24小时热门版块排行榜

南方科技大学公共卫生及应急管理学院2025级博士研究生招生报考通知

返回列表

vista朱

木虫 (著名写手)

应助: 39 (小学生)
金币: 5493
散金: 348
红花: 42
帖子: 1891
在线: 331.9小时
虫号: 1930187
注册: 2012-08-08
性别: GG
专业: 动力学与控制

[求助] 紧急求助已有2人参与

用 ode45解下面非线性方程组，
clc
clear
syms  A1 A2 A3 B1 B2 B3 x2 x3 x4 x5 x6 x7 pho1 pho2 pho3 mdot_in h_in T_a1 T_a2 T_a3 T_r1 T_r2 T_r3 alpha_o1 alpha_o2 alpha_o3 alpha_i1 alpha_i2 alpha_i3;
a=[-2.8e-3,2.8e-3*pho1,0,0,0,0,0
2.8e-3*A1*(x2-x3)-2.8e-3, 2.8e-3*B1*(x2-x3), 2.8e-3*pho2,0,0,0,0
(2.8e-3*A2+2.8e-3*A1)*(x3-x4)-2.8e-3,2.8e-3*B1*(x3-x4),2.8e-3*B2*(x3-x4),2.8e-3*pho3,0,0,0
2.8e-3*(A1+A2+A3),2.8e-3*B1,2.8e-3*B2,2.8e-3*B3,0,0,0
0,0,0,0,2377,0,0
0,0,0,0,0,2377,0
0,0,0,0,0,0,2377];
c=[(mdot_in*(h_in-x2)+1.3*alpha_i1*(x5-T_r1));(mdot_in*(x2-x3)+1.3*alpha_i2*(x6-T_r2));(mdot_in*(x3-x4)+1.3*alpha_i3*(x7-T_r3));0;30.58*alpha_o1*(T_a3-x5)-1.3*alpha_i1*(x5-T_r1);30.58*alpha_o2*(T_a2-x6)-1.3*alpha_i2*(x6-T_r2);30.58*alpha_o3*(T_a1-x7)-1.3*alpha_i3*(x7-T_r3)];
bb=inv(a)*c;

用传统的ode方法编写
即
  dx(1)=  (2500*((13*alpha_i1*(T_r1 - x(5)))/10 - mdot_in*(h_in - x(2)))*(B1*pho2*pho3 - B1*B2*pho3*x(2) + B1*B2*pho3*x(3) - B1*B3*pho2*x(3) + B1*B3*pho2*x(4) - B1*B2*B3*x(3)^2 + B1*B2*B3*x(2)*x(3) - B1*B2*B3*x(2)*x(4) + B1*B2*B3*x(3)*x(4)))/(7*(B1*pho2*pho3 + B2*pho1*pho3 + B3*pho1*pho2 - B1*B2*pho3*x(2) + B1*B2*pho3*x(3) - B1*B3*pho2*x(3) - B2*B3*pho1*x(3) + B1*B3*pho2*x(4) + B2*B3*pho1*x(4) + A1*pho1*pho2*pho3 + A2*pho1*pho2*pho3 + A3*pho1*pho2*pho3 - B1*B2*B3*x(3)^2 + B1*B2*B3*x(2)*x(3) - B1*B2*B3*x(2)*x(4) + B1*B2*B3*x(3)*x(4) - A1*B2*pho1*pho3*x(2) + A1*B2*pho1*pho3*x(3) - A1*B3*pho1*pho2*x(3) + A1*B3*pho1*pho2*x(4) - A2*B3*pho1*pho2*x(3) + A2*B3*pho1*pho2*x(4) - A1*B2*B3*pho1*x(3)^2 + A1*B2*B3*pho1*x(2)*x(3) - A1*B2*B3*pho1*x(2)*x(4) + A1*B2*B3*pho1*x(3)*x(4))) + (2500*((13*alpha_i2*(T_r2 - x(6)))/10 - mdot_in*(x(2) - x(3)))*(B2*pho1*pho3 - B2*B3*pho1*x(3) + B2*B3*pho1*x(4)))/(7*(B1*pho2*pho3 + B2*pho1*pho3 + B3*pho1*pho2 - B1*B2*pho3*x(2) + B1*B2*pho3*x(3) - B1*B3*pho2*x(3)- B2*B3*pho1*x(3) + B1*B3*pho2*x(4) + B2*B3*pho1*x(4) + A1*pho1*pho2*pho3 + A2*pho1*pho2*pho3 + A3*pho1*pho2*pho3 - B1*B2*B3*x(3)^2 + B1*B2*B3*x(2)*x(3) - B1*B2*B3*x(2)*x(4) + B1*B2*B3*x(3)*x(4) - A1*B2*pho1*pho3*x(2) + A1*B2*pho1*pho3*x(3) - A1*B3*pho1*pho2*x(3) + A1*B3*pho1*pho2*x(4) - A2*B3*pho1*pho2*x(3) + A2*B3*pho1*pho2*x(4) - A1*B2*B3*pho1*x(3)^2 + A1*B2*B3*pho1*x(2)*x(3) - A1*B2*B3*pho1*x(2)*x(4) + A1*B2*B3*pho1*x(3)*x(4))) + (2500*B3*pho1*pho2*((13*alpha_i3*(T_r3 - x(7)))/10 - mdot_in*(x(3) - x(4))))/(7*(B1*pho2*pho3 + B2*pho1*pho3 + B3*pho1*pho2 - B1*B2*pho3*x(2) + B1*B2*pho3*x(3) - B1*B3*pho2*x(3) - B2*B3*pho1*x(3) + B1*B3*pho2*x(4) + B2*B3*pho1*x(4) + A1*pho1*pho2*pho3 + A2*pho1*pho2*pho3 + A3*pho1*pho2*pho3 - B1*B2*B3*x(3)^2 + B1*B2*B3*x(2)*x(3) - B1*B2*B3*x(2)*x(4) + B1*B2*B3*x(3)*x(4) - A1*B2*pho1*pho3*x(2) + A1*B2*pho1*pho3*x(3) - A1*B3*pho1*pho2*x(3) + A1*B3*pho1*pho2*x(4) - A2*B3*pho1*pho2*x(3) + A2*B3*pho1*pho2*x(4) - A1*B2*B3*pho1*x(3)^2 + A1*B2*B3*pho1*x(2)*x(3) - A1*B2*B3*pho1*x(2)*x(4) + A1*B2*B3*pho1*x(3)*x(4)));
dx(2)=  (2500*((13*alpha_i2*(T_r2 - x(6)))/10 - mdot_in*(x(2) - x(3)))*(B2*pho3 - B2*B3*x(3) + B2*B3*x(4)))/(7*(B1*pho2*pho3 + B2*pho1*pho3 + B3*pho1*pho2 - B1*B2*pho3*x(2) + B1*B2*pho3*x(3) - B1*B3*pho2*x(3) - B2*B3*pho1*x(3) + B1*B3*pho2*x(4) + B2*B3*pho1*x(4) + A1*pho1*pho2*pho3 + A2*pho1*pho2*pho3 + A3*pho1*pho2*pho3 - B1*B2*B3*x(3)^2 + B1*B2*B3*x(2)*x(3) - B1*B2*B3*x(2)*x(4) + B1*B2*B3*x(3)*x(4) - A1*B2*pho1*pho3*x(2) + A1*B2*pho1*pho3*x(3) - A1*B3*pho1*pho2*x(3) + A1*B3*pho1*pho2*x(4) - A2*B3*pho1*pho2*x(3) + A2*B3*pho1*pho2*x(4) - A1*B2*B3*pho1*x(3)^2 + A1*B2*B3*pho1*x(2)*x(3) - A1*B2*B3*pho1*x(2)*x(4) + A1*B2*B3*pho1*x(3)*x(4))) - (2500*((13*alpha_i1*(T_r1 - x(5)))/10 - mdot_in*(h_in - x(2)))*(B2*pho3 + B3*pho2 - B2*B3*x(3) + B2*B3*x(4) + A1*pho2*pho3 + A2*pho2*pho3 + A3*pho2*pho3 - A1*B2*pho3*x(2) + A1*B2*pho3*x(3) - A1*B3*pho2*x(3) + A1*B3*pho2*x(4) - A2*B3*pho2*x(3) + A2*B3*pho2*x(4) - A1*B2*B3*x(3)^2 + A1*B2*B3*x(2)*x(3) - A1*B2*B3*x(2)*x(4) + A1*B2*B3*x(3)*x(4)))/(7*(B1*pho2*pho3 + B2*pho1*pho3 + B3*pho1*pho2 - B1*B2*pho3*x(2) + B1*B2*pho3*x(3) - B1*B3*pho2*x(3) - B2*B3*pho1*x(3) + B1*B3*pho2*x(4) + B2*B3*pho1*x(4) + A1*pho1*pho2*pho3 + A2*pho1*pho2*pho3 + A3*pho1*pho2*pho3 - B1*B2*B3*x(3)^2 + B1*B2*B3*x(2)*x(3) - B1*B2*B3*x(2)*x(4) + B1*B2*B3*x(3)*x(4) - A1*B2*pho1*pho3*x(2) + A1*B2*pho1*pho3*x(3) - A1*B3*pho1*pho2*x(3) + A1*B3*pho1*pho2*x(4) - A2*B3*pho1*pho2*x(3) + A2*B3*pho1*pho2*x(4) - A1*B2*B3*pho1*x(3)^2 + A1*B2*B3*pho1*x(2)*x(3) - A1*B2*B3*pho1*x(2)*x(4) + A1*B2*B3*pho1*x(3)*x(4))) + (2500*B3*pho2*((13*alpha_i3*(T_r3 - x(7)))/10 - mdot_in*(x(3) - x(4))))/(7*(B1*pho2*pho3 + B2*pho1*pho3 + B3*pho1*pho2 - B1*B2*pho3*x(2) + B1*B2*pho3*x(3) - B1*B3*pho2*x(3) - B2*B3*pho1*x(3) + B1*B3*pho2*x(4) + B2*B3*pho1*x(4) + A1*pho1*pho2*pho3 + A2*pho1*pho2*pho3 + A3*pho1*pho2*pho3 - B1*B2*B3*x(3)^2 + B1*B2*B3*x(2)*x(3) - B1*B2*B3*x(2)*x(4) + B1*B2*B3*x(3)*x(4) - A1*B2*pho1*pho3*x(2) + A1*B2*pho1*pho3*x(3) - A1*B3*pho1*pho2*x(3) + A1*B3*pho1*pho2*x(4) - A2*B3*pho1*pho2*x(3) + A2*B3*pho1*pho2*x(4) - A1*B2*B3*pho1*x(3)^2 + A1*B2*B3*pho1*x(2)*x(3) - A1*B2*B3*pho1*x(2)*x(4) + A1*B2*B3*pho1*x(3)*x(4)));

不仅程序冗长，而且debug太麻烦！！！（主要是自己太笨，工作量太大），没其他好法子，所以希望大家给些好的方法，自己对ode45很多用法不太了解，谢谢大家了。

回复此楼

» 猜你喜欢

三甲基碘化亚砜的氧化反应已经有4人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有5人回复
孩子确诊有中度注意力缺陷已经有12人回复
2025冷门绝学什么时候出结果已经有3人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复
论文投稿，期刊推荐已经有4人回复
硕士和导师闹得不愉快已经有13人回复

dream

1楼2018-10-06 12:08:51

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

天天进步啊

木虫 (著名写手)

应助: 87 (初中生)
金币: 6058.7
散金: 101
红花: 42
沙发: 9
帖子: 1793
在线: 1005.5小时
虫号: 3689648
注册: 2015-02-15
性别: GG
专业: 理论和计算化学

★
jjdg: 金币+1, 感谢参与 2018-10-09 21:15:31

楼主应该直接把m文件贴出来
而且源代码里面要有注释啊

赞一下

回复此楼

分子筛、分子动力学模拟，有问题可咨询我的知乎主页 https://www.zhihu.com/people/rao777

2楼2018-10-06 14:43:09

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

vista朱

木虫 (著名写手)

应助: 39 (小学生)
金币: 5493
散金: 348
红花: 42
帖子: 1891
在线: 331.9小时
虫号: 1930187
注册: 2012-08-08
性别: GG
专业: 动力学与控制

★
jjdg: 金币+1, 感谢参与 2018-10-09 21:15:44

引用回帖:

2楼: Originally posted by 天天进步啊 at 2018-10-06 14:43:09
楼主应该直接把m文件贴出来
而且源代码里面要有注释啊

哦哦不好意思哒
这是代码

CODE:

%%

%%

 clc

 clear

   syms  A1 A2 A3 B1 B2 B3 x2 x3 x4 x5 x6 x7 pho1 pho2 pho3 mdot_in h_in T_a1 T_a2 T_a3 T_r1 T_r2 T_r3 alpha_o1 alpha_o2 alpha_o3 alpha_i1 alpha_i2 alpha_i3;

 a=[-2.8e-3,2.8e-3*pho1,0,0,0,0,0

    2.8e-3*A1*(x2-x3)-2.8e-3, 2.8e-3*B1*(x2-x3), 2.8e-3*pho2,0,0,0,0

    (2.8e-3*A2+2.8e-3*A1)*(x3-x4)-2.8e-3,2.8e-3*B1*(x3-x4),2.8e-3*B2*(x3-x4),2.8e-3*pho3,0,0,0

    2.8e-3*(A1+A2+A3),2.8e-3*B1,2.8e-3*B2,2.8e-3*B3,0,0,0

    0,0,0,0,2377,0,0

    0,0,0,0,0,2377,0

    0,0,0,0,0,0,2377];

c=[(mdot_in*(h_in-x2)+1.3*alpha_i1*(x5-T_r1));(mdot_in*(x2-x3)+1.3*alpha_i2*(x6-T_r2));(mdot_in*(x3-x4)+1.3*alpha_i3*(x7-T_r3));0;30.58*alpha_o1*(T_a3-x5)-1.3*alpha_i1*(x5-T_r1);30.58*alpha_o2*(T_a2-x6)-1.3*alpha_i2*(x6-T_r2);30.58*alpha_o3*(T_a1-x7)-1.3*alpha_i3*(x7-T_r3)];

 bb=inv(a)*c;     %求取 dx(1),dx(2),...dx(6)的值

%%

ode 45 函数

function dx=li8(t,x) 

global mdot_in h_in A1 A2 A3 B1 B2 B3 T_a1 T_a2 T_a3 T_r1  T_r2  T_r3  pho1 pho2 pho3 alpha_i1 alpha_i2 alpha_i3  alpha_o1 alpha_o2 alpha_o3 %A1=diff(Rho1 p),A2=diff(Rho2 p),A3=diff(Rho3 p),B1=diff(Rho1 h1),B2=diff(Rho2 h2),B3=diff(Rho3 h3)

dx= zeros(1,7);                %x(1)=P,x(2)=h1,x(3)=h2,x(4)=h3,x(5)=T_w1,x(6)=T_w2,x(7)=T_w3

dx(1)=  (2500*((13*alpha_i1*(T_r1 - x(5)))/10 - mdot_in*(h_in - x(2)))*(B1*pho2*pho3 - B1*B2*pho3*x(2) + B1*B2*pho3*x(3) - B1*B3*pho2*x(3) + B1*B3*pho2*x(4) - B1*B2*B3*x(3)^2 + B1*B2*B3*x(2)*x(3) - B1*B2*B3*x(2)*x(4) + B1*B2*B3*x(3)*x(4)))/(7*(B1*pho2*pho3 + B2*pho1*pho3 + B3*pho1*pho2 - B1*B2*pho3*x(2) + B1*B2*pho3*x(3) - B1*B3*pho2*x(3) - B2*B3*pho1*x(3) + B1*B3*pho2*x(4) + B2*B3*pho1*x(4) + A1*pho1*pho2*pho3 + A2*pho1*pho2*pho3 + A3*pho1*pho2*pho3 - B1*B2*B3*x(3)^2 + B1*B2*B3*x(2)*x(3) - B1*B2*B3*x(2)*x(4) + B1*B2*B3*x(3)*x(4) - A1*B2*pho1*pho3*x(2) + A1*B2*pho1*pho3*x(3) - A1*B3*pho1*pho2*x(3) + A1*B3*pho1*pho2*x(4) - A2*B3*pho1*pho2*x(3) + A2*B3*pho1*pho2*x(4) - A1*B2*B3*pho1*x(3)^2 + A1*B2*B3*pho1*x(2)*x(3) - A1*B2*B3*pho1*x(2)*x(4) + A1*B2*B3*pho1*x(3)*x(4))) + (2500*((13*alpha_i2*(T_r2 - x(6)))/10 - mdot_in*(x(2) - x(3)))*(B2*pho1*pho3 - B2*B3*pho1*x(3) + B2*B3*pho1*x(4)))/(7*(B1*pho2*pho3 + B2*pho1*pho3 + B3*pho1*pho2 - B1*B2*pho3*x(2) + B1*B2*pho3*x(3) - B1*B3*pho2*x(3)- B2*B3*pho1*x(3) + B1*B3*pho2*x(4) + B2*B3*pho1*x(4) + A1*pho1*pho2*pho3 + A2*pho1*pho2*pho3 + A3*pho1*pho2*pho3 - B1*B2*B3*x(3)^2 + B1*B2*B3*x(2)*x(3) - B1*B2*B3*x(2)*x(4) + B1*B2*B3*x(3)*x(4) - A1*B2*pho1*pho3*x(2) + A1*B2*pho1*pho3*x(3) - A1*B3*pho1*pho2*x(3) + A1*B3*pho1*pho2*x(4) - A2*B3*pho1*pho2*x(3) + A2*B3*pho1*pho2*x(4) - A1*B2*B3*pho1*x(3)^2 + A1*B2*B3*pho1*x(2)*x(3) - A1*B2*B3*pho1*x(2)*x(4) + A1*B2*B3*pho1*x(3)*x(4))) + (2500*B3*pho1*pho2*((13*alpha_i3*(T_r3 - x(7)))/10 - mdot_in*(x(3) - x(4))))/(7*(B1*pho2*pho3 + B2*pho1*pho3 + B3*pho1*pho2 - B1*B2*pho3*x(2) + B1*B2*pho3*x(3) - B1*B3*pho2*x(3) - B2*B3*pho1*x(3) + B1*B3*pho2*x(4) + B2*B3*pho1*x(4) + A1*pho1*pho2*pho3 + A2*pho1*pho2*pho3 + A3*pho1*pho2*pho3 - B1*B2*B3*x(3)^2 + B1*B2*B3*x(2)*x(3) - B1*B2*B3*x(2)*x(4) + B1*B2*B3*x(3)*x(4) - A1*B2*pho1*pho3*x(2) + A1*B2*pho1*pho3*x(3) - A1*B3*pho1*pho2*x(3) + A1*B3*pho1*pho2*x(4) - A2*B3*pho1*pho2*x(3) + A2*B3*pho1*pho2*x(4) - A1*B2*B3*pho1*x(3)^2 + A1*B2*B3*pho1*x(2)*x(3) - A1*B2*B3*pho1*x(2)*x(4) + A1*B2*B3*pho1*x(3)*x(4)));

dx(2)=  (2500*((13*alpha_i2*(T_r2 - x(6)))/10 - mdot_in*(x(2) - x(3)))*(B2*pho3 - B2*B3*x(3) + B2*B3*x(4)))/(7*(B1*pho2*pho3 + B2*pho1*pho3 + B3*pho1*pho2 - B1*B2*pho3*x(2) + B1*B2*pho3*x(3) - B1*B3*pho2*x(3) - B2*B3*pho1*x(3) + B1*B3*pho2*x(4) + B2*B3*pho1*x(4) + A1*pho1*pho2*pho3 + A2*pho1*pho2*pho3 + A3*pho1*pho2*pho3 - B1*B2*B3*x(3)^2 + B1*B2*B3*x(2)*x(3) - B1*B2*B3*x(2)*x(4) + B1*B2*B3*x(3)*x(4) - A1*B2*pho1*pho3*x(2) + A1*B2*pho1*pho3*x(3) - A1*B3*pho1*pho2*x(3) + A1*B3*pho1*pho2*x(4) - A2*B3*pho1*pho2*x(3) + A2*B3*pho1*pho2*x(4) - A1*B2*B3*pho1*x(3)^2 + A1*B2*B3*pho1*x(2)*x(3) - A1*B2*B3*pho1*x(2)*x(4) + A1*B2*B3*pho1*x(3)*x(4))) - (2500*((13*alpha_i1*(T_r1 - x(5)))/10 - mdot_in*(h_in - x(2)))*(B2*pho3 + B3*pho2 - B2*B3*x(3) + B2*B3*x(4) + A1*pho2*pho3 + A2*pho2*pho3 + A3*pho2*pho3 - A1*B2*pho3*x(2) + A1*B2*pho3*x(3) - A1*B3*pho2*x(3) + A1*B3*pho2*x(4) - A2*B3*pho2*x(3) + A2*B3*pho2*x(4) - A1*B2*B3*x(3)^2 + A1*B2*B3*x(2)*x(3) - A1*B2*B3*x(2)*x(4) + A1*B2*B3*x(3)*x(4)))/(7*(B1*pho2*pho3 + B2*pho1*pho3 + B3*pho1*pho2 - B1*B2*pho3*x(2) + B1*B2*pho3*x(3) - B1*B3*pho2*x(3) - B2*B3*pho1*x(3) + B1*B3*pho2*x(4) + B2*B3*pho1*x(4) + A1*pho1*pho2*pho3 + A2*pho1*pho2*pho3 + A3*pho1*pho2*pho3 - B1*B2*B3*x(3)^2 + B1*B2*B3*x(2)*x(3) - B1*B2*B3*x(2)*x(4) + B1*B2*B3*x(3)*x(4) - A1*B2*pho1*pho3*x(2) + A1*B2*pho1*pho3*x(3) - A1*B3*pho1*pho2*x(3) + A1*B3*pho1*pho2*x(4) - A2*B3*pho1*pho2*x(3) + A2*B3*pho1*pho2*x(4) - A1*B2*B3*pho1*x(3)^2 + A1*B2*B3*pho1*x(2)*x(3) - A1*B2*B3*pho1*x(2)*x(4) + A1*B2*B3*pho1*x(3)*x(4))) + (2500*B3*pho2*((13*alpha_i3*(T_r3 - x(7)))/10 - mdot_in*(x(3) - x(4))))/(7*(B1*pho2*pho3 + B2*pho1*pho3 + B3*pho1*pho2 - B1*B2*pho3*x(2) + B1*B2*pho3*x(3) - B1*B3*pho2*x(3) - B2*B3*pho1*x(3) + B1*B3*pho2*x(4) + B2*B3*pho1*x(4) + A1*pho1*pho2*pho3 + A2*pho1*pho2*pho3 + A3*pho1*pho2*pho3 - B1*B2*B3*x(3)^2 + B1*B2*B3*x(2)*x(3) - B1*B2*B3*x(2)*x(4) + B1*B2*B3*x(3)*x(4) - A1*B2*pho1*pho3*x(2) + A1*B2*pho1*pho3*x(3) - A1*B3*pho1*pho2*x(3) + A1*B3*pho1*pho2*x(4) - A2*B3*pho1*pho2*x(3) + A2*B3*pho1*pho2*x(4) - A1*B2*B3*pho1*x(3)^2 + A1*B2*B3*pho1*x(2)*x(3) - A1*B2*B3*pho1*x(2)*x(4) + A1*B2*B3*pho1*x(3)*x(4)));

dx(3)= (2500*((13*alpha_i1*(T_r1 - x(5)))/10 - mdot_in*(h_in - x(2)))*(B1*pho3 + B1*B3*x(2) - 2*B1*B3*x(3) + B1*B3*x(4) + A2*B1*pho3*x(2) - A2*B1*pho3*x(3) + A3*B1*pho3*x(2) - A3*B1*pho3*x(3) + A2*B1*B3*x(3)^2 - A2*B1*B3*x(2)*x(3) + A2*B1*B3*x(2)*x(4) - A2*B1*B3*x(3)*x(4)))/(7*(B1*pho2*pho3 + B2*pho1*pho3 + B3*pho1*pho2 - B1*B2*pho3*x(2) + B1*B2*pho3*x(3) - B1*B3*pho2*x(3) - B2*B3*pho1*x(3) + B1*B3*pho2*x(4) + B2*B3*pho1*x(4) + A1*pho1*pho2*pho3 + A2*pho1*pho2*pho3 + A3*pho1*pho2*pho3 - B1*B2*B3*x(3)^2 + B1*B2*B3*x(2)*x(3) - B1*B2*B3*x(2)*x(4) + B1*B2*B3*x(3)*x(4) - A1*B2*pho1*pho3*x(2) + A1*B2*pho1*pho3*x(3) - A1*B3*pho1*pho2*x(3) + A1*B3*pho1*pho2*x(4) - A2*B3*pho1*pho2*x(3) + A2*B3*pho1*pho2*x(4) - A1*B2*B3*pho1*x(3)^2 + A1*B2*B3*pho1*x(2)*x(3) - A1*B2*B3*pho1*x(2)*x(4) + A1*B2*B3*pho1*x(3)*x(4))) - (2500*((13*alpha_i2*(T_r2 - x(6)))/10 - mdot_in*(x(2) - x(3)))*(B1*pho3 + B3*pho1 - B1*B3*x(3) + B1*B3*x(4) + A1*pho1*pho3 + A2*pho1*pho3 + A3*pho1*pho3 - A1*B3*pho1*x(3) + A1*B3*pho1*x(4) - A2*B3*pho1*x(3) + A2*B3*pho1*x(4)))/(7*(B1*pho2*pho3 + B2*pho1*pho3 + B3*pho1*pho2 - B1*B2*pho3*x(2) + B1*B2*pho3*x(3) - B1*B3*pho2*x(3) - B2*B3*pho1*x(3) + B1*B3*pho2*x(4) + B2*B3*pho1*x(4) + A1*pho1*pho2*pho3 + A2*pho1*pho2*pho3 + A3*pho1*pho2*pho3 - B1*B2*B3*x(3)^2 + B1*B2*B3*x(2)*x(3) - B1*B2*B3*x(2)*x(4) + B1*B2*B3*x(3)*x(4) - A1*B2*pho1*pho3*x(2) + A1*B2*pho1*pho3*x(3) - A1*B3*pho1*pho2*x(3) + A1*B3*pho1*pho2*x(4) - A2*B3*pho1*pho2*x(3) + A2*B3*pho1*pho2*x(4) - A1*B2*B3*pho1*x(3)^2 + A1*B2*B3*pho1*x(2)*x(3) - A1*B2*B3*pho1*x(2)*x(4) + A1*B2*B3*pho1*x(3)*x(4))) + (2500*B3*((13*alpha_i3*(T_r3 - x(7)))/10 - mdot_in*(x(3) - x(4)))*(pho1 - B1*x(2) + B1*x(3) - A1*pho1*x(2) + A1*pho1*x(3)))/(7*(B1*pho2*pho3 + B2*pho1*pho3 + B3*pho1*pho2 - B1*B2*pho3*x(2) + B1*B2*pho3*x(3) - B1*B3*pho2*x(3) - B2*B3*pho1*x(3) + B1*B3*pho2*x(4) + B2*B3*pho1*x(4) + A1*pho1*pho2*pho3 + A2*pho1*pho2*pho3 + A3*pho1*pho2*pho3 - B1*B2*B3*x(3)^2 + B1*B2*B3*x(2)*x(3) - B1*B2*B3*x(2)*x(4) + B1*B2*B3*x(3)*x(4) - A1*B2*pho1*pho3*x(2) + A1*B2*pho1*pho3*x(3) - A1*B3*pho1*pho2*x(3) + A1*B3*pho1*pho2*x(4) - A2*B3*pho1*pho2*x(3) + A2*B3*pho1*pho2*x(4) - A1*B2*B3*pho1*x(3)^2 + A1*B2*B3*pho1*x(2)*x(3) - A1*B2*B3*pho1*x(2)*x(4) + A1*B2*B3*pho1*x(3)*x(4)));

dx(4)= (2500*((13*alpha_i2*(T_r2 - x(6)))/10 - mdot_in*(x(2) - x(3)))*(B2*pho1 + A3*B2*pho1*x(3) - A3*B2*pho1*x(4)))/(7*(B1*pho2*pho3 + B2*pho1*pho3 + B3*pho1*pho2 - B1*B2*pho3*x(2) + B1*B2*pho3*x(3) - B1*B3*pho2*x(3) - B2*B3*pho1*x(3) + B1*B3*pho2*x(4) + B2*B3*pho1*x(4) + A1*pho1*pho2*pho3 + A2*pho1*pho2*pho3 + A3*pho1*pho2*pho3 - B1*B2*B3*x(3)^2 + B1*B2*B3*x(2)*x(3) - B1*B2*B3*x(2)*x(4) + B1*B2*B3*x(3)*x(4) - A1*B2*pho1*pho3*x(2) + A1*B2*pho1*pho3*x(3) - A1*B3*pho1*pho2*x(3) + A1*B3*pho1*pho2*x(4) - A2*B3*pho1*pho2*x(3) + A2*B3*pho1*pho2*x(4) - A1*B2*B3*pho1*x(3)^2 + A1*B2*B3*pho1*x(2)*x(3) - A1*B2*B3*pho1*x(2)*x(4) + A1*B2*B3*pho1*x(3)*x(4))) - (2500*((13*alpha_i3*(T_r3 - x(7)))/10 - mdot_in*(x(3) - x(4)))*(B1*pho2 + B2*pho1 - B1*B2*x(2) + B1*B2*x(3) + A1*pho1*pho2 + A2*pho1*pho2 + A3*pho1*pho2 - A1*B2*pho1*x(2) + A1*B2*pho1*x(3)))/(7*(B1*pho2*pho3 + B2*pho1*pho3 + B3*pho1*pho2 - B1*B2*pho3*x(2) + B1*B2*pho3*x(3) - B1*B3*pho2*x(3) - B2*B3*pho1*x(3) + B1*B3*pho2*x(4) + B2*B3*pho1*x(4) + A1*pho1*pho2*pho3 + A2*pho1*pho2*pho3 + A3*pho1*pho2*pho3 - B1*B2*B3*x(3)^2 + B1*B2*B3*x(2)*x(3) - B1*B2*B3*x(2)*x(4) + B1*B2*B3*x(3)*x(4) - A1*B2*pho1*pho3*x(2) + A1*B2*pho1*pho3*x(3) - A1*B3*pho1*pho2*x(3) + A1*B3*pho1*pho2*x(4) - A2*B3*pho1*pho2*x(3) + A2*B3*pho1*pho2*x(4) - A1*B2*B3*pho1*x(3)^2 + A1*B2*B3*pho1*x(2)*x(3) - A1*B2*B3*pho1*x(2)*x(4) + A1*B2*B3*pho1*x(3)*x(4))) + (2500*((13*alpha_i1*(T_r1 - x(5)))/10 - mdot_in*(h_in - x(2)))*(B1*pho2 - B1*B2*x(2) + B1*B2*x(3) + A3*B1*pho2*x(3) - A3*B1*pho2*x(4) + A3*B1*B2*x(3)^2 - A3*B1*B2*x(2)*x(3) + A3*B1*B2*x(2)*x(4) - A3*B1*B2*x(3)*x(4)))/(7*(B1*pho2*pho3 + B2*pho1*pho3 + B3*pho1*pho2 - B1*B2*pho3*x(2) + B1*B2*pho3*x(3) - B1*B3*pho2*x(3) - B2*B3*pho1*x(3) + B1*B3*pho2*x(4) + B2*B3*pho1*x(4) + A1*pho1*pho2*pho3 + A2*pho1*pho2*pho3 + A3*pho1*pho2*pho3 - B1*B2*B3*x(3)^2 + B1*B2*B3*x(2)*x(3) - B1*B2*B3*x(2)*x(4) + B1*B2*B3*x(3)*x(4) - A1*B2*pho1*pho3*x(2) + A1*B2*pho1*pho3*x(3) - A1*B3*pho1*pho2*x(3) + A1*B3*pho1*pho2*x(4) - A2*B3*pho1*pho2*x(3) + A2*B3*pho1*pho2*x(4) - A1*B2*B3*pho1*x(3)^2 + A1*B2*B3*pho1*x(2)*x(3) - A1*B2*B3*pho1*x(2)*x(4) + A1*B2*B3*pho1*x(3)*x(4)));

dx(5)=(1529*alpha_o1*(T_a3 - x(5)))/118850 + (13*alpha_i1*(T_r1 - x(5)))/23770;

dx(6)=(1529*alpha_o2*(T_a2 - x(6)))/118850 + (13*alpha_i2*(T_r2 - x(6)))/23770;

dx(7)= (1529*alpha_o3*(T_a1 - x(7)))/118850 + (13*alpha_i3*(T_r3 - x(7)))/23770;

dx  =dx(:);

end

其实就是利用前面的 m程序通过矩阵计算求取出dx(1), dx(2)...的值，然后把得出的值带到后面m文件中 dx(1)=..., dx(2)=...中，可以看到特别麻烦，我想问一下可不可以调取ode45直接对矩阵求解，就是不用写出后面dx(1)=..., dx(2)=...中一大堆的式子了，希望大家给出些例子讲解下，不甚感激啊！！！

赞一下

回复此楼

dream

3楼2018-10-06 20:54:46

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

vista朱

木虫 (著名写手)

应助: 39 (小学生)
金币: 5493
散金: 348
红花: 42
帖子: 1891
在线: 331.9小时
虫号: 1930187
注册: 2012-08-08
性别: GG
专业: 动力学与控制

引用回帖:

3楼: Originally posted by vista朱 at 2018-10-06 20:54:46
哦哦不好意思哒
这是代码
%%
%%
clc
clear
syms A1 A2 A3 B1 B2 B3 x2 x3 x4 x5 x6 x7 pho1 pho2 pho3 mdot_in h_in T_a1 T_a2 T_a3 T_r1 T_r2 T_r3 alpha_o1 alpha_o2 alpha_o3 alpha_i1 alpha_i ...

抱歉，最后一行是
dx =dx(

;