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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » MATLAB非线性方程组求解~~~反应釜模拟
101/1返回列表
查看: 2134  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

jiakai003

禁虫 (小有名气)
本帖内容被屏蔽

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-05-14 11:10:45
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

change0618

铁杆木虫 (著名写手)

方丈大师

【答案】应助回帖

★
jiakai003(金币+3): 我还没输入方程组呢，数据输入完酒错误了，貌似是调用函数错误~ 2011-05-14 12:17:26
sunyang1988(金币+1): 谢谢交流 2011-05-14 12:33:34
建立一个完整的m函数，而不是在命令窗口一行一行的输入
赞一下(1人) 回复此楼
2楼2011-05-14 11:45:10
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

change0618

铁杆木虫 (著名写手)

方丈大师

【答案】应助回帖

★
sunyang1988(金币+1): 谢谢交流 2011-05-14 12:33:42
你还是先看看书再搞。不然真的是说什么你都不明白。不能一窍不通就想春暖花开
赞一下(1人) 回复此楼
3楼2011-05-14 12:21:47
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

xianghongw

铜虫 (小有名气)

【答案】应助回帖

★ ★ ★
sunyang1988(金币+3): 谢谢人热心帮助 2011-05-14 12:33:59
jiakai003(金币+6): 貌似我一开始时直接在命令窗口输入的，晕死，我没用过，谢谢了，我正在输方程，谢谢!!!只是丢死人了~~~ 2011-05-14 13:55:33
首先建立函数文件，然后输入：
global  Cin1 Cin2 Cin3 Cin4 Cin5 k1 k2 k3 k4 k5 k6 d1 d2 d3 ...
       d4 m Macoh b2 b4 b6 Cinh2o Cgh2o Clh2o t tol tog tl ...
       cpacoh cph2o cppx cpair Cinacoh Cgacoh Clacoh Fing ...
       Hvacoh Hvh2o Cinpx Hr C10 C20 C30 C40 C50 C60 Finl0 Fout0 L G
%parameter values
Cin1=2.85; Cin2=0; Cin3=0; Cin4=0; Cin5=0; k1=5.50e5; k2=15.88e3; ...
k3=17.09e3; k4=3.28e3; k5=9.81e3; k6=0.54; d1=1.4247; d2=0; d3=0; ...
d4=4.8419; m=0.0146; Macoh=182.1; b2=0.5254; b4=0.8111; b6=0.9302;...
Cinh2o=8.33; Cgh2o=3.124; Clh2o=1.98; t=469; tog=433; tol=412; tl=363;...
cpacoh=140; cph2o=4200; cppx=214; cpair=50.9; Cinacoh=0.0167; Cgacoh=0.0167;...
Clacoh=0.0167; Fing=2.12; Hvacoh=19.9; Hvh2o=35.3; Cinpx=2.85; Hr=1332.2;...
C10=0.025; C20=0.025; C30=0.025; C40=0.005; C50=0.005; C60=0.005; ...
Finl0=1.05e5; Fout0=1.00e5; L=3000; G=1000;
Finli=Finl0:1000:Finl0+5000;%initial conditions
n=length(Finli);
for i=1:n
sol=fsolve('方程',[C10,C20,C30,C40,C50,C60,Finli(i),Fout0]); y{i}=sol;
end
你把方程输进去
赞一下(1人) 回复此楼
4楼2011-05-14 12:30:13
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

jiakai003

禁虫 (小有名气)
本帖内容被屏蔽
5楼2011-05-14 14:36:47
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

xianghongw

铜虫 (小有名气)
我不知道你要干什么，数据输入的恐怖，然后方程没有，没法给你解答啊!能仔细描述一下不？
赞一下 回复此楼
6楼2011-05-14 18:58:32
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

jiakai003

禁虫 (小有名气)
本帖内容被屏蔽
7楼2011-05-14 19:09:35
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

fl1fl

金虫 (小有名气)
这个东西不是一个星期可以搞定的，建议你先看看徐老师的第六章的内容，里面附有相关的程序，改一下方程就可以了，
赞一下 回复此楼
衣带渐宽终不悔，为伊消得人憔悴！
8楼2011-05-14 20:08:20
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

jiakai003

禁虫 (小有名气)
本帖内容被屏蔽
9楼2011-05-14 20:16:00
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

change0618

铁杆木虫 (著名写手)

方丈大师

【答案】应助回帖

jiakai003(金币+51): 真NB,其实我对程序一窍不通，非常感谢你，我以后也要努力学习下这方面的知识，略表心意，希望笑纳，以后还望多多指教！！！ 2011-05-14 20:27:00
function liweimin
clear
clc
format long
global  Cin1 Cin2 Cin3 Cin4 Cin5 k1 k2 k3 k4 k5 k6 d1 d2 d3 d4 m Macoh b2 b4 b6 Cinh2o Cgh2o Clh2o t
global tol tog tl cpacoh cph2o cppx cpair Cinacoh Cgacoh Clacoh Fing Hvacoh Hvh2o Cinpx Hr Finl0 Fout0 L G a b c d
%parameter values
Cin1=2.85; Cin2=0; Cin3=0; Cin4=0; Cin5=0; k1=5.50e5; k2=15.88e3; k3=17.09e3; k4=3.28e3; k5=9.81e3; k6=0.54;
d1=1.4247; d2=0; d3=0; d4=4.8419; m=0.0146; Macoh=182.1; b2=0.5254; b4=0.8111; b6=0.9302;
Cinh2o=8.33; Cgh2o=3.124; Clh2o=1.98; t=469; tog=433; tol=412; tl=363; cpacoh=140; cph2o=4200; cppx=214; cpair=50.9;
Cinacoh=0.0167; Cgacoh=0.0167; Clacoh=0.0167; Fing=2.12; Hvacoh=19.9; Hvh2o=35.3; Cinpx=2.85; Hr=1332.2;
C10=0.025; C20=0.025; C30=0.025; C40=0.005; C50=0.005; C60=0.005; Finl0=1.05e5; Fout0=1.00e5; L=3000; G=1000; a=2; b=0.5; c=1.5; d=2;
Finli=Finl0:1000:Finl0+5000;    %initial conditions
n=length(Finli);
for i=1:n
sol=fsolve(@Equations,[C10,C20,C30,C40,C50,C60,Finli(i),Fout0],optimset('Display','off'));
y(:,i)=real(sol);
end
y
%--------------------------------------------------------------------------
function f= Equations(y)
global Cin1 Cin2 Cin3 Cin4 Cin5 k1 k2 k3 k4 k5 k6 d1 d2 d3 d4 m Macoh b2 b4 b6 Cinh2o Cgh2o Clh2o t
global tol tog tl cpacoh cph2o cppx cpair Cinacoh Cgacoh Clacoh Fing Hvacoh Hvh2o Cinpx Hr L G a b c d
C1=y(1); C2=y(2); C3=y(3); C4=y(4); C5=y(5); C6=y(6); Finl=y(7); Fout=y(8);
f1=Fout*C1+(k1*C1/(d1*C1+d2*C2+d3*C3+d4*C4+m)+k5*C5/(d1*C1+d2*C2+d3*C3+d4*C4+m))*Macoh-Finl*Cin1;
f2=Fout*C2+(k2*C2/(d1*C1+d2*C2+d3*C3+d4*C4+m)^b2-k1*C1/(d1*C1+d2*C2+d3*C3+d4*C4+m)-k6*C6/(d1*C1+d2*C2+d3*C3+d4*C4+m)^b6)*Macoh-Finl*Cin2;
f3=Fout*C3+(k3*C3/(d1*C1+d2*C2+d3*C3+d4*C4+m)-k2*C2/(d1*C1+d2*C2+d3*C3+d4*C4+m)^b2)*Macoh-Finl*Cin3;
f4=Fout*C4+(k4*C4/(d1*C1+d2*C2+d3*C3+d4*C4+m)^b4-k3*C3/(d1*C1+d2*C2+d3*C3+d4*C4+m))*Macoh-Finl*Cin4;
f5=Fout*C5+(k6*C6/(d1*C1+d2*C2+d3*C3+d4*C4+m)^b6-k5*C5/(d1*C1+d2*C2+d3*C3+d4*C4+m))*Macoh-Finl*Cin5;
f6=Fout+G-Finl-L;
f7=Fout*Cinh2o+G*Cgh2o-Finl*Cinh2o-a*(Finl*C1-Fout*C1)+b*Fout*C4+Fout*C3+c*Fout*C2+d*Fout*C6-L*Clh2o;
f8=Finl*(t-tol)*(Cinacoh*cpacoh+Cinh2o*cph2o+Cinpx*cppx)+L*(t-tl)*(Clacoh*cpacoh+Clh2o*cph2o)+Fing*cpair*...
    (t-tog)+G*(Cgacoh*Hvacoh+Cgh2o*Hvh2o)-(Finl*Cinpx-Fout*(Cinh2o+Cinacoh+Cinpx))*Hr;
f=[f1;f2;f3;f4;f5;f6;f7;f8];





y =

  1.0e+003 *

  Columns 1 through 4

   0.000000006802276   0.000000006802276   0.000000006802276   0.000000006802276
  -0.000000135512226  -0.000000135512226  -0.000000135512226  -0.000000135512226
  -0.000000016422428  -0.000000016422428  -0.000000016422428  -0.000000016422428
  -0.000000192485189  -0.000000192485189  -0.000000192485189  -0.000000192485189
  -0.000000409983711  -0.000000409983711  -0.000000409983711  -0.000000409983711
  -0.010049532361599  -0.010049532361599  -0.010049532361599  -0.010049532361599
  -1.311221203993331  -1.311221203993331  -1.311221203993331  -1.311221203993331
   0.688778796006669   0.688778796006670   0.688778796006670   0.688778796006669

  Columns 5 through 6

   0.000000006802276   0.000000006802276
  -0.000000135512226  -0.000000135512226
  -0.000000016422428  -0.000000016422428
  -0.000000192485189  -0.000000192485189
  -0.000000409983711  -0.000000409983711
  -0.010049532361599  -0.010049532361599
  -1.311221203993331  -1.311221203993331
   0.688778796006670   0.688778796006669
赞一下(1人) 回复此楼
10楼2011-05-14 20:19:21
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jiakai003 的主题更新
101/1返回列表
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 大佬帮着分析今年的评审意见？该课题连续申请了三年了，明年要不要换个方向。 100+3 cai-11138 2026-08-31 11/550 2026-09-06 20:49 by archvillain
[考博] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 19:21 by eObJPa4gWmp6
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 jPDp0sc2B7zQ 2026-09-05 4/200 2026-09-06 18:09 by eObJPa4gWmp6
[博后之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 18:09 by eObJPa4gWmp6
[论文投稿] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 18:01 by eObJPa4gWmp6
[考研] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-04 8/400 2026-09-06 17:09 by eObJPa4gWmp6
[教师之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 13:09 by eObJPa4gWmp6
[考博] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 12:49 by eObJPa4gWmp6
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 12:49 by eObJPa4gWmp6
[考博] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 12:29 by eObJPa4gWmp6
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 7/350 2026-09-06 11:20 by eObJPa4gWmp6
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 8/400 2026-09-06 11:09 by eObJPa4gWmp6
[论文投稿] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 jPDp0sc2B7zQ 2026-09-04 8/400 2026-09-06 11:00 by eObJPa4gWmp6
[考研] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 4/200 2026-09-06 07:20 by E1iiBLMU2XnD
[教师之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +8 CRJ0Aq2FtVH0 2026-09-03 15/750 2026-09-06 05:48 by E1iiBLMU2XnD
[基金申请] 面上没中，邀请各位路过的虫友分析一下分数 +12 jackleilei 2026-09-03 13/650 2026-09-05 18:10 by fireguard
[文学芳草园] 初秋的晨风 +4 阿美_Lml888 2026-09-04 6/300 2026-09-05 13:51 by 阿美_Lml888
[硕博家园] Ei源刊怎么投 +4 Fengshun9711 2026-09-04 4/200 2026-09-04 19:30 by 研途知予
[硕博家园] 哈尔滨工业大学韩晓军教授课题组招收2027年硕士推免生及博士研究生 +3 灯灯灯灯DDDD 2026-09-03 3/150 2026-09-03 21:22 by 科研皮皮猪
[基金申请] 要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够 +9 瞬息宇宙 2026-08-31 17/850 2026-09-02 14:04 by anjeeshine
信息提示
关闭
请填处理意见
关闭