plotcdf.m

来自「可以计算电力系统概率潮流」· M 代码 · 共 38 行

clc
obsline=csvread('obsline.csv');
nobl=length(obsline(:,1));

for i=1:nobl
    nfrom=obsline(i,1);
    nto=obsline(i,2);
    ni=locLN(nfrom, nto, line);
    mui_p=PLineX(ni,1);
    smgi_p=sqrt(PLineX(ni,2));
    xi_p=mui_p;
    for j=1:300
        PCL(i,j)=xi_p+(j-150)*1;
        PPL(i,j)=CalCdf1(PCL(i,j), mui_p, smgi_p, PLineC(ni, 3), PLineC(ni,4), PLineC(ni,5), PLineC(ni,6), PLineC(ni,7), PLineC(ni,8), PLineC(ni,9));
    end
        
    figure(2*i-1);
    plot(PCL(i,:)+mui_p, PPL(i,:))
    htitle=title('线路有功功率cdf曲线图 ','color',[0 1 1]);
    ylabel('概率','color',[1 0 1]);
    xlabel('MW','color',[1 1 0]);
    legend('pdf的曲线')
    
    mui_q=PLineX(ni+nline,1);
    smgi_q=sqrt(PLineX(ni,2));
    xi_q=mui_q;
    for j=1:300
        PCL(i,j)=xi_q+(j-150)*1;
        PQL(i,j)=CalCdf1(PCL(i,j), mui_q, smgi_q, PLineC((ni+nline), 3), PLineC((ni+nline),4), PLineC((ni+nline),5), PLineC((ni+nline),6), PLineC((ni+nline),7), PLineC((ni+nline),8), PLineC((ni+nline),9));
    end 
    figure(2*i+1);
    plot(PCL(i,:), PQL(i,:))
    htitle=title('线路无功功率cdf曲线图 ','color',[0 1 1]);
    ylabel('概率','color',[1 0 1]);
    xlabel('Mvar','color',[1 1 0]);
    legend('cdf的曲线')
end