| 4 | 1/1 | 返回列表 |
| 查看: 561 | 回復(fù): 3 | ||
仰望星空7新蟲 (初入文壇)
|
[求助]
幫忙調(diào)試一段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; %燃?xì)鈱α?br /> 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; %內(nèi)壁導(dǎo)熱 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; %外壁導(dǎo)熱 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 得到的結(jié)果應(yīng)該是溫度值,但我卻得到復(fù)數(shù),初值大概為1500 1200 1000 800 500 |
新蟲 (初入文壇)
用戶注銷 (著名寫手)
奔跑路上的小汪
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0854人工智能方向招收調(diào)劑 +3 | 章小魚567 2026-03-24 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 調(diào)劑 +4 | 13853210211 2026-03-24 | 4/200 |
|
|
[考研] 材料292調(diào)劑 +8 | 橘頌思美人 2026-03-23 | 8/400 |
|
|
[考研] 307求調(diào)劑 +5 | 超級伊昂大王 2026-03-24 | 5/250 |
|
|
[考研] 307求調(diào)劑 +3 | 余意卿 2026-03-21 | 6/300 |
|
|
[考研] 材料與化工328分調(diào)劑 +4 | 。,。,。,。i 2026-03-23 | 4/200 |
|
|
[考研]
求調(diào)劑材料學(xué)碩080500,總分289分
5+3
|
@taotao 2026-03-19 | 21/1050 |
|
|
[考研] 070300,一志愿北航320求調(diào)劑 +3 | Jerry0216 2026-03-22 | 5/250 |
|
|
[考研] 323求調(diào)劑 +6 | 洼小桶 2026-03-18 | 6/300 |
|
|
[考研] 293求調(diào)劑 +3 | 濤濤Wjt 2026-03-22 | 5/250 |
|
|
[考研] 求助 +5 | 夢里的無言 2026-03-21 | 6/300 |
|
|
[考研] 313求調(diào)劑 +4 | 肆叁貳壹22 2026-03-19 | 4/200 |
|
|
[考研] 材料 271求調(diào)劑 +5 | 展信悅_ 2026-03-21 | 5/250 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑 +7 | 墨墨漠 2026-03-20 | 7/350 |
|
|
[考研] 085601調(diào)劑 358分 +3 | zzzzggh 2026-03-20 | 4/200 |
|
|
[考研] 295求調(diào)劑 +4 | 一志愿京區(qū)211 2026-03-18 | 6/300 |
|
|
[考研] 329求調(diào)劑 +9 | 想上學(xué)吖吖 2026-03-19 | 9/450 |
|
|
[考研] 求調(diào)劑 +3 | @taotao 2026-03-20 | 3/150 |
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +5 | February_Feb 2026-03-19 | 5/250 |
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|