[数学]Matlab基础2012-2.pdfVIP

MATLAB基础 ———————— 连续时间信号与系统连续时间信号与系统连续时间信号与系统连续时间信号与系统 课程安排课程安排 MATLAB简介 连续时间信号与系统连续时间信号与系统 离散时间信号与系统 MATLAB实验 内容内容 连续时间系统的时域分析 傅里叶变换 、 拉普拉斯变换 s域分析 内容内容 1 连续时间系统的时域分析连续时间系统的时域分析 微分方程式的建立与求解 零输入响应与零状态响应 冲激响应与阶跃响应 微分方程式的建立与求解微分方程式的建立与求解 • 连续时间系统的研究方法包括 输输入－输输出法法 （端端口描描述法法） 系统状态变量分析法 • 输入－输出法 LTILTI 系统可以用系统可以用一元高阶微分方程描述元高阶微分方程描述 a = [1[1,77,00,12]12]; bb = 11; %%定义微分方程的行矢量定义微分方程的行矢量 sys = tf(b,a) %生成系统描述sys （tf=transfer function ） 微分方程式的建立与求解微分方程式的建立与求解 • 状态变量描述法 一个一元高阶微分方程，必然可以化成两个多元一阶微分方程组， 即状态方程和输出方程（观测方程） sys = ss(A,B,C,D) %ss=state space，ABCD为状态方程和输出方 程的矩阵 微分方程式的建立与求解微分方程式的建立与求解 • 微分方程的解包括齐次解和特解两部分 • 齐次解即系统特征方程的根，用roots函数计算 p = [1 7 16 12]; %定义多项式p a = roots(p); %求解多项式p=0的根a a %打印显示 a = -3.0000 -2.0000 + 0.0000i -2.0000 - 0.0000i 微分方程式的建立与求解微分方程式的建立与求解 • 特解即系统（采用微分方程表示）在给定信号激励下的输出 • 用lsim函数进行仿真 a = [1,2,3]; %定义多项式p b = [1,1]; %求解多项式p=0的根a syssys == tftf((bb,aa);); %%建立系统描述建立系统描述syssys t = [0:0.1:10]’; %定义仿真时间为0s到10s，抽样间隔为0.1s e1 = t.^2’; %定义激励信号e1 r1r1 == lsimlsim(sys(sys,e1e1,t);t); %%用用e1e1激励激励syssys，输出为输出为r1r1 e2 = exp(t); … 微分方程式的建立与求解微分方程式的建立与求解 微分方程式的建立与求解微分方程式的建立与求解 • Help polyfun roots - Find polynomial roots. poly - Convert roots to polynomial. pollyvall - EEvalluatte pollynomiiall. polyvalm - Evaluate polynomial with matrix argument. residueresidue - PartialPartial-fractionfraction expansionexpansion (residues).(residues). polyfit - Fi