离散系统变换域分析.

离散系统变换域分析

实验目的

1 熟悉对离散系统的频率响应分析方法；

2 加深对零，极点分布的概念的理解。

实验内容：

程序：

b=[0.036 0.143 0.214 0.143 0.036];

a=[1 -1.035 0.826 0.260 0.040];

[z,p,k]=tf2zp(b,a);

disp('零点');disp(z);

disp('极点');disp(p);

disp('增益系数');disp(k);

sos=zp2sos(z,p,k);

disp('二阶节');disp(real(sos));

zplane(b,a);

k=256;

b=[0.036 0.143 0.214 0.143 0.036];

a=[1 -1.035 0.826 0.260 0.040];

w=0:pi/k:pi;

h=freqz(b,a,w);

subplot(2,2,1);

plot(w/pi,real(h));grid

title('实部')

xlabel('\omega/\pi');

subplot(2,2,2);

plot(w/pi,imag(h));grid

title('虚部')

xlabel('\omega/\pi');

subplot(2,2,3);

plot(w/pi,abs(h));grid

title('幅度谱')

xlabel('\omega/\pi');ylabel('幅值')

subplot(2,2,4)

plot(w/pi,angle(h));grid

title('相位谱')

xlabel('\omega/\pi');ylabel('相位')

实验结果：

零点

-0.9861 + 0.1661i

-0.9861 - 0.1661i

-1.0000

-1.0000

极点

0.6516 + 0.8460i 0.6516 - 0.8460i -0.1341 + 0.1308i -0.1341 - 0.1308i

增益系数 0.0360

二阶节

0.0360 0.0720 0.0360 1.0000 0.2681 0.0351

1.0000 1.9722 1.0000 1.0000 -1.3031 1.1403

-1

-0.5

00.5

1

-1

-0.8-0.6-0.4-0.200.20.40.60.8

1Real Part

I m a g i n a r y P a r t

0.51

-1012

3实部

ω/π0

0.51

-1012

3虚部

ω/π0

0.510

12

3幅度谱

ω/π

幅值

0.51

-4

-20

2相位谱

ω/π

相位