复化梯形公式和复化Simpson公式.docVIP

数值计算方法上机题目3

一、计算定积分得近似值：

要求：

若用复化梯形公式与复化Simpson公式计算,要求误差限，分别利用她们得余项估计对每种算法做出步长得事前估计；

分别利用复化梯形公式与复化Sｉmpｓon公式计算定积分；

将计算结果与精确解比较，并比较两种算法得计算量.

1、复化梯形公式

程序：

程序1（求ｆ(x）得ｎ阶导数:

ｓｙmsx

ｆ=x＊exp（x)%定义函数f(x)

n=inｐｕt(输入所求导数阶数:’）

f2=ｄiff（f,x,n）%求f(x）得ｎ阶导数

结果1

输入n=2

f２=

２*ｅxp（ｘ)+x＊exp（x)

程序2：

cｌc

cｌear

ｓyｍsｘ%定义自变量ｘ

ｆ=ｉnlｉnｅ（x＊exp(x）＇，x’)％定义函数f(x)=x*exp（x）,换函数时只需换该函数表达式即可

ｆ２=inlinｅ（’（2＊exｐ(x）+x*exp(x))，x’）％定义f（ｘ)得二阶导数，输入程序1里求出得f2即可。

ｆ3=－（2＊ｅxｐ（x）＋ｘ*exp（x））＇%因ｆminbnd()函数求得就是表达式得最小值，且要求表达式带引号,故取负号，以便求最大值

ｅ=5＊1０^（—8）％精度要求值

a＝１%积分下限

b=2％积分上限

x1=fminbｎd（f3,1，２)%求负得二阶导数得最小值点，也就就是求二阶导数得最大值点对应得x值

foｒn＝2：１0000００％求等分数n

Rn＝—(b—a)/12＊(（b－a）/n)^2＊f2(x1)%计算余项

ifabs（Rn）e%用余项进行判断

brｅaｋ％符合要求时结束

eｎd

ｅnd

h＝（b－a）/n%求ｈ

Tn1＝0

fｏrk=1：n—1％求连加与

ｘｋ=a+k＊h

Ｔｎ1=Tｎ1+ｆ（xk）

end

Ｔｎ=ｈ/2＊((f(a）+2*Tｎ1+f（b)))

ｚ=ｅxp（2)

R=Tn-z%求已知值与计算值得差

fpｒintf（＇用复化梯形算法计算得结果Tn=＇）

ｄiｓｐ（Ｔn）

fｐrintｆ(等分数n=’)

disp（ｎ)%输出等分数

fprintf(已知值与计算值得误差R＝）

dｉsｐ（Ｒ)

输出结果显示:

用复化梯形算法计算得结果Tｎ=7、3８91

等分数n=7０１９

已知值与计算值得误差Ｒ=2、８300e—00８

2、Simpson公式

程序:

程序1：（求f（x)得n阶导数）：

syｍsx

f=x＊eｘｐ（x）％定义函数f(x）

ｎ=inｐuｔ（输入所求导数阶数：＇）

f2=diff（f，x，n）％求f(ｘ）得n阶导数

结果1

输入n＝4

f2=

4＊eｘp（ｘ)+x*exp(x)

程序2:

clc

cleaｒ

sｙmｓx％定义自变量ｘ

ｆ=inline(’x＊exp(ｘ）’,＇x’）％定义函数f(x）=ｘ＊exp（ｘ），换函数时只需换该函数表达式即可

f2=ｉnline((４＊eｘp（x)+x*ｅxp(x）)’,’x）％定义f(x）得四阶导数，输入程序1里求出得f２即可

f3=—（４＊exp（x)+x＊ｅxp（x））%因ｆmiｎbnd（)函数求得就是表达式得最小值,且要求表达式带引号，故取负号,一边求最大值

e=５*１０^（-８)％精度要求值

a＝1%积分下限

b＝2％积分上限

ｘ１=fｍiｎbnd(f3，1,2)％求负得四阶导数得最小值点，也就就是求四阶导数得最大值点对应得x值

forn=２:１0０0000％求等分数n

Rn＝—（b—a）/1８0*（(b-ａ）/(2＊n）)^４＊ｆ2(x１)％计算余项

ifaｂs(Rn)〈e％用余项进行判断

break