| 查看: 634 | 回复: 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 |
回复此楼» 猜你喜欢
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有2人回复
-
物理学I论文润色/翻译怎么收费?
已经有79人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
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
28