Corsi di Laurea

[Kinn]

Salve, 
Stiamo cercando di svolgere quest'ultimo quesito (allego la prova in calce).
... Abbiamo calcolato la funzione di rete e ricostruito l'ingresso a onda quadra. Il problema è quello di individuare il numero corretto di armoniche necessarie per ricostruire correttamente l'uscita.
 

http://pastebin.com/cBdR6asK

 

Questo è il codice da noi sviluppato.
Allego anche il grafico che ne esce:
 

File allegato

•  etinf_provascritta_19092013.pdf   158,74K   9 Download
•  Schermata 2013-12-19 alle 15.51.28.png   26,04K   0 Download

[LeoVa]

Non ho capito due cose del codice: come ti esce la funzione di rete? Cioè il risultato. Poi perché alfa=0.05? Nella traccia non c'è indizio su questo.

 

Ad ogni modo io l'ho buttata lì così: a me la funz di rete esce un filtro passa basso, quindi attenua molto le armoniche ad alta frequenza (o pulsazione), perciò ho commentato che mi bastano le prime 10 armoniche (quelle a bassa frequenza cioè) per avere una buona ricostruzione del segnale, in quanto quelle a frequenza superiore (dato che la funz di rete mi è uscita un passa basso) verranno attenuate abbastanza da non avere un gran peso sulla ricostruzione dell'uscita.

 

L'uscita che mi trovo è perfettamente uguale alla tua, e questo è il codice che ho usato:

Spoiler

%Punto 5

syms AH AL t1 t2 T n

eq1=AH-AL-1;

eq2=AH*0.33+AL*0.67;

sol2=solve(eq1,eq2,'AH,AL');

AH=sol2.AH; AL=sol2.AL; f=100; T=1/f; wq=2*pi*f; t1=0.33*T; t2=0.67*T;

% %Calcolo funzione di rete

R23=(R2*R3)/(R2+R3)

H=(R3*1/(j*n*wq*C2))/((1/(j*n*wq*C2)+R23)*(R2+R3));

%Dall'analisi della forma della funzione di rete vediamo che è un filtro

%passa basso quindi ci vanno bene 20 armoniche

N=20;

An_sym=(2/T)*(int(AH*cos(n*wq*t),t,0,t1)+int(AL*cos(n*wq*t),t,t1,T));

Bn_sym=(2/T)*(int(AH*sin(n*wq*t),t,0,t1)+int(AL*sin(n*wq*t),t,t1,T));

j_t=0;

iR2=0;

for n=1:N

    An=subs(An_sym);

    Bn=subs(Bn_sym);

    FAN=An;

    FBN=Bn*exp(-j*pi/2);

    FCN=FAN+FBN;

    j_t=j_t+abs(FCN)*cos(n*wq*t+angle(FCN));

    FIRn=subs(H,'n',n)*FCN;

    iR2=iR2+abs(FIRn)*cos(n*wq*t + angle(FIRn));

end

 

figure(3), ezplot(j_t,[0,3*T]),hold on, grid on, axis auto, ezplot(iR2,[0,3*T]), axis auto, grid on