信号语系统实验指导书.docVIP

《信号与系统》实验指导书 电子与信息工程教研室 序 论 Mat lab是Math works公司于1984年正式推出的一套集数值计算、符号运算、图形处理等于一身的科学计算语言。目前，国内绝大多数院校的电子与信息课程的4门主要课程：电路、信号与系统、数字信号处理、自动控制原理的课堂教学、实验均与之关联。其优点在于：克服了物理器件的自身制约；克服了实验室规模的限制等。在我院电子与信息工程专业的教学及实践环节中，我们引入Mat lab这个语言工具，可以使信息课程的教学更加直观，可以加深同学门对课本中的理论的深层次理解。 同时，同学们在实验过程中，可以掌握一门新的工具、新的语言、新的思考问题解决问题的方法。 实验一、RC电路的响应分析 实验的内容 正弦激励的一阶电路（如图1所示）。已知R＝2Ω，UC(0+)=4V,激励的正弦电压US(t)=Umcos(ωt)，其中Um＝10V，ω＝2rad/s。当t＝0时，开关S闭合，求电容电压的全响应，区分其暂态响应和稳态响应，并画出波形图。 建模 电容电压的微分方程为： 令T＝RC，利用电路的三要素法，其解为： UC(t)=UCP(t)+[UC(0+)-UCP(0+)]e-t/T,(t≥0) 其中：UC(0+)为电容的初始电压； UCP(t)为方程的特解；当激励为US(t)=Umcos(ωt)，则特解UCP(t)＝Ucmcos(ωt＋φ)，其中：Ucm＝；φ＝900－arctan(ωCR)-1. UCP(0+)= Ucmcos(φ) 最后得出的电容电压的全响应为：UC(t)= Ucmcos(ωt＋φ)+[UC(0+)-UCP(0+)]e-t/T 暂态响应（固态响应）：Uuctr(t)= [UC(0+)-UCP(0+)]e-t/T,(t≥0) 稳态相应（强迫响应）：Uucsr(t)= Ucmcos(ωt＋φ) Mat lab程序 r=2;c=0.5;T=r*c; %输入参数 uc0=4; um=10; %给定数据 w=2; zc=1/j/w/c; t=0:0.1:10; %设定时间数组 us=um*cos(w*t); %激励信号 ucst=us*zc/(r+zc); %稳态计算 ucp0=ucst(1); %稳态初值 uctr=(uc0-ucp0)*exp(-t/T); % 暂态 uc=uctr+ucst; %全响应 plot(t,uc,-,t,uctr,:,t,ucst,.), %绘图 legend(nc,uctr,ucst) %图例标注 结果： 图2 电容电压波形 2 连续LTI系统的零输入响应 2.1实验内容 n阶线性时不变系统的微分方程为： 式中：n≥m。已知y及其各阶导数的初始值y(0)、y(1)(0)、y(2)(0)、y(3)(0)、…、y(n－1)(0)，求零输入响应 2.2 Mat lab环境下，编程实现。 2.3 三、分析解决问题 建模 零输入响应为微分方程的奇次解。其解的形式为（设特征根均为单根）： y(t)=C1ep1t+ C2ep2t+…+ Cnepnt 特征根可以用Mat lab中的roots(a)语句求得。各系数Ci由y的初值确定。因此有： C1＋C2＋…＋Cn＝y(0) C1 P1＋C2 P2＋…＋Cn Pn＝Dy0= y(1)(0) ┆ ┆ 写成矩阵形式为： 即：VC＝Y0；C＝V\Y0； V为范德蒙矩阵，在Mat lab中有vander函数调用。 Mat lab程序设计： a=input(请输入微分方程的系数 a＝[a1 a2 ...]=); n=length(a)-1; %计算微分方程的阶数 y0=input(请输入初始条件 y=[y0,dy0 d2y0,...]=); p=roots(a); v=rot90(vander(p)); c=v\y0; dt=input(时间步长 dt=); tf=input(时间终值 tf=); t=0:dt:tf; y=zeros(1,length(t)); for k=1:n y=y+c(k)*exp(p(k)*t); end plot(t,y) grid 程序运行结果： 请输入微分方程的系数 a＝[a1 a2 ...]=[3,5,7,1]; 请输入初始