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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » MATLAB非线性方程组求解~~~反应釜模拟
汕头大学海洋科学接受调剂
101/1返回列表
查看: 1871  |  回复: 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返回列表
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 山东省基金2026 +8 jerry681 2026-04-08 11/550 2026-04-15 10:33 by dx设计师
[考研] 290调剂生物0860 +35 哇哈哈，。 2026-04-11 41/2050 2026-04-14 23:47 by Xurambo2014
[考研] 材料299专硕求调剂 +16 +21 2026-04-09 16/800 2026-04-14 17:06 by jiangguiquan11
[考研] 297，工科调剂? +6 河南农业大学-能 2026-04-14 6/300 2026-04-14 16:04 by Art1977
[基金申请] RY：中国产出的科学垃圾论文，绝对数量和比例都世界第一 +6 zju2000 2026-04-14 17/850 2026-04-14 14:34 by jurkat.1640
[考研] 085600材料与化工349分求调剂 +16 李木子啊哈哈 2026-04-12 17/850 2026-04-14 09:11 by fenglj492
[考研] 一志愿085802 323分求调剂 +13 drizzle_9 2026-04-12 14/700 2026-04-13 10:26 by Faiz5552
[考研] 电气专硕320求调剂 +6 小麻子111 2026-04-10 6/300 2026-04-12 10:54 by lemon6009
[考研] 085410 273求调剂 +10 X1999 2026-04-09 10/500 2026-04-12 09:24 by 逆水乘风
[考研] 280求调剂 +7 兮兮夜夜 2026-04-09 10/500 2026-04-12 00:33 by 蓝云思雨
[考研] 270求调剂 +14 杨乐369 2026-04-11 14/700 2026-04-11 20:16 by 蓝云思雨
[考研] 283求调剂 +22 那个噜子 2026-04-09 22/1100 2026-04-11 10:41 by 逆水乘风
[考研] 293求调剂 +6 勇远库爱314 2026-04-08 6/300 2026-04-11 10:08 by zhq0425
[考研] 337求调剂 +4 研s. 2026-04-10 4/200 2026-04-11 08:57 by zhq0425
[考研] 284求调剂 +9 让我上岸吧阿西 2026-04-09 11/550 2026-04-10 19:18 by 靖jing
[考研] 0858求调剂 5+5 Gky09300550， 2026-04-10 8/400 2026-04-10 19:13 by chemisry
[考研] 调剂申请086000一志愿西北农林科技大学生物与医药320分-本科齐鲁工业大学 +3 美美女士 2026-04-09 3/150 2026-04-10 10:31 by liuhuiying09
[考研] 调剂 +19 2261744733 2026-04-08 19/950 2026-04-09 19:11 by vgtyfty
[考研] 332，085601求调剂 +12 ydfyh 2026-04-09 14/700 2026-04-09 17:28 by wp06
[考研] 复试调剂，一志愿郑州大学材料与化工289分 +31 硕星赴 2026-04-08 31/1550 2026-04-09 16:54 by Delta2012
信息提示
关闭
请填处理意见
关闭