Listato

function [ z ] = potenzaSin( fase )

%POTENZASIN Calcolo potenza in regime sinusoidale

 

Fs=1000;

t=[0:1/Fs:1-1/Fs];

vce=sin(2*pi*3*t+fase);

ic=10*sin(2*pi*3*t);

 

N=length(t);

 

%Potenza attiva

Pint=sum(vce.*ic)/N;

 

%FFT

In=fft(ic);

Vn=fft(vce);

 

%Potenza attiva

P=sum(abs(Vn).*abs(In).*cos(abs(angle(Vn)-angle(In))))/N2;

 

%Potenza reattiva secondo Budeanu

Q=sum(abs(Vn).*abs(In).*sin(abs(angle(Vn)-angle(In))))/N2;

 

%Potenza apparente

S=sqrt(sum(abs(Vn).2)*sum(abs(In).2))/N2;

St=sqrt(P2+Q2);

 

%Potenza di distorsione

D=sqrt(S2-St2);

 

%Fasori

Pc=10/sqrt(2)*1/sqrt(2)*cos(fase);

Qc=10/sqrt(2)*1/sqrt(2)*sin(fase);

Sc=10/sqrt(2)*1/sqrt(2);

 

stringa=['Potenza attiva con il metodo dell''integrazione discreta: Pa = ' num2str(Pint) ' mW'];

disp(stringa);

 

stringa=['Potenza attiva con il metodo della FFT: Pa = ' num2str(P) ' mW'];

disp(stringa);

 

stringa=['Potenza reattiva secondo Budeanu: Qb = ' num2str(Q) ' mW'];

disp(stringa);

 

stringa=['Potenza apparente: S = ' num2str(S) ' mW'];

disp(stringa);

 

stringa=['Potenza apparente ottenuta dal triangolo delle potenze: St = ' num2str(St) ' mW'];

disp(stringa);

 

stringa=['Potenza di distorsione secondo Budeanu: D = ' num2str(D) ' mW'];

disp(stringa);

 

stringa=['Potenza attiva con il metodo dei fasori: Pa = ' num2str(Pc) ' mW'];

disp(stringa);

 

stringa=['Potenza reattiva con il metodo dei fasori: Q = ' num2str(Qc) ' mW'];

disp(stringa);

 

stringa=['Potenza apparente con il metodo dei fasori: S = ' num2str(Sc) ' mW'];

disp(stringa);