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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » 紧急求助
61/1返回列表
查看: 1050  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

vista朱

木虫 (著名写手)

[求助] 紧急求助 已有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很多用法不太了解，谢谢大家了。
回复此楼

» 猜你喜欢
dream
1楼 2018-10-06 12:08:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

天天进步啊

木虫 (著名写手)
★
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朱

木虫 (著名写手)
★
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朱

木虫 (著名写手)
引用回帖:
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(;
赞一下 回复此楼
dream
4楼2018-10-06 20:56:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chendequan

铁虫 (小有名气)

【答案】应助回帖

★
感谢参与，应助指数 +1
jjdg: 金币+1, 感谢参与 2018-10-09 21:16:04
内容已删除
赞一下 回复此楼
QQ：516477448，真心帮助解决MATLAB相关问题，提供详细资料，Word文档明确具体问题及要求，尽力而为！
5楼2018-10-07 09:11:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

★
感谢参与，应助指数 +1
jjdg: 金币+1, 感谢参与 2018-10-09 21:15:52
构造微分方程是供ode函数调用的先决条件，而且你的方程里面那么多的符号参数，你怎么用ode计算？
赞一下 回复此楼
数值计算
6楼2018-10-08 22:04:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 vista朱 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 资源与环境 调剂申请(333分) +5 holy J 2026-03-21 5/250 2026-03-21 22:42 by Catalysis25
[考研] 269专硕求调剂 +5 金恩贝 2026-03-21 5/250 2026-03-21 22:37 by zhyzzh
[考研] 材料学硕301分求调剂 +7 Liyouyumairs 2026-03-21 7/350 2026-03-21 22:31 by peike
[考研] 考研调剂 +3 呼呼？~+123456 2026-03-21 3/150 2026-03-21 20:04 by 无际的草原
[考研] 求调剂 +4 要好好无聊 2026-03-21 4/200 2026-03-21 18:57 by 学员8dgXkO
[考研] 306求调剂 +4 chuanzhu川烛 2026-03-18 4/200 2026-03-21 08:25 by laoshidan
[考研] 346求调剂[0856] +4 WayneLim327 2026-03-16 7/350 2026-03-21 04:02 by JourneyLucky
[考研] 一志愿天津大学化学工艺专业（081702）315分求调剂 +12 yangfz 2026-03-17 12/600 2026-03-21 03:30 by JourneyLucky
[考研] 材料工程（专）一志愿985 初试335求调剂 +3 hiloiy 2026-03-17 4/200 2026-03-21 03:04 by JourneyLucky
[考研] 287求调剂 +7 晨昏线与星海 2026-03-19 8/400 2026-03-20 22:19 by JourneyLucky
[考研] 08工学调剂 +5 用户573181 2026-03-20 5/250 2026-03-20 15:47 by xia_2003
[考研] 一志愿中国海洋大学，生物学，301分，求调剂 +5 1孙悟空 2026-03-17 6/300 2026-03-19 23:46 by zcl123
[考研] 一志愿福大288有机化学，求调剂 +3 小木虫200408204 2026-03-18 3/150 2026-03-19 13:31 by houyaoxu
[考研] 本科郑州大学物理学院，一志愿华科070200学硕，346求调剂 +4 我不是一根葱 2026-03-18 4/200 2026-03-19 09:11 by 浮云166
[考研] 生物学071000 329分求调剂 +3 我爱生物生物爱� 2026-03-17 3/150 2026-03-18 10:12 by macy2011
[考研] 0703化学336分求调剂 +6 zbzihdhd 2026-03-15 7/350 2026-03-18 09:53 by zhukairuo
[考研] 材料专硕326求调剂 +6 墨煜姒莘 2026-03-15 7/350 2026-03-17 17:10 by ruiyingmiao
[考研] 290求调剂 +3 p asserby. 2026-03-15 4/200 2026-03-17 16:35 by wangkm
[考研] 有没有道铁/土木的想调剂南林，给自己招师弟中～ +3 TqlXswl 2026-03-16 7/350 2026-03-17 15:23 by TqlXswl
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
信息提示
关闭
请填处理意见
关闭