请问高手下面的程序怎么运行啊?
偶是刚接触matlab,弱的不知道为什么不能运行啊!!!请达人指点!!!k=1:3000
D=0.005;
A=3.1415926/4*D*D;
%h=1;
%t=1e-5;
rouN=1.293;
miu=(4.8e-8)*20+1.71e-5;
G(1)=8.3336;
Deta_G(1)=0;
u(1)=2.5745;
rou(1)=3.553;
Re(1)=2532.7;
lamta(1)=0.044603;
for j=2:3000
%if(lamta(j-1)==0 | u(j-1)==0)
%Deta_G(j)=0;
%else
Deta_G(j)=-A*D/(lamta(j-1)*u(j-1))*DP_t(j)*1e-3*(exp(-lamta(j-1)*u(j-1)/D*(1+j*1e-3))-1)/rouN*60000;
G(j)=G(j-1)+Deta_G(j-1);
u(j)=G(j)*rouN/(rou(j)*A*60000);
Re(j)=rou(j)*u(j)*D/miu;
%ifRe(j)<=0
%lamta(j)=0;
%else if (Re(j)>0 & Re(j)<=2300)
%lamta(j)=64/Re(j);
%else
lamta(j)=0.3164*Re(j)^(-0.25);
%end
end
plot(1+k*1e-3,G(k)); 原帖由 yamalu 于 2007-5-22 15:37 发表 http://www.chinavib.com/forum/images/common/back.gif
偶是刚接触matlab,弱的不知道为什么不能运行啊!!!请达人指点!!!
k=1:3000
D=0.005;
A=3.1415926/4*D*D;
%h=1;
%t=1e-5;
rouN=1.293;
miu ...
请先看看 置顶帖:聚宝盆,然后重新整理你的问题,特别是给出错误信息,不要指望别人帮你调试
页:
[1]