matlab实时显示每次计算记过,第六讲 计算方法在MATLAB语言中的体现

2.4.3 自适应Newton-Cotes 法求积quad8(建议使用quadl 代替quad8)

help quad8

QUAD8 Numerically evaluate integral, higher order method.

QUAD8 is obsolete. QUADL is its recommended replacement. 用Newton-Cotes 公式求以下积分

function f=fun1(x) f=exp(-x.^2/2);

quad8('fun1',0,3,1e-3) ans =1.24993044474109 quad8('fun1',0,3,1e-10) ans =1.24993044474154

2.4.4 自适应Gauss Lobatlo 法求积quadl

MATLAB 6.0以后，建议以新指令quadl 取代之quad8，一般而言，quad 较quadl 计算速度较快，但精确度稍差。

对于一些简单的函数，可直接將函数输入quadl 中求解，不需另外编写函数。 例如，

area1=quadl('1./(x.^3-2*x-5)',0,2)

area1 =

-0.4605

或者，利用inline 的指令输出函数，然后求积分分亦可。

f=inline('1./(x.^3-2*x-5)');

area2=quadl(f,0,2)

area2 =

-0.4605

例 数学函数?

-=

x

t dt e x erf 0

2

2

)(π

,由于其解析不可积，求其数值积分。

f=inline('2/sqrt(pi)*exp(-x.^2)','x'); y=quad(f,0,1.5)

y =0.966105186231726

syms x,y0=vpa(int(2/sqrt(pi)*exp(-x^2),0,1.5),60)

y0 =.966105146475310713936933729949905794996224943257461473285749

y=quad(f,0,1.5,1e-20)

Warning: Maximum function count exceeded; singularity likely. In quad at 106 y = 0.966060260015354

%设置高精度，但该方法失败。 y=quadl(f,0,1.5,1e-20)

?-

30

2

2dx

e

x