MATLAB的数据建模.docVIP

第6章 MATLAB数据建模 本章主要学习数据挖掘和建模的一些基本方法和相关的MATLAB命令，包括插值，拟合，回归分析，函数逼近等。 6.1多项式 在MATLAB中，多项式是利用一个行向量来表示的，它的系数是按照降序方式排列的， 例如，多项式p1(x)= x3+21x2+20x可以表示为： P1=[1 21 20 0] %常数项为0 6.1.1多项式的求值、求根和部分分式展开 1. 多项式求值 函数polyval可以用来计算多项式在给定变量时的值，是按数组运算规则进行计算的。 语法： polyval(p,s) 说明：p为多项式, s为给定矩阵。 【例6. 1】计算p1(x)= x3+21x2+20x多项式的值。 p1=[1 21 20 0]; polyval(p1,2) %计算x=2时多项式的值 ans = 132 x=0:0.5:3; polyval(p1,x) %计算x为向量时多项式的值 ans = 0 15.3750 42.0000 80.6250 132.0000 196.8750 276.0000 2. 多项式求根 roots用来计算多项式的根。 语法： r=roots(p) 说明：p为多项式；r为计算的多项式的根，以列向量的形式保存。 与函数roots相反，根据多项式的根来计算多项式的系数可以用poly函数来实现。 语法： p=poly (r) 【例6.1续】计算多项式p1(x)= x3+21x2+20x的根以及由多项式的根得出系数。 roots(p1) %计算多项式的根 ans = 0 -20 -1 poly([0;-20;-1]) %计算多项式的系数 ans = 1 21 20 0 3. 多项式部分分式展开** 在许多的工程实际应用中，例如傅里叶变换，拉普拉斯变化和Z变换中，都会出现两个多项式的比值，这时就需要对多项式进行部分分式展开运算， 用residue函数来实现将分式表达式进行多项式的部分分式展开成以下形式： 语法： [r,p,k]=residue(b,a) 说明：b和a分别是分子和分母多项式系数行向量；r是[r1 r2 …rn]行向量，表示部分分式展开的常数项；p为[p1 p2 …pn]极点行向量；k为余数。 【例6.2】将表达式进行部分分式展开。 p1=[1 21 20 0]; p3=[100 200]; [r,p,k]=residue(p3,p1) r = -4.7368 -5.2632 10.0000 p = -20 -1 0 k = [] 程序分析：表达式展开结果为+0。 6.1.2多项式的四则运算 1. 加减法 MATLAB没有提供专门进行多项式加减运算的函数，事实上，多项式的加减运算，就是其对应的系数向量的加减运算，加减运算时，向量的大小必须相同，缺项的用零补齐。 【例6.3】求多项式x3-2x2+5x+3和6x-1的和。 clear all; p1=[1 -2 5 3]; p2=[0 0 6 -1]; c=p1+p2 c= 1 -2 11 2 也即c= x3-2x2+11x+2 2. 多项式的乘法和除法 多项式的乘法 语法： p=conv(pl,p2) 说明：p是多项式p1和p2的乘积多项式。 多项式的除法 语法： [q,r]=deconv(pl,p2) 说明：除法不一定会除尽，会有余子式。多项式p1被p2除的商为多项式q，而余子式是r。 【例6.4】计算表达式。 a1=[1 0]; %对应多项式s a2=[1 1]; %对应多项式s+1 a3=[1 20]; %对应多项式s+20 p1=conv(a1,a2) p1 = 1 1 0 p1=conv(p1,a3) %计算s(s+1)(s+20) p1 = 1 21 20 0 [p2,r]=deconv(p1,a3) %计算多项式除法的商和余子式 p2 = 1 1 0 r = 0 0 0 0 conv(p2,a3)+r %用商*除式+余子式验算 ans = 1 21 20 0 3. 多项式的求导** 对多项式求导的函数是polyder，其调用格式为： p=polyder(p1)；求多项式p1的导函数； p=