MATLAB在信号和系统中的应用.pptVIP

例6.21 连续系统状态方程的解 线性连续系统在输入信号为零时的状态方程表示式为 ，其中x为n×1的列向量，A为n×n阶的方阵，求其在初始条件作用下的解。 解：建模 ：按线性方程理论，此齐次方程的解为 其中 此处要把方阵At作为一个整体求指数函数，就需要调用MATLAB矩阵指数函数expm。注意它与exp函数不同。 例6.21 连续系统状态方程的解 如果A是n×n阶，则expm(At)也是n×n阶，这时如果要算一系列的t值所对应的expm(At)，就不可能像标量指数那样用元素群运算方法了，必须用for循环。不仅如此，如果t的长度为nt，则状态转移矩阵F(t)=expm(At)将是一个n×n×nt的三维矩阵，需要用MATLAB中高维矩阵的概念来解决问题。举以下数字例来说明。 设A = [-2，1;-17，-4]，x0 = [3 ;4]; 求上述线性方程的状态转移矩阵F(t)及x的解。 例6.22 离散系统状态方程的响应 离散系统的状态方程为 初始条件 输入信号为阶跃函数，求前10步的解x。 解：已知x(0)和u(0)就可求出k=1时的x(1)，再依次递推求x(2)，…，x(10)。在程序中可用循环语句来完成。离散系统的计算比连续系统简单得多，只要作四则矩阵运算即可。 第六章  MATLAB在信号与系统中的应用 例6.1 连续信号的MATLAB描述 （1）单位冲激函数 （2）单位阶跃函数： （3）复指数函数 例6.4 卷积的计算 根据卷积公式： 因此编程的过程为： （1）写出h(t)的MATLAB表达式； （2）写出u(t)的MATLAB表达式； （3）利用MATLAB的卷积语句y=conv(u，h)求解 （4）画曲线plot(t,y)。 例6.7 方波分解为多次正弦波之和 图示的周期性方波，其傅里叶级数为 分别计算 直到9次谐波，并做图。 图6.7-1 输入周期性方波 - - x 例6.13 离散信号的MATLAB表述 编写MATLAB程序，产生下列基本脉冲序列： （1）单位脉冲序列： 起点n0，终点nf，在ns处有一单位脉冲(n0≤ns≤nf)。 n1=n0:nf; x1=[zeros(1,ns-n0),1,zeros(1,nf-ns)]; 用逻辑关系编程： n1 = n0:nf; x1=[(n1-ns)==0] 例6.13 离散信号的MATLAB表述 编写MATLAB程序，产生下列基本脉冲序列： （2）单位阶跃序列： 起点n0，终点nf，在ns前为0,在ns后为1(n0≤ns≤nf)。 n2=n0:nf; x2=[zeros(1,ns-n0),ones(1,nf-ns+1)]; 用逻辑关系编程： n1 = n0:nf; x1=[(n1-ns)=0] 例6.13 离散信号的MATLAB表述 编写MATLAB程序，产生下列基本脉冲序列： （3）复数指数序列： n3 = n0:nf; x3=exp((-0.2+0.5j)*n3); 例6.11 方波的频谱分析 频谱分析用福利也变实现，傅里叶变换公式： 积分可以用离散求和来近似。将积分区间t分成N等份，用求和代替积分。这样，傅立叶变换式可近似为 求和可以用f(t)行向量乘以e-j?tn列向量来实现。式中Δt是t的增量，在程序中，用dt表示。 频域离散化 将频域取值区间也等分，在各点上的取值为： 利用元素群运算简化表达式： 令 正变化简写为： 同理反变换为： 例6.11 方波的频谱分析(续) 求不同?处的F值，都用同一公式，这就可以利用MATLAB中的元素群运算。将?设为一个行数组，代入上式，则可写为（程序中?用w表示） 其中，F是与w等长的行向量，t ‘是列向量，w是行向量，t ’*w是一矩阵，其行数与t相同，列数与w相同。这样，此式就完成了傅里叶变换。 类似地也可得到傅里叶逆变换表示式为 例6.11 方波的频谱分析(续) 算得的时域信号波形及其频谱图如右。 下图为采样周期较低时的情况，有明显的频率混叠。 例6.12 信号通过滤波器 计算幅度为1，宽度为5s的矩形脉冲（同例6.11）通过下列滤波器的响应。 （1）理想低通滤波器， （2）低通滤波器 解： 滤波器输出的频谱 Y(j?)=F(j?)·H(j?) 其时间响应y(t)是Y(j?)的傅里叶反变换。 例6.12 信号通过滤波器(续) （1）理想低通滤波器的截止角频率?c=10，故只