实验二信号与系统时域分析2.pptVIP

实验二 信号与系统时域分析 Signals and Systems 实验目的 了解信号表示的Matlab函数、多项式运算的相关函数 理解一般微分方程的求解方法、连续积分的数值计算方法 掌握用Matlab表示信号、进行信号运算 掌握用Matlab求解连续系统的零状态响应、冲激响应和阶跃响应的方法 一、基本信号的MATLAB表示 二、信号基本运算的MATLAB实现 已知三角波x(t)，用MATLAB画出x(2t)和x(2-2t) 波形 三角波f(t)微分与积分的波形 实验原理 三、系统微分方程的描述 四、零输入响应，零状态响应的求解 五、连续系统冲激响应和阶跃响应求解 六、卷积的计算 实验内容 一、计算并绘制由如下微分方程表示的系统在输入信号为x(t) = (e-2t - e-3t)u(t)时，零状态响应曲线。 实验内容 二、计算并绘制系统 的冲激响应与阶跃响应曲线。 实验内容 三、计算并用MATLAB实现下列信号的卷积，要求绘出信号及卷积图形。 小 结 Y = diff(X) calculates differences between adjacent elements of X. y=diff(f)/h?d(x)/dt，h取很小?dt趋于0 y = rectpuls(t,width) (t-2T):中心平移(0-2) y = rectpuls(t) returns a continuous, aperiodic, unity-height rectangular pulse at the sample times indicated in array t, centered about t = 0 and with a default width of 1. Note that the interval of non-zero amplitude is defined to be open on the right, that is, rectpuls(-0.5) = 1 while rectpuls(0.5) = 0. y = rectpuls(t,w) generates a rectangle of width w. y = tripuls(t,width,skew) t,w与rectpuls类似，skew?(倾斜角度or程度度量）:2倍的三角波顶点坐标tmax与三角波宽度之比，即skew=2 tmax/width. y = tripuls(T) returns a continuous, aperiodic, symmetric, unity-height triangular pulse at the times indicated in array T, centered about T=0 and with a default width of 1. y = tripuls(T,w) generates a triangular pulse of width w. y = tripuls(T,w,s) generates a triangular pulse with skew s, where -1 s 1. When s is 0, a symmetric triangular pulse is generated. 自编单位脉冲序列impseq函数，利用 关系运算“等于”来实现（k=k0时返回“真”值1）。 同前，自编stepseq函数，利用关系运算“大于等于”来实现（k=k0时返回“真”值1）。 前： t=-3:0.001:3; ft=tripuls(t,4,0.5); plot(t,ft) function yt=f2_2(t) yt=tripuls(t,4,0.5); y2(n)=quad(f2_2, -3,t(n)); 这里就是，t每取一个值，就算一次积分(面积），作为y2的元素，最后构成所需序列 连续信号的微分与积分都是近似实现（序列-plot）,t是对模拟精度的控制 引出 卷积 Yzs(t)=x(t)*h(t) Help conv 利用conv函数（序列卷积和）实现连续时间卷积积分。 冲激响应：初始状态为零，输入激励——单位冲激信号，单独作用而产生的输出响应。（即 单位冲激信号 作用于系统的 零状态响应） 实验书p35,36讲了利用lsim函数仿真冲激、阶跃响应。 y=impulse(sys,t) 冲激响应只取决于系统，而非x(t). * * * * * * 第三军医大学生物医学工程与医学影像学院 电子学与电子工程学教研室 指数信号Aeat y = A*exp(a*t); 指数序列ak 幂运算a.^k实现 正弦型信号 内部