Ch5数值逼近.pptVIP

自适应Simpleson积分法：quad、quadl、quad8 q = quad(fun,a,b) 近似地从a到b计算函数fun的数值积分，误差为10-6。若给fun输入向量x，应返回向量y，即fun是一单值函数。 q = quad(fun,a,b,tol) 用指定的绝对误差tol代替缺省误差。tol越大，函数计算的次数越少，速度越快，但结果精度变小。 q = quad(fun,a,b,tol,trace,p1,p2,…) 将可选参数p1,p2,…等传递给函数fun(x,p1,p2,…)，再作数值积分。若tol=[]或trace=[]，则用缺省值进行计算。 [q,n] = quad(fun,a,b,…) 同时返回函数计算的次数n [q,n] = quadl(fun,a,b,…) 用高精度进行计算，效率可能比quad更好。 [q,n] = quad8(fun,a,b,…) 该命令是将废弃的命令，用quadl代替。 fun = inline(3*x.^2./(x.^3-2*x.^2+3)); Q1 = quad(fun,0,2) Q2 = quadl(fun,0,2) Q1 = 3.7224 Q2 = 3.7224 一元函数的数值积分：方法比较 （1）将(0，pi/2)10等分，步长为pi/20，利用trapz(x)函数按梯形算法计算 （2）使用quad(‘fun’,a,b)函数命令的辛普森公式法： （3）蒙特卡罗均值估计法：蒙特卡罗法是一种随机实验方法。本例使用sum(x)函数，sum(x)函数的功能是输入数组x，输出x的和。 z = 0.99794 m=pi/20; x=0:m:pi/2; y=sin(x); z=trapz(y)*m z=quad(sin,0,pi/2) r=100000; x=rand(1,r); y=sin(x*pi/2); z=sum(y)*pi/2/r z = 1 z = 1.001 上面三种方法得出的结果不同这是由于所选方法各自产生的误差不同所造成的。 应用1：饱和蒸气压的拟合 1、问题描述 饱和蒸气压的拟合 2、理论分析 题目要求是求模型方程ln p=A+B/(t+C)中的三个参数A、B和C的值。因为该模型为非线性关系，所以这是一个非线性拟合问题。 3、数值方法及程序设计分析 对于该问题可以采用MATLAB中的非线性拟合命令lsqcurvefit求解。求解分两步 第一步：建立拟合模型的函数。 function y=funantoine(x,t)，其中x是拟合输出参数，t是拟合数据的，y是拟合模型计算值。 第二步：主程序 先输入拟合数据，然后用x = lsqcurvefit(@funantoine,x0,t,p)进行非线性拟合，其中t，p是输入拟合数据，funantoine是拟合模型方程，x0时拟合参数的初始值，x是输出的拟合参数。那么A=x(1),B=x(2),C=x(3)。 源程序(1) %主程式 function antoine clear clc % %定义参数： % t:温度[oC] % p:温度t对应的纯液体的饱和蒸汽压[mmHg] % A,B,C:要求的参数 % global tp % %给定数据 % t=[8.86 50.0 76.22 20.59 60.78 85.92 16.31 56.95 80.00 26.89 67.13 91.78 26.1 63.33 81.30 35.19 74.03 97.69 40.25 69.57 14.55 49.07 78.89 103.65]; p=[42.6 268.3 671.9 77.28 402.4 906.06 62.4 347.1 756.2 103.64 500.7 1074.6 98.4 437.0 789.2 149.4 627.9 1268.0 183.0 540.9 57.41 261.8 732.1 1489.1]; 源程序(2) x0=[16 -3000 230]; %求出三个参数解 x = lsqcurvefit(@funantoine,x0,t,p) A=x(1); B=x(2); C=x(3); %disply disp(sprintf(\n 参数估计值为： A=%5.2f, B=%7.2f,and C=%6.2f,A,B,C)) t2=linspace(8,104,100) ; p2=funantoine(x,t2); plot(t,p,b*,t2,p2,r-) xlabel(temper