 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 1 | 1/1 | 返回列表 |
| 查看: 293 | 回復(fù): 0 | ||
yangjian911x鐵蟲(chóng) (初入文壇)
|
[求助]
求大神幫忙解決一個(gè)matlab算例問(wèn)題
|
|
下面是我寫(xiě)的一串代碼： k0=0; for h=0.01:0.01:50; %////////cengshu p=101325;%wendu t=300; a=100; N=50; l=100; g=h/a; px=124.5;%xianweimidu pa=110;%qiningjiaomidu p1 = 16.8367680681458; p2 = 351.719716799583; p3 = 170.039755829837; p4 = 23.1916198325517; p5 = 1269.41233622451; p6 = -4.76831022571228; p7 = -124.095938410379; p8 =0.158164420679919; p9 = 3.48022582551722; kgo=1.6241*10^(-3)+8.4798*10^(-5)*t+2.8587*10^(-9)*t^2-3.7581*10^(-11)*t^3+1.6705*10^(-14)*t^4; kg=kgo/(1+sqrt(2)*1.75*1.38*10^(-23)*t/(pi*0.366^(2)*10^(-18)*p*l*10^(-9))); % a2=0.8; s=(324.3/pa+5.03)*10^5;%比表面積 kng=60.22*10^(5)*p*t^(-0.5)/(0.25*s*(1-pa/2200)^(-1)*pa+4.01*10^(9)*p*t^(-1)); d=12*(1-(1-pa/2200))/((2+a2^2)*pa*s); D=(-pi*sqrt(1-a2^2)*d^2/(pa*s)+((pi*sqrt(1-a2^2)*d^2/(pa*s))^2+(-pi*d/pa/s)^3)^0.5)^(1/3)+(-pi*sqrt(1-a2^2)*d^2/(pa*s)-((pi*sqrt(1-a2^2)*d^2/(pa*s))^2+(-pi*d/pa/s)^3)^0.5)^(1/3); a1=d/D; C1=5.2*10^(-6); kSO=0.75264+0.0031286*t-0.0000045242*t^2+3.5253*10^(-9)*t^3; kS0=0.75264+0.0031286*300-0.0000045242*300^2+3.5253*10^(-9)*300^3; kns=C1*pa^1.5*kSO/kS0;%氣凝膠固體熱導(dǎo)率 b=1-kng/kns; n=50; ka=(pi*a1^2*a2^2/(4*(1-b))+(1-a1^2)-pi*a1^2*(1-a2^2)*(b+log(1-b))/(2*b^2)+pi*(1/sqrt(1-a2^2)-a1)*(log((1-b*a1*a2)/(1-b*a1))/(b*a1)-(1-a2))/b)*kng; %//氣凝膠熱導(dǎo)率 l1=0.508742684953864; l2=0.0783770550695127; l3=-0.004081649272934; l4=0.000177301411935905; l5=2.59695488944144*10^(-7); l6=-3.97089287852354*10^(-7); l7=1.22847546022621*10^(-8); l8=-7.92885898906499*10^(-11); q1 = 596118.307345212; q2 = -113950.80374852; q3 = 9006.52284616573; q4 = -372.348258004001; q5 = 8.38526805820327; q6 = -0.0955903338053084; q7 = 0.000414336358839436;%氣凝膠 ks=l1+l2*t^0.5+l3*t+l4*t^1.5+l5*t^2+l6*t^2.5+l7*t^3+l8*t^3.5; kc=0.8*ka+0.2*ks; kcg=2*kc*kg/(kc+kg); %復(fù)合材料單元體密度 KR1=(p1+p3*t^0.5+p5*t+p7*t^1.5+p9*t^2)/(1+p2*t^0.5+p4*t+p6*t^1.5+p8*t^2); KR2=(q1+q2*t^0.5+q3*t+q4*t^1.5+q5*t^2+q6*t^2.5+q7*t^3)/pa; KR=px*KR1+0.8*0.95*pa*KR2; kr=(16*5.67*10^(-8)*t^3)/(3*KR);%fushe ke=kg+4*g*(sqrt(3)-g)*(kc-kg)/3;%danceng if k0<=0; for i=1:N-1; A1=rand(1)*((1.5*sqrt(3)*a^2-6*a*h+2*sqrt(3)*h^2)-(1.5*sqrt(3)*a^2-12*a*h+8*sqrt(3)*h^2))+1.5*sqrt(3)*a^2-12*a*h+8*sqrt(3)*h^2; A2=rand(1)*((6*a*h-2*sqrt(3)*h^2)-(4*sqrt(3)*h^2))+4*sqrt(3)*h^2; A3=1.5*sqrt(3)*a^2-(A1+A2); K=(A1*kg+A2*kc+A3*kcg)/(1.5*sqrt(3)*a^2); k0=k0+1/K; hold on; end; end;%/////// k=N/(k0+1/ke)+kr; hold on; plot(h,k,'b-.*'); end 想求各位大神，能不能保證代碼中（紅色部分）隨機(jī)取樣對(duì)應(yīng)于不同的h只運(yùn)行一次。 [ Last edited by jjdg on 2012-11-17 at 23:55 ] |
回復(fù)此樓» 猜你喜歡
08工科 320總分 求調(diào)劑
已經(jīng)有9人回復(fù)
-
生物學(xué)071000 329分求調(diào)劑
已經(jīng)有5人回復(fù)
-
306求調(diào)劑
已經(jīng)有5人回復(fù)
-
284求調(diào)劑
已經(jīng)有4人回復(fù)
-
尋找調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿中南化學(xué)（0703）總分337求調(diào)劑
已經(jīng)有10人回復(fù)
-
324求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿 西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。
已經(jīng)有3人回復(fù)
-
289求調(diào)劑
已經(jīng)有7人回復(fù)
-
319求調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 yangjian911x 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 生物學(xué)071000 329分求調(diào)劑 +5 | 我愛(ài)生物生物愛(ài)?/a> 2026-03-17 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 尋找調(diào)劑 +4 | 倔強(qiáng)芒? 2026-03-21 | 4/200 |
|
|
[考研] 319求調(diào)劑 +4 | 小力氣珂珂 2026-03-20 | 4/200 |
|
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-22 | 6/300 |
|
|
|
[基金申請(qǐng)] 山東省面上項(xiàng)目限額評(píng)審 +4 | 石瑞0426 2026-03-19 | 4/200 |
|
|
|
[考研] 085600材料與化工306 +4 | z1z2z3879 2026-03-21 | 4/200 |
|
|
|
[考研] 廣西大學(xué)材料導(dǎo)師推薦 +3 | 夏夏夏小正 2026-03-17 | 5/250 |
|
|
|
[考研] 材料工程專(zhuān)碩 348分求調(diào)劑 +3 | 冬辭. 2026-03-17 | 5/250 |
|
|
|
[考研] 一志愿深大，0703化學(xué)，總分302，求調(diào)劑 +4 | 七月-七七 2026-03-21 | 4/200 |
|
|
|
[考研] 313求調(diào)劑 +4 | 肆叁貳壹22 2026-03-19 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +3 | 白QF 2026-03-21 | 3/150 |
|
|
|
[考研] 22408 344分 求調(diào)劑 一志愿 華電計(jì)算機(jī)技術(shù) +4 | solanXXX 2026-03-20 | 4/200 |
|
|
|
[考研] 南京大學(xué)化學(xué)376求調(diào)劑 +3 | hisfailed 2026-03-19 | 6/300 |
|
|
|
[考研] 288求調(diào)劑 +16 | 于海海海海 2026-03-19 | 16/800 |
|
|
|
[考研] 一志愿 西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +4 | 晨昏線與星海 2026-03-19 | 4/200 |
|
|
|
[考研] 085600材料與化工 +8 | 安全上岸！ 2026-03-16 | 8/400 |
|
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 5/250 |
|
|
|
[考研] 295材料求調(diào)劑，一志愿武漢理工085601專(zhuān)碩 +5 | Charlieyq 2026-03-19 | 5/250 |
|
|
|
[考研] 工科材料085601 279求調(diào)劑 +7 | 困于星晨 2026-03-17 | 9/450 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 18889395102 2026-03-18 | 4/200 |
|