哈工大数值分析报告上机实验报告

数值分析2024上机实验报告数值分析是计算数学的一个重要分支，它研究如何用数值方法来解决数学问题。

在数值分析的学习过程中，学生需要通过上机实验来巩固理论知识，并学会使用相应的数值方法来解决实际问题。

本篇报告将详细介绍2024年度数值分析上机实验的内容和结果。

一、实验内容2024年度数值分析上机实验分为四个部分，分别是：方程求根、插值与拟合、数值积分和常微分方程的数值解。

1.方程求根这部分实验要求使用数值方法求解给定的非线性方程的根。

常见的数值方法有二分法、牛顿法、割线法等。

在实验过程中，我们需要熟悉这些数值方法的原理和实现步骤，并对不同方法的收敛性进行分析和比较。

2.插值与拟合这部分实验要求使用插值和拟合方法对给定的一组数据进行拟合。

插值方法包括拉格朗日插值、牛顿插值等；拟合方法包括最小二乘拟合、多项式拟合等。

在实验中，我们需要熟悉插值和拟合方法的原理和实现步骤，并对不同方法的精度和稳定性进行比较。

3.数值积分这部分实验要求使用数值方法计算给定函数的积分。

常见的数值积分方法有梯形法则、辛普森法则、龙贝格积分等。

在实验过程中，我们需要熟悉这些数值积分方法的原理和实现步骤，并对不同方法的精度和效率进行比较。

4.常微分方程的数值解这部分实验要求使用数值方法求解给定的常微分方程初值问题。

常见的数值方法有欧拉法、改进的欧拉法、四阶龙格-库塔法等。

在实验中，我们需要熟悉这些数值解方法的原理和实现步骤，并对不同方法的精度和稳定性进行比较。

二、实验结果在完成2024年度数值分析上机实验后，我们得到了以下实验结果：1.方程求根我们实现了二分法、牛顿法和割线法，并对比了它们的收敛速度和稳定性。

结果表明，割线法的收敛速度最快，但在一些情况下可能会出现振荡；二分法和牛顿法的收敛速度相对较慢，但稳定性较好。

2.插值与拟合我们实现了拉格朗日插值和最小二乘拟合，并对比了它们的拟合效果和精度。

结果表明，拉格朗日插值在小区间上拟合效果较好，但在大区间上可能出现振荡；最小二乘拟合在整体上拟合效果较好，但可能出现过拟合。

指导教师：姓名：学号：专业：联系电话：上海交通大学目录序言 (3)实验课题(一) 雅可比迭代法和高斯-塞得尔迭代法的收敛性和收敛速度 (4)数值分析 (6)实验课题(二) 松弛因子对SOR法收敛速度的影响 (6)数值分析 (12)总结 (13)附录（程序清单） (14)1.雅可比迭代法和高斯-塞得尔迭代法的收敛性和收敛速度 (14)雅可比迭代法： (14)高斯-塞得尔迭代法： (16)2.松弛因子对SOR法收敛速度的影响 (18)松弛法（SOR） (18)序言随着科学技术的发展，提出了大量复杂的数值计算问题，在实际解决这些计算问题的长期过程中，形成了计算方法这门学科，专门研究各种数学问题的数值解法（近似解法），包括方法的构造和求解过程的误差分析，是一门内容丰富，有自身理论体系的实用性很强的学科。

解决工程问题，往往需要处理很多数学模型，这就要花费大量的人力和时间，但是还有不少数学模型无法用解析法得到解。

使用数值方法并利用计算机，就可以克服这些困难。

事实上，科学计算已经与理论分析、科学实验成为平行的研究和解决科技问题的科学手段，经常被科技工作者所采用。

作为科学计算的核心内容——数值分析（数值计算方法），已逐渐成为广大科技工作者必备的基本知识并越来越被人重视。

由于数值方法是解数值问题的系列计算公式，所以数值方法是否有效，不但与方法本身的好坏有关，而且与数值问题本身的好坏也有关，因此，研究数值方法时，不但需要研究数值方法的好坏，即数值稳定性问题，而且还需要研究数值问题本身的好坏，即数值问题的性态，以及它们的判别问题。

数值计算的绝大部分方法都具有近似性，而其理论又具有严密的科学性，方法的近似值正是建立在理论的严密性基础上，根据计算方法的这一特点。

因此不仅要求掌握和使用算法，还要重视必要的误差分析，以保证计算结果的可靠性。

数值计算还具有应用性强的特点，计算方法的绝大部分方法如求微分方程近似解，求积分近似值，求解超越方程，解线性方程组等都具有较强的实用性，而插值法，最小二乘法，样条函数等也都是工程技术领域中常用的，有实际应用价值的方法。

数值分析第一次上机练习实验报告——Lagrange 插值与三次样条插值一、 问题的描述设()2119f x x =+, []1,1x ∈-,取15iix =-+,0,1,2,...,10i =.试求出10次Lagrange 插值多项式()10L x 和三次样条插值函数()S x (采用自然边界条件),并用图画出()f x ,()10L x ,()S x .二、 方法描述——Lagrange 插值与三次样条插值我们取15i ix =-+,0,1,2,...,10i =，通过在i x 点的函数值()2119i i f x x =+来对原函数进行插值，我们记插值函数为()g x ，要求它满足如下条件：()()21,0,1,2,...,1019i i ig x f x i x ===+ (1) 我们在此处要分别通过Lagrange 插值(即多项式插值)与三次样条插值的方法对原函数()2119f x x=+进行插值，看两种方法的插值结果，并进行结果的比较。

10次的Lagrange 插值多项式为：()()10100i i i L x y l x ==∑ (2)其中：()21,0,1,2,...,1019i i i y f x i x ===+以及()()()()()()()()()011011......,0,1,2,...,10......i i n i i i i i i i n x x x x x x x x l x i x x x x x x x x -+-+----==----我们根据(2)进行程序的编写，我们可以通过几个循环很容易实现函数的Lagrange 插值。

理论上我们根据区间[]1,1-上给出的节点做出的插值多项式()n L x 近似于()f x ，而多项式()n L x 的次数n 越高逼近()f x 的精度就越好。

但实际上并非如此，而是对任意的插值节点，当n →+∞的时候()n L x 不一定收敛到()f x ；而是有时会在插值区间的两端点附近会出现严重的()n L x 偏离()f x 的现象，即所谓的Runge 现象。

数值分析第一次上机练习实验报告一、实验目的本次实验旨在通过上机练习，加深对数值分析方法的理解，并掌握实际应用中的数值计算方法。

二、实验内容1. 数值计算的基本概念和方法在本次实验中，我们首先回顾了数值计算的基本概念和方法。

数值计算是一种通过计算机进行数值近似的方法，其包括近似解的计算、误差分析和稳定性分析等内容。

2. 方程求解的数值方法接下来，我们学习了方程求解的数值方法。

方程求解是数值分析中非常重要的一部分，其目的是找到方程的实数或复数解。

我们学习了二分法、牛顿法和割线法等常用的数值求解方法，并对它们的原理和步骤进行了理论学习。

3. 插值和拟合插值和拟合是数值分析中常用的数值逼近方法。

在本次实验中，我们学习了插值和拟合的基本原理，并介绍了常见的插值方法，如拉格朗日插值和牛顿插值。

我们还学习了最小二乘拟合方法，如线性拟合和多项式拟合方法。

4. 数值积分和数值微分数值积分和数值微分是数值分析中的两个重要内容。

在本次实验中，我们学习了数值积分和数值微分的基本原理，并介绍了常用的数值积分方法，如梯形法和辛卜生公式。

我们还学习了数值微分的数值方法，如差商法和牛顿插值法。

5. 常微分方程的数值解法常微分方程是物理和工程问题中常见的数学模型，在本次实验中，我们学习了常微分方程的数值解法，包括欧拉法和四阶龙格-库塔法。

我们学习了这些方法的步骤和原理，并通过具体的实例进行了演示。

三、实验结果及分析通过本次实验，我们深入理解了数值分析的基本原理和方法。

我们通过实际操作，掌握了方程求解、插值和拟合、数值积分和数值微分以及常微分方程的数值解法等数值计算方法。

实验结果表明，在使用数值计算方法时，我们要注意误差的控制和结果的稳定性。

根据实验结果，我们可以对计算结果进行误差分析，并选择适当的数值方法和参数来提高计算的精度和稳定性。

此外，在实际应用中，我们还需要根据具体问题的特点和条件选择合适的数值方法和算法。

四、实验总结通过本次实验，我们对数值分析的基本原理和方法有了更加深入的了解。

实验报告一题目：非线性方程求解摘要：非线性方程的解析解通常很难给出，因此线性方程的数值解法就尤为重要。

本实验采用两种常见的求解方法二分法和Newton法及改进的Newton法。

前言：（目的和意义）掌握二分法与Newton法的基本原理和应用。

数学原理：对于一个非线性方程的数值解法很多。

在此介绍两种最常见的方法：二分法和Newton法。

对于二分法，其数学实质就是说对于给定的待求解的方程f(x)，其在[a,b]上连续，f(a)f(b)<0，且f(x)在[a,b]内仅有一个实根x*，取区间中点c，若，则c恰为其根，否则根据f(a)f(c)<0是否成立判断根在区间[a,c]和[c,b]中的哪一个，从而得出新区间，仍称为[a,b]。

重复运行计算，直至满足精度为止。

这就是二分法的计算思想。

Newton法通常预先要给出一个猜测初值x0，然后根据其迭代公式产生逼近解x*的迭代数列{x k}，这就是Newton法的思想。

当x0接近x*时收敛很快，但是当x0选择不好时，可能会发散，因此初值的选取很重要。

另外，若将该迭代公式改进为其中r为要求的方程的根的重数，这就是改进的Newton法，当求解已知重数的方程的根时，在同种条件下其收敛速度要比Newton法快的多。

程序设计：本实验采用Matlab的M文件编写。

其中待求解的方程写成function的方式，如下function y=f(x);y=-x*x-sin(x);写成如上形式即可，下面给出主程序。

二分法源程序：clear%%%给定求解区间b=1.5;a=0;%%%误差R=1;k=0;%迭代次数初值while (R>5e-6) ;c=(a+b)/2;if f12(a)*f12(c)>0;a=c;elseb=c;endR=b-a;%求出误差k=k+1;endx=c%给出解Newton法及改进的Newton法源程序：clear%%%% 输入函数f=input('请输入需要求解函数>>','s')%%%求解f(x)的导数df=diff(f);%%%改进常数或重根数miu=2;%%%初始值x0x0=input('input initial value x0>>');k=0;%迭代次数max=100;%最大迭代次数R=eval(subs(f,'x0','x'));%求解f(x0)，以确定初值x0时否就是解while (abs(R)>1e-8)x1=x0-miu*eval(subs(f,'x0','x'))/eval(subs(df,'x0','x'));R=x1-x0;x0=x1;k=k+1;if (eval(subs(f,'x0','x'))<1e-10);breakendif k>max;%如果迭代次数大于给定值，认为迭代不收敛，重新输入初值ss=input('maybe result is error,choose a new x0,y/n?>>','s');if strcmp(ss,'y')x0=input('input initial value x0>>');k=0;elsebreakendendendk;%给出迭代次数x=x0；%给出解结果分析和讨论：1.用二分法计算方程在[1，2]内的根。

数值分析上机实验报告摘要：本报告是对数值分析课程上机实验的总结和分析，涵盖了多种算法和数据处理方法，通过对实验结果的分析，探究了数值计算的一般过程和计算的稳定性。

1. 引言数值计算是数学的一个重要分支，广泛应用于物理、金融、工程等领域。

本次实验是对数值分析课程知识的实际应用，通过上机实现算法，探究数值计算的可靠性和误差分析。

2. 实验方法本次实验中，我们实现了多种算法，包括：（1）牛顿迭代法求方程的根；（2）高斯消元法求线性方程组的解；（3）最小二乘法拟合数据点；（4）拉格朗日插值法估计函数值；（5）梯形公式和辛普森公式求积分近似值。

对于每个算法，我们都进行了多组数值和不同参数的实验，并记录了相关数据和误差。

在实验过程中，我们着重考虑了算法的可靠性和计算的稳定性。

3. 实验结果与分析在实验中，我们得到了大量的实验数据和误差分析，通过对数据的展示和分析，我们得到了以下结论：（1）牛顿迭代法求解非线性方程的根能够对算法的初始值和迭代次数进行适当的调整，从而达到更高的稳定性和可靠性。

（2）高斯消元法求解线性方程组的解需要注意到矩阵的奇异性和精度的影响，从而保证计算的准确性。

（3）最小二乘法拟合数据点需要考虑到拟合的函数形式和数据的误差范围，采取适当的数据预处理和拟合函数的选择能够提高计算的准确性。

（4）拉格朗日插值法估计函数值需要考虑到插值点的选择和插值函数的阶数，防止出现龙格现象和插值误差过大的情况。

（5）梯形公式和辛普森公式求积分近似值需要考虑到采样密度和拟合函数的选择，从而保证计算的稳定性和收敛速度。

4. 结论通过本次实验的分析和总结，我们得到了深入的认识和理解数值计算的一般过程和算法的稳定性和可靠性，对于以后的数值计算应用也提供了一定的指导和参考。

数值分析上机实验报告《数值分析》上机实验报告1.用Newton 法求方程 X 7-X 4+14=0在（0.1，1.9）中的近似根（初始近似值取为区间端点，迭代6次或误差小于0.00001）。

1.1 理论依据：设函数在有限区间[a ，b]上二阶导数存在，且满足条件{}αϕ上的惟一解在区间平方收敛于方程所生的迭代序列迭代过程由则对任意初始近似值达到的一个中使是其中上不变号在区间],[0)(3,2,1,0,)(')()(],,[x |))(),((|,|,)(||)(|.4;0)(.3],[)(.20)()(.110......b a x f x k x f x f x x x Newton b a b f a f mir b a c x f ab c f x f b a x f b f x f k k k k k k ==-==∈≤-≠>+令)9.1()9.1(0)8(4233642)(0)16(71127)(0)9.1(,0)1.0(,1428)(3225333647>⋅''<-=-=''<-=-='<>+-=f f x x x x x f x x x x x f f f x x x f故以1.9为起点⎪⎩⎪⎨⎧='-=+9.1)()(01x x f x f x x k k k k 如此一次一次的迭代，逼近x 的真实根。

当前后两个的差<=ε时，就认为求出了近似的根。

本程序用Newton 法求代数方程(最高次数不大于10)在（a,b ）区间的根。

1.2 C 语言程序原代码：#include<stdio.h>#include<math.h> main(){double x2,f,f1;double x1=1.9; //取初值为 1.9 do{x2=x1;f=pow(x2,7)-28*pow(x2,4)+14; f1=7*pow(x2,6)-4*28*pow(x2,3); x1=x2-f/f1;}while(fabs(x1-x2)>=0.00001||x1<0.1); //限制循环次数 printf("计算结果：x=%f\n",x1);}1.3 运行结果：1.4 MATLAB 上机程序function y=Newton(f,df,x0,eps,M) d=0;for k=1:Mif feval(df,x0)==0d=2;breakelsex1=x0-feval(f,x0)/feval(df,x0);ende=abs(x1-x0);x0=x1;if e<=eps&&abs(feval(f,x1))<=epsd=1;breakendendif d==1y=x1;elseif d==0y='迭代M次失败';elsey= '奇异'endfunction y=df(x)y=7*x^6-28*4*x^3;Endfunction y=f(x)y=x^7-28*x^4+14;End>> x0=1.9;>> eps=0.00001;>> M=100;>> x=Newton('f','df',x0,eps,M);>> vpa(x,7)1.5 问题讨论：1.使用此方法求方解，用误差来控制循环迭代次数，可以在误差允许的范围内得到比较理想的计算结果。

实验成绩实验报告
课程名称：数学实验
实验项目：常微分方程数值解
所在院系：
学生姓名：
学生学号：
授课学期： 2020秋
完成时间：
当1)0(=y ，取区间为[0,10]，用改进欧拉法得到的结果与ode23得到的结果如下图所示：
（2）0cos =+''x y y 1)0(=y ，0)0(='y 取区间[0,10]，用改进欧拉法得到的结果与ode23得到的结果如下图所示：
四、实验结果或结论分析
由实验结果可知，改进欧拉法与用ode23输出的结果十分相近，欧拉法在求解微分方程中具有极其重要的作用。

若改进式欧拉法取点越多，迭代次数越多，所得结果越接近真实值。

五、Matlab源程序
问题1的程序：
n=101;
y=zeros(n,1);
x=zeros(n,1);
y(1)=1;
x(1)=0;
long=10;
b=0;
m=0;
for i=1:100
m=y(i+1);
x(i+1)=(long/n)*i;
K1=8;
while k1>0.1
m=b;
b=y(i)+ long/n/2*(fun(x(i+1),m)+fun(x(i),y(i)));
%b=2*y(i+1)-y(i)+(long/n)^2/2*(fun(x(i+1),m)+fun(x(i),y(i))); K1=abs(m-b);
end
y(i+1)=b;
end
[m,n]=ode23(@(x,y)(x^2-y^2),[0,10],0.10);
plot(x,y);。

数学实验报告实验一Matlab的使用1.上机实验各种数据输入方法：程序语句：a=[1 2 3;4 5 6 ;7,8,9] 程序语句：linspace(1,10,5) 等等…………计算结果：a = 计算结果：ans =1 2 34 5 6 1.0000 3.2500 5.5000 7.7500 10.00007 8 92.(1) (a)方法：(b) 方法：程序语句：程序语句：a=[-3 5 0 8;1 -8 2 -1;0 -5 9 3;-7 0 -4 5]; a=[-3 5 0 8;1 -8 2 -1;0 -5 9 3;-7 0 -4 5];b=[0;2;-1;6]; b=[0;2;-1;6];inv(a)*b a\b计算结果：计算结果：ans = ans =-0.6386 -0.6386-0.4210 -0.4210-0.3529 -0.35290.0237 0.0237(2) 4个矩阵的生成语句：矩阵a 的生成语句：e=eye(3,3); a=[e r;o s]r=rand(3,2); 验证语句：o=zeros(2,3); a^2s=diag([1,2]);%此为一个任取的2X2 矩阵b=[e r+r*s; o s^2]计算结果相同：ans =1.0000 0 0 1.9003 1.45790 1.0000 0 0.4623 2.67390 0 1.0000 1.2137 2.28630 0 0 1.0000 00 0 0 0 4.00003.生成多项式的语句：poly ([2,-3,1+2i,1-2i,0,-6])计算结果：ans = 1 5 -9 -1 72 -180 0 计算x＝0.8,-x=-1.2 之值的指令与结果：指令：polyval([1,5,-9,-1,72,-180,0],0.8) 结果：ans= -100.2179指令：polyval([1,5,-9,-1,72,-180,0],-1.2) 结果：ans= 293.29004.求a的指令与结果：指令：a=compan([1,0,-6,3,-8])结果：a =0 6 -3 81 0 0 00 1 0 00 0 1 0求a的特征值的指令与结果：roots(p)的指令与结果为：指令：eig(a) 指令：roots([1,0,-6,3,-8])结果：结果：ans = ans =-2.8374 -2.83742.4692 2.46920.1841 + 1.0526i 0.1841 + 1.0526i0.1841 - 1.0526i 0.1841 - 1.0526i结论：利用友元阵函数a=company(p) 和eig(a) 可以与roots(p)有相同的作用，结果相同。

数值分析上机实验报告导言：本次上机实验主要是针对数值分析课程中的一些基本算法进行实验验证。

实验内容包括迭代法、插值法、数值积分和常微分方程的数值解等。

在实验过程中，我们将会使用MATLAB进行算法的实现，并对结果进行分析。

一、迭代法迭代法是解决函数零点、方程解等问题的常用方法。

我们将选择几个常见的函数进行迭代求根的实验。

（1）二分法二分法是一种简单而有效的迭代求根法。

通过函数在区间两个端点处的函数值异号来确定函数在区间内存在零点，并通过不断缩小区间来逼近零点。

（2）牛顿法牛顿法利用函数的一阶导数和二阶导数的信息来逼近零点。

通过不断迭代更新逼近值，可以较快地求得零点。

实验结果表明，对于简单的函数，这两种迭代法都具有很好的收敛性和稳定性。

但对于一些复杂的函数，可能会出现迭代失效或者收敛速度很慢的情况。

二、插值法插值法是在给定一些离散数据点的情况下，通过构造一个插值函数来逼近未知函数的值。

本实验我们将使用拉格朗日插值和牛顿插值两种方法进行实验。

（1）拉格朗日插值拉格朗日插值通过构造一个多项式函数来逼近未知函数的值。

该多项式经过离散数据点，并且是唯一的。

该方法简单易懂，但插值点越多，多项式次数越高，插值函数的精度也就越高。

（2）牛顿插值牛顿插值利用差商的概念，通过构造一个插值多项式来逼近未知函数的值。

与拉格朗日插值相比，牛顿插值的计算过程更加高效。

但同样要求插值点的选择要合理，否则可能出现插值函数不收敛的情况。

实验结果表明，这两种插值方法都能够很好地逼近未知函数的值。

插值点的选择对插值结果有很大的影响，过多或者过少的插值点都可能导致插值结果偏离真实函数的值。

三、数值积分数值积分是一种将定积分问题转化为数值求和的方法。

本实验我们将使用复合梯形求积法和复合辛普森求积法进行实验。

（1）复合梯形求积法复合梯形求积法将定积分区间等分为若干小区间，然后使用梯形公式对每个小区间进行近似求积，最后将结果相加得到整个定积分的近似值。

数值分析上机实习报告随着现代科学技术的迅猛发展，计算机科学的应用日益广泛，数值分析作为计算机科学中重要的分支之一，其在工程、物理、生物学等领域的应用也越来越受到重视。

本学期，我们在数值分析课程的学习中，进行了多次上机实习，通过实习，我们对数值分析的基本方法和算法有了更深入的理解和掌握。

在实习过程中，我们使用了MATLAB软件作为主要的工具，MATLAB是一种功能强大的数学软件，它提供了丰富的数值计算函数和图形显示功能，使我们能够更加方便地进行数值计算和分析。

第一次实习是线性插值和函数逼近。

我们学习了利用已知数据点构造插值函数的方法，并通过MATLAB软件实现了线性插值和拉格朗日插值。

通过实习，我们了解了插值的基本原理，掌握了插值的计算方法，并能够利用MATLAB软件进行插值计算。

第二次实习是解线性方程组。

我们学习了高斯消元法、列主元高斯消元法和克莱姆法则等解线性方程组的方法，并通过MATLAB软件实现了这些算法。

在实习过程中，我们通过实际例子了解了这些算法的应用，掌握了它们的计算步骤，并能够利用MATLAB软件准确地求解线性方程组。

第三次实习是求解非线性方程和方程组。

我们学习了二分法、牛顿法、弦截法和迭代法等求解非线性方程的方法，以及雅可比法和高斯-赛德尔法等求解非线性方程组的方法。

通过实习，我们了解了非线性方程和方程组的求解方法，掌握了它们的计算步骤，并能够利用MATLAB软件求解实际问题。

通过这次上机实习，我们不仅深入学习了数值分析的基本方法和算法，而且锻炼了利用MATLAB软件进行数值计算和分析的能力。

同时，我们也认识到了数值分析在实际问题中的应用价值，增强了解决实际问题的能力。

总之，这次上机实习使我们受益匪浅，对我们学习数值分析课程起到了很好的辅助作用。

数值分析上机实验报告实验报告插值法与数值积分实验（数值计算方法，3学时）一实验目的1．掌握不等距节点下的牛顿插值公式以及拉格朗日插值公式。

2．掌握复化的梯形公式、辛扑生公式、牛顿-柯特斯公式计算积分。

3. 会用龙贝格公式和高斯公式计算积分。

二实验内容用拉格朗日插值公式计算01.54.1==y x 以及所对应的近似值。

用牛顿插值公式求)102(y 的近似值。

三实验步骤（算法）与结果1拉格朗日插值法：（C 语言版）#include "Stdio.h" #include "Conio.h"int main(void) {float X[20],Y[20],x; int n;void input(float *,float *,float *,int *); float F(float *,float *,float,int); input(X,Y,&x,&n);printf("F(%f)=%f",x,F(X,Y,x,n));getch(); return 0; }void input(float *X,float *Y,float *x,int *n) {int i;printf("Please input the number of the data:");scanf("%d",n);printf("\nPlease input the locate of each num:\n");for(i=0;i<*n;i++){scanf("%f,%f",X+i,Y+i);}printf("\nPlease input the chazhi:"); scanf("%f",x);}float F(float *X,float *Y,float x,int n){int i,j;float Lx,Fx=0;for(i=0;i<n;i++)< p="">{Lx=1;for(j=0;j<n;j++)< p="">{if(j!=i) Lx=Lx*((x-*(X+j))/(*(X+i)-*(X+j))); } Fx=Fx+Lx*(*(Y+i));}return Fx;}得出结果如图：所以Y(1.4)=3.7295252#include#define N 10double X[N], Y[N], A[N][N];int n;double Newton(double x);double f(double x);void main() {printf("请输入已知x与对应y=f(x)的个数: n = "); scanf("%d", &n);getchar();if(n>N||n<=0) {printf("由于该维数过于犀利, 导致程序退出!"); return;}printf("\n请输入X[%d]: ", n);for (int i=0; i<="" p="">scanf("%lf", &X[i]);getchar();printf("\n请输入Y[%d]: ", n);for (i=0; i<="" p="">scanf("%lf", &Y[i]);getchar();double x;printf("\n请输入所求结点坐标x = ");scanf("%lf", &x);getchar();printf("\nf(%.4lf)≈%lf\n\n", x, Newton(x));}double Newton(double x) {int i, j;// 求均差for (i =0; i<="" p="">A[i][0] = Y[i];for (i=1; i<="" p="">for (j =1; j<=i; j++)A[i][j] = (A[i][j-1] - A[i-1][j-1]) / (X[i] - X[i-j]); // 求结点double result = A[0][0];for (i=1; i<="">double tmp = 1.0;for (int j=0; j<="" p="">tmp *= (x - X[j]);result += tmp * A[i][i];}return result;}四实验收获与教师评语</n;j++)<></n;i++)<>。

数值分析上机实验报告《数值分析》上机实验报告1.用Newton 法求方程 X 7-X 4+14=0在（0.1，1.9）中的近似根（初始近似值取为区间端点，迭代6次或误差小于0.00001）。

1.1 理论依据：设函数在有限区间[a ，b]上二阶导数存在，且满足条件{}αϕ上的惟一解在区间平方收敛于方程所生的迭代序列迭代过程由则对任意初始近似值达到的一个中使是其中上不变号在区间],[0)(3,2,1,0,)(')()(],,[x |))(),((|,|,)(||)(|.4;0)(.3],[)(.20)()(.110......b a x f x k x f x f x x x Newton b a b f a f mir b a c x f ab c f x f b a x f b f x f k k k k k k ==-==∈≤-≠>+令)9.1()9.1(0)8(4233642)(0)16(71127)(0)9.1(,0)1.0(,1428)(3225333647>⋅''<-=-=''<-=-='<>+-=f f x x x x x f x x x x x f f f x x x f故以1.9为起点⎪⎩⎪⎨⎧='-=+9.1)()(01x x f x f x x k k k k 如此一次一次的迭代，逼近x 的真实根。

当前后两个的差<=ε时，就认为求出了近似的根。

本程序用Newton 法求代数方程(最高次数不大于10)在（a,b ）区间的根。

1.2 C语言程序原代码：#include<stdio.h>#include<math.h>main(){double x2,f,f1;double x1=1.9; //取初值为1.9do{x2=x1;f=pow(x2,7)-28*pow(x2,4)+14;f1=7*pow(x2,6)-4*28*pow(x2,3);x1=x2-f/f1;}while(fabs(x1-x2)>=0.00001||x1<0.1); //限制循环次数printf("计算结果：x=%f\n",x1);}1.3 运行结果：1.4 MATLAB上机程序function y=Newton(f,df,x0,eps,M)d=0;for k=1:Mif feval(df,x0)==0d=2;breakelsex1=x0-feval(f,x0)/feval(df,x0);ende=abs(x1-x0);x0=x1;if e<=eps&&abs(feval(f,x1))<=epsd=1;breakendendif d==1y=x1;elseif d==0y='迭代M次失败';elsey= '奇异'endfunction y=df(x)y=7*x^6-28*4*x^3;Endfunction y=f(x)y=x^7-28*x^4+14;End>> x0=1.9;>> eps=0.00001;>> M=100;>> x=Newton('f','df',x0,eps,M);>> vpa(x,7)1.5 问题讨论：1.使用此方法求方解，用误差来控制循环迭代次数，可以在误差允许的范围内得到比较理想的计算结果。

材料力学I上机实验设计报告院系：机电学院班级： 1308***姓名： ***学号： 11308*****指导教师：张桂莲时间： 2015年6月一、问题描述1、应力状态分析对于空间或者是平面应力状态的相关计算，如果采用人工计算的方式比较繁琐而且容易出错，对于这种简单的重复计算，编制相应的程序则可以大大提高计算准确度和人工计算强度。

对于平面应力状态，输入量应为（,,x y xy σστ），以及某截面的方位角α，其输出数据应为该单元体所受主应力（123,,σσσ），所受最大剪应力（13max 132σσττ-==），以及方位角为α的斜截面上的应力（,ααστ）以及主方向角σα，同时还要画出其应力圆示意图，以直观的显示其应力状态。

对于空间应力状态，输入量则应该为各应力（,,,,,x y z xy yz xz σσστττ），其输出数据应该为该单元体所受主应力（123,,σσσ），所受最大剪应力（13max 132σσττ-==），同时还要画出其应力圆示意图，以直观的显示其应力状态。

这样，应力状态分析的基本任务就可以完成。

2、常用截面图形几何性质的分析在生活中，有各种各样的几何形状，但是对于工程实际中经常用到的构件，其截面的几何形状则非常有限。

对于不同的截面，其形心位置、对于形心轴的惯性矩也就有所不同，这样在进行如弯曲、扭转等的应力分析时就会到来不便，因此编制相应的程序来计算相关截面的几何性质也就具有了实际应用价值和可行性。

在这部分程序中，截面几何形状分为三角形、矩形、椭圆形、梯形、圆形、扇形等多种形式，对于不同的截面形状，输入量也就不同。

例如，对于扇形应输入直径和圆心角（,d α）；对于梯形则应输入上底、下底和高（,,a b h ）；对于椭圆形，则要输入长轴长和短轴长（,a b ）等等，在此不一一列举，具体输入数据请参看程序运行。

不过对于不同的截面，其输出的量都是相同的，即截面形心的位置、面积、对于形心轴的惯性矩（,,,,C C C C y z y z S I I ），这些输出量就是这些截面的基本几何性质参数，有了这些参数之后则可以对其进行进一步的计算和接下来的分析等问题。

哈工大数值分析实验报告标题：哈工大数值分析实验报告一、实验目的：本实验的目的是探究在数值分析中使用的各种数值方法，对于解决实际问题的有效性和可靠性进行评估。

二、实验内容：本实验主要包括以下几个方面的内容：1. 熟悉数值分析中常用的数值方法，如数值积分、数值微分、迭代法等；2. 在MATLAB等数学软件平台上，编写程序实现所学的数值方法；3. 使用所编写的程序，对给定的实际问题进行求解，并分析其结果的有效性和可靠性；4. 根据实际问题的特点，评估不同数值方法的适用性，并给出相应的结论和建议。

三、实验步骤：1. 阅读相关的理论知识，熟悉数值分析中常用的数值方法；2. 编写数值分析实验的程序代码，包括数值积分、数值微分和迭代法等；3. 使用编写的程序，对所给的实际问题进行求解，记录并分析结果；4. 根据实际问题的特点，评估所使用的数值方法的可靠性和有效性；5. 根据实验结果，撰写实验报告，包括实验目的、实验内容、实验步骤和实验结果的分析等。

四、实验结果：根据实际问题的不同，实验结果也会有所差异。

在实验报告中，可以详细叙述对所给实际问题的求解过程，并对结果进行分析和解释。

同时，还可以比较不同数值方法的结果，评估其优劣和适用性。

五、实验结论：根据实验结果的分析，可以得出结论，总结不同数值方法的优缺点，并对其在实际问题中的应用进行评价。

同时，还可以给出相应的建议，为以后的数值分析工作提供参考。

六、实验总结：通过本次实验，进一步加深了对数值分析中常用数值方法的理解和掌握。

通过实际问题的求解，对于这些数值方法的应用和效果有了更深入的认识。

同时，也提高了编程和科研报告撰写的能力，为以后的学习和工作打下了坚实的基础。

以上是关于哈工大数值分析实验报告的基本内容，具体实验细节和结果请根据实际情况进行补充。

实验报告哈尔滨工程大学教务处制实验三 数值积分一．数值积分的基本思想1.复合梯形公式：Tn=++)()([2b f a f h2∑-=11)](n k xk f ；2.复合辛普森公式：Sn=6h[f(a)+f(b)+2∑-=11)](n k xk f +4∑-=+1)2/1(n k x f ]；以上两种算法都是将a-b 之间分成多个小区间（n ）,则h=(b-a)/n,x k =a+kh,x k+1/2=a+(k+1/2)h,利用梯形求积根据两公式便可。

3.龙贝格算法：在指定区间内将步长依次二分的过程中运用如下公式（1）Sn=34T2n-31Tn（2）Cn=1516S2n-151Sn（3）Rn=6364C2n-631Cn4T)(k m=144-m m T )1(1+-k m - 141-mT )(1k m -，k = 1，2，… 二.实验题目及实验目的(第4章计算实习题第1题)用不同数值方法计算积分xdx x ln 1⎰= -94。

（1）取不同的步长h 。

分别用复合梯形及复合辛普森求积计算积分，给出误差中关于h 的函数，并与积分精确值比较两个公式的精度，是否存在一个最小的h ，使得精度不能再被改善？（2）用龙贝格求积计算完成问题（1）。

（3）用自适应辛普森积分，使其精度达到104-。

三．实验手段：指操作环境和平台:win7系统下MATLAB R2009a程序语言：一种类似C 语言的程序语言，但比C 语言要宽松得多，非常方便。

四.程序①复合梯形求积程序function t=TiXing_quad(a,b,.h) format longx=a:h:b;y=sqrt(x).*log(x);y(1)=0;t=0;for k=1:(b-a)/h,t= t+y(k)+y(k+1);endt=t*h/2;②复合辛普森求积程序function s=Simpson_quad(a,b,h) format longx=a:h:b;y=sqrt(x).*log(x);z=sqrt(x+h/2).*log(x+h/2);y(1)=0;s=0;for k=1:(b-a)/h,s= s+y(k)+y(k+1)+4*z(k);ends=s*h./6;③龙贝格求积程序function [q,R]=Romberg(a,b,eps) h=b-a;R(1,1)=h*(0+sqrt(b).*log(b))/2; M=1;J=0;err=1;while err>epsJ=J+1;h=h/2;S=0;for p=1:Mx=a+h*(2*p-1);S=S+sqrt(x).*log(x);endR(J+1,1)=R(J,1)/2+h*S;M=2*M;for k=1:JR(J+1,k+1)=R(J+1,k)+(R(J+1,k)-R(J,k))/(4^k-1);enderr=abs(R(J+1,J)-R(J+1,J+1));endq=R(J+1,J+1);控制台输入代码：（1）>> a=0;>> b=1;>> h=0.1;>> t=TiXing_quad(a,b,h)>> s=Simpson_quad(a,b,h)>> h=0.01;>> t=TiXing_quad(a,b,h)>> s=Simpson_quad(a,b,h)>> h=0.001;>> t=TiXing_quad(a,b,h)>> s=Simpson_quad(a,b,h)（2）>> a=0;>> b=1;>> eps=10^-8;>> [quad,R]=Romberg(a,b,eps)（3）>> a=0;>> b=1;>> eps=10^-4;>> q=ZiShiYingSimpson('sqrt(x).*log(x)',a,b,eps) 五．实验结果比较与分析（1）h = 0.1时h = 0.01时h = 0.001时由结果（1）可知对于同一步长h，复合辛普森法求积分精度明显比复合梯形法求积的精度要高，且当步长取不同值时即h越小时，积分精度越高。