南京邮电大学DSP实验报告

c) x(n)?3cos(0.125?n?0.2?)?2sin(0.25?n?0.1?) 0?n?15 解： n=0:1:15;

x3=3*cos(0.125*pi*n+0.2*pi) + 2*sin(0.25*pi*n+0.1*pi); stem(n,x3) xlabel('n') ylabel('x(n)')

(4) 绘出下列时间函数的图形，对x轴、y轴以及图形上方均须加上适当的标注： a) x(t)?sin(2?t) 0?t?10s 解：

t=0:0.001:10; x=sin(2*pi*t); plot(t,x,'r-')

xlabel('t'),ylabel('x(t)'),title('sin(2\\pit)')

b) x(t)?cos(100?t)sin(?t) 0?t?4s 解：

t=0:0.001:4;

x=cos(100*pi*t).*sin(pi*t); plot(t,x,'r-')

xlabel('t'),ylabel('x(t)'),title('cos(100\\pi*t).*sin(\\pi*t)')

(6) 给定一因果系统H(z)?(1?2z?1?z?2)/(1?0.67z?1?0.9z?2)，求出并绘制H(z)的幅

频响应和相频响应。 解： k=256;

num=[1 1.414 1]; den=[1 -0.67 0.9]; 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 title(' 虚部')

xlabel('\\omega/\\pi');ylabel('幅度') subplot(2,2,3);

plot(w/pi,abs(h));grid title(' 幅度谱')

xlabel('\\omega/\\pi');ylabel('幅度') subplot(2,2,4);

plot(w/pi,angle(h));grid title(' 相位谱')

xlabel('\\omega/\\pi');ylabel('幅度')

(7) 计算序列{8 -2 -1 2 3}和序列{2 3 -1 -3}的离散卷积，并作图表示卷积结果。 解：

x=[8,-2,-1,2,3]; kx=0:4;

h=[2,3,-1,-3]; kh=0:3;

y=conv(x,h);

k=kx(1)+kh(1):kx(end)+kh(end); stem(k,y);

xlabel('n');ylabel('y(n)');

(8) 求以下差分方程所描述系统的单位脉冲响应h(n), 0?n?50

y(n)?0.1y(n?1)?0.06y(n?2)?x(n)?2x(n?1)

解： N=51; a=[1 -2];

b=[1 0.1 -0.06]; x=[1 zeros(1,N-1)]; k=0:1:N-1; y=filter(a,b,x); stem(k,y); xlabel('n'); ylabel('幅度');