| 3 | 1/1 | 返回列表 |
| 查看: 541 | 回復(fù): 2 | |||
langzitianya銀蟲(chóng) (正式寫(xiě)手)
|
[交流]
金幣請(qǐng)教一個(gè)matlab畫(huà)三維圖的問(wèn)題
|
|
我想用matlab畫(huà)一個(gè)三維曲面圖。操作如下: >> clear; >> tt=[1700:1:2200]; >> pp=[0:5:25000]; >> [t,p]=meshgrid(tt,pp); >> c=exp(0.53*((7783/t+5.779+0.5*log(p/101325))-((-632.2666*0.53*(1+0.04*(3.65-4.2)*(1-0.53))*(1-0.53)*(1+0.07*(4.2-3.65)*0.53)/(4.8*0.53*(1+0.04*(1-0.53)*(3.65-4.2))+(1-0.53)*4.6*(1+0.07*0.53*(4.2-3.65))))/(2*8.314*t)-(8.314/4-6.9504*10^(-6)*8.314^2*t)*(1-4*(-632.2666*0.53*(1+0.04*(3.65-4.2)*(1-0.53))*(1-0.53)*(1+0.07*(4.2-3.65)*0.53)/(4.8*0.53*(1+0.04*(1-0.53)*(3.65-4.2))+(1-0.53)*4.6*(1+0.07*0.53*(4.2-3.65))))/(8.314*t)).^(1/2)-(3/2)*(-0.28007)*(1-0.53).^2+(8.314/4-6.9504*10^(-6)*8.314^2*t)/8.314+(1-0.53)*(0.5+2*(8.314/4-6.9504*10^(-6)*8.314^2*t)*(1-4*(-632.2666*0.53*(1+0.04*(3.65-4.2)*(1-0.53))*(1-0.53)*(1+0.07*(4.2-3.65)*0.53)/(4.8*0.53*(1+0.04*(1-0.53)*(3.65-4.2))+(1-0.53)*4.6*(1+0.07*0.53*(4.2-3.65))))/(8.314*t)).^(-1/2)/8.314)*(1/(8.314*t))*((-632.2666*0.53*(1+0.04*(3.65-4.2)*(1-0.53))*(1-0.53)*(1+0.07*(4.2-3.65)*0.53)/(4.8*0.53*(1+0.04*(1-0.53)*(3.65-4.2))+(1-0.53)*4.6*(1+0.07*0.53*(4.2-3.65))))*(1/0.53-1/(1-0.53)-0.04*(3.65-4.2)/(1+0.04*(1-0.53)*(3.65-4.2))+0.07*(4.2-3.65)/(1+0.07*0.53*(4.2-3.65))-(4.8*(1+0.04*(1-2*0.53)*(3.65-4.2))+4.6*(-1+0.07*(1-2*0.53)*(4.2-3.65)))/(0.53*4.8*(1+0.04*(1-0.53)*(3.65-4.2))+(1-0.53)*4.6*(1+0.07*0.53*(4.2-3.65)))))))+(1-0.53)*((-6159/t+6.247+0.5*log(p/101325))-(((-632.2666*0.53*(1+0.04*(3.65-4.2)*(1-0.53))*(1-0.53)*(1+0.07*(4.2-3.65)*0.53)/(4.8*0.53*(1+0.04*(1-0.53)*(3.65-4.2))+(1-0.53)*4.6*(1+0.07*0.53*(4.2-3.65))))-t*((-0.4201*8.314*0.53*(1-0.53)*+6.9504*10^(-6)*8.314^2*t*(1-(1-4*(-632.2666*0.53*(1+0.04*(3.65-4.2)*(1-0.53))*(1-0.53)*(1+0.07*(4.2-3.65)*0.53)/(4.8*0.53*(1+0.04*(1-0.53)*(3.65-4.2))+(1-0.53)*4.6*(1+0.07*0.53*(4.2-3.65))))/(8.314*t))^0.5))+((-632.2666*0.53*(1+0.04*(3.65-4.2)*(1-0.53))*(1-0.53)*(1+0.07*(4.2-3.65)*0.53)/(4.8*0.53*(1+0.04*(1-0.53)*(3.65-4.2))+(1-0.53)*4.6*(1+0.07*0.53*(4.2-3.65))))/(2*t)+8.314*(1-4*(-632.2666*0.53*(1+0.04*(3.65-4.2)*(1-0.53))*(1-0.53)*(1+0.07*(4.2-3.65)*0.53)/(4.8*0.53*(1+0.04*(1-0.53)*(3.65-4.2))+(1-0.53)*4.6*(1+0.07*0.53*(4.2-3.65))))/(8.314*t)).^0.5/4-8.314/4))-0.53*8.314*t*((-632.2666*0.53*(1+0.04*(3.65-4.2)*(1-0.53))*(1-0.53)*(1+0.07*(4.2-3.65)*0.53)/(4.8*0.53*(1+0.04*(1-0.53)*(3.65-4.2))+(1-0.53)*4.6*(1+0.07*0.53*(4.2-3.65))))/(2*8.314*t)-(8.314/4-6.9504*10^(-6)*8.314^2*t)*(1-4*(-632.2666*0.53*(1+0.04*(3.65-4.2)*(1-0.53))*(1-0.53)*(1+0.07*(4.2-3.65)*0.53)/(4.8*0.53*(1+0.04*(1-0.53)*(3.65-4.2))+(1-0.53)*4.6*(1+0.07*0.53*(4.2-3.65))))/(8.314*t)).^(1/2)-(3/2)*(-0.28007)*(1-0.53).^2+(8.314/4-6.9504*10^(-6)*8.314^2*t)/8.314+(1-0.53)*(0.5+2*(8.314/4-6.9504*10^(-6)*8.314^2*t)*(1-4*(-632.2666*0.53*(1+0.04*(3.65-4.2)*(1-0.53))*(1-0.53)*(1+0.07*(4.2-3.65)*0.53)/(4.8*0.53*(1+0.04*(1-0.53)*(3.65-4.2))+(1-0.53)*4.6*(1+0.07*0.53*(4.2-3.65))))/(8.314*t)).^(-1/2)/8.314)*(1/(8.314*t))*((-632.2666*0.53*(1+0.04*(3.65-4.2)*(1-0.53))*(1-0.53)*(1+0.07*(4.2-3.65)*0.53)/(4.8*0.53*(1+0.04*(1-0.53)*(3.65-4.2))+(1-0.53)*4.6*(1+0.07*0.53*(4.2-3.65))))*(1/0.53-1/(1-0.53)-0.04*(3.65-4.2)/(1+0.04*(1-0.53)*(3.65-4.2))+0.07*(4.2-3.65)/(1+0.07*0.53*(4.2-3.65))-(4.8*(1+0.04*(1-2*0.53)*(3.65-4.2))+4.6*(-1+0.07*(1-2*0.53)*(4.2-3.65)))/(0.53*4.8*(1+0.04*(1-0.53)*(3.65-4.2))+(1-0.53)*4.6*(1+0.07*0.53*(4.2-3.65)))))))/((1-0.53)*8.314*t)))); 輸入這一步,得到如下結(jié)果: ??? Error using ==> mldivide Matrix dimensions must agree. 公式有點(diǎn)長(zhǎng),不過(guò)我驗(yàn)證了,公式輸入沒(méi)有問(wèn)題。不知道為什么會(huì)出錯(cuò),請(qǐng)高手指點(diǎn)。 |
» 搶金幣啦!回帖就可以得到:
+1/184
+3/128
+1/91
+1/82
+1/74
+1/63
+2/40
+1/40
+1/33
+1/26
+1/18
+1/14
+3/13
+1/10
+1/7
+1/6
+1/6
+1/4
+1/4
+1/3
金蟲(chóng) (正式寫(xiě)手)
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 309求調(diào)劑 +3 | 誰(shuí)不是少年 2026-03-29 | 3/150 |
|
|---|---|---|---|---|
|
[碩博家園] 招收生物學(xué)/細(xì)胞生物學(xué)調(diào)劑 +4 | IceGuo 2026-03-26 | 5/250 |
|
|
[碩博家園] 求調(diào)劑 有機(jī)化學(xué)考研356分 +3 | Nadiums 2026-03-25 | 4/200 |
|
|
[考研] 332求調(diào)劑 +6 | 蕉蕉123 2026-03-28 | 6/300 |
|
|
[考研]
|
nnnnnnn5 2026-03-25 | 7/350 |
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-27 | 6/300 |
|
|
[考研] 299求調(diào)劑 +7 | 嗯嗯嗯嗯2 2026-03-27 | 7/350 |
|
|
[考研] 286求調(diào)劑 +12 | PolarBear11 2026-03-26 | 12/600 |
|
|
[考研] 一志愿南京航空航天大學(xué)材料學(xué)碩求調(diào)劑 +3 | @taotao 2026-03-28 | 3/150 |
|
|
[考研] 材料求調(diào)劑 一志愿哈工大總分298分,前三科223分 +5 | dongfang59 2026-03-27 | 5/250 |
|
|
[考研] 340求調(diào)劑 +5 | jhx777 2026-03-27 | 5/250 |
|
|
[考研] 085601 材料工程 313分 求調(diào)劑 +5 | Ong3 2026-03-27 | 5/250 |
|
|
[考研] 324求調(diào)劑 +5 | hanamiko 2026-03-26 | 5/250 |
|
|
[考研] 317求調(diào)劑 +7 | 蛋黃咸肉粽 2026-03-26 | 7/350 |
|
|
[考研] 求調(diào)劑 +8 | Auroracx 2026-03-22 | 8/400 |
|
|
[考研] 一志愿中南大學(xué)化學(xué)學(xué)碩0703總分337求調(diào)劑 +7 | niko- 2026-03-22 | 7/350 |
|
|
[考研] 347求調(diào)劑 +4 | L when 2026-03-25 | 4/200 |
|
|
[考研] 一志愿武理085500機(jī)械專業(yè)總分300求調(diào)劑 +3 | an10101 2026-03-24 | 7/350 |
|
|
[考研] 341求調(diào)劑(一志愿湖南大學(xué)070300) +5 | 番茄頭--- 2026-03-22 | 6/300 |
|
|
[考研] 361求調(diào)劑 +3 | Glack 2026-03-22 | 3/150 |
|