24小时热门版块排行榜

返回列表

langzitianya

银虫 (正式写手)

应助: 0 (幼儿园)
金币: 1442.9
帖子: 553
在线: 73.5小时
虫号: 226434

[交流] 金币请教一个matlab画三维图的问题

我想用matlab画一个三维曲面图。操作如下：
>> 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))));
输入这一步，得到如下结果：
??? Error using ==> mldivide
Matrix dimensions must agree.

公式有点长，不过我验证了，公式输入没有问题。不知道为什么会出错，请高手指点。

回复此楼