南京邮电大学实验报告实验名称熟悉MATLAB环境快速傅里叶变换(FFT)及其应用IIR数字滤波器的设计FIR数字滤波器的设计课程名称数字信号处理A班级学号___14110___姓名_____________开课时间2015/2016学年，第二学期实验一熟悉MATLAB环境实验目的（1）熟悉MATLAB的主要操作命令。（2）学会简单的矩阵输入和数据读写。（3）掌握简单的绘图命令。（4）用MATLAB编程并学会创建函数。（5）观察离散系统的频率响应。实验内容(1)数组的加、减、乘、除和乘方运算。输入A=[1234]，B=[3,4,5,6]，求C=A+B，D=A-B，E=A.*B,F=A./B,G=A.^B。并用stem语句画出A、B、C、D、E、F、G。实验程序：n=0:1:3;A=[1234];subplot(3,1,1)stem(n,A)xlabel('n')ylabel('A')B=[3,4,5,6];subplot(3,1,2)stem(n,B)xlabel('n')ylabel('B')C=A+B;subplot(3,1,3)stem(n,C)xlabel('n')ylabel('C')figuresubplot(4,1,1)D=A-Bstem(n,D)xlabel('n')ylabel('D')subplot(4,1,2)E=A.*Bstem(n,E)xlabel('n')ylabel('E')subplot(4,1,3)F=A./Bstem(n,F)xlabel('n')ylabel('F')subplot(4,1,4)G=A.^Bstem(n,G)xlabel('n')ylabel('G')图像：(2)用MATLAB实现下列序列：a)实验程序：n=0:1:15;x1=0.8.^n;stem(n,x1)xlabel('n')ylabel('x(n)')title('2(a)')图像：b)实验程序：n=0:1:15;i=sqrt(-1);a=0.2+3*i;x2=exp(a*n);figuresubplot(1,2,1)stem(n,real(x2))xlabel('n')ylabel('x(n)实部')subplot(1,2,2)stem(n,imag(x2))xlabel('n')ylabel('x(n)虚部')图像：c)实验程序：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)实验程序：t=0:0.001:10;x=sin(2*pi*t);plot(t,x,'r-')xlabel('t'),ylabel('x(t)'),title('sin(2\pit)')图像：b)实验程序：t=0:0.001:4;x=cos(100*pi*t).*sin(pi*t);plot(t,x,'r-')xlabel('t'),ylabel('x(t)'),title('cos(100pit)sin(pit)')图像：(6)给定一因果系统，求出并绘制H(z)的幅频响应和相频响应。实验程序：num=[1sqrt(2)1];den=[1-0.670.9];w=0:pi/512:pi;h=freqz(num,den,w);subplot(2,2,1);plot(w/pi,real(h));gridtitle('实部')xlabel('\omega/\pi');ylabel('幅度')subplot(2,2,2);plot(w/pi,imag(h));gridtitle('虚部')xlabel('\omega/\pi');ylabel('幅度')subplot(2,2,3);plot(w/pi,abs(h));gridtitle('幅度谱')xlabel('\omega/\pi');ylabel('幅度')subplot(2,2,4);plot(w/pi,angle(h));gridtitle('相位谱')xlabel('\omega/\pi');ylabel('相位（rad）');图像：(7)计算序列{8-2-123}和序列{23-1-3}的离散卷积，并作图表示卷积结果。实验程序：%输入x（n）及其下标x=