信号与系统上机试验11级试用.DOC

《信号与系统》上机实验 实验指导（11级试用） 说明：、MATLAB学习要点 t=0:0.1:30; t = 起点：间隔：终点，强调横轴的间隔。 t=linspace(0,30,300); t有300点，从起点0到终点30,等间隔自动计算，强调横轴的总点数。 （2）纵轴赋值 f=sin(pi*t.*sign(t)); 根据所给函数写出MATLAB表达式 （3）波形绘制 figure(1)，plot(t,f,linewidth,2) rectpuls（） 调用格式 y=rectpuls(t) 产生高度为1,宽度为1的门函数 调用格式 y=rectpuls(t,W) 产生高度为1,宽度为W的门函数 tripuls（） 调用格式 y=tripuls(t) 产生高度为1,宽度为1的三角脉冲函数 调用格式 y=tripuls(t,w) 产生高度为1,宽度为w的三角脉冲函数 调用格式 y=tripuls(t,w,s) 产生高度为1,宽度为w的三角脉冲函数, -1s1， s=0时为对称三角形; s=-1, 三角形顶点左边。 sign（）, 指数函数exp（）, 正弦函数sin（），余弦函数cos（） （2）外部函数：根据编程的需要，自己编写的函数。介绍几个自编函数 阶跃函数u( ), 文件名u.m，经常要用到，调用：y=u(t) function y=u(t) y=(t=0); 画图函数myplot( ),函数名为myplot.m，可以自动将波形放在图的中间，使其美观。 Matlab程序为 function myplot(x,y) % x为横坐标数组,y为纵坐标数组. x0=x(1);xe=x(end); max_y=max(y);min_y=min(y);dy=(max_y-min_y)/10; plot(x,y,linewidth,2);grid; axis([x0,xe,min_y-dy,max_y+dy]) set(gca,FontSize,8) （1）自定义外部函数 zdyf.m 见例2.2-3程序 function y=zdyf(t) y=t.*(u(t)-u(t-1))-(u(t-1)-u(t-2)); （）定义函数见例.2程序 zdyf=inline(t.*(u(t)-u(t-1))-(u(t-1)-u(t-2))); 5、 6、思考题 若已知的波形如图所示。试画出的波形。 实验2 微分方程的符号计算和波形绘制 1、实验目的 掌握MATLAB解微分方程的方法，并能绘制零输入响应、零状态响应和全响应的波形。 2、预习阅读内容 《信号与系统实验教程》实验4。《信号与系统分析》第2章。 3、参考程序 《信号与系统实验教程》第2章 例2.4-3程序（如何用MATLAB的符号计算方法解微分方程并画波形）， 《信号与系统分析》第2章 计算机例题C2.1（如何用MATLAB的符号计算方法解微分方程）， 计算机例题C2.2（如何用MATLAB的符号计算方法求零输入响应并画波形）， 4、MATLAB学习要点 dsolve函数可以计算微分方程的解析式。 其调用格式 r=dsolve(eq1,eq2,..., cond1,cond2,..., v) r=dsolve(eq1,eq2,...,cond1,cond2,...,v) 式中，eq1,eq2,...表示常系数微分方程（组）。cond1,cond2,...表示初始条件。V表示求解的变量。在缺省情况下，默认变量为t。 符号量与数值量的转换 函数subs（） 参见参考程序，或《信号与系统实验教程》第33页 在一幅图上画多条曲线，并标记图例。 plot(t,y_n,t,yzi_n,m:,t,yzs_n,r-.,linewidth,2) legend(全响应,零输入响应,零状态响应,0)参见参考程序，或《信号与系统实验教程》第20页 、 6、思考题 描述某线性时不变系统的方程为 已知初始条件为，输入，能否用MATLAB的符号计算方法计算系统全响应的解析解。 实验3 绘制连续系统的四幅图 1、实验目的 用MATLAB绘制连续系统的零极点图、冲激响应和阶跃响应图、频率响应和波特图。 四幅图为：零极点图，冲激响应和阶跃响应图（两个响应竖排），频率响应（幅频响应和相频响应竖排），波特图。 2、预习阅读内容 《信号与系统实验教程》实验9。《信号与系统分析》第3章。 3、参考程序 《信号与系统实验教程》实验9 例4.2-1程序（如何用MATLAB画零极点图、冲激响应和阶跃响应波形图）， 例4.2-2程序（如何用MATL