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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » MATLAB非线性方程组求解~~~反应釜模拟
101/1返回列表
查看: 2089  |  回复: 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返回列表
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 今天系统多次维护，明天很可能放榜！ +5 zju2000 2026-08-16 5/250 2026-08-16 23:23 by Haru815
[基金申请] 时间戳又变了8-15 +13 archvillain 2026-08-15 25/1250 2026-08-16 20:13 by zhaosm1982
[基金申请] 快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。 +3 Tide man 2026-08-14 3/150 2026-08-16 17:47 by jurkat.1640
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +9 Ldrop2023 2026-08-13 9/450 2026-08-16 13:24 by wukun_007
[基金申请] 我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人 +11 家与远方 2026-08-10 16/800 2026-08-16 10:28 by ray43
[基金申请] 有时候，自然基金真的不能太认真 （我的申报经验） +10 majunge000 2026-08-11 12/600 2026-08-16 08:18 by xli1984
[精细化工] 招聘 金属平磨液，抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] filecode +8 cratir 2026-08-14 12/600 2026-08-15 18:08 by zyfgau
[基金申请] 小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊 +12 Tide man 2026-08-10 13/650 2026-08-15 16:34 by 氺木
[基金申请] 各位道友，我要去昆明玩几天，回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果？ +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] FileCode能看出啥？ +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 重要来源：本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能，周易传统文化，filecode打分制来了，3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见，大家看看几个A，几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 什么时候出结果，有咨询渠道？？？ +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 确定了，国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
信息提示
关闭
请填处理意见
关闭