yamalu 发表于 2007-5-22 15:37

请问高手下面的程序怎么运行啊?

偶是刚接触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));

eight 发表于 2007-5-22 15:39

原帖由 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]
查看完整版本: 请问高手下面的程序怎么运行啊?