| 3 | 1/1 | 返回列表 |
| 查看: 542 | 回復(fù): 2 | |||
[交流]
金幣請(qǐng)教一個(gè)matlab畫三維圖的問題
|
|
我想用matlab畫一個(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)長,不過我驗(yàn)證了,公式輸入沒有問題。不知道為什么會(huì)出錯(cuò),請(qǐng)高手指點(diǎn)。 |
» 搶金幣啦!回帖就可以得到:
+1/192
+1/89
+1/38
+2/18
+1/18
+1/13
+1/10
+1/10
+1/9
+1/8
+1/6
+1/6
+1/6
+1/4
+1/3
+1/2
+1/2
+1/2
+1/1
+1/1
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 352分-085602-一志愿985 +5 | 海納百川Ly 2026-03-29 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 0703化學(xué) +11 | 妮妮ninicgb 2026-03-27 | 11/550 |
|
|
[考研] 0703 化學(xué) 求調(diào)劑,一志愿山東大學(xué) 342 分 +4 | Shern—- 2026-03-28 | 4/200 |
|
|
[考研] 316求調(diào)劑 +7 | 江辭666 2026-03-26 | 7/350 |
|
|
[考研] 復(fù)試調(diào)劑 +3 | raojunqi0129 2026-03-28 | 3/150 |
|
|
[考研] 0703本科鄭州大學(xué)求調(diào)劑 +3 | nhj_ 2026-03-25 | 3/150 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +8 | 邱gl 2026-03-27 | 8/400 |
|
|
[考研] 材料277求調(diào)劑 +7 | min3 2026-03-24 | 7/350 |
|
|
[考研] 藥學(xué)105500求調(diào)劑 +3 | Ssun。。 2026-03-28 | 3/150 |
|
|
[考研] 352分 化工與材料 +5 | 海納百川Ly 2026-03-27 | 5/250 |
|
|
[考研] 305求調(diào)劑 +5 | 哇盧卡庫 2026-03-26 | 5/250 |
|
|
[考研] 314求調(diào)劑 +3 | 溪云珂 2026-03-26 | 3/150 |
|
|
[考研] 359求調(diào)劑 +4 | 王了個(gè)楠 2026-03-25 | 4/200 |
|
|
[考研]
材料調(diào)劑
5+4
|
想要一壺桃花水 2026-03-25 | 10/500 |
|
|
[考研] 334分 一志愿武理 材料求調(diào)劑 +4 | 李李不服輸 2026-03-26 | 4/200 |
|
|
[考研] 332求調(diào)劑 +6 | 032500 2026-03-25 | 6/300 |
|
|
[考研] 網(wǎng)絡(luò)空間安全0839招調(diào)劑 +4 | w320357296 2026-03-25 | 6/300 |
|
|
[考研] 一志愿吉林大學(xué)材料與化工303分求調(diào)劑 +4 | 為學(xué)666 2026-03-24 | 4/200 |
|
|
[考研] 材料專碩找調(diào)劑 +5 | 哈哈哈吼吼吼哈 2026-03-23 | 5/250 |
|
|
[考研] 一志愿國科過程所081700,274求調(diào)劑 +3 | 三水研0水立方 2026-03-23 | 3/150 |
|