|
|
【答案】应助回帖
★ ★ ★
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
你把方程输进去 |
|
建立一个完整的m函数，而不是在命令窗口一行一行的输入
回复此楼