MATLAB数值积分与拟合.pptVIP

谢谢 数值积分 数值积分原理： 原函数不存在，采用数值积分： 函数: quad 功能: 数值定积分， Quad : 自适应Simpleson积分法。 格式 q = quad(fun,a,b) %近似地从a到b计算函数fun的数值积分，误差为10-6。若给fun输入向量x，应返回向量y，即fun是一单值函数。 q = quad(fun,a,b,tol) %用指定的绝对误差tol代替缺省误差。tol越大，函数计算的次数越少，速度越快，但结果精度变小。 程序： fun = inline(’3*x.^2./(x.^3-2*x.^2+3)’); Q1 = quad(fun,0,2) 计算结果为： Q1 = 3.7224 例1. 计算[0,2]上如下函数的积分 梯形法数值积分 T = trapz(Y) %用等距梯形法近似计算Y的积分。若Y是一向量，则trapz(Y)为Y的积分；若Y是一矩阵，则trapz(Y)为Y的每一列的积分；若Y是一多维阵列，则trapz(Y)沿着Y的第一个非单元集的方向进行计算。 T = trapz(X,Y) %用梯形法计算Y在X点上的积分。若X为一列向量，Y为矩阵，且size(Y,1) = length(X)，则trapz(X,Y)通过Y的第一个非单元集方向进行计算。 T = trapz(…,dim) %沿着dim指定的方向对Y进行积分。若参量中包含X，则应有length(X)=size(Y,dim)。 例2. 用梯形公式计算[-1,1]上1/(1+25*X^2)的积分。 X = -1:.1:1; Y = 1./(1+25*X.^2); T = trapz(X,Y) 计算结果为： T = 0.5492 二元函数重积分的数值计算 1.矩形区域上的二重积分的数值计算 格式 q = dblquad(fun,xmin,xmax,ymin,ymax) 调用函数quad在区域[xmin,xmax, ymin,ymax]上计算二元函数z=f(x,y)的二重积分。输入向量x，标量y，则f(x,y)必须返回一用于积分的向量。 q = dblquad(fun,xmin,xmax,ymin,ymax,tol) %用指定的精度tol代替缺省精度10-6，再进行计算。 q = dblquad(fun,xmin,xmax,ymin, ymax, tol, method) %用指定的算法method代替缺省算法quad。method的取值有@quadl或用户指定的、与命令quad与quadl有相同调用次序的函数句柄。 例如： fun = inline(’y./sin(x)+x.*exp(y)’); Q = dblquad(fun,1,3,5,7) 输出结果： Q = 3.8319e+003 插值与拟合 Matlab中Lagrange 插值函数命令： Lagrange(x,y,x0) 例1：已知数据如下表，试用Lagrange插值多项式求x分别为0.5626,0.5635,0.5645函数的近似值。 xi 0.56160 0.56280 0.56401 0.56521 Yi 0.82741 0.82659 0.82577 0.81495 在命令窗口输入如下执行命令： x=[0.56160;0.56280;0.56401;0.56521]; y=[0.82741;0.82659;0.82577;0.81495]; x0=[0.5626;0.5635;0.5645]; y0=lagrang(x,y,x0) Y0=0.8265 0.8268 0.8231 plot(x,y,’o’,x0,y0,’k*’) 2.分段三次埃尔米特插值 Matlab 中所用命令 pchip(x,y,x0) x=[0.56160;0.56280;0.56401;0.56521]; y=[0.82741;0.82659;0.82577;0.81495]; x0=[0.5626;0.5635;0.5645]; y0=pchip(x,y,x0) Y0=0.8267 0.8262 0.8232 2.分段三次埃尔米特插值 Matlab 中所用命令 pchip(x,y,x0) x=[0.56160;0.56280;0.56401;0.56521]; y=[0.82741;0.82659;0.82577;0.81495]; x0=[0.5626;0.5635;0.5645]; y0=spline(x,y,x0) Y0=0.8265 0.8268 0.8231 拟 合 问 题 已知热敏电阻数据：