[所有分类]第6章matlab在信号与系统中的应用.pptVIP

第 6章 在信号与系统中的应用 1.连续信号和系统 连续信号和线性时不变(LTI)系统特点 连续信号,在计算机中不能用数值计算的方法来表示，因为给出的是各个样本点的数据，应当采用符号推理的方法来分析. 只有当样本点取的很密时才可看成连续信号，即假设对采样点密度而言，信号变化足够慢。 例6.1 连续信号的MATLAB描述 （1）单位冲激函数 看作宽度为Δ(dt),幅度为1/Δ的矩形脉冲. clear, t0=0;tf=5;dt=0.05;t1=1; t=[t0:dt:tf]; % (1)单位脉冲信号, % 在t1(t0≤t1≤tf)处有一持续时间为dt,面积为1的脉冲信号,其余时间均为零。 t=[t0:dt:tf];st=length(t); n1=floor((t1-t0)/dt); % 求t1对应的样本序号 x1=zeros(1,st); % 把全部信号先初始化为零 x1(n1)=1/dt; % 给出t1处的脉冲信号 %x1=[zeros(1,n1-1),1/dt,zeros(1,st-n1)]; subplot(2,2,1), stairs(t,x1),grid % 绘图，注意为何用stairs命令 axis([0,5,0,1.1/dt]) 或 冲击函数 t=-20:0.01:20; y=(t==0); plot(t,y) axis([-20 20 -1 2]) %(2)单位阶跃信号， %信号从t0到tf，在t1(t0≤t1≤tf) 前为0，到t1处有一跃变，以后为1. % 程序前几句即求t,st,n1的语句与上同，只把x1处改为x2 x2 = [zeros(1,n1-1),ones(1,st-n1+1)]; %产生阶跃信号 subplot(2,2,3),stairs(t,x2),grid % 绘图 axis([0,5,0,1.1]) 或 单位阶跃函数 t=-20:0.1:20; y=(t=0); plot(t,y) axis([-30 30 -1 2]) %(3)复数指数信号 u=-0.5;w=10;x3=exp((u+j*w)*t); subplot(2,2,2),plot(t,real(x3)),grid % 绘图, subplot(2,2,4),plot(t,imag(x3)),grid % 绘图, 1.3 连续信号的时域运算与时域变换 信号的尺度变换、翻转、平移 已知信号 用matlab求f(t+2),f(t-2),f(-t),f(2t),-f(t),并绘出时域波形。 syms t f=sym((1+t/2)*(heaviside(t+2)-heaviside(t-2))) subplot(2,3,1),ezplot(f,[-3,3]); y1=subs(f,t,t+2) subplot(2,3,2),ezplot(y1,[-5,1]) y2=subs(f,t,t-2) subplot(2,3,3),ezplot(y2,[-1,5]) y3=subs(f,t,-t) subplot(2,3,4),ezplot(y3,[-3,3]) y4=subs(f,t,2*t) subplot(2,3,5),ezplot(y4,[-2,2]) y5=-f subplot(2,3,6),ezplot(y5,[-3,3]) 例6.2 LTI系统的零输入响应 n阶线性时不变连续系统的微分方程 已知y及其各阶导数的初始值为y(0),y(1)(0)，…，y(n-1)(0)，求系统的零输入响应。 解：当LTI系统的输入为0时,其输入响应为微分方程的齐次解.(右端为0), 方程的解为 p1, p2…,pn是方程a1?n+a2?n-1+…+ an?+ an+1 =0的根, 由roots求得，而C1，…，Cn由y及其各阶导数的初始值来确定。 例6.2 LTI系统的零输入响应(续) C1+ C2+…+Cn = y0 y0 = y(0) p1C1+ p2C2+…+ pnCn=Dy0 例6.2 LTI系统的零输入响应(续) 即 V·C = Y0 其解为 C =V \ Y0 式中 V为范德蒙矩阵，在MATLAB的特殊矩阵库中有vander。调用方法: V=vander(p) a=input(输入分母系数向量a=[a1,a2,...]=(书上取 [3,5,7,1] ) ); n=length(a)-1; disp(输入初始条件向量(长度应比a小1) Y0=[y0,Dy0,D2y0,...] ) Y0=input(Y0=(书上分别取[1,0,0],[0,1,0],[0