熟悉MATLAB环境,快速傅里叶变换(FFT)及其应用,IIR数字滤波器的设计,FIR数字滤波器的设计实验报告

b) x(n)?e(0.2?3j)n 0?n?15

输入代码如下： a=(0.2+3*i)*n; x2=exp(a) stem(x2) 运行图片：

c) x(n)?3cos(0.125?n?0.2?)?2sin(0.25?n?0.1?) 0?n?15 输入代码如下：

x3=3*cos(0.125*pi*n+0.2*pi)+2*sin(0.25*pi*n+0.1*pi) stem(x3) 运行图片：

d) 将c)中的x(n)扩展为以16为周期的函数x16(n)?x(n?16)，绘出四个周期。 输入代码如下：

x3=3*cos(0.125*pi*n+0.2*pi)+2*sin(0.25*pi*n+0.1*pi) x4=[x3 x3 x3 x3]; stem(x4) 运行图片：

(4) 绘出下列时间函数的图形，对x轴、y轴以及图形上方均须加上适当的标注： a) x(t)?sin(2?t) 0?t?10s 输入代码如下： t1=0:0.001:10; x1=sin(2*pi*t1); subplot(2,1,1)

plot(t1,x1,'r-')

xlabel('t'),ylabel('x(t)'),title('sin') b) x(t)?cos(100?t)sin(?t) 0?t?4s 输入代码如下： t2=0:0.001:4;

x2=cos(100*pi*t2).*sin(pi*t2); subplot(2,1,2) plot(t2,x2)

xlabel('t'),ylabel('x(t)'),title('cossin') 运行图片：

(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,abs(h));grid; title('幅度谱');

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

plot(w/pi,angle(h));grid

title('相位谱');

xlabel('\\omega/ \\pi');ylabel('幅度'); 运行图片： 幅度谱252015幅度相位谱10-1幅度1050-2-3-400.20.40.60.8100.20.40.60.81?/ ??/ ? (7) 计算序列{8 -2 -1 2 3}和序列{2 3 -1 -3}的离散卷积，并作图表示卷积结果。 输入代码如下： x=[8,-2,-1,2,3]; kx=-2;2; h=[2,3,-1,-3]; kh=-2;1; y=conv(x,h); k=kx(1)+kh(1):kx(end)+kh(end); stem(k,y); xlabel('n');ylabel('y(n)'); 运行图片： 20151050y(n)-5-10-15-20-25-3-2-10n1234 (8) 求以下差分方程所描述系统的单位脉冲响应h(n), 0?n?50 y(n)?0.1y(n?1)?0.06y(n?2)?x(n)?2x(n?1)

输入代码如下：