文档库

最新最全的文档下载
当前位置:文档库 > 实验4离散系统地分析报告

实验4离散系统地分析报告

实验四 离散系统分析

一、 实验目的

深刻理解离散时间系统的系统函数在分析离散系统的时域特性、频域特性以及稳定性中的重要作用及意义,熟练掌握利用MATLAB 分析离散系统的时域响应、频响特性和零极点的方法。掌握利用DTFT 和DFT 确定系统特性的原理和方法。

二、实验原理

可以在时域、复频域(Z 域)及频域分析系统,在以上三种域表征系统固有特性的量分别为:

① 单位冲激响应 h (n )(时域表征);

② 系统函数 H (z ) ( Z 域表征);

③ 频率响应 H (e j ω)(频域表征)。

MATLAB 主要从以上三方面提供了许多可用于分析线性时不变系统的函数,包含系统时域响应、系统函数、系统频域响应等分析函数。

本实验通过调用各种系统预置函数来求系统的以上几个表征量以及零极点图。

三、实验内容

1.已知某LTI 系统的差分方程为:

(1)初始状态 ,输入 ,

计算系统的全响应。 程序段:

N=40;

b=[0.0675,0.1349,0.0675];

a=[1,-1.143,0.412];

x=ones(1,N);

zi=filtic(b,a,[1,2]);

y=filter(b,a,x,zi);

stem(y)

xlabel('k');

title('y[k]');

]

2[0675.0]1[1349.0][0675.0]2[412.0]1[143.1][-+-+=-+--k x k x k x k y k y k y 2]2[,1]1[=-=-y y ][][k u k x =

结果:

实验4离散系统地分析报告

2)当以下三个信号分别通过系统时,分别计算离散系统 的零状态响应: 程序

N=30;

k=0:N;

b=[0.0675,0.1394,0.0675];

a=[1,-1.143,0.412];

x1=cos(pi*0.1.*k);

x2=cos(pi*0.2*k);

x3=cos(pi*0.7*k);

y1=filter(b,a,x1);

y2=filter(b,a,x2);

y3=filter(b,a,x3);

subplot(3,1,1);

stem(y1)

subplot(3,1,2);

stem(y2)

][)107cos(][];[)5cos(][];[)10cos(][321k u k k x k u k k x k u k k x πππ===

subplot(3,1,3);

stem(y3):

结果:

实验4离散系统地分析报告

(3)该系统具有什么特性?

答:因果稳定。

2.已知某因果LTI 系统的系统函数为:

(1)计算系统的单位冲激响应。

程序:

N=50;

k=1:N;

b=[0.0357,0.1428,0.2143,0.1428,0.0357];

4

3214

32104033.02605.08264.0035.1103571.01428.02143.01428.003571.0)(--------+-+-++++=z z z z z z z z z H

a=[1,-1.035,0.8264,-0.2605,0.04033];

y1=impz(b,a,N);

stem(y1)

结果为:

实验4离散系统地分析报告

(2)当信号 通过系统时,计算系统的零状态响应。

程序为:

N=50;

k=1:N;

b=[0.0357,0.1428,0.2143,0.1428,0.0357];

a=[1,-1.035,0.8264,-0.2605,0.04033];

x=ones(1,N)+cos(pi*0.25*k)+cos(0.5*pi*k);

y1=impz(b,a,N);

y2=filter(b,a,x);

subplot(2,1,1);

stem(y1)

subplot(2,1,2);

stem(y2)

结果为:

][)2cos(][)4cos(][][k u k k u k k u k x ππ++=

实验4离散系统地分析报告

3.已知LTI 系统的输入输出序列分别为

(1)利用解析方法分别求解系统的单位取样响应。

程序:

N=20;

a1=[1,-0.75];

b1=[1.25,-19/16];

h1=impz(b1,a1,N);

subplot(2,1,1);

stem(h1)

b2=[1,-1];

a2=[1];

h2=impz(b2,a2,N);

subplot(2,1,2);

stem(h2)

结果为:

]1[)41(][)41(][],[)41(][)(][)4

1(][)21(41][],[)21(][)(1--==+==-k u k u k y k u k x b k u k u k y k u k x a k k k k k k ][)2cos(][)4cos(][][k u k k u k k u k x ππ++

=

实验4离散系统地分析报告

(2)利用系统辨识原理确定并求出系统的单位脉冲响应h[k]。比较解析方法与系统辨识方法得到的系统单位冲激响应,分析误差原因。

程序:

N=50;

k=0:N;

x1=0.5.^k;

y1=0.25*0.5.^k+0.25.^k;

X1=fft(x1);

Y1=fft(y1);

H1=Y1/X1;

h1=ifft(H1)

n=1:50;

x2=0.25.^n;

y2=0.25.^n-4*0.25.^n;

X2=fft(x2);

Y2=fft(y2);

H2=Y2/X2;

h2=ifft(H2)

结果:

H e

()j

误差原因:取样点数不够多引起的误差,而且不易与k对应。

4. 已知某离散系统的输入输出序列。

输入序列:2,0.8333,0.3611,0.162,0.0748,0.0354,0.017,0.0083,0.0041,0.002,0.001,0.0005,0.0002,0.0001,0.0001,后面的数值均趋于0;

输出序列:0.0056,-0.0259,0.073,-0.1593,0.297,

-0.4974,0.7711,-1.1267,1.5702,-2.1037,2.724,-3.4207,4.174,-4.9528,5.7117,-6.3889,6.9034,-7.1528,7.012,-6.3322,4.9416,-2.648,-0.7564,5.4872,-11.7557,19.7533,-29.6298,41.4666,-55.2433,70.7979,-87.7810 (1)绘出输入输出信号的波形。

程序:

x=[2,0.8333,0.3611,0.162,0.0748,0.0354,0.017,0.0083,0.0041,0.002,

0.001,0.0005,0.0002,0.0001,0.0001,zeros(1,16)];

y=[0.0056,-0.0259,0.073,-0.1593,0.297,-0.4974,0.7711,-1.1267,

1.5702,-

2.1037,2.724,-

3.4207,

4.174,-4.9528,

5.7117,-

6.3889,6.9034,

-7.1528,7.012,-6.3322,4.9416,-2.648,-0.7564,5.4872,-11.7557,19.7533, -29.6298,41.4666,-55.2433,70.7979,-87.7810];

subplot(2,1,1);

stem(x)

subplot(2,1,2);

stem(y)

结果:

实验4离散系统地分析报告

(2)计算该系统的频率响应,并绘出其幅频特性。

程序:

x=[2,0.8333,0.3611,0.162,0.0748,0.0354,0.017,0.0083,0.0041,0.002,

0.001,0.0005,0.0002,0.0001,0.0001,zeros(1,16)];

y=[0.0056,-0.0259,0.073,-0.1593,0.297,-0.4974,0.7711,-1.1267,

1.5702,-

2.1037,2.724,-

3.4207,

4.174,-4.9528,

5.7117,-

6.3889,6.9034,

-7.1528,7.012,-6.3322,4.9416,-2.648,-0.7564,5.4872,-11.7557,19.7533,

-29.6298,41.4666,-55.2433,70.7979,-87.7810];

X=fft(x,128);

Y=fft(y,128);

H=Y./X;

H0=abs(H);

plot(H0)

结果为:

实验4离散系统地分析报告

3)计算该系统的单位冲激响应,并绘出其波形。

程序:

x=[2,0.8333,0.3611,0.162,0.0748,0.0354,0.017,0.0083,0.0041,0.002, 0.001,0.0005,0.0002,0.0001,0.0001,zeros(1,16)];

y=[0.0056,-0.0259,0.073,-0.1593,0.297,-0.4974,0.7711,-1.1267,

1.5702,-

2.1037,2.724,-

3.4207,

4.174,-4.9528,

5.7117,-

6.3889,6.9034,

-7.1528,7.012,-6.3322,4.9416,-2.648,-0.7564,5.4872,-11.7557,19.7533, -29.6298,41.4666,-55.2433,70.7979,-87.7810];

X=fft(x,64);

Y=fft(y,64);

H=Y./X;

H0=abs(H);

h0=ifft(H0)

stem(h0)

结果为:

实验4离散系统地分析报告

实验4离散系统地分析报告

实验4离散系统地分析报告

图形为:

实验4离散系统地分析报告

5. 利用load mtlb 命令读入一段语音信号得到序列 ,然后在该段语音信号中加入500Hz 的正弦型干扰信号得到信号,利用FFT 分析其频谱。

(1)下列数字滤波器能够滤除信号中500Hz 正弦型干扰信号。

利用zplane 命令做出其零极点分布图,利用freqz 命令分析该滤波器的幅频特性和相频特性,比较零极点分布与滤波器频率特性的关系。

程序:

a=[1,-3.594,5.17,-3.494,0.945];

b=[0.6877,-2.509,3.664,-2.509,0.6877];

z=roots(b);

p=roots(a);

subplot(2,2,1);

zplane(b,a)

title('零极点分布

')

4

3214

321945.0494.3172.5594.316877.0509.2664.3509.26877.0)(--------+-+-+-+-=z z z z z z z z z H

[H,w]=freqz(b,a);

subplot(2,2,2);

plot(w,abs(H))

xlabel('\omega(rad)');

ylabel('H0');

title('幅频特性');

subplot(2,2,3);

plot(w,angle(H))

xlabel('\omega(rad)');

ylabel('\phi ');

title('相频特性');

结果:

实验4离散系统地分析报告

关系:系统零点越接近1,则幅频响应的波谷越低,即滤波器的滤波衰减最低点越低;系统极点越接近1,则幅频响应的波峰越高,即滤波器的滤波通带最高点越高。极点主要影响频率响应的峰值,极点愈靠近单位圆,峰值愈尖锐;零点主要影响频率特性的谷值,零点愈靠近单位圆,谷值愈深,当零点在单位圆上时,频率特性为零,一个传递函数有几个极点幅度响应就有几个峰值,对应出现一些谷值。

(2)利用该数字滤波器滤除信号中的噪声,利用FFT观察其频谱,利用sound 函数播放处理前后的信号,比较处理前后的效果。

程序;

load mtlb

N=150;

x=mtlb(1:N);

k=1:N;

subplot(1,2,1);

plot(k,x);

title('处理前');

X=fft(x,150);

h=ifft(H);

H=fft(H,150);

Y=X.*H; y=ifft(Y);

subplot(1,2,2);

plot((1:N),y);

title('处理后');

结果:

实验4离散系统地分析报告

实验思考题

1.系统函数的零极点对系统频率特性有何影响?

ANS:系统零点越接近1,则幅频响应的波谷越低,即滤波器的滤波衰减最低点越低;系统极点越接近1,则幅频响应的波峰越高,即滤波器的滤波通带最高点越高。极点主要影响频率响应的峰值,极点愈靠近单位圆,峰值愈尖锐;零点主要影响频率特性的谷值,零点愈靠近单位圆,谷值愈深,当零点在单位圆上时,频率特性为零,一个传递函数有几个极点幅度响应就有几个峰值,对应出现一些谷值。

2.对于因果稳定实系数的低通、高通、带通、带阻数字滤波器,零极点分布有何特点?

ANS:因为是因果稳定系统,所以极点都在单位圆内。若为最小相位系统,其零点也在单位圆内。

3.离散系统的系统函数的零极点对系统脉冲响应有何影响?

ANS:系统函数的极点位置决定序列包络的变化趋势和变化频率,极点的半径决定了序列包络的变化趋势,而极点的幅角决定序列包络的变化频率,而零点位置只影响冲激响应的幅度大小和相位。

4.若某因果系统不稳定,有哪些主要措施可使之稳定?

ANS:改变参数,让极点在单位圆内。

5.从频域利用DFT确定离散LTI系统的特性,一般会产生哪些误差,如何改善?

ANS: 频谱混叠,对于带限连续信号,只要提高抽样频率使之满足时域抽样定理;对于非带限信号,更具实际情况对其进行低通滤波,使之成为带限信号;频谱泄露,时域加窗使之成为有限长序列;栅栏现象,在序列后补零,构成新序列后再求频谱

6.若使用DFT对连续LTI系统进行辨识,需要解决哪些问题?

ANS:循环卷积长度的确定。