数字信号处理实验20151127.pdf

《数字信号处理》教材的有关matlab 代码 第二章 离散时间信号及产生（教材P17-19） 2 任意位置上的序列计算线性卷积函数（P25 ） 3 系统时域响应计算 6 DTFT （P82 ） 7 Z 逆变换（P84 ） 8 第三章 实现DFT 的函数（P138 ） 9 实现IDFT 的函数（P138 ） 9 循环移位函数（P142 ） 12 循环卷积函数circonvt （P143 ） 13 第四章 冲激响应不变法设计巴特沃斯数字滤波器（P182 ） 16 窗函数法设计FIR 数字滤波器（P233 ） 18 1 离散时间信号及产生 （教材P17-19） %单位取样序列的Matlab 函数 %n0是单位脉冲产生的位置，n1为显示时的下界， %n2为显示时的上界P17 function [x,n]=impseq(n0,n1,n2) % n0=3; % n1=1; % n2=5; n=[n1:n2]; x=[(n-n0)==0]; %单位阶跃序列的Matlab 函数 %n0信号由0变1的位置，n1显示的下界， %n2显示的上界P17 function [x,n] = stepseq(n0,n1,n2) % n0=1; % n1=-1; % n2=4; n=[n1:n2]; x=[(n-n0)= 0]; 例2.2 （1）p18 %example_2_2_1,p18 2 n=[-5 : 5]; x=3 * impseq(-3,-5,5)-impseq(3,-5,5); subplot(2,2,1) ; stem(n,x) ; xlabel(n); ylabel(xl(n)); 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; subplot(2,2,2); stem(n,x); xlabel(n); ylabel(x2(n)); 任意位置上的序列计算线性卷积函数 （P25 ） %self defined function conv_m, P25 function[y, ny] = conv_m(x, nx, h, nh) x=[1, 4, 3, 5, 1, 2, 3, 2]; nx=[-4:3]; h=[3, 2, 4, 1, 3, 2]; 3 nh=[-3:2]; nyb=nx(1) + nh(1) ; nye = nx(length(x)) +nh(length(h)) ; ny=[nyb : nye]; y=conv(x,h); 例2.8 P26 MATLAB 自带卷积函数-两序列均要求从n=0 开始定义 %example_2_8.m P26 n = 0 : 4; x=[1,1,1,1,1]; n = 0 : 6; h=(1.10).^n; [y] =conv(x,h); ny=length(y) ; n=0 : ny-1; stem(n,y) ; xlabel(n) ; ylabel(y(n)); 4 例2.9 P26 用自定义函数conv_m 计算任意 %P26 example_2_9.m x=[1,4,3,5,1,2,3,2]; nx=[-4:3]; h=[3,2,4,1,3,2]; nh=[-3:2]; [y, ny] = conv_m(x, nx, h, nh)%ny为y在横坐标轴上的位 置取值 %扩展x(n)和h(n)的长度 n=length(ny); x1 = zeros(1,n); h1=zeros(1,n); x1(find((ny=min(nx))(ny=max(nx))==1))=x; h1(find(ny=min(nh)(ny= max(nh))==1))=h; %画图 subplot(3,1,1); stem(n