资源描述:
《信号与系统上机实验(11级试用).doc》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、《信号与系统》上机实验实验指导(11级试用)说明:上机16学时,每次2学时,共8次。共7个实验项目,实验7为综合性实验项目,需4学时。上机时间:由任课老师、学习委员和实验室管理员共同决定。实验1信号变换与系统非时变性质的波形绘制1、实验目的掌握MATLAB画连续信号波形的方法,并能绘制信号变换和系统线性非时变性质的输出波形。2、预习阅读内容《信号与系统实验教程》实验1、实验2。《信号与系统分析》第1章。3、参考程序《信号与系统实验教程》第2章例2.1-2程序(如何用MATLAB画波形),例2.2-3程序(如何进行波形变换),例2.2-4程序
2、(如何进行波形分解)。《信号与系统分析》第1章计算机例题C1.1(如何用MATLAB画波形),计算机例题C1.2(如何进行波形变换),4、MATLAB学习要点l用MATLAB画波形的三要素(1)横轴赋值,两种方法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'
3、,2)lMATLAB的内部函数和外部函数(1)内部函数:MATLAB软件内部就存在,调用就可以了,本次介绍几个常用函数门函数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的三角脉冲函数,-1
4、为对称三角形;s=-1,三角形顶点左边。符号函数sign(),指数函数exp(),正弦函数sin(),余弦函数cos()(2)外部函数:根据编程的需要,自己编写的函数。介绍几个自编函数阶跃函数u(),文件名u.m,经常要用到,调用:y=u(t)functiony=u(t)y=(t>=0);画图函数myplot(),函数名为myplot.m,可以自动将波形放在图的中间,使其美观。Matlab程序为functionmyplot(x,y)%x为横坐标数组,y为纵坐标数组.x0=x(1);xe=x(end);max_y=max(y);min_y=m
5、in(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)l建立简单函数的两种方法如(1)自定义外部函数zdyf.m见例2.2-3程序functiony=zdyf(t)y=t.*(u(t)-u(t-1))-(u(t-1)-u(t-2));(2)定义在线函数,见例C1.2程序zdyf=inline('t.*(u(t)-u(t-1))-(u(t-1)-u(t-2))');5、上机内容l用MA
6、TLAB画出习题1-8的波形。l用MATLAB画出习题1-10的波形。6、思考题若已知的波形如图所示。试画出的波形。实验2微分方程的符号计算和波形绘制1、实验目的掌握MATLAB解微分方程的方法,并能绘制零输入响应、零状态响应和全响应的波形。2、预习阅读内容《信号与系统实验教程》实验4。《信号与系统分析》第2章。3、参考程序《信号与系统实验教程》第2章例2.4-3程序(如何用MATLAB的符号计算方法解微分方程并画波形),《信号与系统分析》第2章计算机例题C2.1(如何用MATLAB的符号计算方法解微分方程),计算机例题C2.2(如何用MA
7、TLAB的符号计算方法求零输入响应并画波形),4、MATLAB学习要点l微分方程的符号计算用Matlab的符号计算方法dsolve函数可以计算微分方程的解析式。其调用格式r=dsolve('eq1,eq2,...','cond1,cond2,...','v')r=dsolve('eq1','eq2',...,'cond1','cond2',...,'v')式中,eq1,eq2,...表示常系数微分方程(组)。cond1,cond2,...表示初始条件。V表示求解的变量。在缺省情况下,默认变量为t。参见参考程序,或《信号与系统实验教程》第40
8、页l符号量与数值量的转换函数subs()参见参考程序,或《信号与系统实验教程》第33页l在一幅图上画多条曲线,并标记图例。plot(t,y_n,t,yzi_n,'m:',t,yz