离散信号matlab.docVIP

《数字信号处理》课程实验报告 实验一 离散时间信号与系统 一、实验目的 掌握运用MATLAB表示常用离散时间信号 掌握运用MATLAB实现离散信号的基本运算 掌握运用MATLAB求解离散时间系统的零状态响应 掌握运用MATLAB求解离散时间系统的单位取样响应 掌握运用MATLAB求解离散时间系统的卷积和 二、实验内容与步骤 离散时间信号在MATLAB中的表示 单位抽样序列 通过编写impDTS.m文件来实现单位抽样序列的生成函数。 利用MATLAB的impDTS函数绘出单位抽样序列的波形图。 代码： function y=impDTS(n) %%%脉冲 y=(n==0); 单位阶跃序列 通过编写uDTS.m文件来实现单位抽样序列的生成函数。 利用MATLAB的uDTS函数绘出单位抽样序列的波形图。 代码： function y=uDST(n) y=(n=0); %======阶跃 正弦序列 用MATLAB命令绘制下方正弦序列的波形图。 验证正弦序列的周期性。 复指数序列 用MATLAB命令画出复指数序列的实部、虚部、模及相角随时间变化的曲线，并观察其时域特性。 代码main: %%%=======yinxiaogen %%%=======2013.11.22 %%%======= clear all close all clc n=-5:0.1:5; y=impDTS(n); %figure subplot(311); stem(n,y); title(抽样函数); xlabel(n); ylabel(y); x=uDST(n); subplot(312); stem(n,x); title(阶跃函数); xlabel(n); ylabel(x); z=sin(n*pi/6); subplot(313); stem(n,z); title(正弦函数); xlabel(n); ylabel(z); 离散时间信号的基本运算 用MATLAB命令画出下列离散时间信号的波形图。 代码： %%%%%%========yinxiaogen %%%%%%========2013.11.24 clear all close all clc a=2; N=4; n=-10:0.1:10; x=heaviside(n); x1=(a.^n).*[heaviside(n)-heaviside(n-N)] ; %%u(n)表示为heaviside() x2=(a.^(n-2)).*[heaviside(n-2)-heaviside(n-N-2)] ;%.支持向量运算 x3=-x1; figure (1); subplot(221);stem(n,x);title(阶跃函数); subplot(222);stem(n,x1);title(x1); subplot(223);stem(n,x2);title(x2); subplot(224);stem(n,x3);title(x3); 3、离散时间系统的响应 已知某LSI系统的差分方程为 试用MATLAB命令绘出当激励信号为时，该系统的零状态响应。 离散时间系统的单位取样响应 上例中系统的单位抽样响应。 代码： %%%=======yinxiaogen %%%=======2013.11.25 clear all close all clc n=-1:0.1:10; a=[3,-4,2]; b=[1,2]; x=((1/2).^n).*(n=0); %heaviside(n); x1=impDTS(n); %冲激函数 y=filter(b,a,x); y1=filter(b,a,x1); %冲激函数响应 figure (1); stem(n,y); title(零状态响应); xlabel(n); ylabel(y); figure(2); stem(n,y1); title( 单位取样响应); xlabel(n); ylabel(y1); 离散时间信号的卷积和运算 已知某系统的单位取样响应为 试用MATLAB求当激励信号为时，系统的零状态响应。 代码： %%%=======yinxiaogen %%%=======2013.11.26 %%%======= clear all close all clc p=0.1; %抽样间隔 n1=-5:0.1:10; h=(0.8.^n1).*[(n1=0)(n18)]; n2=-10:0.1:5; x=(n2=0)(n24); y=conv(x,h); ny1=n1+n2; %n1, n2的长度应相等 ny=ny1:p:ny1+(len