Octave/Matlab-oppgave: Frekvensresponsanalyse

L�sning

• Octave: Definerer transferfunksjonen som system-variabel av typen transferfunksjon med tf2sys:

  H=tf2sys([1],[1,1]);

• Matlab: Definerer en LTI-modell (Linear Time Invariant) av typen transferfunksjon med tf:

  H=tf([1],[1,1]);

Plotter frekvensresponsen i Bodediagram med bode:

wmin=0.5;
wmaks=5;
dw=0.1;
bode(H,[wmin:dw:wmaks])

Fra amplitude(-forsterknings)diagrammet i Bodediagrammet ser vi at amplituden er -3dB ved frekvensen w=1rad/s, som derfor er b�ndbredden.

Nedenfor brukes lsim til � simulere og plotte responsen i utgangsvariabelen y n�r inngangssignalet er en sinus med frekvens lik 1rad/s. Windowsversjonen av Octave st�tter visstnok ikke plotting i flere figurvinduer samtidig. Derfor m� evt. tidligere genererte plott slettes f�r resultatet av lsim-simuleringen plottes. Dette gj�res med kommandoene

clearplot (som lukker eksisterende figur)

etterfulgt av

closeplot (som avslutter og nullstiller Gnu-plott"maskinen" som jobber i bakgrunnen i Octave)

Etter ovennevnte to kommandoer skal alt v�re klart for � generere nye plott.

t=[0:0.1:30]';
w=1;%rad/s
u=sin(w*t);
y=lsim(H,u,t);
plot(t,u,t,y);

Hvis du vil ha rutenett:

• Octave:

  grid;
  replot

• Matlab:

  grid;

Amplitudeforholdet mellom den resulterende sinusformede utgangssignalet og sinusinngangen b�r v�re -3dB = 0,71. Ser du dette fra simuleringen?

Oppdatert 7.10.04 av Finn Haugen, fagl�rer. E-postadresse: Finn.Haugen@hit.no.