Listato

function [ z ] = es6(  )

%ES6 Calcolo potenza

 

%Caricamento dati

prn=load('s00.prn');

t=prn(:,1);

vce=prn(:,2);

vr=prn(:,4);

dati=load('data.txt');

vcemax=dati(1);

R1=dati(2);

vbeSat=dati(3);

ib=dati(4);

 

ic=vr/R1;

 

%Correzione dell'errore di offset

eo=max(vce)-vcemax;

vce=vce-eo;

 

%Selezione dell'intervallo

j=find(t>-0.2e-3 & t<0.6e-3);

vce=vce(j);

ic=ic(j);

t=t(j);

 

j=find(vce>0);

vce=vce(j);

ic=ic(j);

t=t(j);

 

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);

 

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);