数字信处理实验指导书.docVIP

数字信号处理实验指导书 Matlab在数字信号处理中的应用 聂小燕 编著 电子科技大学成都学院 2012年3月 数字信号处理上机实验 本课程上机实验8课时，共4个实验。 学时 实验内容 备注 2 实验1 离散时间信号与系统 2 实验2 FFT及其应用 2 实验3 IIR滤波器的设计 2 实验4 FIR滤波器的设计 实验一 离散时间信号与系统 【实验目的】 熟悉用matlab分析离散时间信号与系统 【实验内容】 1．用MATLAB实现函数impseq(n0,n1,n2)，使函数实现。 （1）该函数的MATLAB程序为 function [x,n]=impseq(n0,n1,n2) n=[n1:n2]; x =[(n-n0)==0]; （2）函数调用举例 % ＝2－, -5=n=5 n = [-5:5]; x = 2*impseq(-2,-5,5)-impseq(4,-5,5); stem(n,x); title(调用自定义的单位脉冲序列函数) xlabel(n); ylabel(x(n)); axis([-5,5,-2,3]) 2．用MATLAB实现函数stepseq(n0,n1,n2)，使函数实现u(n-n0)， 。 （1）该函数的MATLAB程序为： function [x,n]=stepseq(n0,n1,n2) n=[n1:n2]; x=[(n-n0)=0]; （2）函数调用举例 % x(n) = n[u(n)-u(n-10)]+10*exp(-0.3(n-10))(u(n-10)-u(n-20)); 0=n=20 n = [0:20]; x1 = n.*(stepseq(0,0,20)-stepseq(10,0,20)); x2 = 10*exp(-0.3*(n-10)).*(stepseq(10,0,20)-stepseq(20,0,20)); x = x1+x2; stem(n,x); title(调用单位脉冲序列实例) xlabel(n);ylabel(x(n));axis([0,20,-1,11]) 3．用MATLAB实现常见的离散时间序列。 （1） 参考程序：n=[0:10]; x=(0.8).^n; stem(n,x,’ro’) （2）x(n) = exp((2+3j)*n) n=[-10:10]; x=exp((2+3j)*n); stem(n,x,ro) title(复指数序列) xlabel(n);ylabel(x(n)); （3）x(n) = cos(0.04*pi*n) + 0.2*w(n) n = [0:50]; x = cos(0.04*pi*n)+0.2*randn(size(n)); stem(n,x);title(正弦序列＋随机序列) xlabel(n);ylabel(x(n));axis([0,50,-1.4,1.4]) （4）x(n) = {...,5,4,3,2,1,5,4,3,2,1,...}; -10=n=9 序列的周期化 n=[-10:9]; x=[5,4,3,2,1]; xtilde=x * ones(1,4); xtilde=(xtilde(:)); stem(n,xtilde);title(有限长序列的周期化) xlabel(n);ylabel(xtilde(n));axis([-10,9,-1,6]) 4．计算两个序列的卷积。教材P39页（7） a=[8 -2 -1 2 3] b=[2 3 -1 -3] c=conv(a,b) M=length(c)-1 n=0:1:M stem(n,c) xlabel(n) ylabel(幅度); 5．计算幅频响应与相频响应.教材P39页（6） b = [1 -sqrt(2) 1]; % H(Z)分子多项式系数 a = [1 -0.67 0.9]; % H(Z)分母多项式系数 [h,w] = freqz(b,a); % 计算频率响应 am = 20*log10(abs(h)); % 计算幅频响应并绘图 subplot(2,1,1); plot(w,am); ph = angle(h); % 计算相频响应并绘图 subplot(2,1,2); plot(w,ph); 6．计算离散系统的单位脉冲响应.教材P39页（8） N=50; 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(幅度) 实验二 快速傅立叶变换FFT及其应用 【实验目的】 熟悉matlab中的有关函数 应用FFT对典型信号进行频谱分析，并了解可能出现的问题，以便在实际中正确应用FFT 应用FF