数值分析作业：逐次二分梯形公式与Romberg求积法程序设计及应用.pdfVIP

作业十四(第七章)：

2.编写使用逐次二分梯形求数值积分的程序，并计算下式

的定积分。

3.编写Romberg求积法的程序，并计算下式的定积分。

4.编写使用Richardson外推法求数值微分的程序，并计算f(x)=

sinx/x在x=2处的导数f(2)。

11

dx

01+x2

（

2）

程序

%复合梯形求函数f在区间[a,b]上的定积分

%函数名：f%积分下限：a%积分上限：b%积分

精度：eps%积分值：I%积分划分的子区间个数：

step

=00e‑4;%默认精度为

clear;clc;f=@(x)1/(1+x^2);a;b=1;eps=1.

0.1=0

0001n=1;h=(b‑a)/2;I;I2=

ssm(f),findsym(sym(f)),a)+

(sub(y

subs(sym(f),findsym(sym(f)),b))/h;whileabs(I2‑I1)epsn=n+1;

h=(b‑a)/n;

作业十四(第七章)：

2.编写用逐次二分梯形求数值积分的程序，并计算下式的定

积分。

3.编写Romberg求积法的程序，并计算下式的定积分。

4.编写用Richardson外推法求数值微分的程序，并计算f(x)=sinx/x

在x=2处的导数f(2)。

11

dx

01+x2

（2）

程序

%复合梯形求函数f在区间[a,b]上的定积分

%函数名：f

%积分下限：a

%积分上限：b

%积分精度：eps

%积分值：I

%积分划分的子区间个数：step

clear;clc;

f=@(x)1/(1+x^2);

a=0;

b=1;

eps=1.0e-4;

%默认精度为0.0001

n=1;

h=(b-a)/2;

I1=0;

I2=(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b))/h;

whileabs(I2-I1)eps

n=n+1;

h=(b-a)/n;

2=0=0‑1%第n次的复合梯形积分x=a+h*

I1=I2;I;fori:ni;%i=0和n‑1时，

分别代表积分区间的左右端点x1=x+h;I2=I2+

(h/2)*(subs(sym(f),findsym(sym(f)),x)+subs(sym(f),findsym(sym(f)),x1));

endendstep=nI=I2

程序结果

（3）

程序

clear;clc;%被积函数为f(x)=1/(1+x^2);积分区

间为[0,1]%误差界eps(‑5)b=1;a;h=

=10^=0

b‑a;eps(‑5);%误差界epskmax;%最

=10^=10

大递推次数m;fx;f=@(x)1/(1+x^2);

=0=0

I1=I2;

I2=0;

fori=0:n-1

%第n次的复合梯形积分

x=a+h*i;

%i=