MATLAB在电类专业课程中的应用第8章.pptxVIP

第八章 MATLAB在信号与系统中的应用;8.1 信号8.1.1 各种信号的MATLAB实现;【例8_1】产生宽度为1的矩形波和宽度为2π的矩形信号，在如图8_1中，左子图显示的是宽度为1的矩形信号，右子图显示的是宽度为2π。 t=-2:0.01:2; y1=rectpuls(t); %产生矩形信号，省略宽度默认为1 subplot(1,2,1) %一行两列子图的左图 plot(t,y1) axis([-2,2,0,4]) %修改坐标轴范围 t=-5:0.01:5; T=2*pi; y2=rectpuls(t,T); %产生宽度为2π的矩形信号 subplot(1,2,2) plot(t,y2) axis([-5,5,0,4]) ;;?;;8.1.2 信号的基本运算;2. 信号的反褶、尺度变换和平移 【例8_6】将三角波信号f(t)进行运算，变换成f1(t)=f(3t-2)，曲线转换过程如图8_7所示。 t=-3:0.001:3; f=tripuls(t,4,0.5); %产生三角波f(t) hold on subplot(3,1,1) plot(t,f); subplot(3,1,2) f0=tripuls(3*t,4,0.5); %产生尺度变换三角波f(3t) plot(t,f0) subplot(3,1,3) f1=tripuls(3*t-2,4,0.5); %产生尺度变换三角波f(3t-2) plot(t,f1) ;8.2线性系统的时域分析8.2.1 线性时不变系统的零输入响应;【例8_7】 已知y(0)=2，y’(0)=-1，描述某线性系统的微分方程为： ，求解零输入响应。 a=input(Input coefficient=); %输入系统参数 Input denominator coefficient=[1 5 6] n=length(a)-1; y0=input(Input Initial State=); %输入初始状态 Input Initial State=[2 -1] p=roots(a) %求特征根 v=rot90(vander(p)); %得出范德孟德矩阵 v = 1 1 -3 -2 A=v\y0 %计算待定系数 A = -3.000 5．000 y=0; t=0:0.2:10; for k=1:n %绘制零输入响应曲线 y=A(k)*exp(p(k)*t)+y end ;8.2.2 线性时不变系统的全响应;【例8_8】 已知y(0)=2，y’(0)=-1，f(t)=e-t，描述某线性系统的微分方程为： 。求该系统的零输入响应、零状态响应、自由响应、强迫响应和全响应。 yt=dsolve(D2y+5*Dy+6*y=2*exp(-t)) %非齐次方程通解 yt = exp(-t) + C2*exp(-2*t) + C3*exp(-3*t) C2=0;C3=0; ypt=subs(yt) %特解即强迫响应 yt = exp(-t) y=dsolve(D2y+5*Dy+6*y=2*exp(-t),y(0)=2,Dy(0)=-1)%系统的全响应 y = exp(-t) + 3*exp(-2*t) - 2*exp(-3*t) yht=y-ypt yht = 3*exp(-2*t) - 2*exp(-3*t) ;8.2.3 线性时不变系统的冲激响应;8.2.4用卷积计算线性时不变系统的时域响应;【例8_10】计算两个函数f1(t)和f2(t)的卷积，其中f1(t)和f2(t)的波形图如图所示。 clear dt=0.01; t=-1:dt:5; %横坐标t范围 L=length(t); f1=2*rectpuls(t-0.25,3/2); %产生矩形波f1(t) subplot(3,1,1) plot(t,f1) axis([-1 5 0 3]) subplot(3,1,2) f2=t.*rectpuls(t-0.5,2).*(t0t2); %产生f2(t) plot(t,f2) f=conv(f1,f2); %计算卷积 subplot(3,1,3) tf=[2*t(1):dt:2*t(L)] %横坐标范围扩大一倍 plot(t