[工学]dsp第6章.pptVIP

第 6章 MATLAB 在信号与系统中的应用 6.1　用MATLAB实现系统的时域分析 1　连续信号的MATLAB表示　　严格地说,MATLAB不能处理连续信号,它是用连续信号在等间隔点的样值来近似表示连续信号的。当采样间隔足够小时,这些样值就能较好地近似表示连续信号。　　MATLAB提供了大量的生成基本信号的函数。最常用的指数信号、正弦信号是MATLAB的内部函数,即不安装任何工具箱就可调用的函数。 例6.1 连续信号的MATLAB描述 %《MATLAB在电子信息课程中的应用》第六章例6.1程序q601 % 连续信号的MATLAB表述 clear, t0=0;tf=5;dt=0.05;t1=2; 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处的脉冲信号 subplot(2,2,1),stairs(t,x1),grid % 绘图，注意为何用stairs命令 axis([0,5,0,1.1/dt]) %(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]) %(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 % 绘图, 例6.2 LTI系统的零输入响应 n阶线性时不变连续系统的微分方程 已知y及其各阶导数的初始值为y(0),y(1)(0)，…，y(n-1)(0)，求系统的零输入响应。 解：方程的解为 p1, p2…,pn是方程a1?n+a2?n-1+…+ an?+ an+1 =0的根, 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) %《MATLAB在电子信息课程中的应用》第六章例6.2程序q602% 求任意线性系统的零输入分量： % 电子工业出版社出版 陈怀琛 吴大正 高西全合著 clear all;close all; a=[3,5,7,1]; n=length(a)-1; disp(输入初始条件向量(长度应比a小一) Y0=[y0,Dy0,D2y0,...] ) for m=1:3; Y0=input(Y0=(书上分别取[1,0,0],[0,1,0],[0,0,1]) ); p=roots(a); % 求系统的极点 V=rot90(vander(p));c= V\Y0; % 用范德蒙特矩阵求对应于各极点的系数 dt=0.1; %input(dt=(书上取0.1) ); tf=5; %input(tf=(书上取 5) ) 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), hold on end grid on text(3.5,0.7,Y0=[1 0 0]); text(1.5,0.7,Y0=[0 1 0]); text(2,0.35,Y0=[0 0 1]); 　　3 连续系统