| 查看: 574 | 回复: 3 | |||
[求助]
帮忙调试一段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 |
回复此楼» 猜你喜欢
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有2人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有12人回复
-
物理学I论文润色/翻译怎么收费?
已经有104人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有13人回复
-
桂林理工大学物理学专业招收调剂
已经有17人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
-
此贴已删除
已经有7人回复
-
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有0人回复
-
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树
已经有0人回复
-
【方案分享】单环磁场+轴心控制+偏转导出电子束约束系统（可行性实验）
已经有0人回复
2楼2016-03-12 20:14:49
3楼2016-03-12 20:15:01
匿名
用户注销 (著名写手)
奔跑路上的小汪
- 应助: 15 (小学生)
- 金币: 1858.1
- 散金: 346
- 红花: 12
- 沙发: 1
- 帖子: 1199
- 在线: 155.1小时
- 虫号: 0
- 注册: 2016-03-14
- 性别: GG
- 专业: 传热传质学
4楼2017-04-05 17:11:20
