o_callcheck.m

来自「电力能源预测」· M 代码 · 共 46 行

function varargout=O_callcheck(O_h2text,str1,o,I,G,K,N,O,Gc)


global o_num;
cla;

set(O_h2text,'string',[str1,sprintf('%1.4g\',o)]);


O(o_num)=o;

for i=26:42
    I(i)=I(i-1)*K(i);
end

for i=26:42
    G(i)=G(i-1)+Gc(i-25);
end

for i=1:42
     L(i)=G(i)-I(i)/(87.6*(1-0.479195-0.005*N(i)+0.005*O(i)));
 end
 
 for i=1:42
     J(i)=(I(i)/(G(i)*87.6)+0.479195-0.005*N(i)+0.005*O(i))/(I(i)/((G(i)-L(i))*87.6)+0.479195-0.005*N(i)+0.005*O(i));
 end
 if nargout>0
 varargout{1}=I;
 varargout{2}=J;
 varargout{3}=K;
 varargout{4}=N;
 varargout{5}=O;
 varargout{6}=Gc;
 varargout{7}=L;
 end
  h1_axes=axes('position',[0.05,0.1,0.6,0.35]);
 t=1979:1:2020;
hbar=bar(t,L);
axis([1979 2020 -40 40]);
 h2_axes=axes('position',[0.05,0.55,0.6,0.35]);
t=1979:1:2020;
hline=plot(t,J);
axis([1979 2020 0.9  1.1]);