24小时热门版块排行榜

>论坛更新日志 (3160)
>虫友互识 (669)
>休闲灌水 (108)
>导师招生 (98)
>博后之家 (92)
>基金申请 (92)
>论文投稿 (85)
>文献求助 (76)
>教师之家 (46)
>硕博家园 (41)
>考博 (25)
>药学 (20)
>招聘信息布告栏 (14)
>医学 (11)
>公派出国 (11)
>考研 (10)

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

仰望星空7

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 1237.7
帖子: 42
在线: 21.3小时
虫号: 4127268
注册: 2015-10-08
专业: 工程热物理相关交叉领域

[求助] 帮忙调试一段MATLAB程序已有1人参与

程序如下
function F=myfun(X)
T1=X(1);T2=X(2);T3=X(3);T4=X(4);Tf=X(5);q=X(6);
c1=0.003;
c2=0.003;
h=0.002;
w1=0.002;
w2=0.002;
v=5;
T0=3287;
%燃气对流
syms sigma At A Pr hc1
dt=0.013;d=0.03;Pc=3.5*10^6;c=2000;C0=9.18*10^5;Cp=1500;u=1.56*10^(-5);
m=0.022;k=1.26;Ma=3;
sigma=(T1*(1+(k-1)*Ma*Ma/2)/(2*T0)+0.5)^(-0.68)*(1+(k-1)*Ma*Ma/2)^(-0.12);
At=pi*dt*dt/4;A=pi*d*d/4;
Pr=u*Cp/m;
hc1=(0.026/dt^0.2)*(u^0.2*Cp/Pr^0.6)*(Pc/c)^0.8*(At/A)^0.9*sigma;
F(1)=hc1*(T0-T1)+C0-q;
%内壁导热
syms k2
m1=387;
k2=m1/c1;
F(2)=k2*(T1-T2)-q;
%冷却剂对流
syms Cpl ml lf lt pf pt uf ut Ref Prf hc20 hc2 yita
Cpl=4.1*Tf+1892.2;
ml=-0.0001*Tf+0.1383;lf=10^(-8)*(-5.3*Tf+683.7);lt=10^(-8)*(-5.3*T2+683.7);
pf=-0.8102*Tf+831.4661;pt=-0.8102*T2+831.4661;uf=pf*lf;ut=pt*lt;
de=2*h*w1/(h+w1);Ref=v*w1/lf;Prf=uf*Cpl/ml;
hc20=0.023*Ref^0.8*Prf^0.34*(uf/ut)*ml/de;
m2=202.4;m=(2*hc20/(m2*w2))^0.5;yita=tanh(m*h)/(m*h);
hc2=hc20*yita;
F(3)=hc2*(T2-Tf)-q;
F(4)=T2+T3-2*Tf;
%外壁导热
m3=387;
F(5)=(T3-T4)/(c2/m3)-q;
%空气对流
syms av Gra hc3
g=9.8;Ta=293;av=1/((T4+Ta)/2);d=0.2;pa=1.06;Pra=0.7;ma=0.04;la=4*10^(-5);
Gra=g*av*(T4-Ta)*d^3/la^2;hc3=0.48*(Gra*Pra)^0.25*ma/d;
F(6)=hc3*(T4-Ta)-q;
end

得到的结果应该是温度值，但我却得到复数，初值大概为1500 1200 1000 800 500

回复此楼