信号描述与信号变换中MATLAB应用.docVIP

PAGE PAGE 11 实验一 信号描述与信号变换中的MATLAB应用 实验目的 熟悉MATLAB软件的使用； 用MATLAB语句实现信号描述与信号变换。 实验原理 信号描述（连续时间信号） 向量表示法 例：，可采用下列两个向量表示： 若要绘出该连续时间信号的波形，可加入下列语句。 plot(t,f) grid title(‘f(t)=sin(t)’) xlable(‘t’) axis(-10,10,-1.1,1.1) 符号运算表示法 例：，采用符号表示式表述为 syms t 或 f=sym(‘sin(pi/4*t)’) f=’sin(pi/4*t)’ f= sin(pi/4*t) ezplot(f,[-16,16]) 注意：画图命令的不同， 画出离散时间序列的波形用stem； 画出连续时间信号的波形用plot，如用符号运算表示法表示画图命令是ezplot。 信号运算与变换(连续时间信号) 信号的相加、相乘，信号的导数、积分 例：对于，若，则采用下列语句实现： syms t f1=sin(t) f2=t f=f1+f2 ezplot(f) 例：对于，若，则采用下列语句实现： f=f1*f2 例：对于，若，则采用下列语句实现： syms t f=sin(t) df=diff(f) ezplot(df) 例：对于，若，则采用下列语句实现： s=int(f) ezplot(s) 信号时移、折叠、尺度变换 时移，可用下列语句实现： y=subs(f,t,t-t0) 折叠，可用下列语句实现： y=subs(f,t,-t) 尺度变换，可用下列语句实现： y=subs(f,t,at) 三．实验所用参考程序 例：设信号，试用MATLAB求解，，，画出相应的波形。 如用符号运算表示法，用ezplot实现阶跃函数的可视化时，首先需要在自己的工作目录work下创建Heaviside的M文件（实现阶跃函数，注意函数符号的大小写），该文件如下： Function f=Heaviside(t) f=(t0); 正确定义出该函数并保存为Heaviside.m文件后，就可以调用该函数了。 用符号运算表示法，用ezplot实现函数的可视化，所用MATLAB命令如下： syms t f=sym((1+t/2)*(Heaviside(t+2)-Heaviside(t-2))) subplot(2,2,1) ezplot(f,[-3,3]) y1=subs(f,t,t+2) subplot(2,2,2) ezplot(y1,[-5,1]) y2=subs(f,t,-t) subplot(2,2,3) ezplot(y2,[-3,3]) y3=subs(f,t,2*t) subplot(2,2,4) ezplot(y3,[-2,2]) 或： 用向量表示法，用plot实现函数的可视化，所用MATLAB命令如下： p=0.1; t= -2:p:2; f=1+t/2; subplot(2,2,1),plot(t,f) grid ylabel(f(t)) subplot(2,2,2),plot(t+2,f) grid ylabel(f(t+2)) subplot(2,2,3),plot(-t,f) grid ylabel(f(-t)) subplot(2,2,4),plot(2.*t,f) grid ylabel(f(2t)) 四．实验要求 上机实现下列信号的分析 1．已知信号，画出信号 ，，，的波形。 2．已知，，分别画出下列信号波形：（1）；（2）。 实验二 连续时间信号卷积运算的MATLAB实现 实验目的 掌握卷积的概念及计算方法； 熟悉通过调用conv()函数求解连续时间信号卷积的数值分析方法。 实验原理 在线性系统的时域分析方法中，卷积是个极其重要的概念，占有重要地位，只要知道了系统在单位冲激信号作用下的零状态响应即系统的单位冲激响应，就可以利用卷积积分求出系统在任何激励作用下的零状态响应。 连续信号和的卷积定义为： 卷积计算可以通过信号分段求和来实现，即 数值计算时近似为： 记作： 式中、和分别为对以、和为时间间隔进行采样所得的离散序列。 用MATLAB实现连续信号和卷积的过程如下： 将连续信号和以等时间间隔Δ进行取样，得到离散序列和和； 构造与和相对应的时间向量k1