第六matlab在信号与系统中的应用精要.ppt

无穷级数 第六章 matlab在信号与系统中的应用 6.1 连续信号和系统的MATLAB描述 6.1.2 LTI连续系统分析 C1+ C2+…+Cn = y0 y0 = y(0) p1C1+ p2C2+…+ pnCn=Dy0 即 V·C = Y0 其解为 C =V \ Y0 式中 V为范德蒙矩阵，在MATLAB的特殊矩阵库中有vander。调用方法: V=vander(p) n阶LTI系统的冲激响应 n阶微分方程，写成系统函数为： 冲击响应就是H(s)的拉普拉斯反变换，可以把H(s)展开为极点留数式。 其反变换为 课后练习 6.2 傅里叶分析 图示的周期性方波，其傅里叶级数为 分别计算 直到9次谐波，并做图。 6.2.2 信号的频谱分析-以方波为例 将积分上下限定为0～10s，并将t分成N等份，用求和代替积分。这样，傅立叶变换式可写为 求和可以用f(t)行向量乘以e-j?tn列向量来实现。式中Δt是t的增量，在程序中，用dt表示。 求不同?处的F值，都用同一公式，这就可以利用MATLAB中的元素群运算。将?设为一个行数组，代入上式，则可写为（程序中?用w表示） 其中，F是与w等长的行向量，t ‘是列向量，w是行向量，t ’*w是一矩阵，其行数与t相同，列数与w相同。这样，此式就完成了傅里叶变换。 类似地也可得到傅里叶逆变换表示式为 算得的时域信号波形及其频谱图如右。 下图为采样周期较低时的情况，有明显的频率泄漏。 6.2.3 信号滤波 计算幅度为1，宽度为5s的矩形脉冲（同例6.11）通过下列滤波器的响应。 （1）理想低通滤波器， （2）低通滤波器 解： 滤波器输出的频谱 Y(j?)=F(j?)·H(j?) 其时间响应y(t)是Y(j?)的傅里叶反变换。 （1）理想低通滤波器的截止角频率?c=10，故只取F(j?)中? = 0～10的部分，用MATLAB语言表述，输出频率分量对应的?的下标数组为 n2 = find ((w=- wc) (w= wc) 其对应的频率数组为w2 = w(n2) 频段内的频谱数组为F2 = F(n2:)，它就是滤波后的频谱数组Y2，其逆变换即 y2 = F2*exp(j*n2*t)/pi*dw （2）三阶低通滤波器的频率响应 滤波器的输出为 6.4 线性时不变系统的模型 1．连续系统 ? 状态空间型 设x为状态变量，u为输入，y为输出，系统的状态方程为： 如果系统是n阶的，输入有nu个，输出有ny个，则A为n×n阶，B为n×nu阶，C为ny×n阶，而D为ny×nu阶矩阵，对单输入单输出(SISO)系统，ny=nu=1。 ? 传递函数型 单输入单输出(SISO)n阶系统的传递函数为 知道分子系数矢量f =［f (1)，f (2)，…，f (m + 1)］，分母系数矢量g =［g(1)，g(2)，…，g (n+1)］，就惟一地确定了系统的模型（注意系统的阶次n）。而对物理可实现的系统，必有n≥m。 ? 零极增益型 对传递函数分子和分母进行因式分解，可得 令z = [z(1)，z(2)，…，z(m)]为系统的零点矢量，p = [p(1)，p(2)，…，p(n)]为系统的极点矢量，k为系统增益，它是一个标量。可以看出，H(s)有m个零点、n个极点。物理可实现系统的n≥m，系统的模型将由矢量z，p及增益k惟一确定，故称为零极增益模型。零极增益模型通常用于描述SISO系统。并可以推广到MISO系统。 ? 极点留数型 如果零极增益模型中的极点都是单极点，将它分解为部分分式，可得 其中p =［p(1)，p(2)，…，p(n)］仍为极点矢量，而r =［r(1)，r(2)，…，r(n)］为对应于各极点的留数矢量，p，r两个矢量及常数h惟一地决定了系统的模型。 典型模型表达式的比较 比较一下这四种情况下模型系数的总个数。假定都是SISO系统，阶数为n，则状态空间型有n2+2n+1个系数；传递函数型为m+n+1个(不含g(1)) (注意，由于m≤n，系数的数目小于等于2n+1)；零极增益型的系数个数为n+m+1；而极点留数型为2n+1。 因此，传递函数法的待定系数最少，而状态空间法的待定系数最多。这说明了状态空间法中有许多冗余的系数。事实上，同一个系统可以有无数个状态空间矩阵A，B，C，D的组合来描述，其他描述方法则都是惟一的。 2．离散系统 以上四种表示模型的方法可以全部推广至离散系统。只是将系数矩阵后面加小写字母d，便有： ? 状态空间型 x[n + 1] = Adx[n] + Bdu[n]