计算机上机实验报告
验
1掌握用MATLA 计算拉格朗日、分段线性、三次样条三种插值
目 2、 掌握用MATAB 作线性最小二乘的方法。 的方法，改变节点的数目，对三种插值结果进行初步分析。
3、
通过实例学习如何
用插值方法与拟合方法解决实际问题， 注意
二者的联系和区别。 实验的主要内容
1编制拉格朗日、牛顿插值程序，并运行一个简单的实例。
(1)拉格朗日插值程序：
fun cti on v 二pol yin terp(x, y, u) n=len gth(x); v=zeros (size(u));
for k=l: n w=on es(size(u));
for j=[l:k~l k+1:n] w=(u-x(j)). /(x (k) -x(j)). *w;
end v 二v+w*y(k); end
实例:当x=144, 169, 225时，y=12, 13, 15,用拉格朗日差值法求根
号如下:
的 和 要
求 实 验 内 容 和 步 骤
*Ut 丫2 二l “匚二Q2UL3P t-JiLJW |K4i
。莎tl马c＜安殳厂□ B ? □
(2)牛顿插值程序：
fun ctio n y=n ewi nterp (X, Y, x)% 牛顿插值函数m=le ngth(X);
for k=2:m
for j=l:k~l
Y(k)= (Y(k)- Y(j))/(X(k)-X(j));
end
end
y=Y(m);
for j=m-l：-1:1
y 二y?*(x-X(j))+Y(j);
end
实例：当x=144,169,225时，y=42,13,15,用牛顿差值法求根号175。如下：
2、给定函数f(x) \己知:
f(2. 0) . 1.414214 f (2. 1)? 1.449138
f (2.2) . 1.483240 f (2. 3) . 1.516575 f (2. 4)
用牛顿插值法求4次Newton插值多项式在2. 15处的值, Comma nd Window
? 口2, 0 2. 1 2. 2 2. 3 2
Y二[L414214 L449138 叨IS; ?
1.483240 LE1S5F5 1.549193]
? y二ne"interp (Xj Yj
1.4S?3
? z=[l. 9;0. 0tl?:2. 6]? y=nevrirrt exp 傍了」莖)i ? pltrt % y) ? 1. 549193以此作为函
3?选择函数y=exp (-x 2) (-2 作初步分析。 程序： %不同插值方法是否会出现震荡runge现象 %文件 fun cti on run gelO ［X,Y］ =fenduan(10, 1) ;%将［-1,］区间分成10等份，返回对应的(x, y)五组数据 x=linspace (-2, 2,100);%将［-1,1］划分成100等份，以便作岀样条插值多项式的图形。 for i=l: le ngth(x)% 绘制原函数曲线图 y (i)=exp (-x(i)A2); end hold on plot (x, y); text (0, 1,' \ left arrow原函数')％对曲线添加标注y=n ewi nterp(X, Y, x) ;%多项式插值中的牛顿插值法hold on Plot (x, y); titlef插值函数中的runge现象，区间等分为10段');濟加标题 xlabel (' X 轴')； ylabel (J Y 轴')； text (-0. 9, 1. 5,' \ left arrow牛顿插值')％对曲线添加标注 y=i nterpl (X, Y, x); plot (x, y); text (-0. 4, 0. 8521, J \ left arrow 分段线性插值') cs=spline(X, [0 Y 0]);%调用spline函数插值，y比x多两个元素。 plot (X, Y, ' o', x, ppval (cs, x),' 做样条多项式的图形 text (-1. 2, 0. 2369,J \ left arrow 样条插值') function [X, Y]=fenduan(n, b)%将区间等分成n份,并求对应点上的函数值for i=l: n+1 X(i)=-2+(4*(i-l))/n; Y(i)=exp(-X(i)A2); end fun ctio n y=n ewi nterp (X, Y, x)% 牛顿插值函数 m=le ngth (X); for k=2:m for j=l:k~l Y(k) = (Y(k)- Y(j))/(X(k)-X(j)); end end y=Y(m); for j=m-l:~1:1 y 二y?*(x-X(j))+Y(j); end 将三种插值结果相比较，显然分段线性插值法在节点处不光滑，拉格朗日值出现较大的振荡，样条差值的结果是最好的，改变n 的值，运行程序，得到的图形如右图所示，比较这两个图可发现，节点增加后，三种插值方法结果的准确度均有所提高，因此可近似地认为：增加节点的个数可以提高插值结果的准确程度。 学生实验报告 了解插值与拟合的基本原理和方法；掌握用MATLAB计算插值与作最小二乘多项式拟合和曲线拟合的方法；通过范例展现求解实际问题的初步建模过程； 通过动手作实验学习如何用插值与拟合方法解决实际问题，提高探索和解决问题的能力。这对于学生深入理解数学概念，掌握数学的思维方法，熟悉处理大量的工程计算问题的方法具有十分重要的意义。 二、实验仪器、设备或软件：电脑,MATLAB软件 三、实验内容 1．编写插值方法的函数M文件； 2．用MATLAB中的函数作函数的拟合图形； 3．针对实际问题，试建立数学模型，并求解。 四、实验步骤 1．开启软件平台——MATLAB，开启MATLAB编辑窗口； 2．根据各种数值解法步骤编写M文件； 3．保存文件并运行； 4．观察运行结果(数值或图形)； 5．写出实验报告，并浅谈学习心得体会。 五、实验要求与任务 根据实验内容和步骤，完成以下具体实验，要求写出实验报告（实验目的→问题→数学模型→算法与编程→计算结果→分析、检验和结论→心得体会）。 1．天文学家在1914年8月的7次观测中，测得地球与金星之间距离（单位：米），并取得常用对数值，与日期的一组历史数据如下表： 由此推断何时金星与地球的距离（米）的对数值为9.93518？ 解：输入命令 days=[18 20 22 24 26 28 30]; distancelogs=[9.96177 9.95436 9.94681 9.93910 9.93122 9.92319 9.91499]; t1=interp1(distancelogs,days,9.93518) %线性插值 t2=interp1(distancelogs,days,9.93518,'nearest') %最近邻点插值 t3=interp1(distancelogs,days,9.93518,'spline') %三次样条插值 t4=interp1(distancelogs,days,9.93518,'cubic') %三次插值 计算结果： t1 = 24.9949 t2 = 24 t3 = 25.0000 t4 = 大学物理仿真实验报告－牛顿环法测曲率半径 ————————————————————————————————作者: ————————————————————————————————日期: 大学物理仿真实验报告 实验名称 牛顿环法测曲率半径 班级: 姓名: 学号: 日期: 牛顿环法测曲率半径 实验目的 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０) 插值法和拟合实验报告 一、 实验目的 1.通过进行不同类型的插值，比较各种插值的效果，明确各种插值的优越性； 2.通过比较不同次数的多项式拟合效果，了解多项式拟合的原理； 3.利用matlab 编程，学会matlab 命令； 4.掌握拉格朗日插值法； 5.掌握多项式拟合的特点和方法。 二、 实验题目 1.、插值法实验 将区间[-5,5]10等分，对下列函数分别计算插值节点 k x 的值，进行不同类型 的插值，作出插值函数的图形并与)(x f y =的图形进行比较： ;11)(2x x f += ;a r c t a n )(x x f = .1)(42 x x x f += (1) 做拉格朗日插值； (2) 做分段线性插值； (3) 做三次样条插值. 2、拟合实验 给定数据点如下表所示： 分别对上述数据作三次多项式和五次多项式拟合，并求平方误差，作出离散函数 ),(i i y x 和拟合函数的图形。 三、 实验原理 1.、插值法实验 ∏∑∏∏∏∑∑≠==≠=≠=≠=+-==--= =-= ==-=-=----==++==j i j j i i i i i n i i n n j i j j n j i j j i i n j i j j n i i i n i i n n n o i n i i n x x x x x y x l x L x x c n i x x c x x x c x x x x x x x x c y x l x L y x l y x l y x l x L ,00 ,0,0,01100 00 )(l )()() (1 ,1,0, 1)()(l ) ()())(()()()()()()()(， 故， 得 再由，设 2、拟合实验 等厚干涉——牛顿环 等厚干涉是薄膜干涉的一种。薄膜层的上下表面有一很小的倾角是，从光源发出的光经上下表面反射后在上表面附近相遇时产生干涉，并且厚度相同的地方形成同一干涉条纹，这种干涉就叫等厚干涉。其中牛顿环是等厚干涉的一个最典型的例子，最早为牛顿所发现，但由于他主张微粒子学说而并未能对他做出正确的解释。光的等厚干涉原理在生产实践中育有广泛的应用，它可用于检测透镜的曲率，测量光波波长，精确地测量微笑长度、厚度和角度，检验物体表面的光洁度、平整度等。 一.实验目的 (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 -=- 《数值分析》课程实验一：插值与拟合 一、实验目的 1. 理解插值的基本原理，掌握多项式插值的概念、存在唯一性； 2. 编写MA TLAB 程序实现Lagrange 插值和Newton 插值，验证Runge 现象； 3. 通过比较不同次数的多项式拟合效果，理解多项式拟合的基本原理； 4. 编写MA TLAB 程序实现最小二乘多项式曲线拟合。 二、实验内容 1. 用Lagrange 插值和Newton 插值找经过点(-3, -1), (0, 2), (3, -2), (6, 10)的三次插值公式，并编写MATLAB 程序绘制出三次插值公式的图形。 2. 设 ]5,5[,11 )(2 -∈+= x x x f 如果用等距节点x i = -5 + 10i /n (i = 0, 1, 2, …, n )上的Lagrange 插值多项式L n (x )去逼近它。不妨取n = 5和n = 10，编写MATLAB 程序绘制出L 5(x )和L 10(x )的图像。 (2) 编写MA TLAB 程序绘制出曲线拟合图。 三、实验步骤 1. (1) Lagrange 插值法：在线性空间P n 中找到满足条件： ?? ?≠===j i j i x l ij j i , 0,, 1)(δ 的一组基函数{}n i i x l 0)(=，l i (x )的表达式为 ∏ ≠==--= n i j j j i j i n i x x x x x l ,0),,1,0()( 有了基函数{}n i i x l 0)(=，n 次插值多项式就可表示为 ∑==n i i i n x l y x L 0 )()( (2) Newton 插值法：设x 0, x 1, …, x n 是一组互异的节点，y i = f (x i ) (i = 0, 1, 2, …, n )，f (x )在处的n 阶差商定义为 等厚干涉牛顿环实验报告 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.实验目的 1．学习和掌握拉格朗日插值多项式。 2.运用拉格朗日插值多项式进行计算。 2.实验过程 作出插值点（1.00，0.00），（-1.00，-3.00），（2.00，4.00）二、算法步骤 已知：某些点的坐标以及点数。 输入：条件点数以及这些点的坐标。 输出：根据给定的点求出其对应的拉格朗日插值多项式的值。 3.程序流程： （1）输入已知点的个数； （2）分别输入已知点的X坐标； （3）分别输入已知点的Y坐标； 程序如下： #include 插值算法*/ { int i,j; float *a,yy=0.0; /*a作为临时变量，记录拉格朗日插值多项*/ a=(float*)malloc(n*sizeof(float)); for(i=0;i<=n-1;i++) { a[i]=y[i]; for(j=0;j<=n-1;j++) if(j!=i) a[i]*=(xx-x[j])/(x[i]-x[j]); yy+=a[i]; } free(a); return yy; } int main() { int i; int n; float x[20],y[20],xx,yy; printf("Input n:"); scanf("%d",&n); if(n<=0) { printf("Error! The value of n must in (0,20)."); getch();return 1; } for(i=0;i<=n-1;i++) { printf("x[%d]:",i); scanf("%f",&x[i]); } printf("\n"); for(i=0;i<=n-1;i++) { printf("y[%d]:",i);scanf("%f",&y[i]); } printf("\n"); printf("Input xx:"); scanf("%f",&xx); yy=lagrange(x,y,xx,n); printf("x=%f,y=%f\n",xx,yy); getch(); } 举例如下：已知当x=1，-1,2时f(x)=0，-3,4，求f(1.5)的值。 北京师范大学珠海分校大学物理实验报告 实验名称：牛顿环实验测量 学院工程技术学院 专业测控技术与仪器 学号 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 一、给定函数y=sinx的函数表如下表，用拉格朗日插值求sin0.57891的近似 值 M文件： function yh=lagrange2(x0,y0,xh) n = length(x0); m = length(xh); yh=zeros(1,m); for k = 1:m for i = 1:n xp = x0([1:i-1 i+1:n]); yp = prod((xh(k)-xp)./(x0(i)-xp)); yh(k) = yh(k) + yp*y0(i); end end 执行：>> x0=[0.4,0.5,0.6,0.7] x0 = 0.4000 0.5000 0.6000 0.7000 >> y0=[0.38942,0.47943,0.56464,0.64422] y0 = 0.3894 0.4794 0.5646 0.6442 >> lagrange2(x0,y0,0.57891) 执行结果： ans = 0.5471 二、 1. 给定sin110.190809,sin120.207912,sin130.224951,o o o ===构造牛顿 插值函数计算'sin1130o 。 M 文件： function fp = newpoly(x,y,p) n = length(x); a(1) = y(1); for k = 1 : n - 1 d(k, 1) = (y(k+1) - y(k))/(x(k+1) - x(k)); end for j = 2 : n - 1 for k = 1 : n - j d(k, j) = (d(k+1, j - 1) - d(k, j - 1))/(x(k+j) - x(k)); end end d for j = 2 : n a(j) = d(1, j-1); end Df(1) = 1; c(1) = a(1); for j = 2 : n Df(j)=(p - x(j-1)) .* Df(j-1); c(j) = a(j) .* Df(j); 实验三 插值法与拟合实验 一、实验目的 1. 通过本实验学会利用程序画出插值函数，并和原图形相比较 2. 通过本实验学会拟合函数图形的画法，并会求平方误差 二、实验题目 1. 插值效果的比较 实验题目：区间[]5,5-10等分，对下列函数分别计算插值节点k x 的值，进行不同类型的插值，作出插值函数的图形并与)(x f y =的图形进行比较： 2 11)(x x f +=； x x f arctan )(=； 4 41)(x x x f += （1） 做拉格朗日插值； （2） 做三次样条插值. 2. 拟合多项式实验 实验题目：给定数据点如下表所示： 分别对上述数据作三次多项式和五次多项式拟合，并求平方误差，作出离散函数),(i i y x 和拟合函数的图形. 三、实验原理 本实验应用了拉格朗日插值程序、三次样条插值程序、多项式拟合程序等实验原理. 四、实验内容 1（1） figure x=-5:0.2:5; y=1./(1+x.^2); plot(x,y,'r'); hold on %拉格朗日插值 x1=-5:1:5; y1=1./(1+x1.^2); xx=-4.5:0.5:4.5; yy=malagr(x1,y1,xx); plot(xx,yy,'+') %三次样条插值 dy0=1./(1+25); dyn=1./(1+25); m=maspline(x1,y1,dy0,dyn,xx); plot(xx,m,'ok') 1（2） x=-5:0.2:5; y=atan(x); plot(x,y,'r'); hold on %拉格朗日插值 x1=-5:1:5; y1=atan(x1); xx=-4.5:0.5:4.5; yy=malagr(x1,y1,xx); plot(xx,yy,'+') %三次样条插值 dy0=1./(1+25); dyn=1./(1+25); m=maspline(x1,y1,dy0,dyn,xx); plot(xx,m,'ok') 1（3） x=-5:0.2:5; y=x.^2./(1+x.^4); plot(x,y,'r'); hold on %拉格朗日插值 x1=-5:1:5; y1=x1.^2./(1+x1.^4); xx=-4.5:0.5:4.5; yy=malagr(x1,y1,xx); plot(xx,yy,'+') %三次样条插值 dy0=1./(1+25); dyn=1./(1+25); m=maspline(x1,y1,dy0,dyn,xx); plot(xx,m,'ok') 2. x=[-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5]'; y=[-4.45 -0.45 0.55 0.05 -0.44 0.54 4.55]'; plot(x,y,'or'); hold on %三次多项式拟合 p1=mafit(x,y,3); 实验二插值法 1、实验目的： 1、掌握直接利用拉格郎日插值多项式计算函数在已知点的函数值；观察拉格郎日插值的龙格现象。 2、了解Hermite插值法、三次样条插值法原理，结合计算公式，确定函数值。 2、实验要求: 1)认真分析题目的条件和要求，复习相关的理论知识，选择适当的解决方案和算法； 2)编写上机实验程序，作好上机前的准备工作； 3)上机调试程序，并试算各种方案，记录计算的结果（包括必要的中间结果）； 4)分析和解释计算结果； 5)按照要求书写实验报告； 3、实验内容： 1) 用拉格郎日插值公式确定函数值；对函数f(x)进行拉格郎日插值，并对f(x)与插值多项式的曲线作比较。 已知函数表：（0.56160,0.82741）、（0.56280,0.82659）、（0.56401,0.82577）、（0.56521,0.82495）用三次拉格朗日插值多项式求x=0.5635时函数近似值。 2) 求满足插值条件的插值多项式及余项 1) 4、题目：插值法 5、原理： 拉格郎日插值原理： n次拉格朗日插值多项式为：L n (x)=y l (x)+y 1 l 1 (x)+y 2 l 2 (x)+…+y n l n (x) n=1时，称为线性插值， L 1(x)=y (x-x 1 )/(x -x 1 )+y 1 (x-x )/(x 1 -x )=y +(y 1 -x )(x-x )/(x 1 -x ) n=2时，称为二次插值或抛物线插值， L 2(x)=y (x-x 1 )(x-x 2 )/(x -x 1 )/(x -x 2 )+y 1 (x-x )(x-x 2 )/(x 1 -x )/(x 1 -x 2 )+y 2 (x -x 0)(x-x 1 )/(x 2 -x )/(x 2 -x 1 ) n=i时， Li= （X-X0）……(X-X i-1)(x-x i+1) ……(x-x n) （X-X0）……(X-X i-1)(x-x i+1) ……(x-x n) 6、设计思想： 拉格朗日插值法是根据n + 1个点x0, x1, ... x n(x0 < x1 < ... x n)的函数值f (x0), f (x1) , ... , f (x n)推出n次多項式p(x)，然后n次多項式p (x)求出任意的点x对应的函数值f (x)的算法。 7、对应程序： 1 ) 三次拉格朗日插值多项式求x=0.5635时函数近似值 #include"stdio.h" #define n 5 void main() { int i,j; float x[n],y[n]; float x1; float a=1; float b=1; float lx=0; printf("\n请输入想要求解的X：\n x="); scanf("%f",&x1); printf("请输入所有点的横纵坐标：\n"); for(i=1;i 计算方法实验报告 专业班级：医学信息工程一班姓名：陈小芳学号：201612203501002 实验成绩： 1．【实验题目】 插值法与曲线拟合 2．【实验目的】 3．【实验内容】 4. 【实验要求】 5. 【源程序（带注释）】 （1）拉格朗日插值 #include 广州大学学生实验报告 开课学院及实验室： 2014年 月 日 学院 数学与信息科学学院 年级、专业、班 姓名 学号 实验课程名称 数学实验 成绩 实验项目名称 实验2 插值与拟合 指导老师 一、实验目的 1、掌握用MATLAB 计算拉格朗日、分段线性、三次样条三种插值的方法，改变节点的数目，对三种插值结果进行初步分析。 2、掌握用MATLAB 作线性最小二乘拟合的方法。 3、通过实例学习如何用插值方法与拟合方法解决实际问题，注意二者的联系和区别。 二、实验设备 电脑、MATLAB 三、实验要求 1..选择一些函数，在n 个节点上（n ）不要太大，如5～11）用拉格朗日，分段线性，三次样条三种插值方法，，计算m 各插值点的函数值（m 要适中，如50～100）.通过数值和图形的输出，将三种插值结果与精确值进行比较.适当增加n ，再作比较，由此作初步分析.下列函数供选择参考： a. y=sin x ,0≦x ≦2π； 2.用 1 2 y x =在x=0,1,4,9,16产生5个节点15,...,P P .用不同的节点构造插值公式来计算x=5处的插值（如用 15,...,P P ；14,...,P P ；24,...,P P 等）与精确值比较进行分析。 5.对于实验1中的录像机计数器，自己实测一组数据（或利用给出的数据），确定模型2 t an bn =+中的系数a,b. 6.用电压V=10伏的电池给电容器充电，电容器上t 时刻的电压为 0()()t v t V V V e -τ =--,其中 0V 是电容器的初始 电压，τ是充电常数。试由下面一组t ，V 数据确定0V 和τ. t/s 0.5 1 2 3 4 5 7 9 V/V 6.36 6.48 7.26 8.22 8.66 8.99 9.43 9.63 8. 弹簧在力F 的作用下伸长x ，一定范围内服从胡克定律：F 与x 成正比，即F=kx,k 为弹性系数.现在得到下面一组x ，F 数据，并在（x,F ）坐标下作图(图13).可以看出，当Ｆ大到一定数值(如x=9以后)后，就不服从这个定律了。试由数据拟合直线F=kx,并给出不服从胡克定律时的近似公式（曲线）。 1）要求直线与曲线在x=9处相连接。 2）要求直线与曲线在x=9处光滑连接. 四、实验程序 预备： function y=lagr1(x0,y0,x) n=length(x0);m=length(x); for i=1:m z=x(i); s=0.0; for k=1:n p=1.0; for j=1:n if j~=k p=p*(z-x0(j))/(x0(k)-x0(j)); end end s=s+p*y0(k); end y(i)=s; end 五、实验操作过程 当n=5时 clear; n=5; %在n 个节点上进行插值 m=75; %产生m 个插值点，计算函数在插值点处的精确值，将来进行对比 x=0:4/(m-1):2*pi; y=sin(x); z=0*x; x0=0:4/(n-1):2*pi; y0=sin(x0); y1=lagr1(x0,y0,x); % y1为拉格朗日插值 y2=interp1(x0,y0,x); % y2为分段线性插值 y3=spline(x0,y0,x); % y3为三次样条插值 [x' y' y1' y2' y3'] plot(x,z,'k',x,y,'r:',x,y1,'g-.',x,y2,'b',x,y3,'y--') gtext('Lagr.'), gtext('Pieces. linear'), gtext('Spline'), gtext('y=sin(x)') hold off; %比较插值所得结果与函数在插值点处的精确值 s = ' x y y1 y2 y3' [x' y' y1' y2' y3'] 结果 ans = 0 0 0 0 0 0.0541 0.0540 0.0495 0.0455 0.0611 0.1081 0.1079 0.0999 0.0910 0.1207 0.1622 0.1615 0.1510 0.1365 0.1787 0.2162 0.2145 0.2025 0.1819 0.2350 0.2703 0.2670 0.2541 0.2274 0.2896 0.3243 0.3187 0.3054 0.2729 0.3425 0.3784 0.3694 0.3563 0.3184 0.3936 0.4324 0.4191 0.4066 0.3639 0.4429 0.4865 0.4675 0.4559 0.4094 0.4904 0.5405 0.5146 0.5040 0.4548 0.5359 0.5946 0.5602 0.5508 0.5003 0.5796 0.6486 0.6041 0.5961 0.5458 0.6212 0.7027 0.6463 0.6396 0.5913 0.6609 0.7568 0.6866 0.6812 0.6368 0.6985 0.8108 0.7248 0.7208 0.6823 0.7341 0.8649 0.7610 0.7583 0.7278 0.7675 大学物理仿真实验报告 实验名称：牛顿环法测曲率半径实验日期： 专业班级： 姓名：学号： 教师签字：________________ 一、实验目的 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） 实验题目： 用多项式模型进行数据拟合实验 1 实验目的 本实验使用多项式模型对数据进行拟合，目的在于： （1）掌握数据拟合的基本原理，学会使用数学的方法来判定数据拟合的情况； （2）掌握最小二乘法的基本原理及计算方法； （3）熟悉使用matlab 进行算法的实现。 2 实验步骤 2.1 算法原理 所谓拟合是指寻找一条平滑的曲线，最不失真地去表现测量数据。反过来说，对测量 的实验数据，要对其进行公式化处理，用计算方法构造函数来近似表达数据的函数关系。由于函数构造方法的不同，有许多的逼近方法，工程中常用最小平方逼近（最小二乘法理论）来实现曲线的拟合。 最小二乘拟合利用已知的数据得出一条直线或曲线，使之在坐标系上与已知数据之间的距离的平方和最小。模型主要有：1.直线型2.多项式型3.分数函数型4.指数函数型5.对数线性型6.高斯函数型等，根据应用情况，选用不同的拟合模型。其中多项式型拟合模型应用比较广泛。 给定一组测量数据()i i y x ,，其中m i ,,3,2,1,0Λ=，共m+1个数据点，取多项式P （x ），使得 min )]([020 2=-=∑∑==m i i i m i i y x p r ，则称函数P （x ）为拟合函数或最小二乘解，此时，令 ∑==n k k k n x a x p 0 )(，使得min ])([02 002=??? ? ??-=-=∑∑∑===m i n k i k i k m i i i n y x a y x p I ，其中 n a a a a ,,,,210Λ为待求的未知数，n 为多项式的最高次幂，由此该问题化为求),,,(210n a a a a I I Λ=的极值问题。 由多元函数求极值的必要条件：0)(200 =-=??∑∑==m i j i n k i k i k i x y x a a I ，其中n j ,,2,1,0Λ= 得到： ∑∑∑===+=n k m i i j i k m i k j i y x a x )(，其中n j ,,2,1,0Λ=，这是一个关于n a a a a ,,,,210Λ的线 性方程组，用矩阵表示如下所示： 插值与多项式逼近的数组计算方法实验 郑发进 2012042020022 【摘要】计算机软件中经常要用到库函数，如） cos，x e，它们 （x （x sin，） 是用多项式逼近来计算的。虽然目前最先进的逼近方法是有理函数（即多项式的商），但多项式逼近理论更适于作为数值分析的入门课程。在已知数据具有高精度的情况下，通常用组合多项式来构造过给定数据点的多项式。构造组合多项式的方法有许多种，如线性方程求解、拉格朗日系数多项式以及构造牛顿多项式的方分和系数表。 关键字泰勒级数、拉格朗日插值法、牛顿插值法、帕德逼近 一、实验目的 1.通过具体实验，掌握泰勒级数、拉格朗日插值法、牛顿插值法、帕德逼近的编程技巧。 2.比较各插值方法的优劣并掌握。 二、实验原理 1.泰勒级数 在数学中，泰勒级数（英语：Taylor series）用无限项连加式——级数来表示一个函数，这些相加的项由函数在某一点的导数求得。 如果在点x=x 具有任意阶导数，则幂级数 称为在点x 处的泰勒级数。 =0，得到的级数 在泰勒公式中，取x 称为麦克劳林级数。函数的麦克劳林级数是x的幂级数，那么这种展开 是唯一的，且必然与的麦克劳林级数一致。 2.拉格朗日插值法 如对实践中的某个物理量进行观测，在若干个不同的地方得到相应的观测值，拉格朗日插值法可以找到一个多项式，其恰好在各个观测的点取到观测到的值。这样的多项式称为拉格朗日（插值）多项式。数学上来说，拉格朗日插值法可以给出一个恰好穿过二维平面上若干个已知点的多项式函数。 在平面上有(x 1,y 1)(x 2,y 2)...(x n ,y n )共n 个点，现作一条函数f （x ）使其图像经过这n 个点。 作n 个多项式p i (x),i=1,2,3...,n,使得 最后可得 3.牛顿插值法 插值法利用函数f (x)在某区间中若干点的函数值，作出适当的特定函数，在这些点上取已知值，在区间的其他点上用这特定函数的值作为函数f (x)的近似值。如果这特定函数是多项式，就称它为插值多项式。利用插值基函数很容易得到拉格朗日插值多项式，公式结构紧凑，在理论分析中甚为方便，但当插值节点增减时全部插值基函数均要随之变化，整个公式也将发生变化， 这在实际计算中是很不方便的，为了克服这一缺点，提出了牛顿插值。 牛顿插值通过求各阶差商，递推得到的一个公式： 10121()()()()()()N N N N P x P x a x x x x x x x x --=+---- 牛顿插值与拉格朗日插值具有唯一性。 4.帕德逼近 它不仅与逼近论中其他许多方法有着密切的关系，而且在实际问题特别是许多物理问题中有着广泛的应用。设是在原点某邻域内收敛的、具有复系数的麦克劳林级数。欲确定一个有理函数，式中,使得前次方的系数为0，即使得 此处约定qk ＝0（k>n ）。虽然所求得的Pm(z)和Qn(z)不惟一，但是比式却总是惟一的。有理函数称为F(z)的(m,n)级帕德逼近,记为(m/n)。由(m/n)所形成的阵列称为帕德表。 实验报告 学生姓名：学号：指导教师： 实验地点：实验时间： 一、实验室名称： 二、实验项目名称：牛顿环测曲面半径和劈尖干涉 三、实验学时： 四、实验原理： 1、等厚干涉 如图1所示，在C点产生干涉，光线11`和22`的光程差为△=2d+λ/2 式中λ/2是因为光由光疏媒质入射到光密媒质上反射时，有一相位 突变引起的附加光程差。 当光程差△=2d+λ/2=(2k+1)λ 即d=k λ/2时产生暗条纹; 当光程差△=2d+λ/2=2kλ/2, 即d=(k－1/2)λ/2时产生明条纹 图1 因此,在空气薄膜厚度相同处产生同一级的干涉条纹,叫等厚干涉条纹。 2、用牛顿环测透镜的曲率半径 将一个曲率半径较大的平凸透镜的凸面置于一块光学平板玻璃上则 可组成牛顿环装置。如图2所示。 这两束反射光在AOB 表面上的某一点E 相遇，从而产生E 点的干涉。由于AOB 表面是球面，所产生的条纹是明暗相间 的圆环，所以称为牛顿环，如图3所示。 将两块光学平玻璃重叠在一起，在一端插入一薄纸片，则在两玻璃板间形成一空气劈尖，如图4所示。K 级干涉暗条纹对应的薄膜厚度为d=k λ/2 k=0时，d=0, 即在两玻璃板接触处为零级暗条纹；若在薄纸处呈现k=N 级条纹，则薄纸片厚度为 d ’=N λ/2 若劈尖总长为Ｌ，再测出相邻两条纹之间的距离为△x,则暗条纹总数为Ｎ＝Ｌ／△x ， 即 d ’=L λ/2 △x 。 五、实验目的： 深入理解光的等厚干涉及其应用，学会使用移测显微镜。 六、实验内容： １、用牛顿环测透镜的曲率半径 ２、用劈尖干涉法测薄纸片的厚度 七、实验器材（设备、元器件）： 牛顿环装置，移测显微镜，两块光学平玻璃板，薄纸片，钠光灯及电 图２ L d插值与拟合实验报告
大学物理仿真实验报告牛顿环法测曲率半径
插值法和拟合实验报告(数值计算)
等厚干涉牛顿环实验报告材料97459
数值分析实验插值与拟合
等厚干涉牛顿环实验报告
数值分析拉格朗日插值法上机实验报告
牛顿环实验报告
插值与拟合实验报告
计算方法--插值法与拟合实验
插值法实验报告
实验四 插值法与曲线拟合
数学实验-实验2 插值与拟合
大学物理仿真实验报告 牛顿环
用多项式模型进行数据拟合实验报告(附代码)
插值与多项式逼近的数组计算方法实验讲解
牛顿环等厚干涉标准实验报告
相关主题
文本预览