MATLAB基础——数值计算.ppt

Matlab及其应用 Matlab及其应用 魏立秋 哈工大先进动力研究所 weiliqiu@hcms.hit.edu.cn 基于Matlab的数值分析方法 Matlab多项式处理 Matlab的多项式函数位于\toolbox\matlab\polyfun 1、多项式的表示方法： 采用行向量表示多项式，将多项式系数按照降幂排列； 鉴于MATLAB无零下标，故把多项式的一般形式表达为： 系数行向量P=[a1,a2,a3….,an,an+1] poly2str函数 多项式加、减：对于次数相同的若干个多项式，可直接对多项式系数向量进行加、减的运算。如果多项式的次数不同，则应该把低次的多项式系数不足的高次项用零补足，使同式中的各多项式具有相同的次数。 2、多项式求根、由根创建多项式（exp31.m） 函数格式：R=roots(C)。 这里C为多项式的系数C(1),C(2),…,C(N),C(N+1)；解得的根赋值给列数组R，即R(1),R(2), …,R(N)。 函数格式：C=poly(R) R为N*N矩阵时，返回值C是一个N+1个元素的行向量，是该矩阵特征多项式的系数 R为向量时，返回值是多项式系数行向量 3、多项式的乘除 多项式的乘法和除法运算实质上就是多项式系数向量的卷积和解卷运算。 长度为m的向量a和长度为n的向量b的卷积c定义为： c的长度为（m+n-1） 解卷是卷积的逆运算，向量a对向量c解卷将得到商向量q和余量r，并满足： 3、多项式的乘除 乘法函数格式：c=conv(a,b)。 其中a，b为两个多项式函数的系数向量，返回值c为乘积后多项式的系数向量 除法函数格式：[q,r]=deconv(c,a) 其中，c为作为被除数的多项式矩阵系数向量，a为除数多项式矩阵系数，q为商系数，r为余量 3、多项式求值 函数格式：y=polyval(p,x)。 其中，p为多项式的系数向量，x是指定的自变量的值（标量/向量/矩阵），对每一个元素计算多项式的值返回给y 函数格式：y=polyvalm(p,x)。 计算矩阵多项式的值，即以矩阵为整体作为多项式自变量进行计算 4、多项式导数 函数格式： polyder(p) 计算系数向量为p的多项式的导数 polyder(a,b) 计算多项式a*b的导数 [q,d]=polyder(a,b) 计算多项式的商a/b的导数，并用q/d表示 5、多项式的部分分式展开 展开或合并分式多项式对系统传递函数特别有用 函数格式： [r,p,k]=residue(b,a) [b,a]=residue(r,p,k) 向量r、p的长度和向量a、b之间有如下关系： length(r)=length(p)=length(a)-1 当向量b的长度小于a时，向量k中没有元素，否则满足如下关系：length(k)=length(b)-length(a)+1 例： 1、已知： 求： 2、求极限： 3、求积分： Matlab 符号计算 符号计算是在运算时，无需事先对变量赋值，而将所得到的结果以标准的符号形式来表示——功能：代替手工的推公式。 以加拿大滑铁卢大学研制的一种数学软件Maple为内核，1993年MathWorks公司购买了Maple的使用权，依赖Maple已有的函数库，开发实现符号计算的工具箱。 1. 符号变量的定义 参与符号运算的对象是符号变量、符号表达式或是符号矩阵，符号变量要先定义、后引用，可用sym函数和syms函数将运算量定义为符号型数据 sym(‘x’); sym(‘a’) syms x y 符号常数和符号变量 sin(at) y=ax+b 例：解一元二次方程solve(a*x^2+b*x+c) findsym(f,n) 2、符号表达式 符号表达式有符号变量、函数、算术运算符等组成，符号表达式的书写格式与数值表达式相同。 例：sym(‘x’) f=1+sqrt(5*x)/2 sym(‘a*x^2+b*x+c’) Matlab不会自动创建表达式中的符号运算变量，符号变量或是符号表达式必须经过定义后才能进行符号运算 ezplot(f,[a,b]) 3、极限运算 从高等数学的发展来看，极限概念占有极其重要的地位，在实际运算中，由于需要很多技巧，因而还是比较复杂。 函数格式： limit(f) x?0 limit(f,a) x?a limit(f,x,a) x?a findsym limit(f,x,a,’left’) x?a- limit(f,x,a,’right’)