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('幅度');
