牛顿插值法试验报告

.

牛顿插值法一、实验目的：学会牛顿插值法，并应用算法于实际问题。

x?x)f(二、实验内容：给定函数，已知：

4832401.2)?.?1449138f(2.f.f(20)?1.414214(2.1)

549193.)?1f(2.4516575(f2.3)?1.

三、实验要求：以此作为函数2.15插值多项式在处的值，用牛顿插值法求4

次Newton（ 1）2.15?N(2.15)。在MATLAB中用内部函数ezplot绘制出的近似值4次Newton插值多项式的函数图形。

（2）在MATLAB中用内部函数ezplot可直接绘制出以上函数的图形，并与作出的4次Newton插值多项式的图形进行比较。

四、实验过程：

1、编写主函数。打开Editor编辑器，输入Newton插值法主程序语句：

function [y,L]=newdscg(X,Y,x)

n=length(X); z=x; A=zeros(n,n);A(:,1)=Y';s=0.0; p=1.0;

for j=2:n

for i=j:n

A(i,j)=(A(i,j-1)- A(i-1,j-1))/(X(i)-X(i-j+1));

end

end

C=A(n,n);

for k=(n-1):-1:1

C=conv(C,poly(X(k)));

d=length(C);C(d)=C(d)+A(k,k);

end

y(k)= polyval(C, z);

L(k,:)=poly2sym(C);

0 / 3

.

%%%%%%%%%%%%%%%%%%

t=[2,2.1,2.2,2.3,2.4];

fx=sqrt(t);

wucha=fx-Y;

以文件名newdscg.m保存。

2、运行程序。

（1）在MATLAB命令窗口输入：

>> X=[2,2.1,2.2,2.3,2.4]; Y

=[1.414214,1.449138,1.483240,1.516575,1.549193];

x=2.15;[y,P]=newdscg(X,Y,x)

回车得到：

y =1.4663

wucha =1.0e-06 *

-0.4376 -0.3254 -0.3026 0.0888 0.3385

P = - (4803839603609061*x^4)/2305843009213693952 + (7806239355294329*x^3)/288230376151711744 - (176292469178709*x^2)/1125899906842624 + (1624739243112817*x)/2251799813685248 + 1865116246031207/4503599627370496

（2）在MATLAB命令窗口输入：

>> v=[0,6,-1,3];

>> ezplot(P),axis(v),grid

>> hold on

>> x=0:0.1:6;

>> yt=sqrt(x);plot(x,yt,':')

>> legend('插值效果','原函数')

>> xlabel('X')

>> ylabel('Y')

>>title('Newton插值与原函数比较')

回车即可得到图像1-1。

1 / 3

.

牛顿插值效果图1-1五、实验结果分析：2.15Newton插值多项式在用牛顿插值法求4次的程序可得，由上运行（1）2.15 N(2.15)=1.4663。处函数的近似值由在MATLAB中用内部函数ezplot直接绘制出出的4次Newton插值图形与原函数的图形知，4次Newton插值图形在区间[0，1]与区间[4，5]内与原函数存在一定的偏差，而在区间[1,4]内误差在10的-6次方，这个精度是非常高的。因此，在计算区间[1,4]内的值时结果是比较准确的。

2 / 3

牛顿插值法的C语言编程

Newton 插值 Newton 插值函数 Newton 插值函数是用差商作为系数，对于01,,,n x x x …这1n +个点，其一般形式为： 00100120101011()[][,]()[,,]()()[,,,]()()() n n n N x f x f x x x x f x x x x x x x f x x x x x x x x x ?=+?+??++???…………对于011,,,n x x x ?…这n 个点， 100100120101012()[][,]()[,,]()()[,,,]()()() n n n N x f x f x x x x f x x x x x x x f x x x x x x x x x ??=+?+??++???…………差商的定义 若已知函数()f x 在点(0,1,2,,)i x i n =???处的函数值()i f x 。则称： 00[]()f x f x =为函数()f x 在点0x 的0阶差商； 100110 [][] [,]f x f x f x x x x ?= ?为函数()f x 关于01,x x 的1阶差商； 120101220 [,][,] [,,]f x x f x x f x x x x x ?= ?为函数()f x 过点012,,x x x 的2阶差商； 依此类推，一般地称 121012101210 [,,,,][,,,,] [,,,,,]k k k k k k k f x x x x f x x x x f x x x x x x x ??????????????= ?为函数()f x 关于01,,,k x x x ???的 k 阶差商。 表1 差商表 i x ()i f x 1阶差商 2阶差商 3阶差商 4阶差商 0x 1x 2x 3x 4x …… 0()f x 1()f x 2()f x 3()f x 4() f x …… 01[,]f x x 12[,]f x x 23[,]f x x 34[,]f x x …… 012[,,]f x x x 123[,,]f x x x 234[,,] f x x x …… 0123[,,,]f x x x x 1234[,,,] f x x x x …… 01234[,,,,]f x x x x x …… 根据Newton 插值函数编写的C 语言编程 根据Newton 插值函数并对照上面的差商表，可编写出Newton 插值法的C 语言程序如下： #include #include #include double NewtonInterpolation(double *x,double *y,int n,double xx,double *pyy) {

大学物理仿真实验报告牛顿环法测曲率半径

大学物理仿真实验报告－牛顿环法测曲率半径

————————————————————————————————作者: ————————————————————————————————日期:

大学物理仿真实验报告 实验名称 牛顿环法测曲率半径 班级: 姓名: 学号: 日期:

牛顿环法测曲率半径 实验目的 1．学会用牛顿环测定透镜曲率半径。 2.正确使用读书显微镜，学习用逐差法处理数据。 实验原理 如下图所示，在平板玻璃面ＤＣF上放一个曲率半径很大的平凸透镜ACB，C点为接触点,这样在ACB和DCF之间,形成一层厚度不均匀的空气薄膜,单色光从上方垂直入射到透镜上，透过透镜,近似垂直地入射于空气膜。分别从膜的上下表面反射的两条光线来自同一条入射光线，它们满足相干条件并在膜的上表面相遇而产生干涉，干涉后的强度由相遇的两条光线的光程差决定,由图可见,二者的光程差等于膜厚度ｅ的两倍。此外，当光在空气膜的上表面反射时，是从光密媒质射向光疏媒质，反射光不发生相位突变，而在下表面反射时，则会发生相位突变，即在反射点处,反射光的相位与入射光的相位之间相差π，与之对应的光程差为λ/2 ,所以相干的两条光线还具有λ/2的附加光程差，总的光程差为 (1) 当?满足条件（2） 时，发生相长干涉，出现第Ｋ级亮纹，而当 （k = 0,１,２…) (３) 时,发生相消干涉,出现第k级暗纹。因为同一级条纹对应着相同的膜厚,所以干涉条纹是一组等厚度线。可以想见，干涉条纹是一组以C点为中心的同心圆，这就是所谓的牛顿环。 如图所示,设第ｋ级条纹的半径为，对应的膜厚度为，则

（４） 在实验中，R的大小为几米到十几米,而的数量级为毫米,所以Ｒ＞> eｋ,ｅk ２相对于2Re 是一个小量，可以忽略,所以上式可以简化为 k （５） 如果rｋ是第ｋ级暗条纹的半径,由式（１)和（３）可得 （6) 代入式（5）得透镜曲率半径的计算公式 （７） 对给定的装置，R为常数,暗纹半径 （８） 和级数ｋ的平方根成正比,即随着k的增大,条纹越来越细。 同理,如果r k是第k级明纹，则由式（１)和(2)得 （９） 代入式(5)，可以算出 （1０)

插值法数值上机实验报告

插值法数值上机实验报告 实验题目： 利用下列条件做插值逼近，并与R (x) 的图像比较 考虑函数：R x y=1 1+x2 （1）用等距节点X i=?5+i,i=0,1,...,10.给出它的10次Newton插值多项式的图像； π),i=0,1,...,20.给出它的20次Lagrange插值多项式（2）用节点X i=5cos(2i+1 42 的图像； （3）用等距节点X i=?5+i,i=0,1,...,10.给出它的分段线性插值函数的图像；（4）用等距节点X i=?5+i,i=0,1,...,10.给出它的三次自然样条插值函数的图像； （5）用等距节点X i=?5+i,i=0,1,...,10.给出它的分段三次Hermite插值函数的图像； 实验图像结果：

实验结果分析： 1.为了验证Range现象，我还特意做了10次牛顿插值多项式和20次牛顿插值多项式的对比图像，结果如下图（图对称，只截取一半） 可以看出，Range现象在高次时变得更加明显。这也是由于高次多项式在端点处的最值随次数的变大很明显。可以料定高次多项式在两侧端点处剧烈震荡，在更小的间距内急剧上升然后下降，Range现象非常明显。

2.分析实验（2）的结果，我们会惊讶地发现，由于取21个点逼近，原本预料的Range现象会很明显，但这里却和f(x)拟合的很好。（即下图中Lagrange p(x)的图像）。可是上图中取均匀节点的20次牛顿多项式逼近的效果在端点处却很差。料想是由于节点X i=5cos2i+1 42 π ,i=0,1,...,20 取得很好。由书上第五章的 知识，对于函数y=1 1+x ，y 1 2对应的cherbyshev多项式的根恰好为X i= 5cos2i+1 42 π ,i=0,1,...,20 。由于所学限制，未能深入分析。 （3）比较三次样条插值图像和Hermit插值图像对原函数图像的逼近情形。见下图：

计算方法第二章方程求根上机报告

实验报告名称 班级：学号：姓名：成绩： 1实验目的 1）通过对二分法与牛顿迭代法作编程练习与上级运算，进一步体会二分法与牛顿迭代法的不同特点。 2）编写割线迭代法的程序，求非线性迭代法的解，并与牛顿迭代法。 2 实验内容 用牛顿法和割线法求下列方程的根 x^2-e^x=0; x*e^x-1=0; lgx+x-2=0; 3实验步骤 1)根据二分法和牛顿迭代法，割线法的算法编写相应的求根函数； 2)将题中所给参数带入二分法函数，确定大致区间； 3)用牛顿迭代法和割线法分别对方程进行求解； 3 程序设计 牛顿迭代法x0=1.0; N=100; k=0; eps=5e-6; delta=1e-6; while(1) x1=x0-fc1(x0)/fc2(x0); k=k+1; if k>N disp('Newmethod failed')

break end if(abs(x1-x0)=delta) c=x1; x1=cutnext(x0,x1); x0=c; %x0 x1μYí?μ?μ?x1 x2 è?è?±￡′??úx0 x1 end k=k+1; if k>N disp('Cutline method failed') break; end if(abs(x1-x0)

牛顿插值法原理及应用

牛顿插值法 插值法是利用函数f (x)在某区间中若干点的函数值，作出适当的特定函数，在这些点上取已知值，在区间的其他点上用这特定函数的值作为函数f (x)的近似值。如果这特定函数是多项式，就称它为插值多项式。当插值节点增减时全部插值基函数均要随之变化，这在实际计算中很不方便。为了克服这一缺点，提出了牛顿插值。牛顿插值通过求各阶差商，递推得到的一个公式： f(x)=f[x0]+f[x0,x1](x-x0)+f[x0,x1,x2](x-x0)(x-x1)+...f[x0,...xn](x-x0 )...(x-xn-1)+Rn(x)。 插值函数 插值函数的概念及相关性质[1] 定义：设连续函数y-f(x) 在区间[a,b]上有定义，已知在n+1个互异的点 x0,x1,…xn上取值分别为y0,y1,…yn （设a≤ x1≤x2……≤xn≤b)。若在函数类中存在以简单函数P(x) ，使得P(xi)=yi,则称P(x) 为f(x)的插值函数. 称x1,x2,…xn 为插值节点，称[a,b]为插值区间。 定理：n次代数插值问题的解存在且唯一。

牛顿插值法C程序 程序框图#include void main() { float x[11],y[11][11],xx,temp,newton; int i,j,n; printf("Newton插值:\n请输入要运算的值:x="); scanf("%f",&xx); printf("请输入插值的次数(n<11):n="); scanf("%d",&n); printf("请输入%d组值:\n",n+1); for(i=0;i

等厚干涉牛顿环实验报告材料97459

等厚干涉——牛顿环 等厚干涉是薄膜干涉的一种。薄膜层的上下表面有一很小的倾角是，从光源发出的光经上下表面反射后在上表面附近相遇时产生干涉，并且厚度相同的地方形成同一干涉条纹，这种干涉就叫等厚干涉。其中牛顿环是等厚干涉的一个最典型的例子，最早为牛顿所发现，但由于他主张微粒子学说而并未能对他做出正确的解释。光的等厚干涉原理在生产实践中育有广泛的应用，它可用于检测透镜的曲率，测量光波波长，精确地测量微笑长度、厚度和角度，检验物体表面的光洁度、平整度等。 一.实验目的 (1)用牛顿环观察和分析等厚干涉现象； (2)学习利用干涉现象测量透镜的曲率半径； 二.实验仪器 读数显微镜钠光灯牛顿环仪

三. 实验原理 牛顿环装置是由一块曲率半径较大的平凸面放在一块光学玻璃平板（平镜）上构成的，如图。平凸透镜的凸面与玻璃平板之间的空气层厚度从中心到边缘逐渐增加，若以平行单光垂直照射到牛顿环上，则经空气层上、下表面反射的两光束存在光程差，他们在平凸透镜的凸面相遇后，将发生干涉。从透镜上看到的干涉花样是以玻璃接触点为中心的一系列明暗相间的圆环，称为牛顿环。同一干涉环上各处的空气层厚度是相同的，因此他属于等厚干涉。 图2 图3 由图2可见，若设透镜的曲率半径为R ，与接触点O 相距为r 处空气层的厚度为d ，其几何关系式为 2222222)(r d Rd R r d R R ++-=+-= 由于r R >>，可以略去d 2得

R r d 22 = （1） 光线应是垂直入射的，计算光程差时还要考虑光波在平玻璃上反射会有半波损失，，从而带来2λ的附加程差，所以总光程差为 2 2λ + =?d （2） 所以暗环的条件是 2 ) 12(λ +=?k （3） 其中 3,2,1,0=k 为干涉暗条纹的级数。综合（1）（2）（3）式可得第可k 级暗环的半径为 λkR r k =2 （4） 由式（4）可知，如果单色光源的波长λ已知，测出第m 级的暗环半径r m,,即可得出平图透镜的曲率半径R ；反之，如果R 已知，测出r m 后，就可计算出入射单色光波的波长λ。但是用此测量关系式往往误差很大，原因在于凸面和平面不可能是理想的点接触；接触压力会引起局部形变，使接触处成为一个圆形平面，干涉环中心为一暗斑。或者空气间隙层有了灰尘，附加了光程差，干涉环中心为一亮（或暗）斑，均无法确定环的几何中心。实际测量时，我们可以通过测量距中心较远的两个暗环半径r m 和r n 的平方差来计算曲率半径R 。因为 λMR r m =2 λnR r n =2 两式相减可得 λ)(22n m R r r n m -=-

牛顿插值法的分析与应用

牛顿插值法的分析与应用 学生： 班级： 学号： ： 指导教师： 成绩：

一．定义 )(x f 关于i x 的零阶差商 )(][i i x f x f = )(x f 关于i x ，j x 的一阶差商 i j i j j i x x x f x f x x f --= ][][],[ 依次类推，)(x f 关于i x ，1+i x ，……，k i x +的k 阶差商 i k i k i i k i i k i i i x x x x f x x f x x x f --= +-+++++] ,,[],,[],,,[111 二. 牛顿插值多项式 设给定的n+1个互异点))(,(k k x f x ，n k ,,1,0 =，j i x x ≠，j i ≠， 称满足条件 )()(k k n x f x N =，n k ,,1,0 = 的n 次多项式 )()](,,,[)](,[][)(10100100---++-+=n n n x x x x x x x f x x x x f x f x N 为Newton 插值多项式，称 ],[,)(],,,[)()()(0 10b a x x x x x x f x N x f x E n j j n n ∈-=-=∏= 为插值余项。 三.算法 步骤1：输入节点（xj ，yj ），精度ξ，计值点xx ，f0→p ，1→T ，1→i ； 步骤2：对k=1，2，……，i 依次计算k 阶均差 f[xi-k,xi-k+1,…,xi] = (f[xi-k+1,…,xi]- f[xi-k,…,xi])/( xi -xi-k ) 步骤3：（1）、若| f[x1,…,xi]- f[x0,…,xi-1]|< ξ，则p 为最终结果Ni-1(x)，余项Ri-1= f[x0,…,xi](xx-xi-1)T 。 （2）、否则(xx-xi-1)*T →T ，p+ f[x0,…,xi]*T →p ，转步骤4。 步骤4：若i

等厚干涉牛顿环实验报告

等厚干涉牛顿环实验报告 Prepared on 22 November 2020

等厚干涉——牛顿环等厚干涉是薄膜干涉的一种。薄膜层的上下表面有一很小的倾角是，从光源发出的光经上下表面反射后在上表面附近相遇时产生干涉，并且厚度相同的地方形成同一干涉条纹，这种干涉就叫等厚干涉。其中牛顿环是等厚干涉的一个最典型的例子，最早为牛顿所发现，但由于他主张微粒子学说而并未能对他做出正确的解释。光的等厚干涉原理在生产实践中育有广泛的应用，它可用于检测透镜的曲率，测量光波波长，精确地测量微笑长度、厚度和角度，检验物体表面的光洁度、平整度等。 一. 实验目的 (1)用牛顿环观察和分析等厚干涉现象； (2)学习利用干涉现象测量透镜的曲率半径； 二. 实验仪器 读数显微镜钠光灯牛顿环仪 三. 实验原理 牛顿环装置是由一块曲率半径较大的平凸面放在一块光 学玻璃平板（平镜）上构成的，如图。平凸透镜的凸面与玻 璃平板之间的空气层厚度从中心到边缘逐渐增加，若以平行单光垂直照射到牛顿环上，则经空气层上、下表面反射的两光束存在光程差，他们在平

凸透镜的凸面相遇后，将发生干涉。从透镜上看到的干涉花样是以玻璃接触点为中心的一系列明暗相间的圆环，称为牛顿环。同一干涉环上各处的空气层厚度是相同的，因此他属于等厚干涉。 图2 图3 由图2可见，若设透镜的曲率半径为R ，与接触点O 相距为r 处空气层的厚度为d ，其几何关系式为 由于r R >>，可以略去d 2得 R r d 22 = （1） 光线应是垂直入射的，计算光程差时还要考虑光波在平玻璃上反射会有半波损失，，从而带来2λ的附加程差，所以总光程差为 2 2λ + =?d （2） 所以暗环的条件是 2 ) 12(λ +=?k （3） 其中 3,2,1, 0=k 为干涉暗条纹的级数。综合（1）（2）（3）式可得第可k 级暗环的半径为 λkR r k =2 （4） 由式（4）可知，如果单色光源的波长λ已知，测出第m 级的暗环半径r m,,即可得出平图透镜的曲率半径R ；反之，如果R 已知，测出r m 后，就可计算出入射单色光波的波长λ。但是用此测量关系式往往误差很大，原因在于凸面和平面不可能是理想的点接触；接触压力会引起局部形变，使接触处成为一个圆形平面，干涉环中心为一暗斑。或者空气间隙层有了灰尘，附加了光程差，干涉环中心为一亮（或暗）斑，均无法确定环的几何

数值分析实验报告记录

数值分析实验报告记录

————————————————————————————————作者：————————————————————————————————日期：

数值分析实验报告 （第二章） 实验题目： 分别用二分法、牛顿迭代法、割线法、史蒂芬森迭代法求方程 的根，观察不同初始值下的收敛性，并给出结论。 问题分析： 题目有以下几点要求： 1.不同的迭代法计算根，并比较收敛性。 2.选定不同的初始值，比较收敛性。 实验原理： 各个迭代法简述 二分法：取有根区间的重点，确定新的有根区间的区间长度仅为区间长度的一版。对压缩了的有根区间重复以上过程，又得到新的有根区间，其区间长度为的一半，如此反复，……，可得一系列有根区间，区间收敛到一个点即为根。 牛顿迭代法：不动点迭代法的一种特例，具有局部二次收敛的特性。迭代格式为 割线法：是牛顿法的改进，具有超线性收敛的特性，收敛阶为1.618. 迭代格式为 史蒂芬森迭代法：采用不动点迭代进行预估校正。至少是平方收敛的。迭代格式为 这里可采用牛顿迭代法的迭代函数。 实验内容：

1.写出该问题的函数代码如下： function py= f(x) syms k; y=(k^2+1)*(k-1)^5; yy=diff(y,k); py(1)=subs(y,k,x); py(2)=subs(yy,k,x); end 2.分别写出各个迭代法的迭代函数代码如下： 二分法： function y=dichotomie(a,b,e) i=2; m(1)=a; while abs(a-b)>e t=(a+b)/2; s1=f(a); s2=f(b); s3=f(t); if s1(1)*s3(1)<=0 b=t; else a=t; end m(i)=t; i=i+1; end y=[t,i+1,m]; end 牛顿迭代法： function y=NewtonIterative(x,e) i=2; en=2*e;m(1)=x; while abs(en)>=e s=f(x); t=x-s(1)/s(2); en=t-x; x=t; m(i)=t; i=i+1; end y=[x,i+1,m]; end 牛顿割线法： function y=Secant(x1,x2,e) i=3; m(1)=x1,m(2)=x2; while abs(x2-x1)>=e s1=f(x1); s2=f(x2); t=x2-(x2-x1)*s2(1)/(s2(1)-s1( 1)); x1=x2; x2=t; m(i)=t; i=i+1; end

数值计算(二分法、简单迭代法、Newton迭代法、弦截法(割线法、双点弦法))

本科生实验报告 实验课程数值计算方法 学院名称信息科学与技术学院 专业名称计算机科学与技术 学生姓名 学生学号 指导教师 实验地点 实验成绩 二〇一六年五月二〇一六年五月

实验一非线性方程求根 1.1问题描述 实验目的：掌握非线性方程求根的基本步骤及方法，。 实验内容：试分别用二分法、简单迭代法、Newton迭代法、弦截法(割线法、双点弦法)，求x5-3x3+x-1= 0 在区间 [-8,8]上的全部实根，误差限为10-6。 要求：讨论求解的全过程，对所用算法的局部收敛性，优缺点等作分析及比较， 第2章算法思想 2.1二分法 思想：在函数的单调有根区间内，将有根区间不断的二分，寻找方程的解。 步骤： 1.取中点mid=(x0+x1)/2 2.若f(mid)=0,则mid为方程的根，否则比较与两端的符号，若与 f(x0) 异号，则根在[x0,mid]之间，否则在[mid,x1]之间。 3并重复上述步骤，直达达到精度要求，则mid为方程的近似解。

2.2 简单迭代法 思想：迭代法是一种逐次逼近的方法，它是固定公式反复校正跟的近似值，使之逐步精确，最后得到精度要求的结果。 步骤：1.构造迭代公式f(x)，迭代公式必须是收敛的。 2.计算x1,x1=f(x0). 3.判断|x1-x0|是否满足精度要求，如不满足则重复上述步骤。 4．输出x1,即为方程的近似解。

开始 输入x0,e X1=f(x0)|x1-x0|

牛顿环实验报告

北京师范大学珠海分校大学物理实验报告 实验名称：牛顿环实验测量 学院工程技术学院 专业测控技术与仪器 学号 1218060075 姓名钟建洲 同组实验者 1218060067余浪威 1218010100杨孟雄 2013 年 1 月 17日

实验名称 牛顿环实验测量 一、实验目的 1.观察牛顿环干涉现象条纹特征； 2.学习用光的干涉做微小长度的测量； 3.利用牛顿环干涉测量平凸透镜的曲率半径； 4.通过实验掌握移测显微镜的使用方法 二、实验原理 在一块平面玻璃上安放上一焦距很大的平凸透镜，使其凸面与平面相接触，在接触点 o 附近就形成一层空 气膜。当用一平行的准单色光垂直照射时，在空气膜上表面反射的光束和下表面反射的光束在膜上表面相遇相干，形成以 o 为圆心的明暗相间的环状干涉图样，称为牛顿环。如果已知入射光波长，并测得第 k 级 暗环的半径 r k ，则可求得透镜的曲率半径 R 。但 实际测量时，由于透镜和平面玻璃接触时，接触点有压力产生形变或有微尘产生附加光程差，使得干涉条纹的圆心和环级确定困难。第m 环与第n 环 用直径 D m 、 D n 。 () λ n m n D m D R +-= 42 2此为计算 R 用的公式，它与附加厚度、

圆心位置、绝对级次无关，克服了由这些因素带来的系统误差，并且D m 、 D n 可以是弦长。 三、实验内容与步骤 用牛顿环测量透镜曲率半径 (1).按图布置好实验器材，使用单色扩展光源，将牛顿环装置放在读数显微镜工作台毛玻璃中央，并使显微镜筒正对牛顿环装置中心。 (2).调节读数显微镜。 1.调节目镜，使分划板上的十字刻度线清晰可见，并转动目镜，使十字刻度线的横线与显微镜筒的移动方向平行。 2.调节45度反射镜，使显微镜视觉中亮度最大，这时基本上满足入射光垂直于待测量透镜的要求。 1.转动手轮A，使显微镜平移到标尺中部，并调节调焦手轮B，使物镜接近牛顿环装置表面。 2.对显微镜调焦。缓慢地转动调焦手轮B，使显微镜筒由下而上移动进行调焦，直到从目镜中清楚地看到牛顿环干涉条纹且无视差为止；然后移动牛顿环装置，使目镜中十字刻度线交点与牛顿环中心重合 (1).观察条纹的特征。 观察各级条纹的粗细是否一致，其间距有无差异，并做出解释。观察牛顿环中心是亮斑还是暗斑？ (2).测量暗环的直径 转动读数显微镜的读数鼓轮，同时在目镜中观察，使十字刻度线由牛顿环中心缓慢地向一侧移动到43环；然后再回到第42环。自42环起，单方向移动十字刻度，每移3环读数一——直到测量完成另一侧的第42环。并将所测量的第42环到第15环各直径的左右两边的读数记录在表格内。 四、数据处理与结果 1.求透镜的曲率半径。 测出第15环到第42环暗环的直径，取m-n=15,用逐差法求出暗环的直径平方 差的平均值，按算出透镜的曲率半径的平均值R。 R1＝(d422-d272)/[4(42-27]λ= 895.85 mm R2＝(d392-d242)/[4(39-24]λ= 896.97 mm R3＝(d362-d212)/(4(36-21)λ= 887.94mm R4＝(d332-d182)/(4(33-18)λ= 893.30mm

牛顿插值法试验报告

. 牛顿插值法一、实验目的：学会牛顿插值法，并应用算法于实际问题。 x?x)f(二、实验内容：给定函数，已知： 4832401.2)?.?1449138f(2.f.f(20)?1.414214(2.1) 549193.)?1f(2.4516575(f2.3)?1. 三、实验要求：以此作为函数2.15插值多项式在处的值，用牛顿插值法求4 次Newton（ 1）2.15?N(2.15)。在MATLAB中用内部函数ezplot绘制出的近似值4次Newton插值多项式的函数图形。 （2）在MATLAB中用内部函数ezplot可直接绘制出以上函数的图形，并与作出的4次Newton插值多项式的图形进行比较。 四、实验过程： 1、编写主函数。打开Editor编辑器，输入Newton插值法主程序语句： function [y,L]=newdscg(X,Y,x) n=length(X); z=x; A=zeros(n,n);A(:,1)=Y';s=0.0; p=1.0; for j=2:n for i=j:n A(i,j)=(A(i,j-1)- A(i-1,j-1))/(X(i)-X(i-j+1)); end end C=A(n,n); for k=(n-1):-1:1 C=conv(C,poly(X(k))); d=length(C);C(d)=C(d)+A(k,k); end y(k)= polyval(C, z); L(k,:)=poly2sym(C); 0 / 3 . %%%%%%%%%%%%%%%%%% t=[2,2.1,2.2,2.3,2.4]; fx=sqrt(t); wucha=fx-Y; 以文件名newdscg.m保存。 2、运行程序。 （1）在MATLAB命令窗口输入： >> X=[2,2.1,2.2,2.3,2.4]; Y =[1.414214,1.449138,1.483240,1.516575,1.549193]; x=2.15;[y,P]=newdscg(X,Y,x) 回车得到：

用牛顿迭代法求解非线性方程

数值分析实验报告（一） 实验 名称 用牛顿迭代法求解非线性方程实验时间2011年11 月19日姓名班级学号成绩 一、实验目的 1．了解求解非线性方程的解的常见方法。 2．编写牛顿迭代法程序求解非线性方程。 二、实验内容 分别用初值 0.01 x=， 10 x=和 300 x=求113，要求精度为5 10-。 三、实验原理 设113 x=，则21130 x-=，记f(x)= 2113 x-，问题便成为了求2x -113=0的正根； 用牛顿迭代公式得 2 1 113 2 k k k k x x x x + - =-，即 1 1113 () 2 k k k x x x + =+（其中k=0，1，2，3，…,） 简单推导 假设f(x)是关于X的函数: 求出f(x)的一阶导，即斜率: 简化等式得到: 然后利用得到的最终式进行迭代运算直至求到一个比较精确的满意值。 如果f函数在闭区间[a,b]内连续，必存在一点x使得f(x) = c，c是函数f在闭区间[a,b]内的一点 我们先猜测一X初始值，然后代入初始值，通过迭代运算不断推进，逐步靠近精确值，直到得到我们主观认为比较满意的值为止。 回到我们最开始的那个”莫名其妙”的公式，我们要求的是N的平方根，令x2 = n，假设一关

于X的函数f(x)为: f(X) = X2 - n 求f(X)的一阶导为: f'(X) = 2X 代入前面求到的最终式中: X k+1 = X k - (X k 2 - n)/2X k 化简即得到我们最初提到求平方根的迭代公式: 四、实验步骤 1．根据实验题目，给出题目的C程序。 当初值为0.01、10、300时，即x=0.01,10,300 分别应用程序： #include "stdio.h" int main() { float number; printf("Please input the number:"); scanf("%f", &number); float x=1; int i; for (i=0;i<1000;i++) { x = (x + number/x)/2; } printf("The square root of %f is %8.5f\n", number ,x); } 得出结果 2．上机输入和调试自己所编的程序。 当x=0.01时，结果为：10.63015 x=10时，结果为：10.63015 x=300时，结果也为：10.63015 3．实验结果分析。 当初值取0.01、10、300时取不同的初值得到同样的结果10.63015。 五、程序

数值分析实验报告-插值、三次样条(教育教学)

实验报告：牛顿差值多项式&三次样条 问题：在区间[-1,1]上分别取n=10、20用两组等距节点对龙格函数2 1()25f x x 作多项式插值及三次样条插值，对每个n 值，分别画出插值函数及()f x 的图形。 实验目的：通过编程实现牛顿插值方法和三次样条方法，加深对多项式插值的理解。应用所编程序解决实际算例。 实验要求： 1． 认真分析问题，深刻理解相关理论知识并能熟练应用； 2． 编写相关程序并进行实验； 3． 调试程序，得到最终结果； 4． 分析解释实验结果； 5． 按照要求完成实验报告。 实验原理： 详见《数值分析 第5版》第二章相关内容。 实验内容： （1）牛顿插值多项式 1.1 当n=10时： 在Matlab 下编写代码完成计算和画图。结果如下： 代码： clear all clc x1=-1:0.2:1; y1=1./(1+25.*x1.^2); n=length(x1); f=y1(:); for j=2:n for i=n:-1:j f(i)=(f(i)-f(i-1))/(x1(i)-x1(i-j+1)); end end syms F x p ; F(1)=1;p(1)=y1(1); for i=2:n F(i)=F(i-1)*(x-x1(i-1)); p(i)=f(i)*F(i);

end syms P P=sum(p); P10=vpa(expand(P),5); x0=-1:0.001:1; y0=subs(P,x,x0); y2=subs(1/(1+25*x^2),x,x0); plot(x0,y0,x0,y2) grid on xlabel('x') ylabel('y') P10即我们所求的牛顿插值多项式，其结果为：P10（x）=-220.94*x^10+494.91*x^8-9.5065e-14*x^7-381.43*x^6-8.504e-14*x^5+123.36*x^4+2.0202e-1 4*x^3-16.855*x^2-6.6594e-16*x+1.0 并且这里也能得到该牛顿插值多项式的在[-1，1]上的图形，并和原函数进行对比（见Fig.1）。 Fig.1 牛顿插值多项式（n=10）函数和原函数图形 从图形中我们可以明显的观察出插值函数在两端点处发生了剧烈的波动，产生了极大的误差，即龙格现象，当n=20时，这一现象将更加明显。 1.2 当n=20时： 对n=10的代码进行修改就可以得到n=20时的代码。将“x1=-1:0.2:1;”改为“x1=-1:0.1:1;”即可。运行程序，我们得到n=20时的牛顿插值多项式，结果为：P20(x)= 260188.0*x^20 - 1.0121e6*x^18 + 2.6193e-12*x^17 + 1.6392e6*x^16 + 2.248e-11*x^15 - 1.4429e6*x^14 - 4.6331e-11*x^13 + 757299.0*x^12 + 1.7687e-11*x^11 - 245255.0*x^10 + 2.1019e-11*x^9 + 49318.0*x^8 + 3.5903e-12*x^7 - 6119.2*x^6 - 1.5935e-12*x^5 + 470.85*x^4 + 1.3597e-14*x^3 - 24.143*x^2 - 1.738e-14*x + 1.0 同样的，这里得到了该牛顿插值多项式的在[-1，1]上的图形，并和原函数进行对比（见Fig.2）。

大学物理仿真实验报告 牛顿环

大学物理仿真实验报告 实验名称：牛顿环法测曲率半径实验日期： 专业班级： 姓名：学号： 教师签字：________________ 一、实验目的 1.学会用牛顿环测定透镜曲率半径。 2.正确使用读书显微镜，学习用逐差法处理数据。 二、实验仪器 牛顿环仪，读数显微镜，钠光灯，入射光调节架。 三、实验原理 如图所示，在平板玻璃面DCF上放一个曲率半径很大的平 凸透镜ACB，C点为接触点，这样在ACB和DCF之间，形 成一层厚度不均匀的空气薄膜，单色光从上方垂直入射到 透镜上，透过透镜，近似垂直地入射于空气膜。分别从膜 的上下表面反射的两条光线来自同一条入射光线，它们满 足相干条件并在膜的上表面相遇而产生干涉，干涉后的强 度由相遇的两条光线的光程差决定，由图可见，二者的光 程差等于膜厚度e的两倍，即 此外，当光在空气膜的上表面反射时，是从光密媒质射向光疏媒质，反射光不发生相位突变，而在下表面反射时，则会发生相位突变，即在反射点处，反射光的相位与入射光的相位之间相差π，与之对应的光程差为λ/2 ，所以相干的两条光线还具有λ/2的附加光程差，总的光程差为（1） 当?满足条件（2）时，发生相长干涉，出现第K级亮纹，而当 （k = 0,1,2…）（3）时，发生相消干涉，出现第k级暗纹。因为

同一级条纹对应着相同的膜厚，所以干涉条纹是一组等厚度线。可以想见，干涉条纹是一组以C点为中心的同心圆，这就是所谓的牛顿环。 如图所示，设第k级条纹的半径为，对应的膜厚度为，则 （4） 在实验中，R的大小为几米到十几米，而的数量级为毫米，所以R >> e k，e k2相对于2Re k是一个小量，可以忽略，所以上式可以简化为 （5） 如果r k是第k级暗条纹的半径，由式（1）和（3）可得 （6） 代入式（5）得透镜曲率半径的计算公式 （7） 对给定的装置，R为常数，暗纹半径 （8） 和级数k的平方根成正比，即随着k的增大，条纹越来越细。 同理，如果r k是第k级明纹，则由式（1）和（2）得 （9） 代入式（5），可以算出（10）

牛顿迭代法实验报告

用牛顿迭代法求非线性方程的根 一、 实验题目 求方程()013=--=x x x f 在5.1附近的根。 二、 实验引言 （1）实验目的 1. 用牛顿迭代法求解方程的根 2. 了解迭代法的原理 3. 改进和修缮迭代法 （2）实验意义 牛顿迭代法就是众多解非线性方程迭代法中比较普遍的一种，求解方便实用。 三、 算法设计 （1）基本原理 给定初始值0x ，ε为根的容许误差，η为()x f 的容许误差，N 为迭代次数的容许值。 1．如果()0='x f 或迭带次数大于N ，则算法失败，结束；否则执行2. 2．计算()() 0001x f x f x x '-=. 3．若ε<-21x x 或()η<1x f ，则输出1x ，程序结束；否则执行4. 4．令10x x =，转向1. （2）流程图

四、程序设计program nndd01 implicit none real,parameter::e=0.005 real,parameter::n=9 real::x1 real::x0=1.5 integer::k real,external::f,y do k=1,9 if (y(x0)==0) then write(*,*)"失败" else x1=x0-f(x0)/y(x0) if (abs(x1-x0)

else x0=x1 end if end if end do end function f(x) implicit none real::f real::x f=x*x*x-x-1 return end function function y(x) implicit none real::y real::x y=3*x*x-1 return end function 五、求解结果 3 1.324718 4 1.324718 5 1.324718 6 1.324718 7 1.324718 8 1.324718 9 1.324718 六、算法评价及讨论 1.在求解在1.5处附近的根，不难发现在输入区间左端值为1时 需要迭代6次，而输入区间左端值为1.5时，却只要4次。初

牛顿插值法的应用

牛顿插值法在处理磁化曲线和铁损曲线 中的应用 指导老师：李国霞 院系：物理工程学院 专业：物理电子学 姓名：夏委委 学号：201112131526

一、牛顿插值法简介 在科学研究与其他领域中所遇到的许多实际问题中，经常会出现函数不便于处理或计算的情形。有时候函数关系没有明显的解析表达式，需要根据实验数据或其他方法来确定与自变量的某些值相对应的函数值；有时候函数虽有明显的解析表达式，但是使用很不方便。因此，在实际应用中，往往需要对实际使用的函数建立一个简单的便于处理和计算的近似表达式，即用一个简单的函数表达式来近似替代原来复杂的函数。与用近似数代替准确值一样，这也是计算法中最基本的概念和方法之一。近似代替又称为逼近。用多项式逼近列表函数的问题即为多项式插值问题。根据函数)(x f 已有的数据表格来计算函数)(x f 在一些新的点x 处的函数值，这就是插值法所要解决的问题。因此，所谓的插值法就是在所给定的函数表格中间在插入一些所需要的新的点上的函数值。 插值法的基本思想：首先设法根据表格中已有的函数值来构造一个简单的函数)(x y 作为)(x f 的近似表达式，然后再用)(x y 来计算新的点上的函数值作为 )(x f 的近似值。通常可以选用多项式函数作为近似函数)(x y ，因为多项式具有 各阶的导数，求值比较方便。用代数多项式作为工具研究插值问题，通常称为代数插值。 代数插值法问题的完整提法如下：设函数)(x f y =在区间[]b a ,上是连续的，且已知)(x f 在区间[]b a ,上1+n 个互异点处的函数值，即n i x f y i i ,......1,0),(== 其中，)(j i x x j i ≠≠。寻找一个次数不高于 n 的多项式 0111)(a x a x a x a x P n n n n n +++=-- 使满足条件n i x f x P i i n ,,1,0),()( ==称)(x P n 为)(x f 的插值多项式，),,1,0(n i x i =称为插值结点，[]b a ,称为插值区间。 牛顿(Newton)插值是数值逼近中的一个重要部分,它向前继承了拉格朗日(Lagrange)插值,向后引出了埃尔米特(Hermite)插值,可以看作对多项式插值作了一个简单的统一。牛顿插值公式具有形式简单,便于计算等优点。因此,在插值中得到广泛的应用。牛顿插值公式为)()()(x R x P x f n n +=，其中)(x P n 是牛顿插值多项式，)(x R n 为牛顿插值余项，)(x P n 和)(x R n 的表达式如下式所示：

matlab(迭代法-牛顿插值)Word版

实验报告内容： 一：不动点迭代法解方程 二：牛顿插值法的MATLAB实现 完成日期：2012年6月21日星期四 数学实验报告一 日期：2012-6-21

所以，确定初值为x0=1 二：不断迭代 算法： 第一步：将f(x0)赋值给x1 第二步：确定x1-x0的绝对值大小，若小于给定的误差值，则将x1当做方程的解，否则回到第一步 编写计算机程序： clear f=inline('0.5*sin(x)+0.4'); x0=1; x1=f(x0); k=1; while abs(x1-x0)>=1.0e-6 x0=x1; x1=f(x0); k=k+1; fprintf('k=%.0f,x0=%.9f,x1=%.9f\n',k,x0,x1) end 显示结果如下： k=2,x0=0.820735492,x1=0.765823700 k=3,x0=0.765823700,x1=0.746565483 k=4,x0=0.746565483,x1=0.739560873

k=6,x0=0.736981783,x1=0.736027993 k=7,x0=0.736027993,x1=0.735674699 k=8,x0=0.735674699,x1=0.735543758 k=9,x0=0.735543758,x1=0.735495216 k=10,x0=0.735495216,x1=0.735477220 k=11,x0=0.735477220,x1=0.735470548 k=12,x0=0.735470548,x1=0.735468074 k=13,x0=0.735468074,x1=0.735467157 >>。。。 以下是程序运行截图：