Dus ik heb de overdrachtsfunctie:

$$ H [z] = 1 + \ sqrt {2} z ^ {- 1} + z ^ {- 2} $$

En ik moet $ H (e ^ {j \ omega}) $ voor $ \ omega = 0, \ pi / 4, \ pi / 4 \ ldots $

Ik heb de berekeningen handmatig gedaan met de formule van Euler, maar nu is de toewijzing ik vraag me om deze plots te vergelijken met de plots met freqz in MATLAB. Ik kan “geen instructies vinden over hoe ik dat kan doen met dit type overdrachtsfunctie.

Opmerkingen

  • Ik kan ' zelfs niet: D Dus, hint: elk getal is $ x $ wordt weergegeven door $ \ frac xy $ voor een specifiek getal $ y $. Altijd. Wat ' is dat $ y $?
  • Van wat ik kan zien, heb je de teller (b) van uw filter. Sluit het dus gewoon aan op freqz en voila.

Answer

U specificeert eenvoudig a = 1 (omdat de noemer gelijk is aan $ 1 $). Dus je krijgt

 b = [1,sqrt(2),1]; a = 1; N = 512; [H,w] = freqz(b,a,N); 

Je kunt dit vergelijken met de analytische oplossing:

 H2 = 1 + sqrt(2)*exp(-1i*w) + exp(-1i*2*w); max(abs(H2-H)) % 8.0825e-16 

Opmerkingen

  • Sorry, ik ' ben hier echt nieuw in, maar wat vertegenwoordigt N hier?
  • @Freddie: Het ' is het aantal (equidistante) frequentiepunten waarop de frequentierespons wordt geëvalueerd. Bekijk de Matlab-documentatie van freqz .

Antwoord

Voor de evaluatie alleen bij specifieke frequenties, moet u de frequentievector specificeren met ten minste twee frequenties erin (zie MATLAB “s freqz ). Hieronder staat de MATLAB-code voor de evaluatie bij de frequenties $ \ omega = 0, \ pi / 4, \ pi / 2, 3 \ pi / 4, \ text {and} \ \ pi $ .

>> [h, w] = freqz([1, sqrt(2), 1], 1, [0 , pi/4, pi/2, 3*pi/4 pi]) h = 3.4142 + 0.0000i 2.0000 - 2.0000i 0.0000 - 1.4142i -0.0000 - 0.0000i 0.5858 + 0.0000i w = 0 0.7854 1.5708 2.3562 3.1416 >> 

Voor de visualisatie van de resultaten hierboven, zie de omvang antwoord, dwz $ 20 \ log_ {10} \ left (\ lvert H \ left (\ omega \ right) \ rvert \ right) $ , hieronder uitgezet met de vijf frequenties rood gemarkeerd.

voer hier een afbeeldingsbeschrijving in

Merk op dat voor $ \ pm 3 \ pi / 4 $ je dat hebt (zie coderesultaten hierboven) $$ H \ left (\ p m \ frac {3 \ pi} {4} \ right) = 0 \ impliceert 20 \ log_ {10} \ left (\ bigg \ lvert H \ left (\ pm \ frac {3 \ pi} {4} \ right) \ bigg \ rvert \ right) = – \ infty $$ Ook vanwege het feit dat de nullen op $$ z = – \ frac {\ sqrt {2}} staan {2} \ pm j \ frac {\ sqrt {2}} {2} \ quad \ text {with} \ quad z = e ^ {j \ omega} $$ De corresponderende magnitude voor $ \ omega = 3 \ pi / 4 $ wordt niet weergegeven in de enkelzijdige magnitude-responsplot hierboven, maar je kunt de asymptotische trend zien op $ 3 \ pi / 4 $ .

Geef een reactie

Het e-mailadres wordt niet gepubliceerd. Vereiste velden zijn gemarkeerd met *