Hadi Saadat

Diagrama de carga
data=[0,2,6;2,6,5;6,9,10;9,12,15;12,14,12;14,16,14;16,18,16;18,20,18;20,22,16;22,23,12;23,24,6];
P=data(:,3);
Dt=data(:,2)-data(:,1);
W=P'*Dt;
Pavg=W/sum(Dt);
Peak=max(P);LF=Pavg/Peak*100;
bar(data);
xlabel('time,hr'),ylabel('P,MW')

Poncia 2.1
Vm=100; thetav=0; %amplitud y +angulo del voltaje
Z=1.25; gama=60;
thetai=thetav-gama;theta=(thetav-thetai)*pi/180;
Im=Vm/Z;
wt=0:.05:2*pi; %wt desde 0-2pi
v=Vm*cos(wt);
i=Im*cos(wt+thetai*pi/180);
p=v.*i;
V=Vm/sqrt(2); I=Im/sqrt(2);
P=V*I*cos(theta);
Q=V*I*sin(theta);
S=P+1i*Q;pr=P*(1+cos(2*(wt+thetav)));
px=Q*sin(2*(wt+thetav));
PP=P*ones(1,length(wt));
xline=zeros(1,length(wt));
wt=180/pi*wt;
subplot(2,2,1), plot(wt, v, wt, i, wt, xline), grid
title(['v(t)=Vmcoswt, i(t)=Im cos(wt+',num2str(thetai),')'])
xlabel('wt, degree')
subplot(2,2,2), plot(wt, p, wt, xline), grid
title('p(t)=v(t) i(t)'), xlabel('wt,degree')
subplot(2,2,3), plot(wt, pr, wt, PP,wt,xline), grid
title('pr(t) P*(1+cos(2*(wt+thetav)))'), xlabel('wt, degree')
subplot(2,2,4), plot(wt, px, wt, xline), grid
title('px(t) Q*sin(2*(wt+thetav))'), xlabel('wt, degree'), subplot(111)Pag39
Potencia Real 2.6
E1 = input('source #1 Voltage Mag. = ');
a1 = input('source #1 Phase Angle = ');
E2 = input('source #2 Voltage Mag. = ');
a2 = input('source #2 Phase Angle = ');
R =input('Line Resistance = ');
X = input('Line Reactance = ');
Z = R +1i*X;
a1 = (-30+a1:5:30+a1);
a1r = a1*pi/180;
k=length(a1);
a2=ones(k,1)*a2;
a2r=transpose(a2*pi/180);
V1 = E1.*cos(a1r) +1i*E1.*sin(a1r);
V2 = E2.*cos(a2r) + 1i*E2.*sin(a2r);
I12 = (V1 - V2)./Z;
I21=-I12;
S1=V1.*(conj(I12)); P1=real(S1); Q1=imag(S1);
S2=V2.*conj(I21); P2=real(S2); Q2=imag(S2);
SL=S1+S2;PL=real(SL); QL=imag(SL);
Result1=[a1,P1,P2,PL];
disp(' Delta 1 P-1 P-2 P-L' )
disp(Result1)
plot(a1,P1,a1,P2,a1,PL)
xlabel('Source #1 Voltafe Phase Angle')
ylabel(' P, Watts'),...