Por favor el pzmap código

-que símbolos y nomencaltura usa

[John D'Errico]

Rough translation: Please, the pzmap code. What symbols and nomenclature does it use?

All I can say is, what about the help docs for pzmap was unclear?

https://www.mathworks.com/help/ident/ref/lti.pzmap.html

Answer 1

[William Rose]

@Mario Gustavo,

(Traducido con Google.)

Aquí se muestra un ejemplo del uso de pzmap para un sistema de tiempo continuo: un filtro antialiasing con 8 polos y 8 ceros.

Here is an example of the use of pzmap for a continuous time system: an anti-aliasing filter with 8 poles and 8 zeros.

% Design 8th order anti-aliasiong filter (continuous time).
n=8;                % number of poles
fco=1.00;           % cutoff frequency (Hz)
Wp=2*pi*fco;        % cutoff frequency (rad/s)
Rp=1;               % passband ripple (dB)
Rs=80;              % stopband attenuation (dB)
[z,p,k] = ellip(n,Rp,Rs,Wp,'s');
% Convert the zeros, poles, and gain (k) to a system
sys=zpk(z,p,k);
% Plot the poles and zeros.
figure
pzmap(sys)
title('Pole-Zero Map: 8th Order Lowpass Filter, Continuous-Time')

Este ejemplo muestra el uso de pzmap con un sistema de tiempo discreto: un filtro de muesca para eliminar el ruido de 50 Hz de una señal. Este filtro tiene una atenuación de al menos 40 dB entre 48 Hz y 52 Hz. También agregué un gráfico que muestra el círculo de estabilidad para un sistema discreto.

This example shows the use of pzmap with a discrete-time system: a notch filter to remove 50 Hz noise from a signal. This filter has at least 40 dB attenuation from 48 Hz to 52 Hz. I also added a plot showing the circle of stability for a discrete system.

% Design 4+4 order 50 Hz notch filter, discrete time
n=4;                % number of poles (bandstop filter will have 2n poles)
fs=1000;            % sampling frequency (Hz)
fco=[48,52];        % stop band frequencies (Hz)
Wp=fco/(fs/2);      % stop band frequencies (normalized)
Rp=1;               % passband ripple (dB)
Rs=40;              % stopband attenuation (dB)
[z,p,k] = ellip(n,Rp,Rs,Wp,'stop');
% Convert the zeros, poles, and gain (k) to a system
sys=zpk(z,p,k);
% Plot the poles and zeros.
figure
pzmap(sys)
hold on
% Plot stability boundary
plot(cos(2*pi*(0:100)/100),sin(2*pi*(0:100)/100),'-r')
xlim([-1.2,1.2]); ylim([-1.2,1.2]);
title('Pole-Zero Map: 4+4 Order Notch Filter, Discrete-Time')
legend('Zeros and Poles','Stability Boundary')