----------------------个人作业，如果有后辈的作业习题一致，可以参考学习，一起交流，请勿直接copy

··复合抛物线公式：

··龙贝格公式：

四、实验内容

------1

实验题目1中所用到的三种算法的matlab实现代码具体如下：

%复合梯形公式

function y=funct(f,n,a,b)

fi=f(a)+f(b);

h=(b-a)/n;

d=1;

for i=1:n-1

x=a+i*h;

fi=fi+2*f(x);

d=d+1;

end

f4=h/2*fi,d

%复合抛物线公式

function y=funcp(f,n,a,b)

h=(b-a)/n;

f1=0;

f2=0;

dd=1;

for i=1:n-1

dd=dd+1;

if  rem(i,2)~=0;

x1=a+i*h;

f1=f1+f(x1);

else

rem(i,2)==0;

x2=a+i*h;

f2=f2+f(x2)

end

end

f3=(h/3)*(f(a)+4*f1+2*f2+f(b)),dd

%龙贝格公式

function y=romberg(f,n,a,b)

z=zeros(n,n);

h=b-a;

z(1,1)=(h/2)*(f(a)+f(b));

f1=0;

for i=2:n

for  k=1:2^(i-2)

f1=f1+f(a+(k-0.5)*h);

end

z(i,1)=0.5*z(i-1,1)+0.5*h*f1;

h=h/2;

f1=0;

for  j=2:i

z(i,j)=z(i,j-1)+(z(i,j-1)-z(i-1,j-1))/(4^(j-1)-1);

end

end

z,n

------2

实验题目2中所用到的三种算法的matlab实现代码具体如下：

%采用三点公式法求函数func在x0处的二阶导数

function df= ThreePoint2(func,x0,type,h)

%函数名: func

%求导点: x0

%三点公式法的形式： type

%步长：h

%导数值： df

if nargin == 3

h = 0.1;

else if (nargin == 4 && h == 0.0)

disp(‘ h不能为0！‘);

return;

end

end

for i=1:5

y(i) = subs(sym(func), findsym(sym(func)), x0-2*h+i*h-h);

end

hd = 1/h/h;

switch type

case 1

df = (y(3) - 2*y(4) + y(5))*hd;      %用第1个公式求二阶导数

case 2

df = (y(2) - 2*y(3) + y(4))*hd;      %用第2个公式求二阶导数

case 3

df = (y(1) - 2*y(2) + y(3))*hd;      %用第3个公式求二阶导数

end

%采用四点公式法求函数func在x0处的二阶导数

function df= FourPoint2(func,x0,type,h)

%函数名: func

%求导点: x0

%四点公式法的形式： type

%步长：h

%导数值： df

if nargin == 3

h = 0.1;

else if (nargin == 4 && h == 0.0)

disp(‘ h不能为0！‘);

return;

end

end

for i=1:7

y(i) = subs(sym(func), findsym(sym(func)), x0-3*h+i*h-h);

end

hd = 1/6/h/h;

switch type

case 1

df = (12*y(4) - 30*y(5) + 24*y(6) - 6*y(7))*hd; %用第1个公式求二阶导数

case 2

df = (6*y(3) - 12*y(4) + 6*y(5))*hd; %用第2个公式求二阶导数

case 3

df = (12*y(4) - 30*y(3) + 24*y(2) - 6*y(1))*hd; %用第3个公式求二阶导数

end

%采用五点公式法求函数func在x0处的二阶导数

function df= FivePoint2(func,x0,type,h)

%函数名: func

%求导点: x0

%五点公式法的形式： type

%步长：h

%导数值： df

if nargin == 3

h = 0.1;

else if (nargin == 4 && h == 0.0)

disp(‘h不能为0！‘);

return;

end

end

for i=1:9

y(i) = subs(sym(func), findsym(sym(func)), x0-4*h+i*h-h);

end

hd = 1/12/h/h;

switch type

case 1

df = (35*y(5)-104* y(6)+114* y(7)-56* y(8)+11* y(9))*hd; %用第1个公式求二阶导数

case 2

df = (11*y(4)-20* y(5)+6* y(6)+4* y(7)- y(8))*hd; %用第1个公式求二阶导数

case 3

df = (-y(3)+16* y(4)-30* y(5)+16* y(6)- y(7))*hd; %用第1个公式求二阶导数

case 4

df = (-y(2)+4*y(3)-6*y(4)+20*y(5)+11*y(6))*hd; %用第1个公式求二阶导数

case 5

df = (35*y(5)-104* y(4)+114* y(3)-56* y(2)+11* y(1))*hd; %用第1个公式求二阶导数

end

------3

%采用五点公式法求函数func在x0处的二阶导数

function df= FivePoint2(func,x0,type,h)

%函数名: func

%求导点: x0

%五点公式法的形式： type

%步长：h

%导数值： df

if nargin == 3

h = 0.1;

else if (nargin == 4 && h == 0.0)

disp(‘h不能为0！‘);

return;

end

end

for i=1:9

y(i) = subs(sym(func), findsym(sym(func)), x0-4*h+i*h-h);

end

hd = 1/12/h/h;

switch type

case 1

df = (35*y(5)-104* y(6)+114* y(7)-56* y(8)+11* y(9))*hd; %用第1个公式求二阶导数

case 2

df = (11*y(4)-20* y(5)+6* y(6)+4* y(7)- y(8))*hd; %用第1个公式求二阶导数

case 3

df = (-y(3)+16* y(4)-30* y(5)+16* y(6)- y(7))*hd; %用第1个公式求二阶导数

case 4

df = (-y(2)+4*y(3)-6*y(4)+20*y(5)+11*y(6))*hd; %用第1个公式求二阶导数

case 5

df = (35*y(5)-104* y(4)+114* y(3)-56* y(2)+11* y(1))*hd; %用第1个公式求二阶导数

end

六、实验结果分析

在题目一中分别使用了复合梯形公式、复合抛物线公式、以及龙贝格公式对两个题目的积分进行了求解。就结果而言，通过实验结果中的各项数据对比可以明显看出：准确度为 ----龙贝格公式>复合抛物线公式>复合梯形公式；其中龙贝格公式近乎接近正确值，精度达到了极高的标准，相比较之下，复合梯度公式的精度就比较低了。故而，在不考虑计算的时间和空间效率时，使用龙贝格公式的精度最高，效果最好；如果加上时间空间来考虑，可以采用复合抛物线公式来求解。

2.

在题目二中分别使用了三点公式、四点公式、以及五点公式对两个题目的函数进行了求导，在题目三种采用了五点公式，在不同的步长设定下进行了求导。就结果而言，通过实验结果中的各项数据对比的表格可以明显看出： ·在相同的步长下，求解精度：五点公式>四点公式>三点公式;  ·在同一公式下，步长越长，求解的精度会越高；  ·在步长较低(h=0.1)时，三点公式所求得的导数值摆动较大，几乎无法使用； ·在取步长较高(h=0.001)时，采用五点公式，精度极高，几乎可以替代精确值使用。

---------------------------------------------------

排版依旧很别扭，截图更方便阅读，怎么说呢、、、、、凑合还能看吧

原文：http://www.cnblogs.com/nanashi/p/6690337.html