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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » MATLAB非线性方程组求解~~~反应釜模拟
汕头大学海洋科学接受调剂
51/1返回列表
查看: 1875  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

jiakai003

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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-05-14 11:10:45
已阅   关注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的回帖
查看全部 10 个回答

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的回帖
查看全部 10 个回答
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 药学求调剂 +11 喽哈加油 2026-04-14 13/650 2026-04-14 21:14 by qingfeng258
[考研] 一志愿鲁东大学071000生物学学硕初试分数276求调剂 +26 慕绝cc 2026-04-09 30/1500 2026-04-14 18:50 by 蔡苏阳
[考研] 367求调剂 +11 hffQAQ 2026-04-09 11/550 2026-04-14 17:48 by lhj2009
[考研] 化工学硕294分，求导师收留 +32 yzyzx 2026-04-12 36/1800 2026-04-14 17:45 by lhj2009
[考研] 335求调剂 +19 想上岸呀！！ 2026-04-12 21/1050 2026-04-14 16:23 by Art1977
[考研] 290求调剂 +21 luoziheng 2026-04-10 23/1150 2026-04-14 15:49 by zs92450
[考研] 279学硕食品专业求调剂院校 20+6 孤独的狼爱吃羊 2026-04-12 28/1400 2026-04-14 15:44 by zs92450
[基金申请] RY：中国产出的科学垃圾论文，绝对数量和比例都世界第一 +6 zju2000 2026-04-14 17/850 2026-04-14 14:34 by jurkat.1640
[考研] 0854调剂 +10 长弓傲 2026-04-11 11/550 2026-04-13 10:38 by wp06
[考研] 291求调剂 +11 关忆北. 2026-04-09 12/600 2026-04-12 10:32 by 逆水乘风
[考研] 284求调剂 +11 archer.. 2026-04-09 12/600 2026-04-11 20:23 by 蓝云思雨
[考研] 农学0904 312求调剂 +6 Say Never 2026-04-10 6/300 2026-04-11 10:33 by wwj2530616
[考研] 广东省 085601 329分求调剂 +14 Eddieddd 2026-04-10 14/700 2026-04-11 09:58 by bljnqdcc
[考研] 考研调剂 +26 硕星赴 2026-04-09 27/1350 2026-04-10 22:24 by 猪会飞
[考研] 368求调剂 +3 17385968172 2026-04-10 3/150 2026-04-10 20:12 by 电子奥德彪
[考研] 314求调剂 +18 xhhdjdjsjks 2026-04-09 19/950 2026-04-10 18:53 by HPUCZ
[考博] 博士自荐 +7 可可小胖 2026-04-08 7/350 2026-04-10 08:28 by kimhero
[考研] 348求调剂 +3 candyyyi 2026-04-09 3/150 2026-04-09 17:20 by 段伟艳
[考研] 085501机械英二77总分294求调剂，接受跨专业学习 +6 守法公民亓纪 2026-04-08 6/300 2026-04-09 15:55 by wp06
[考研] 招收有机化学、化工，药学，食品灯专业学生 +3 yrfhjgdj 2026-04-08 3/150 2026-04-09 10:15 by QYQX_123
信息提示
关闭
请填处理意见
关闭