第5章积分变换与复变函数问题的计算机求解.ppt

*星期六, 2008-9- 6, 13:09:02 复变函数中一类重要的问题是奇点与留数的概念与求取，本章介绍了利用 MATLAB 语言进行留数运算的方法，并以其为基础编写了有理函数部分分式展开的求解函数，比 MATLAB 自己提供的 residue() 函数更进了一步，本章还讨论了有理函数 Laplace 反变换的求解方法和化简方法，基于留数定理还探讨了封闭曲线积分的求解方法 *星期六, 2008-9- 6, 13:09:02 *星期六, 2008-9- 6, 13:09:02 例 5.37 重新求解：求部分分式展开 MATLAB求解语句： *星期六, 2008-9- 6, 13:09:02 例5-40 求部分分数展开 试图用MATLAB求解 理论上不能展开 *星期六, 2008-9- 6, 13:09:02 如果分母多项式 D(x) = 0 的解析解无法获得，则应使用高精度数值解 数学表示 其中， 是数值解 MATLAB求解新函数 变通方法 *星期六, 2008-9- 6, 13:09:02 函数清单 *星期六, 2008-9- 6, 13:09:02 例 5.41 求下式的部分分式展开 MATLAB近似求解 *星期六, 2008-9- 6, 13:09:02 5.5.6 基于部分分式展开的Laplace反变换 实系数有理函数 如果有 则必有 可以推导出 其中 *星期六, 2008-9- 6, 13:09:02 算法的实现 MATLAB函数 *星期六, 2008-9- 6, 13:09:02 例 5.42 有理函数 计算出部分分式展开 MATLAB求解命令： *星期六, 2008-9- 6, 13:09:02 5.5.7 封闭曲线积分问题计算 封闭曲线积分的数学表达式 其中，G 是一个逆时针方向的闭曲线，那么 该封闭曲线内包围 m 个奇点 *星期六, 2008-9- 6, 13:09:02 例 5.34 给定 计算在 | z | = 6 上的封闭曲线积分 *星期六, 2008-9- 6, 13:09:02 直接积分方法 圆的方程 封闭曲线积分直接计算 *星期六, 2008-9- 6, 13:09:02 例 5.44 计算 其中G 是 | z | = 2（逆时针圆周封闭曲线） 曲线内奇点 MATLAB求解语句： *星期六, 2008-9- 6, 13:09:02 与 进行比较 直接曲线积分方法 改变积分路径 G 为 | z | = 4 *星期六, 2008-9- 6, 13:09:02 5.6 差分方程迭代求解与复平面映射分形 一般差分方程的解析求解方法 线性时变差分方程的数值解法 线性时不变差分方程的解法 一般非线性差分方程的数值求解方法 *星期六, 2008-9- 6, 13:09:02 差分方程求解 常系数线性差分方程的一般形式为 差分方程简单记号 T 采样周期，常系数 *星期六, 2008-9- 6, 13:09:02 5.6.1一般差分方程的解析求解方法 差分方程求 z 变换 *星期六, 2008-9- 6, 13:09:02 得出 其中 *星期六, 2008-9- 6, 13:09:02 算法的实现 MATLAB函数 求解语句 *星期六, 2008-9- 6, 13:09:02 例5.45 差分方程 初值 输入信号 求解 *星期六, 2008-9- 6, 13:09:02 5.6.2 线性时变差分方程的数值解法 线性时变差分状态方程 其中， 递推 *星期六, 2008-9- 6, 13:09:02 继续 最终得出 可通过递推方法求解 *星期六, 2008-9- 6, 13:09:02 例 5.46 求解离散线性时变差分方程 其中 可以通过递推方法直接求解 *星期六, 2008-9- 6, 13:09:02 采用迭代方法 *星期六, 2008-9- 6, 13:09:02 5.6.3 线性时不变差分方程的解法 若线性时不变系统有 则有 线性时不变状态方程 *星期六, 2008-9- 6, 13:09:02 两边取 z 变换 由此得出解析解 由 iztrans 函数可以直接求解 *星期六, 2008-9- 6, 13:09:02 例 5.47 离散系统的状态方程如下 零初始状态 求解 *星期六, 2008-9- 6, 13:09:02