文档介绍:数字信号处理
实验三
离散时间系统的时域分析
学院:信息与通信学院
专业:电子信息工程
学号:0900220418
姓名:梁芝铭
实验目的
(1)理解离散时间信号的系统及其特性。
(2)对简单的离散时间系统进行分析,研,den,x1,ic);
y2 = filter(num,den,x2,ic); y = filter(num,den,x,ic);
yt = a*y1 + b*y2; d = y - yt; subplot(3,1,1);stem(n,y); ylabel('Amplitude'); title('Output Due to Weighted Input: a \cdot x_{1}[n] + b \cdot x_{2}[n]'); subplot(3,1,2);stem(n,yt); ylabel('Amplitude');
title('Weighted Output: a \cdot y_{1}[n] + b \cdot y_{2}[n]'); subplot(3,1,3);stem(n,d); xlabel('Time index n'); ylabel('Amplitude');title('Difference Signal');
结果如下:
Fig^e 1 □冋冈
Output Due to Weighted Input a • x〔[n] + b • x2[n]
QC|l:+:lr?i
口GG':+:G「-,
◎:+:◎■;=
File Edit View Insert Tools Desktop Window Kelp d已h◎「氐⑥祖銘”鳳[□因「r口
-sn=-dE<
-s3-dE<
yjn] + b-y2[n]
、919 S、S,、 G 伞川名?
G W 名"W
6U ®
66 000
1 1
-sn=-dE<
时不变与时变系统的仿真
在 MATLAB 中输入:
% Generate the input sequences
clf; n = 0:40; D = 10; a = ; b = -2;
x = a*cos(2*pi**n) + b*cos(2*pi**n);
xd = [zeros(1,D) x]; num = [ ]; den = [1 - ]; ic = [0 0]; % Set initial conditions
% Compute the output y[n]
y = filter(num,den,x,ic);
% Compute the output yd[n]
yd = filter(num,den,xd,ic);
% Compute the difference output d[n]
d = y - yd(1+D:41+D);
grid;
subplot(3,1,1); stem(n,y); ylabel('Amplitude'); title('Output y[n]'); subplot(3,1,2); stem(n,yd(1:41)); ylabel('Amplitude');
title(['Output