在MATLAB中绘制一个输入和三个输出的子函数

我试图绘制三个不同的图。 h对T，h对P，h对rho。当我运行程序时，我得到三个不同的图表，但没有数据。我不知道我做错了什么。在MATLAB中绘制一个输入和三个输出的子函数

这里是我的阴谋代码

for h=0:1:1000 
[ T,P,rho ] = stdatm(h); 
hold all 
subplot(1,3,1) 
plot(T,h,'r') 

subplot(1,3,2) 
plot(P,h) 

subplot(1,3,3) 
plot(rho,h) 
end

这里是我的功能代码

function [ T,P,rho ] = stdatm(h) 
T0=288.16; 
P0=101.325; 
rho0=1.225; 
a=-6.5*10^-3; 
b=3*10^-3; 
c=-4.5*10^-3; 
d=4.0*10^-3; 
R=286.9; 
g=9.81; 

T1=T0+a*11000; 
P1=P0*(T1/T0)^(-g/(a*R)); 
rho1=rho0*(T1/T0)^((-g/(a*R))-1); 

T2=T1; 
P2=P1*exp((-g/(R*T2))*(25000-11000)); 
rho2=rho1*exp((-g/(R*T2))*(25000-11000)); 

T3=T2+b*(47000-25000); 
P3=P2*(T3/T2)^(-g/(b*R)); 
rho3=rho2*(T3/T2)^((-g/(b*R))-1); 

T4=T3; 
P4=P3*exp((-g/(R*T4))*(53000-47000)); 
rho4=rho3*exp((-g/(R*T4))*(53000-47000)); 

T5=T4+c*(79000-53000); 
P5=P4*(T5/T4)^(-g/(c*R)); 
rho5=rho4*(T5/T4)^((-g/(c*R))-1); 

T6=T5; 
P6=P5*exp((-g/(R*T6))*(90000-79000)); 
rho6=rho5*exp((-g/(R*T6))*(90000-79000)); 

T7=T6+d*(100000-90000); 
P7=P6*(T7/T6)^(-g/(d*R)); 
rho7=rho6*(T7/T6)^((-g/(d*R))-1); 

if h<=11000 
T=T0+a*h; 
P=P0*(T/T0)^(-g/(a*R)); 
rho=rho0*(T/T0)^((-g/(a*R))-1); 
elseif h>11000 && h<=25000 
T=T1; 
P=P1*exp((-g/(R*T))*(h-11000)); 
rho=rho1*exp((-g/(R*T))*(h-11000)); 
elseif h>25000 && h<=47000 
T=T2+b*(h-25000); 
P=P2*(T/T2)^(-g/(b*R)); 
rho=rho2*(T/T2)^((-g/(b*R))-1); 
elseif h>47000 && h<=53000 
T=T3; 
P=P3*exp((-g/(R*T3))*(h-47000)); 
rho=rho3*exp((-g/(R*T3))*(h-47000)); 
elseif h>53000 && h<=79000 
T=T4+c*(h-53000); 
P=P4*(T/T4)^(-g/(c*R)); 
rho=rho4*(T/T4)^((-g/(c*R))-1); 
elseif h>79000 && h<=90000 
T=T5; 
P=P5*exp((-g/(R*T))*(h-79000)); 
rho=rho5*exp((-g/(R*T))*(h-79000)); 
elseif h>90000 && h<=100000 
T=T6+d*(h-90000); 
P=P6*(T/T6)^(-g/(d*R)); 
rho=rho6*(T/T6)^((-g/(d*R))-1); 

end

来源

2017-01-30 Henry Macanas

绘制每个点一个接一个是超效率的。为什么不建立T，P和rho向量（每个都有1000个元素），然后做3个地块...... – Justin

同意@Justin。首先在矢量中生成所有点，然后每个图只调用一次“plot/subplot”。循环执行效率不高，更不用说它会对图形系统造成影响。 – rayryeng

当我运行我的代码时，我得到的vecors h，T，P和rho的维数为100。端保持上副区（1,3,1）情节（T，H）副区（1,3,2）积（P，H）副区（1,3,3）情节（rho，h） –

尝试hold每每一次要情节：

for h=0:1:1000 
[ T,P,rho ] = stdatm(h); 
subplot(1,3,1) 
plot(T,h,'r') 
hold on 
subplot(1,3,2) 
plot(P,h) 
hold on 
subplot(1,3,3) 
plot(rho,h) 
hold on 
end 
hold off

尽管如此，向量化函数stdatm将是更好的解决方案。

来源

2017-01-30 00:40:37 ibancg

当我运行我的代码时，我只得到h，T，P，rho的一维向量。因此每个图上有一个数据点。

for h=0:1:100 
[ T,P,rho ] = stdatm(h); 
end 
hold on 
subplot(1,3,1) 
plot(T,h) 

subplot(1,3,2) 
plot(P,h) 

subplot(1,3,3) 
plot(rho,h)

来源

2017-01-30 01:21:51

您可以在每个循环中覆盖T，P，rho的值。你可能想要做'[T（h + 1），P（h + 1），rho（h + 1）] = stdatm（h）;''另外，使用预分配àla'T = zeros（1,101）。 ......在效率循环之前。 – Florian

在MATLAB中绘制一个输入和三个输出的子函数

回答

相关问题