DSP实验完整版.doc

实验名称： 熟悉MATLAB环境 快速傅里叶变换(FFT)及其应用 IIR数字滤波器的设计 FIR数字滤波器的设计 实验报告内容包括： 实验一 实验名称（注意：9题选4~5题） 实验目的 （1）熟悉MATLAB的主要操作命令。 （2）学会简单的矩阵输入和数据读写。 （3）掌握简单的绘图命令。 （4）用MATLAB编程并学会创建函数。 （5）观察离散系统的频率响应。 实验内容 (2) 用MATLAB实现下列序列： ######a) n=0:1:15; x1=1.2.^n; a=(0.2+3*i)*n; stem(x1) ############b) n=0:1:15; x2=exp(a); a=(0.3+2*i)*n; stem(x2) ################b) t=0:0.01:4; x=cos(200*pi*t).*sin(pi*t); plot(t,x, r-); xlabel(t),ylabel(x(t)),title(cos) ##################(8) 求以下差分方程所描述系统的单位脉冲响应h(n), 实验过程与结果（含实验程序、运行的数据结果和图形）； clear all; N=50; a=[2 -1]; b=[2 0.3 -0.04]; x1=[1 zeros(1,N-1)]; n=0:1:N-1; h=filter(a,b,x1); stem(n,h) axis([-1 53 -2.5 1.2]) 实验二 实验名称（注：共四个，也可以自己在老师给的“实验二的内容”中自己选做） 实验目的 (1) 在理论学习的基础上，通过本实验，加深对FFT的理解，熟悉MATLAB中的有关函数。 (2) 应用FFT对典型信号进行频谱分析。 (3) 了解应用FFT进行信号频谱分析过程中可能出现的问题，以便在实际中正确应用FFT。 (4) 应用FFT实现序列的线性卷积和相关。 实验内容 1、P131修改 x=[2 2 -2 3 1 2 1 2];y=[1 3 -3 2 3 1 -2 2];k=length(x);xk=fft(x,2*k);yk=fft(y,2*k);rm=real(ifft(conj(xk).*yk));rm=[rm(k+2:2*k) rm(1:k)];m=(-k+1):(k-1);stem(m,rm);xlable(m);ylable(幅值); 2、例3。12修改 x=[1 2 1 3 4 -2 -3 0 3 4 4 1 5 6 7 1 0 4];y=[1 2 2 2 4 5 1 1 4 -3 -1 0 4 1 6 4 5 7 6 1 2 2]; k=length(y);e=rand(1,k)-0.5;y=y+e;xk=fft(x,2*k);yk=fft(y,2*k);rm=real(ifft(conj(xk).*yk));rm=[rm(k+2:2*k) rm(1:k)];m=(-k+1):(k-1);stem(m,rm);xlable(m);ylable(幅值); 三角波 xc=[1,3,2,4,3,1,2,6]; n=0:1:7; stem(n,xc) Xc=fft(xc,8); %%改变fft点数 k=0:1:7; plot(k,abs(Xc)) %%%%%%%%% for n=1:4 xc(n)=n-1; end; for n=5:8; xc(n)=8-(n-1); end; xc 画出的DTFT频谱（） 先建高斯序列 建立一个M文件，作用是产生高斯序列，并绘出时域波形和幅频特性。 M文件代码如下： function [xa, N]= gauss(p,q) n=0:1:15; N=16; xa=exp(-(n-p).^2/q); 再输入下面程序 [xa,N]=gauss(10,6); M=10000; w=2*pi/M*(0:1:M-1); Xa=zeros(1,M); for k=1:M Xa(k)=sum(xa*(exp(-j*w(k)*(0:N-1)))); end n=0:1:N-1; subplot(2,1,1); stem(n,xa); subplot(2,1,2); plot(w,abs(Xa)) 实验过程与结果（含实验程序、运行的数据结果和图形）； 实验三 实验名称 实验目的 (1) 掌握双线性变换法及脉冲响应不变法设计IIR数字滤波器的具体设计方法及其原理，熟悉用双线性变换法及脉冲响应不变法设计低通、高通和带通IIR数字滤波器的计算机编程。 (2) 观察双线性变换及脉冲响应不变法设计的滤波器的频域特性，了解双线性变换法及脉冲响应不变法的特点。 (3) 熟悉巴特沃思滤波器、切比雪夫滤波器和椭圆滤波器