第一章 数字信号处理基本概念
输入66输入20x(n-n0)012n输出344420x(n)012n输出3464642020012n输出34Y(n)y(n)012n346y(n-n0)420012n34
题1-10图(a) 时变性验证
b: 线性验证：令x1(n)??(n)?2?(n?1)?3?(n?2)?2?(n?3)，x2(n)?3?(n)
?2?(n?1)??(n?2)??(n?3)，a=1,b?2,p?2,q?1
程序如下： x1=[1 2 3 2];
x2=[3 2 1 1];
a=1;b=2;p=2;q=1;n=0:3; y1=a*x1+b;
y2=a*x2+b;
Y1=a*(x1*p+q*x2)+b;
Y2=p*y1+q*y2;
subplot(1,2,1);stem(n,Y1);
xlabel('n');ylabel('Y1(n)');axis([0,3,0,14]); subplot(1,2,2);stem(n,Y2); xlabel('n');ylabel('Y2(n)');
1-11
第一章 数字信号处理基本概念
14141212101088Y1(n)6Y2(n)64422001n23001n23
题1-10图(b) 线性性验证
1－11 已知系统函数H(z)?1?z?N，试用MATLAB画出该系统的幅频特性。
解： 利用MATLAB中的freqz()函数可以画出该系统的幅频特性曲线，如题1-11图所示。N
取10。 MATLAB程序如下：
N=10;
b=[1 zeros(1,N-1) 1]; a=[1 zeros(1,N)]; OMEGA=0:pi/150:2*pi; H=freqz(b,a,OMEGA); plot(OMEGA,abs(H));
题1-11图 幅频响应特性
1-12 一般的滑动平均由下列方程定义
1-12
第一章 数字信号处理基本概念
y(n)?1M1?M2?1k??M?x(n?k)?
1M2
1M1?M2?1[x(n?M1)?x(n?M1?1)???
[x(n)?x(n?1)???x(n?M2)]
该系统计算输出序列的第n个样本时是将其作为输入序列第n个样本前后的(M1?
M2?1)个样本的平均。
求：（1）该系统的冲激响应h(n)； （2）求该系统的频率响应； （3）对M1?0，M2?4，求H(e形。
解： （1）h(n)?1M1?M2j?)和argH(ej?)，并用MATLAB画出其图
?1k??M??(n?k)
1M21????M1?M2?1?0?1??（2）因为 h(n)??M1?M2?1?0?,?M1?n?M2,其他
,?M1?n?M2,其他
因此频率响应就是
H(ej?)?1M1?M?a2?1?M?e1M2?j?n
利用等比级数求和公式 可以得到：
H(ej??n?N1N2ak?aN1N2?11?a
)?1M1?M2e?1j?M1?e?j?(M2?1)?j??e?j?(M2?M1)/22sin[?(M1?M2?1)/2]1?eM1?M?1sin(?/2)
（3）当M1?0，M2?4时，
j?H(e)?1sin(5?/2)5sin(?/2)，argH(ej?)??2?
1-13
第一章 数字信号处理基本概念
利用MATLAB画出其频率响应图： 由 H(ej?)?1M1?M2?11M1?M2?1ej?M1?e?j?(M2?