代码有什么问题，请帮助我。我想显示从j = 1到j = 24的结果。

lambda0 = 9.6; 
gamma0 = 0.0016; 
Pload = 900; 
nj=1; 
AvW = 32e6; 
t = 12; 
%thermal plant%%% 
ath = 0.002; 
bth = 8; 
gamth = 1200; 
cth = 0; 
Pth0 = (lambda0 - (nj*bth))/((2*nj*ath)-(2*cth*lambda0)); 
%%%%%hydro plant%%%%% 
ah = 6000; 
bh = 100000; 
ch = 100; 
Ph0 = (lambda0 - (nj*ah*gamma0))/(2*lambda0*ch); 
Tl = 0.0000045; 



%%%%%%%%%%%lambda%%%%%%%%%%%%%%%%%%%%%% 
j=1; 
for j=1:24 

    x = Pload - Ph0 - Pth0; 

    e1 = 0.1; 

    while x > e1 

     lambda = lambda0 + 0.25; 
     Pth = (lambda - (nj*bth))/((2*nj*ath)-(2*cth*lambda)); 
     Ph  = (lambda - (nj*ah*gamma0))/(2*lambda*ch); 
     Ploss = 0.0000045 * Ph^2; 
     x = Pload + Ploss - Ph - Pth; 

     if x < e1 
     lambda = lambda0; 
     Pth = (lambda - (nj*bth))/((2*nj*ath)-(2*cth*lambda)); 
     Ph  = (lambda - (nj*ah*gamma0))/(2*lambda*ch); 
     Ploss = 0.0000045 * Ph^2; 
     end 

    W = bth + ath*Ph; 
    fprintf('W=%.3f',W); 
    x = Pload + Ploss - Ph - Pth; 
    e2 = 50 ; 
    dW = W - AvW; 
    end 
end 
while dW > e2 

lambda = (nj *(gamma0*ah)-(2*ath*bth))/((2*cth*Pth)-(2*Tl*Ph)); 
x = Pload-Ph-Pth; 
while x(1)> e1 

    lambda = lambda + 0.25; 
    gamma = gamma0 + 0.00001; 
    Pth = (lambda - (nj*bth))/((2*nj*ath)-(2*cth*lambda)); 
    Ph  = (lambda - (nj*ah*gam))/(2*lambda*ch); 
    Ploss = 0.0000045 * Ph^2; 
    x  = Pload + Ploss - Ph - Pth; 

    if x(1)< e1 
    lambda = lambda0; 
    Pth = (lambda - (nj*bth))/((2*nj*ath)-(2*cth*lambda)); 
    Ph  = (lambda - (nj*ah*gamma))/(2*lambda*ch); 
    Ploss = 0.0000045 * Ph^2; 
    x  = Pload + Ploss - Ph - Pth; 
    W  = bth + ath*Ph; 
    fprintf('w=%.3f',W); 
    end 
end 
    tot = W * t ; 

if e2 > tot 
    fprintf('\t%.3f',lambda); 
    fprintf('\t%.3f',Pth); 
    fprintf('\t%.3f',Ph); 
    fprintf('\t%.3f',gamma); 
    fprintf('\n'); 
end 
end

来源

2011-04-08 user694666

解释你所得到的错误会有所帮助。 – gnovice 2011-04-08 17:04:16