数字信号处理实验报告MATLAB第二章(2)

% y[n]和y2[n]之间的差

d = y - y2;

% 画出输出和差值信号

subplot(3,1,1);

stem(n,y);

ylabel('振幅');

title('四阶实现的输出'); grid;

subplot(3,1,2);

stem(n,y2)

ylabel('振幅');

title('级联实现的输出'); grid;

subplot(3,1,3);

stem(n,d)

xlabel('时间序号 n');ylabel('振幅');

title('差值信号'); grid;

从结果看出y[n]和y2[n]相等。

例2.7 卷积

2.28 运行程序P2.7，对序列h[n]和x[n]求卷积，生成y[n]，并用滤波器h[n]对输入x[n]滤波，求得y1[n]。y[n]和y1[n]有差别吗？为什么要使用对x[n]补零后得到的x1[n]作为输入来产生y1[n]？

%程序P2_7

clf;

h = [3 2 1 -2 1 0 -4 0 3];%冲激

数字信号处理MATLAB

x = [1 -2 3 -4 3 2 1];%输入序列

y = conv(h,x);

n = 0:14;

subplot(2,1,1);

stem(n,y);

xlabel('时间序号 n'); ylabel('振幅');

title('用卷积得到的输出'); grid;

x1 = [x zeros(1,8)];

y1 = filter(h,1,x1);

subplot(2,1,2);

stem(n,y1);

xlabel('时间序号 n'); ylabel('振幅');

title('由滤波生成的输出'); grid;

根据结果可知y[n]和y1[n]没有差别。

2.29修改程序P2.7，计算长度为15的序列h[n]和长度为10的序列x[n]的卷积，重做问题Q2.28。h[n]和x[n]的样本值你自己给定。

%程序P2_7

clf;

h = [3 2 1 -2 1 0 -4 0 3 2 1 5 7 6 2];%冲激

x = [1 -2 3 -4 3 2 1 -5 3 2];%输入序列

y = conv(h,x);

n = 0:23;

数字信号处理MATLAB

subplot(2,1,1);

stem(n,y);

xlabel('时间序号 n'); ylabel('振幅'); title('用卷积得到的输出'); grid;

x1 = [x zeros(1,14)];

y1 = filter(h,1,x1);

subplot(2,1,2);

stem(n,y1);

xlabel('时间序号 n'); ylabel('振幅'); title('由滤波生成的输出'); grid;

例2.8 线性时不变系统的稳定性

%程序P2_8

%基于冲激响应样本的绝对值

%的和的稳定性测试

clf;

num = [1 -0.8];den=[1 1.5 0.9];

N = 200;

h = impz(num,den,N+1)

parsum = 0;

for k =1:N+1;

parsum=parsum+abs(h(k));

if abs(h(k))<10^(-6),break,end end

数字信号处理MATLAB

%画出冲激响应

n=0:N;

stem(n,h)

xlabel('时间序号 n'); ylabel('振幅');

disp('值 =');disp(abs(h(k)));%显示abs(h(k))的值

2.30程序中为什么要使用命令for和end？

为了让程序达到循环的目的。