| 4 | 1/1 | 返回列表 |
| 查看: 562 | 回復(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; %燃氣對流 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ā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿中南大學(xué)化學(xué)學(xué)碩0703總分337求調(diào)劑 +7 | niko- 2026-03-22 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 311求調(diào)劑 +4 | 勇敢的小吳 2026-03-20 | 4/200 |
|
|
[考研] 【2026考研調(diào)劑】制藥工程 284分 求相關(guān)專業(yè)調(diào)劑名額 +4 | 袁奐奐 2026-03-25 | 8/400 |
|
|
[考研] 318求調(diào)劑 +5 | plum李子 2026-03-21 | 8/400 |
|
|
[考研] 材料學(xué)求調(diào)劑 +6 | Stella_Yao 2026-03-20 | 6/300 |
|
|
[考研] 0703化學(xué)調(diào)劑,求導(dǎo)師收 +7 | 天天好運來上岸?/a> 2026-03-24 | 7/350 |
|
|
[考研] 306求0703調(diào)劑一志愿華中師范 +10 | 紙魚ly 2026-03-21 | 11/550 |
|
|
[考研] 291求調(diào)劑 +3 | HanBeiNingZC 2026-03-24 | 3/150 |
|
|
[考研] 一志愿北京化工大學(xué)材料與化工 264分各科過A區(qū)國家線 +3 | 哈哈157349 2026-03-21 | 3/150 |
|
|
[考研] 279分求調(diào)劑 一志愿211 +18 | chaojifeixia 2026-03-19 | 20/1000 |
|
|
[考研] 335求調(diào)劑 +4 | yuyu宇 2026-03-23 | 5/250 |
|
|
[考研] 328求調(diào)劑 +4 | LHHL66 2026-03-23 | 4/200 |
|
|
[考研] 275求調(diào)劑 +6 | shansx 2026-03-22 | 8/400 |
|
|
[考研] 考研調(diào)劑 +3 | 呼呼?~+123456 2026-03-21 | 3/150 |
|
|
[考研] 求調(diào)劑 +4 | 要好好無聊 2026-03-21 | 4/200 |
|
|
[考研] 求助 +5 | 夢里的無言 2026-03-21 | 6/300 |
|
|
[考研] 22 350 本科985求調(diào)劑,求老登收留 +3 | 李軼男003 2026-03-20 | 3/150 |
|
|
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑,總分315(英一) +5 | sbdksD 2026-03-19 | 5/250 |
|
|
[考研] 353求調(diào)劑 +3 | 拉鉤不許變 2026-03-20 | 3/150 |
|
|
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3 | 葵梓衛(wèi)隊 2026-03-20 | 3/150 |
|