[工学]数字信号处理指导书matlab版.docVIP

实验1 时域离散信号的产生 一、实验目的 MATLAB源程序为 1) function [x,n] = impuls (n0,n1,n2) % Generates x(n) = delta(n-n0); n=n0 处建立一个单位抽样序列 % [x,n] = impuls (n0,n1,n2) if ((n0 n1) | (n0 n2) | (n1 n2)) error(arguments must satisfy n1 = n0 = n2) end n = [n1:n2]; x = [zeros(1,(n0-n1)), 1, zeros(1,(n2-n0))]; 将上述文件存为：impuls.m，在命令窗口输入 n0=0,n1=-10,n2=11; [x,n]=impuls (n0,n1,n2); stem(n,x,’filled’) 2) n1=-5;n2=5;n0=0; n=n1:n2; x=[n==n0]; stem(n,x,filled,k); axis([n1,n2,1.1*min(x),1.1*max(x)]); title(单位脉冲序列); xlabel(时间(n)); ylabel(幅度x(n)); 3) n1=-5;n2=5;k=0; n=n1:n2; nt=length(n); %求n点的个数 nk=abs(k-n1)+1; %确定k在n序列中的位置 x=zeros(1,nt); %对所有样点置0 x(nk)=1; %对抽样点置1 stem(n,x,filled,k); axis([n1,n2,0,1.1*max(x)]); title(单位脉冲序列); xlabel(时间(n)); Ylabel(幅度x(n)); 2、单位阶跃序列 MATLAB源程序为： 1) n1=-2;n2=8;n0=0; n=n1:n2; %生成离散信号的时间序列 x=[n=n0]; %生成离散信号x(n) stem(n, x,filled,k); %绘制脉冲杆图，且圆点处用实芯圆表示 axis([n1,n2,0,1.1*max(x)]); title(单位阶跃序列); xlabel(时间(n)); Ylabel(幅度x(n)); 2) n1=-2;n2=8;k=0; n=n1:n2; nt=length(n); %求n点的个数 nk=abs(k-n1)+1; %确定k在n序列中的位置 x=[zeros(1,nk-1),ones(1,nt-nk+1)]; %对所有样点置0 stem(n,x,filled,k); axis([n1,n2,0,1.1*max(x)]); title(单位阶跃序列); xlabel(时间(n)); ylabel(幅度x(n)); 3、正弦序列 x（n）=Um 例、已知一时域周期性正弦信号的频率为1HZ，振幅幅度为1V，在窗口中显示两个周期的信号波形，并对该信号的一个周期进行32点采样获得离散信号。显示连续信号和采样获得离散信号波形。 MATLAB源程序为： f=1;Um=1;nt=2;; %输入信号频率、振幅和显示周期数 N=32; T=1/f; %N为采样点数，T为窗口显示时间 dt=T/N; %采样时间间隔 n=0:nt*N-1; %离散信号的时间序列 tn=n*dt; %时间序列样点在时间轴上的位置 x=Um*sin(2*f*pi*tn); subplot(2,1,1);plot(tn,x); %显示原信号 axis([0,nt*T,1.1*min(x) 1.1*max(x)]); ylabel(x(t)); subplot(2,1,2);stem(tn,x); %显示经采样的信号 axis([0,nt*T,1.1*min(x),1.1*max(x)]); ylabel(x(n)); 4、矩形序列 将square表示式中的t换成n，且n取整数，则可以获得矩形序列。 例、一个连续的周期性矩形信号频率为5kHZ，信号幅度为0-2V之间，脉冲宽