基于MATLAB的数值分析(6)(精品·公开课件).ppt

第六章 数值逼近 6.2 Lageange 插值多项式 6.3 插值多项式的误差及 Chebyshev点的插值 ??????? 6.4 分 段 插 值 ??????? interp1: 一维分段插值 interp2: 二维分片插值 interp3: 三维分块插值 spline: 三次样条插值 6.5 曲线拟合 ??????? 6.6 数值积分 【例 】 已知x=0:0.1:1.0 y=1.0,2.5,3.0,2.0,1.5,0.9,0,-1,-2,-1.5,-0.8 拟合模型s(x)=c1+c2sin(x)+c3xsin(x)。 x=(0:0.1:1); y=[1.0,2.5,3.0,2.0,1.5,0.9,0,-1,-2,-1.5,-0.8 ]; A=[ones(size(x)),sin(x),x.*sin(x)]; C=regress(y,A) C = 2.2262 -0.9358 -4.1223 一、一元函数的数值积分 常见的一元函数的数值积分指令 quad: 自适应Simpson求积公式 quad8：自适应高阶Newton-Cotes求积公式 quadl： 自适应Lobatto求积公式 trapz： 梯形法 sum: 矩形法 fnint: 利用样条函数求不定积分，与spline,ppval 配合使用 y=quad(‘fun’,a,b,tol,trace):自适应Simpson求积公式 y=quad8(‘fun’,a,b,tol,trace)：自适应高阶Newton- Cotes求积公式 如果被积函数带有参变量，如f(x,t1,t2) 调用格式：y=quad(‘fun’,a,b,tol,trace,t1,t2,…) 【例】求 ，其精确值为 。 （1）符号解析法 syms x;IS=int(exp(-x*x),x,0,1) vpa(IS) IS = 1/2*erf(1)*pi^(1/2) ans = .74682413281242702539946743613185 ?（2）MATLAB指令quad和quadl求积 fun=inline(exp(-x.*x),x); Isim=quad(fun,0,1),IL=quadl(fun,0,1) Isim = 0.7468 IL = 0.7468 * * 6.1 关于多项式的MATLAB命令 一、 多项式表达方式的约定 多项式 用行向量表示 用比较习惯的方式显示多项式：pp=poly2str(p,x) 【例】多项式 可表示为 p=[2 1 4 5] pp=poly2str(p,x) 二、多项式运算函数 r=roots(p)：求多项式的零点 p=poly(r) : 以r为零点的多项式 p=poly(A): A的特征多项式 PA=polyval(p,S):按数组运算规则，计算多项式的值 其中S，PA为矩阵 PM=polyvalm(p,S):按矩阵运算规则，计算多项式的值, 其中S，PM为矩阵 p=conv(p1,p2)：多项式的乘积 [q,r]=deconv(p1,p2)：多项式的除法，p1/p2 p1(x)=p2(x)q(x)+r(x) 【例】由给定根向量求多项式系数向量。 R=[-0.5,-0.3+0.4*i,-0.3-0.4*i]; P=poly(R) PPR=poly2str(P,x) P = 1.0000 1.1000 0.5500 0.1250 PPR = x^3 + 1.1 x^2 + 0.55 x + 0.125 【例】求多项式 的零点。 r=roots([1 -6 15 -20 15 -6 1]) r = 1.0042 + 0.0025