MATLAB多项式-符号运算.pptVIP

MATLAB的数值计算;matlab语言把多项式表达成一个行量， 该向量中的元素是按多项式降幂排的。 f(x)=anxn+an-1xn-1+……+loa0 可用行向量 p=[an an-1 …… a1 +a0]表示 如：a(x)=2x3+4x2+6x+8; 对应a=[2,4,6,8]; b(x)= 2x3+6x+8; b=? ;函数1：ploy2str() 作用：生成多项式表达式的字符串 p =[1.00,-6.00,-72.00,-27.00]; p是多项式p(x)=x3-6x2-72x-27的matlab描述方法，我们可用： a=poly2str(p,x)显示数学多项式的形式 a=x^3 - 6 x^2 - 72 x – 27 并将表达式的结果以字符串的形式赋值给变量a;函数2：roots（） 作用： 求多项式的根; 多项式加减运算;函数文件： function y=polyadd(x1,x2) n1=length(x1); n2=length(x2); if n1n2 x2=[zeros(1,n1-n2),x2]; elseif n1n2 x1=[zeros(1,n2-n1),x1]; end; y=x1+x2 ;多项式乘运算conv();多项式除运算deconv();多项式求导;代数方程组求解;1.恰定方程组的解;方程ax=b a=[1 2;2 3];b=[8;13]; ?x=inv(a)*b ? x=a\b x = x = 2.00 2.00 3.00 3.00;2.超定方程组的解;例: x1+2x2=1 2x1+3x2=2 3x1+4x2=3 a=[1 2;2 3;3 4];b=[1;2;3]; 解1 x=a\b 解2 x=inv(a?a) ? a ? b x = x = 1.00 1.00 0 0.00 ;3.欠定方程组的解; x1+2x2+3x3=1 2x1+3x2+4x3=2 a=[1 2 3;2 3 4];b=[1;2]; x=a\b x=pinv(a)?b x = x = 1.00 0.83 0 0.33 0 -0.17;多项式拟合;拟合程序例1：;插值拟合 （1）插值的定义——是在给定的原始数据点之间用某些特定的数学算法插入一些数据点，当原始数据点和插入数据点连线后得到穿过原始数据点的光滑曲线。 （2）当不能很快地求出所需中间点的函数时，插值是一个非常有价值的工具。 （3）Matlab提供了线性插值、三次插值和三次样条插值3种选择。 直线插值：在原始数据点之间用直线拟合 三次插值及三次样条插值：在原始数据点之间用三次多项式曲线分段连接;三次插值和三次样条插值的区别;插值函数：interp1 调用格式：interp1(x,y,xi,’method’) 其中： (1)x、y为离散数据的坐标向量，要求长度必须相同 (2)xi为插入点的数据向量 (3)Method为插值方法，为字符串变量 有三种方法供选择： ’linear’、’cubic’、’spline’ 分别表示线性插值、三次插值和三次样条插值 ;曲线拟合例题2： 将下列离散点用三次样条插值拟合，并绘制： (0,11.2)、(0.2,16.5)、(0.4,20.4)、(0.6,26.3)、(0.8,30.5)、(1,28.7) 程序： x=0:0.2:1; y=[11.2,16.5,20.4,26.3,30.5,28