数字信号DSP实验课件.ppt

实验一 信号、系统及系统响应;理想采样信号 的频谱; 单位圆上序列的Z变换为序列的傅里叶变换（DTFT)，也称为数字序列的频谱。; 序列傅里叶变换是频率w的周期函数，周期为2π。所以通常对 在 上进行M点采样来观察分析。对长度为N的有限长序列;一个时域离散LTI系统的输入/输出关系为;Matlab基本知识;Matlab基本知识;Matlab基本知识;Matlab基本知识;Matlab基本知识;Matlab基本知识;Matlab基本知识;Matlab基本知识;单位阶跃序列的生成函数，stepseq函数如下： function [x,n]=stepseq(n0,n1,n2); %使函数实现u(n-n0)，n1=nn2 n=[n1:n2]; x=[(n-n0)=0]; 保存上面所输入的m文件，命名为stepseq.m Xc=stepseq(1,1,5); stem(n,x,.);%“stem”绘制针状图最简单 xlabel(n); ylabel(x_c(n)); title(initial signal); text((0.3*N),(0.8*t),str);%在图形指定的位置上显示字符串;信号处理工具箱;信号处理工具箱;;2.滤波器分析函数 1）conv 求两个序列的卷积 clear;clf; x(1:5)=1;%注意：MATLAB 中数组下标从1 开始 tx=-2:2; th=0:6; h(1:7)=0.5; ty=-2:8; y=conv(x,h); figure(1)%建立第1幅图形 subplot(2,2,1)%一幅图显示4个图形 stem(tx,x); subplot(2,2,2) stem(th,h); axis([0,6,0,1]); subplot(2,1,2) stem(ty,y); 查看结果;利用这些方法，我们可以产生任意我们需要的离散信号，这些产生离散信号的命令，可以用函数的方式给出，也可以以脚本文件的方式给出，或者直接在Matlab主命令窗口中给出也可以。;示例一 求 参考程序： n=-10:10; x=(-0.9).^n; %实指数序列 k=-200:200;w=(pi/100)*k;%0-2π上去200个采样点 Xw=x*exp(-j*pi/100*n*k);%复指数序列 magXw=abs(Xw); %绘制x(n)的幅度谱 angXw=angle(Xw); %绘制x(n)的相位谱 subplot(2,1,1) plot(w,magXw);grid xlabel(frequency in units of pi);ylabel(|X|); subplot(2,1,2) plot(w,angXw);grid xlabel(frequency in units of pi);ylabel(angle/pi) ;如果需要分析序列的频谱，则要根据序列傅立叶变换公式进行，不过需要注意的是，在DTFT中，频率分量w是连续的，而在计算机实现过程中，需要对其作离散化。