数字信号处理matlab版答案
【篇一：数字信号处理matlab实例】
txt>例1-1 用matlab计算序列{-2 0 1 –1 3}和序列{1 2 0 -1}的
离散卷积。
解 matlab程序如下：
a=[-2 0 1 -1 3];
b=[1 2 0 -1];
c=conv(a,b);
m=length(c)-1;
n=0:1:m;
stem(n,c);
xlabel(n); ylabel(幅度);
图1.1给出了卷积结果的图形，求得的结果存放在数组c中为：{-2
5 1 -3}。
例1-2 用matlab计算差分方程 -4 1 31
当输入序列为
解 matlab程序如下：时的输出结果。
脉冲响应。
n=41; a=[0.8 -0.44 0.36 0.22]; b=[1 0.7 -0.45 -0.6]; x=[1
zeros(1,n-1)]; k=0:1:n-1; y=filter(a,b,x); stem(k,y)
xlabel(n);ylabel(幅度) 1.2 给出了该差分方程的前41个样点的输出，即该系统的单位图
例1-3 用matlab计算例1-2差分方程
所对应的系统函数的dtft。
解例1-2差分方程所对应的系统函数为：
0.8?0.44z?1?0.36z?2?0.02z?3
h(z)?1?0.7z?1?0.45z?2?0.6z?3
其dtft为 0.8?0.44e?j??0.36e?j2??0.02e?j3?
)?1?0.7e?j??0.45e?j2??0.6e?j3? h(e
?j?
用matlab计算的程序如下：
k=256; num=[0.8 -0.44 0.36 0.02]; den=[1 0.7 -0.45 -0.6];
w=0:pi/k:pi; h=freqz(num,den,w); subplot(2,2,1);
plot(w/pi,real(h));grid title(实部) xlabel(\omega/\pi);ylabel(幅度) subplot(2,2,2); plot(w/pi,imag(h));grid
