数值计算方法实验报告(含所有)

数值计算方法实验报告

数值计算方法实验报告一、实验介绍本次实验是关于数值计算方法的实验，旨在通过计算机模拟的方法，实现对于数值计算方法的掌握。

本次实验主要涉及到的内容包括数值微积分、线性方程组的求解、插值与拟合、常微分方程的数值解等。

二、实验内容1. 数值微积分数值微积分是通过计算机模拟的方法，实现对于微积分中的积分运算的近似求解。

本次实验中，我们将会使用梯形公式和辛普森公式对于一定区间上的函数进行积分求解，并比较不同公式的计算误差。

2. 线性方程组的求解线性方程组求解是数值计算领域中的重要内容。

本次实验中，我们将会使用高斯消元法、LU分解法等方法对于给定的线性方程组进行求解，并通过比较不同方法的计算效率和精度，进一步了解不同方法的优缺点。

3. 插值与拟合插值与拟合是数值计算中的另一个重要内容。

本次实验中，我们将会使用拉格朗日插值法和牛顿插值法对于给定的数据进行插值求解，并使用最小二乘法对于给定的函数进行拟合求解。

4. 常微分方程的数值解常微分方程的数值解是数值计算中的难点之一。

本次实验中，我们将会使用欧拉法和龙格-库塔法等方法对于给定的常微分方程进行数值解的求解，并比较不同方法的计算精度和效率。

三、实验结果通过本次实验，我们进一步加深了对于数值计算方法的理解和掌握。

在数值微积分方面，我们发现梯形公式和辛普森公式都能够有效地求解积分，但是辛普森公式的计算精度更高。

在线性方程组求解方面，我们发现LU分解法相对于高斯消元法具有更高的计算效率和更好的数值精度。

在插值与拟合方面，我们发现拉格朗日插值法和牛顿插值法都能够有效地进行插值求解，而最小二乘法则可以更好地进行函数拟合求解。

在常微分方程的数值解方面，我们发现欧拉法和龙格-库塔法都能够有效地进行数值解的求解，但是龙格-库塔法的数值精度更高。

四、实验总结本次实验通过对于数值计算方法的模拟实现，进一步加深了我们对于数值计算方法的理解和掌握。

在实验过程中，我们了解了数值微积分、线性方程组的求解、插值与拟合、常微分方程的数值解等多个方面的内容，在实践中进一步明确了不同方法的特点和优缺点，并可以通过比较不同方法的计算效率和数值精度来选择合适的数值计算方法。

数值计算方法实验报告

数值分析实验报告实验一、解线性方程组的直接方法——梯形电阻电路问题利用追赶法求解三对角方程组的方法，解决梯形电阻电路问题：电路中的各个电流{1i ，2i ，…，8i }须满足下列线性方程组：R V i i =- 22 210 252321=-+-i i i 0 252 432=-+-i i i 0 252 543=-+-i i i 0 252 654=-+-i i i 0 252 765=-+-i i i 0 252 876=-+-i i i 052 87=+-i i设V 220=V ，Ω=27R ，运用追赶法，求各段电路的电流量。

问题分析：上述方程组可用矩阵表示为：⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡--------------00000001481.8522520000002520000002520000002520000002520000002520000002287654321i i i i i i i i问题转化为求解A x b =，8阶方阵A 满足顺序主子式(1,2...7)0i A i =≠，因此矩阵A存在唯一的Doolittle 分解，可以采用解三对角矩阵的追赶法！追赶法a=[0 -2 -2 -2 -2 -2 -2 -2]; b=[2 5 5 5 5 5 5 5]；c=[-2 -2 -2 -2 -2 -2 -2 0]; d=[220/27 0 0 0 0 0 0 0];Matlab 程序function x= zhuiganfa( a,b,c,d )%追赶法实现要求：|b1|>|C1|>0,|bi|>=|ai|+|ci| n=length(b); u=ones(1,n); L=ones(1,n); y=ones(1,n); u(1)=b(1); y(1)=d(1); for i=2:nL(i)=a(i)/u(i-1);u(i)=b(i)-c(i-1)*L(i); y(i)=d(i)-y(i-1)*L(i); endx(n)=y(n)/u(n); for k=n-1:-1:1x(k)=(y(k)-c(k)*x(k+1))/u(k); end endMATLAB 命令窗口输入：a=[0 -2 -2 -2 -2 -2 -2 -2]; b=[2 5 5 5 5 5 5 5];c=[-2 -2 -2 -2 -2 -2 -2 0] d=[220/27 0 0 0 0 0 0 0];x= zhuiganfa(a,b,c,d )运行结果为：x =8.1478 4.0737 2.0365 1.0175 0.5073 0.2506 0.1194 0.0477存在问题根据电路分析中的所讲到的回路电流法，可以列出8个以回路电流为独立变量的方程，课本上给出的第八个回路电流方程存在问题，正确的应该是78240i i -+=；或者可以根据电路并联分流的知识，同样可以确定78240i i -+=。

数值计算基础实验报告(3篇)

第1篇一、实验目的1. 理解数值计算的基本概念和常用算法；2. 掌握Python编程语言进行数值计算的基本操作；3. 熟悉科学计算库NumPy和SciPy的使用；4. 分析算法的数值稳定性和误差分析。

二、实验内容1. 实验环境操作系统：Windows 10编程语言：Python 3.8科学计算库：NumPy 1.19.2，SciPy 1.5.02. 实验步骤（1）Python编程基础1）变量与数据类型2）运算符与表达式3）控制流4）函数与模块（2）NumPy库1）数组的创建与操作2）数组运算3）矩阵运算（3）SciPy库1）求解线性方程组2）插值与拟合3）数值积分（4）误差分析1）舍入误差2）截断误差3）数值稳定性三、实验结果与分析1. 实验一：Python编程基础（1）变量与数据类型通过实验，掌握了Python中变量与数据类型的定义方法，包括整数、浮点数、字符串、列表、元组、字典和集合等。

（2）运算符与表达式实验验证了Python中的算术运算、关系运算、逻辑运算等运算符，并学习了如何使用表达式进行计算。

（3）控制流实验学习了if-else、for、while等控制流语句，掌握了条件判断、循环控制等编程技巧。

（4）函数与模块实验介绍了Python中函数的定义、调用、参数传递和返回值，并学习了如何使用模块进行代码复用。

2. 实验二：NumPy库（1）数组的创建与操作通过实验，掌握了NumPy数组的基本操作，包括创建数组、索引、切片、排序等。

（2）数组运算实验验证了NumPy数组在数学运算方面的优势，包括加、减、乘、除、幂运算等。

（3）矩阵运算实验学习了NumPy中矩阵的创建、操作和运算，包括矩阵乘法、求逆、行列式等。

3. 实验三：SciPy库（1）求解线性方程组实验使用了SciPy库中的线性代数模块，通过高斯消元法、LU分解等方法求解线性方程组。

（2）插值与拟合实验使用了SciPy库中的插值和拟合模块，实现了对数据的插值和拟合，并分析了拟合效果。

数值计算方法实验报告

数值计算方法实验报告实验目的：通过实验验证不同数值计算方法在求解数学问题时的精度和效率，并分析其优缺点。

实验原理：实验内容：本实验选取了三个典型的数值计算问题，并分别采用了二分法、牛顿迭代法和梯度下降法进行求解。

具体问题和求解方法如下：1. 问题一：求解方程sin(x)=0的解。

-二分法：利用函数值的符号变化将解空间不断缩小，直到找到满足精度要求的解。

-牛顿迭代法：通过使用函数的斜率来逼近方程的解，并不断逼近真实解。

-梯度下降法：将方程转化为一个极小化问题，并利用梯度下降的方式逼近极小值点，进而找到方程的解。

2.问题二：求解函数f(x)=x^2-3x+2的极小值点。

-二分法：通过确定函数在一个区间内的变化趋势，将极小值所在的区间不断缩小，从而找到极小值点。

-牛顿迭代法：通过使用函数的导数和二阶导数来逼近极小值点，并不断逼近真实解。

-梯度下降法：将函数转化为一个极小化问题，并利用梯度下降的方式逼近极小值点，进而找到函数的极小值点。

3. 问题三：求解微分方程dy/dx = -0.1*y的解。

-二分法：通过离散化微分方程，将微分方程转化为一个差分方程，然后通过迭代计算不同点的函数值，从而得到函数的近似解。

-牛顿迭代法：将微分方程转化为一个积分方程，并通过迭代计算得到不同点的函数值，从而得到函数的近似解。

-梯度下降法：将微分方程转化为一个极小化问题，并利用梯度下降的方式逼近极小值点，从而得到函数的近似解。

实验步骤：1.编写代码实现各个数值计算方法的求解过程。

2.对每个数值计算问题，设置合适的初始值和终止条件。

3.运行程序，记录求解过程中的迭代次数和每次迭代的结果。

4.比较不同数值计算方法的精度和效率，并分析其优缺点。

实验结果：经过实验测试，得到了如下结果：-问题一的二分法迭代次数为10次，求解结果为x=0；牛顿迭代法迭代次数为4次，求解结果为x=0；梯度下降法迭代次数为6次，求解结果为x=0。

-问题二的二分法迭代次数为10次，求解结果为x=1;牛顿迭代法迭代次数为3次，求解结果为x=1;梯度下降法迭代次数为4次，求解结果为x=1-问题三的二分法迭代次数为100次，求解结果为y=e^(-0.1x);牛顿迭代法迭代次数为5次，求解结果为y=e^(-0.1x);梯度下降法迭代次数为10次，求解结果为y=e^(-0.1x)。

数值实验报告

北京XX大学计算机与通信工程学院实验报告实验名称：《数值计算方法》课程实验学生姓名：_____________XX_______________专业：________计算机科学与技术________班级：_____________XXXX____________学号：_____________XXXXXX___________指导教师：_____________XXXXX_____________实验成绩：________________________________实验地点：__________XXXXXXX____________实验时间：____2017____年___6___月___6___日一、实验目的与实验要求1、实验目的实验1：探究非线性方程的解法，比较不同解法的优劣性，针对具体问题对解法进行实践，并迭代达到指定的精确度。

实验2：探究一般化曲线拟合的方法，采用最小二乘法对含有多项式，指数和对数多种形式的函数进行拟合，确定拟合结果中的系数。

实验3：探究多种积分的数值解法，根据给定的精度，选择恰当的数值积分方法，确定迭代步数与步长，分析积分数值解法的优劣性。

2、实验要求实验1：采用两种算法求解非线性方程，比较两种算法的性能。

实验2：用最小二乘法拟合由不超过三阶多项式和指数、对数函数线性组合成的符合函数，确定各个项的系数实验3：自选一种数值积分方法求解积分值，根据要求的精度给出最终的迭代步数和步长，分析优缺点。

二、实验设备（环境）及要求实验1，3采用C语言编写，编译环境为DEV C++。

实验2涉及矩阵操作，故采用MATLAB编写。

三、实验内容与步骤1、实验1（1）实验内容·采用至少两种不同的算法求解ex+3*x3-x2-2=0在[0,1]范围内的一个根，要求两次迭代误差小于10-4。

·根据实验结果，比较分析不同算法的性能。

（2）主要步骤本实验由C 语言实现。

在两种不同的算法上选用简单迭代法和牛顿迭代法。

数值计算方法实验报告

2如果f[(a+b)/2]<0，则区间（(a+b)/2，b)内存在零点，(a+b)/2≥a；
3如果f[(a+b)/2]>0，则区间(a,(a+b)/2)内存在零点，(a+b)/2≤b；
返回①重新循环，不断接近零点。通过每次把f(x)的零点所在区间收缩一半的方法，使区间内的两个端点逐步逼近函数零点，最终求得零点近似值。
{
int z[10];
int maxi,maxj;
initdata();
for(int i=1;i<=N;i++)
z[i]=i;
for(int k=1;k<N;k++)
{
maxi=k;maxj=k;float maxv=abs(a[k][k]);
for(i=k;i<=N;i++)
for(int j=k;j<=N;j++)
34;请输入矩阵阶数:"<<endl;
cin>>N;
cout<<"请输入矩阵各项:"<<endl;
for(int i=1;i<=N;i++)
for(int j=1;j<=N+1;j++)
{
cin>>a[i][j];
}
cout<<endl;
}
void main()
{
for(i=1;i<=N;i++)
{
float t=a[i][k];a[i][k]=a[i][maxj];a[i][maxj]=t;

数值计算方法实验报告

数值计算⽅法实验报告《数值计算⽅法》实验报告实验题⽬⼆分法求⾮线性⽅程的根专业班级11级数学师范⼆班姓名李洪学号201102024056指导⽼师李梦联系电话188********⼀、实验⽬的熟悉⼆分法求⽅程近似根的数值⽅法，与⽤计算器解出的值进⾏⽐较，并学会误差分析。

⼆、实验原理⼆分法的基本思路是通过计算隔根区间的中点，逐步将隔根区间缩⼩，从⽽可得⽅程的近似根数列}{n x 。

（≤-+1*k x x ?）三、实验内容已知0)()3(3=-=-e x x f 在[]1,0上有⼀个实根*x ，0)1(0)0(>本实验中的⽤到的求根⽅法有①⼆分法，②计算器求根。

四、实验步骤1.输⼊：a ，b 值及精度控制?量；2.if 0)()(>b f a f then 返回第1步，重新输⼊a ，b 值else 转第3步；3.while ?>-b a 时做（1）)(21b a x +=，计算)(x f ；if )(x f =0 then 输出x ,停机。

（2）if0)()(4.输出)(21b a x +=。

五、 Matlab 源程序1.erfen.m:function [c,err,yc]=erfen(f,a,b,delta)ya = feval(f,a);yb = feval(f,b);if ya * yb > 0 ,return,endmax1 = 1+round((log(b-a)-log(delta))/log(2));for k=1:max1c=(a+b)/2;yc=feval(f,c);if yc==0a=c;b=c;elseif yb * yc > 0b=c;yb=yc;elsea=c;ya=yc;endif b-aendc=(a+b)/2;err=abs(b-a);yc=feval(f,c);2.f.m:function f=f(x);f=x^3-exp(-x);六、运⾏结果七、计算机计算结果⼋、实验分析1、⼆分法和计算器均能解出⽅程的根。

江苏大学数值计算方法实习报告

数值计算方法实习报告学院：电气信息工程学院班级：电气xxxx学号：xxxxxxxxxxxx姓名：xxxxxxxxxx1、 cos 2x x π=方程在[0,]上有唯一实根，试用快速弦截法求此根，要求精确到小数点后第3位。

#include<stdio.h>#include<math.h>#define PI 3.1415float f(float x){float y;y=x-cos(x);return (y);}main(){float x0=0,x1=PI/2,e=0.001,f0,f1,x2,f2;clrscr();f0=f(x0);f1=f(x1);while(fabs(f1)>e){ x2=x1-(x1-x0)*f1/(f1-f0);f2=f(x2);x0=x1;x1=x2;f0=f1;f1=f2;}printf("%f",x2);}运行结果：0.7395672、设A XB = 3112123312-1A ⎡⎤⎢⎥⎢⎥=⎢⎥⎢⎥⎣⎦ ，85-2B ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦，用Gauss-Jordan 消元法求方程组的解。

、#include<stdio.h>#include<math.h>main( ){floata[4][4]={0,0,0,0,0,3,1,1,0,2,2.0/3,1.0/3,0,1,2,-1 },b[4]={0,8,5,-2},s,t,max;int i,j,l,k;clrscr();for(k=1;k<=3;k++){max=fabs(a[k][k]);l=k;for(i=k;i<=3;i++)if(fabs(a[i][k])>max){ max=fabs(a[i][k]); l=i; }if(a[l][k]==0){ printf("qi yi ju zhen."); break; }if(l!=k){ for(j=k;j<=3;j++){ t=a[k][j];a[k][j]=a[l][j];a[l][j]=t; }s=b[l]; b[l]=b[k]; b[k]=s;}for(j=k+1;j<=3;j++)a[k][j]=a[k][j]/a[k][k];b[k]=b[k]/a[k][k];for(i=1;i<=3;i++){if(i==k) continue;b[i]=b[i]-a[i][k]*b[k];for(j=k+1;j<=3;j++)a[i][j]=a[i][j]-a[i][k]*a[k][j];}}printf("x1=%f\n",b[1]);printf("x2=%f\n",b[2]);printf("x3=%f\n",b[3]);}运行结果：x1=3.000000X2=-2.000000X3=1.0000003、已知函数表试用牛顿插值法求0.5635x 时的函数值。

数值计算方法实验报告

一、实验目的1. 熟悉数值计算的基本概念和方法；2. 掌握数值计算的基本原理和算法；3. 提高编程能力和数值计算能力；4. 通过实验，加深对数值计算方法的理解和应用。

二、实验内容1. 矩阵运算2. 线性方程组求解3. 函数求值4. 微分方程求解三、实验步骤1. 矩阵运算（1）编写程序实现矩阵的加法、减法、乘法运算；（2）编写程序实现矩阵的转置运算；（3）编写程序实现矩阵的逆运算。

2. 线性方程组求解（1）编写程序实现高斯消元法求解线性方程组；（2）编写程序实现雅可比迭代法求解线性方程组；（3）编写程序实现高斯-赛德尔迭代法求解线性方程组。

3. 函数求值（1）编写程序实现牛顿迭代法求函数的零点；（2）编写程序实现二分法求函数的零点；（3）编写程序实现割线法求函数的零点。

4. 微分方程求解（1）编写程序实现欧拉法求解一阶微分方程；（2）编写程序实现龙格-库塔法求解一阶微分方程；（3）编写程序实现龙格-库塔-法求解二阶微分方程。

四、实验结果与分析1. 矩阵运算（1）矩阵加法、减法、乘法运算结果正确；（2）矩阵转置运算结果正确；（3）矩阵逆运算结果正确。

2. 线性方程组求解（1）高斯消元法求解线性方程组，结果正确；（2）雅可比迭代法求解线性方程组，结果正确；（3）高斯-赛德尔迭代法求解线性方程组，结果正确。

3. 函数求值（1）牛顿迭代法求函数的零点，结果正确；（2）二分法求函数的零点，结果正确；（3）割线法求函数的零点，结果正确。

4. 微分方程求解（1）欧拉法求解一阶微分方程，结果正确；（2）龙格-库塔法求解一阶微分方程，结果正确；（3）龙格-库塔-法求解二阶微分方程，结果正确。

五、实验总结本次实验通过对数值计算方法的学习和实践，使我对数值计算有了更深入的了解。

以下是我对本次实验的总结：1. 矩阵运算是数值计算的基础，熟练掌握矩阵运算对于解决实际问题具有重要意义；2. 线性方程组求解是数值计算中常见的问题，高斯消元法、雅可比迭代法和高斯-赛德尔迭代法是常用的求解方法；3. 函数求值是数值计算中另一个常见问题，牛顿迭代法、二分法和割线法是常用的求解方法；4. 微分方程求解是数值计算中的难点，欧拉法、龙格-库塔法和龙格-库塔-法是常用的求解方法。

数值计算方法实验报告

数值分析实验报告实验一、解线性方程组的直接方法——梯形电阻电路问题利用追赶法求解三对角方程组的方法，解决梯形电阻电路问题：电路中的各个电流{1i ，2i ，…，8i }须满足下列线性方程组：R V i i =- 22 210 252321=-+-i i i 0 252 432=-+-i i i 0 252 543=-+-i i i 0 252 654=-+-i i i 0 252 765=-+-i i i 0 252 876=-+-i i i 052 87=+-i i设V 220=V ，Ω=27R ，运用追赶法，求各段电路的电流量。

问题分析：上述方程组可用矩阵表示为：⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡--------------00000001481.8522520000002520000002520000002520000002520000002520000002287654321i i i i i i i i问题转化为求解A x b =，8阶方阵A 满足顺序主子式(1,2...7)0i A i =≠，因此矩阵A存在唯一的Doolittle 分解，可以采用解三对角矩阵的追赶法！追赶法a=[0 -2 -2 -2 -2 -2 -2 -2]; b=[2 5 5 5 5 5 5 5]；c=[-2 -2 -2 -2 -2 -2 -2 0]; d=[220/27 0 0 0 0 0 0 0];Matlab 程序function x= zhuiganfa( a,b,c,d )%追赶法实现要求：|b1|>|C1|>0,|bi|>=|ai|+|ci| n=length(b); u=ones(1,n); L=ones(1,n); y=ones(1,n); u(1)=b(1); y(1)=d(1); for i=2:nL(i)=a(i)/u(i-1);u(i)=b(i)-c(i-1)*L(i); y(i)=d(i)-y(i-1)*L(i); endx(n)=y(n)/u(n); for k=n-1:-1:1x(k)=(y(k)-c(k)*x(k+1))/u(k); end endMATLAB 命令窗口输入：a=[0 -2 -2 -2 -2 -2 -2 -2]; b=[2 5 5 5 5 5 5 5];c=[-2 -2 -2 -2 -2 -2 -2 0] d=[220/27 0 0 0 0 0 0 0];x= zhuiganfa(a,b,c,d )运行结果为：x =8.1478 4.0737 2.0365 1.0175 0.5073 0.2506 0.1194 0.0477存在问题根据电路分析中的所讲到的回路电流法，可以列出8个以回路电流为独立变量的方程，课本上给出的第八个回路电流方程存在问题，正确的应该是78240i i -+=；或者可以根据电路并联分流的知识，同样可以确定78240i i -+=。

数值计算方法实验报告

数值分析实验报告实验一、解线性方程组的直接方法——梯形电阻电路问题利用追赶法求解三对角方程组的方法，解决梯形电阻电路问题：电路中的各个电流{1i ，2i ，…，8i }须满足下列线性方程组：R V i i =- 22 210 252321=-+-i i i 0 252 432=-+-i i i 0 252 543=-+-i i i 0 252 654=-+-i i i 0 252 765=-+-i i i 0 252 876=-+-i i i 052 87=+-i i设V 220=V ，Ω=27R ，运用追赶法，求各段电路的电流量。

问题分析：上述方程组可用矩阵表示为：⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡--------------00000001481.8522520000002520000002520000002520000002520000002520000002287654321i i i i i i i i问题转化为求解A x b =，8阶方阵A 满足顺序主子式(1,2...7)0i A i =≠，因此矩阵A存在唯一的Doolittle 分解，可以采用解三对角矩阵的追赶法！追赶法a=[0 -2 -2 -2 -2 -2 -2 -2]; b=[2 5 5 5 5 5 5 5]；c=[-2 -2 -2 -2 -2 -2 -2 0]; d=[220/27 0 0 0 0 0 0 0];Matlab 程序function x= zhuiganfa( a,b,c,d )%追赶法实现要求：|b1|>|C1|>0,|bi|>=|ai|+|ci| n=length(b); u=ones(1,n); L=ones(1,n); y=ones(1,n); u(1)=b(1); y(1)=d(1); for i=2:nL(i)=a(i)/u(i-1);u(i)=b(i)-c(i-1)*L(i); y(i)=d(i)-y(i-1)*L(i); endx(n)=y(n)/u(n); for k=n-1:-1:1x(k)=(y(k)-c(k)*x(k+1))/u(k); end endMATLAB 命令窗口输入：a=[0 -2 -2 -2 -2 -2 -2 -2]; b=[2 5 5 5 5 5 5 5];c=[-2 -2 -2 -2 -2 -2 -2 0] d=[220/27 0 0 0 0 0 0 0];x= zhuiganfa(a,b,c,d )运行结果为：x =8.1478 4.0737 2.0365 1.0175 0.5073 0.2506 0.1194 0.0477存在问题根据电路分析中的所讲到的回路电流法，可以列出8个以回路电流为独立变量的方程，课本上给出的第八个回路电流方程存在问题，正确的应该是78240i i -+=；或者可以根据电路并联分流的知识，同样可以确定78240i i -+=。

数值计算方法实验报告(含所有)

本科实验报告课程名称：计算机数值方法实验项目：计算机数值方法实验实验地点：虎峪校区致远楼B401专业班级：软件学院1217班学号：******xxxx 学生姓名：xxx指导教师：xxx2014 年 5 月21 日太原理工大学学生实验报告五、实验结果与分析二分法割线法分析：由程序知，使用二分法和割线法均能计算出方程的根，但利用割线法要比二分法计算的次数少，并且能够较早的达到精度要求。

相比之下，割线法程序代码量较少，精简明了。

六、讨论、心得本次数值计算方法程序设计实验从习题练习中跳脱出来，直接面对实用性较强的程序代码编写。

效果很好，不仅加深对二分法、割线法的理解，还加强了实际用运能力。

将理论知识成功地转化成实践结果。

实验地点虎峪校区致远楼B401指导教师xx太原理工大学学生实验报告l[i][k]=a[i][k];for(r=1;r<k;++r){l[i][k]-=l[i][r]*u[r][k];}l[i][k]/= u[k][k];}l[k][k]=1.0;}for(i=1;i<=n;++i){y[i] = b[i];for(j=1;j<i;++j){y[i]-=l[i][j]*y[j];}}for(i=n;i>0;--i){x[i] = y[i];for(j=i+1;j<=n;++j){x[i]-=u[i][j]*x[j];}x[i]/= u[i][i];}for(i=1;i<=n;++i){printf("%0.2lf\n",x[i]);}return 0;}五、实验结果与分析完全主元素消元法：列主元素消元法：LU分解法：分析：对于两种高斯解方程，完全主元素跟列主元素都是先消元、再回代，由程序段可以发现，始终消去对角线下方的元素。

即，为了节约内存及时效，可以不必计算出主元素下方数据。

列主元素消元法的算法设计上优于完全主元素消元法，它只需依次按列选主元素然后换行使之变到主元素位置，再进行消元即可。

数值计算方法实验报告

（实验报告的首页）本科实验报告课程名称：计算机数值方法实验项目：实验地点：多学科楼专业班级：力学1101 学号：2011005860 学生姓名：王亚博指导教师：刘晓燕2013年6月27日学生姓名王亚博实验成绩实验名称实验一：方程组求根1，用高斯消元法求解下面的方程组：⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-----8118344108318311231224321x x x x#include "stdio.h"double a[15][15],a0[15][15]; double b[15],b0[15],l[15]; int n; int i,j ;void displayA() {printf("\n");for( j=1;j<=n;j++) {for( i=1;i<=n;i++)printf("a[%d][%d]=%f",j,i,a[j][i]); printf("b[%d]=%f\n",j,b[j]); }for(j=1;j<=n;j++)printf("l[%d]=%f ",j,l[j]); printf("\n"); }void main() { int i,j,k;scanf("%d",&n); for(i=1;i<=n;i++) {for(j=1;j<=n;j++) {scanf("%lf",&a[i][j]); a0[i][j]=a[i][j]; }scanf("%lf",&b[i]); b0[i]=b[i]; }displayA(); k=1; do {for(i=1;i<=n;i++){if(i==k) continue;l[i]=a0[i][k]/a0[k][k];}for (j=k+1;j<=n;j++) a[k][j]=a0[k][j]/a0[k][k];b[k]=b0[k]/a0[k][k];for(i=1;i<=n;i++){if(i==k) continue;for(j=k+1;j<=n;j++)a[i][j]=a0[i][j]-l[i]*a0[k][j];b[i]=b0[i]-l[i]*b0[k];}displayA();for(i=1;i<=n;i++){for(j=k+1;j<=n;j++)a0[i][j]=a[i][j];b0[i]=b[i];}if(k==n) break;k++;}while(1);for(i=1;i<=n;i++)printf("b[%2d]=%lf\n",i,b[i]); getch();}实验名称实验二线性方程组的直接求解实验目的和要求合理选择利用Gauss 消元法、LU 分解法、追赶法求解下列方程组：①⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡13814142210321321x x x ②⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡--⨯-2178.4617.5911212592.1121130.6291.51314.59103.0432115x x x x ③⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡----=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡3772201161263841027851244321x x x x ④ ⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡----=⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡-55572112112112121 n n x x x x实验内容高斯消元法：找到与原方程组等价的系数矩阵为三角形方正的方程组：l ik =a ik /a kka ij = a ij - l ik * a kj k=1,2,…,n-1i=k+1,k+2, …,n j=k+1,k+2, …,n+1 由回代过程求得原方程组的解：x n = a nn+1/ a nnx k =( a kn+1-∑a kj x j )/ a kk (k=n-1,n-2, …,2,1)LU 分解：如果A 的各界顺序主子式不为0，则存在唯一的LU 分解。

数值计算方法实验报告

数值计算方法实验报告一、实验目的本实验旨在通过Python语言编写数值计算方法程序，掌握常见数值计算方法的实现原理及应用。

具体包括：插值法、最小二乘法、数值微积分、数值解方程、数值解微分方程等。

二、实验环境Python编程语言、Jupyter Notebook环境三、实验内容1.插值法（1）代码实现：在Python中使用Scipy库中的Interpolate模块实现拉格朗日插值法和牛顿插值法，并通过数据可视化展示其效果。

（2）实验步骤：- 导入所需库，准备所需数据；- 定义拉格朗日插值法函数；- 定义牛顿插值法函数；- 测试函数并可视化结果。

（3）实验结果：2.最小二乘法（1）代码实现：在Python中使用Numpy库实现最小二乘法，并通过数据可视化展示其效果。

（2）实验步骤：- 导入所需库，准备所需数据；- 定义最小二乘法函数；- 测试函数并可视化结果。

（3）实验结果：3.数值微积分（1）代码实现：在Python中实现梯形法和辛普森法，并通过数据可视化展示其效果。

（2）实验步骤：- 导入所需库，准备所需数据；- 定义梯形法函数和辛普森法函数；- 测试函数并可视化结果。

（3）实验结果：4.数值解方程（1）代码实现：在Python中实现二分法、牛顿法和割线法，并通过数据可视化展示其效果。

（2）实验步骤：- 导入所需库，准备所需数据；- 定义二分法函数、牛顿法函数和割线法函数；- 测试函数并可视化结果。

（3）实验结果：5.数值解微分方程（1）代码实现：在Python中实现欧拉法和龙格-库塔法，并通过数据可视化展示其效果。

（2）实验步骤：- 导入所需库，准备所需数据；- 定义欧拉法函数和龙格-库塔法函数；- 测试函数并可视化结果。

（3）实验结果：四、实验总结通过本次实验，我学习了数值计算方法的常用算法和实现原理，掌握了Python 语言实现数值计算方法的方法，加深了对数值计算方法的理解和应用。

实验中遇到的问题，我通过查找资料和与同学的讨论得到了解决，也更加熟练地掌握了Python语言的使用。

数值计算方法实验报告

数值计算方法实验报告实验目的：本次实验的目的是通过对数值计算方法的实践操作，加深对该方法的理解和掌握。

具体来说，本次实验旨在通过使用 MATLAB 软件对一些常见的数值计算问题进行求解，从而掌握和熟练运用一些数值计算方法，如插值、数值微积分、常微分方程数值解等。

实验过程：1.插值(1) Lagrange 插值法(2) Newton 插值法2.数值微积分(1) 梯形公式(2) Simpson 公式3.常微分方程数值解(1) 古典四步 Runge-Kutta 法(2) 改进四步 Runge-Kutta 法实验结果：本次实验中，我们使用 MATLAB 软件对以上数值计算问题进行了求解，成功得到了相应的数值解，并且通过分析和比较不同的数值计算方法的结果，得出了以下结论：1.在插值问题中，Lagrange 插值法和 Newton 插值法的结果相对较为接近，但是 Newton 插值法的计算速度更快。

2.在数值微积分问题中，梯形公式的结果较为精确，但是 Simpson 公式的精度更高。

3.在常微分方程数值解问题中，古典四步 Runge-Kutta 法和改进四步 Runge-Kutta 法均能得到较为准确的结果，但是改进四步Runge-Kutta 法的精度更高，尤其对于复杂的常微分方程求解有更好的效果。

实验结论：本次实验通过对数值计算方法的实践操作，深入理解了该方法的原理和运用，掌握了一些重要的数值计算方法，如插值、数值微积分、常微分方程数值解等，并且通过实验结果的分析比较，得出了相应的结论。

这些知识和技能对于我们在科研和工程实践中的数值计算问题具有非常重要的意义，具有广泛的应用前景。

大学计算实验实验报告

实验名称：数值计算方法实验日期：2021年10月15日实验地点：计算机实验室实验目的：1. 理解数值计算的基本原理和方法。

2. 掌握常用数值计算算法的编程实现。

3. 培养分析和解决实际问题的能力。

实验内容：本次实验主要涉及以下内容：1. 线性方程组的求解2. 函数求值与数值微分3. 函数求根4. 数据插值实验原理：数值计算是计算机科学和工程领域中非常重要的一个分支，它涉及到将数学问题转化为计算机可以处理的数值问题。

本实验主要探讨了以下数值计算方法：1. 高斯消元法：用于求解线性方程组。

2. 牛顿法：用于求解函数的根。

3. 二分法：用于求解函数的根。

4. 拉格朗日插值法：用于数据插值。

实验步骤：1. 线性方程组的求解：- 编写程序实现高斯消元法，用于求解线性方程组。

- 输入方程组的系数和常数项。

- 输出方程组的解。

2. 函数求值与数值微分：- 编写程序实现中点法和辛普森法，用于求函数的近似值。

- 编写程序实现中心差分法，用于求函数的导数近似值。

- 输入函数表达式、求值点和微分点。

- 输出函数的近似值和导数的近似值。

3. 函数求根：- 编写程序实现牛顿法和二分法，用于求解函数的根。

- 输入函数表达式、初始猜测值和误差容忍度。

- 输出函数的根。

4. 数据插值：- 编写程序实现拉格朗日插值法，用于数据插值。

- 输入数据点和待插值点。

- 输出插值结果。

实验结果：1. 线性方程组的求解：- 输入方程组系数和常数项：`2x + 3y = 6`，`x - y = 1`。

- 输出解：`x = 2`，`y = 1`。

2. 函数求值与数值微分：- 输入函数表达式：`f(x) = x^2`，求值点：`x = 2`。

- 输出函数的近似值：`f(2) ≈ 4.000000`。

- 输入微分点：`x = 2`。

- 输出导数的近似值：`f'(2) ≈ 4.000000`。

3. 函数求根：- 输入函数表达式：`f(x) = x^2 - 2`，初始猜测值：`x = 1`，误差容忍度：`10^-6`。

数值计算方法实验报告

数值计算方法实验报告数值计算方法实验报告引言：数值计算方法是一种通过数学模型和计算机算法来解决实际问题的方法。

在科学研究和工程应用中，数值计算方法被广泛应用于求解方程、优化问题、模拟仿真等领域。

本实验报告将介绍数值计算方法的基本原理和实验结果。

一、二分法求根二分法是一种通过不断折半缩小搜索区间来求解方程根的方法。

在实验中，我们选取了一个简单的方程f(x) = x^2 - 4 = 0来进行求根实验。

通过不断将搜索区间进行二分，我们可以逐步逼近方程的根。

实验结果表明，通过二分法，我们可以得到方程的根为x = 2。

二、牛顿迭代法求根牛顿迭代法是一种通过不断逼近方程根的方法。

在实验中，我们同样选取了方程f(x) = x^2 - 4 = 0进行求根实验。

牛顿迭代法的基本思想是通过对方程进行线性近似，求得近似解，并不断迭代逼近方程的根。

实验结果表明，通过牛顿迭代法，我们可以得到方程的根为x = 2。

三、高斯消元法求解线性方程组高斯消元法是一种通过变换线性方程组的系数矩阵，将其化为上三角矩阵的方法。

在实验中，我们选取了一个简单的线性方程组进行求解实验。

通过对系数矩阵进行行变换，我们可以将其化为上三角矩阵，并通过回代求解得到方程组的解。

实验结果表明，通过高斯消元法，我们可以得到线性方程组的解为x = 1，y = 2，z = 3。

四、插值与拟合插值与拟合是一种通过已知数据点来构造函数模型的方法。

在实验中，我们选取了一组数据点进行插值与拟合实验。

通过拉格朗日插值多项式和最小二乘法拟合，我们可以得到数据点之间的函数模型。

实验结果表明，通过插值与拟合，我们可以得到数据点之间的函数关系，并可以通过该函数模型来进行预测和拟合。

结论：数值计算方法是一种通过数学模型和计算机算法来解决实际问题的方法。

通过本次实验，我们学习了二分法求根、牛顿迭代法求根、高斯消元法求解线性方程组以及插值与拟合的基本原理和应用。

这些方法在科学研究和工程应用中具有广泛的应用前景。

数值计算方法实验报告

#include<stdio.h>
#include<math.h>
double f(double x)
{
double s;
s=x*x*x/3-x;
return fabs(s);
}
void main()
{double x=-0.99,y;
int k=0;
printf("%d ,%lf\n",k,x);
{if(r>=x[i]&&r<=x[i+1])
{s=m[i]*pow(x[i+1]-r,3)/6*h[i]+m[i+1]*pow(r-x[i],3)/6*h[i]+(y[i]-m[i]*pow(h[i],2)/6)*(x[i+1]-r)/h[i]+(y[i+1]-m[i+1]*pow(h[i],2)/6)*(r-x[i])/h[i];
28.65
39.62
50.65
5.28794
9.4
13.84
20.2
24.9
28.44
31.1
k
7
8
9
10
11
12
78
104.6
156.6
208.6
260.7
312.5
35
36.5
36.6
34.6
31.6
31.0
k
13
14
15
16
17
18
364.4
416.3
468
494
507
520
20.9
14.8
7.8
do
{y=x;

数值计算方法实验报告

数值计算方法实验报告一、实验目的本实验旨在通过数值计算方法的实验操作，深入理解数值计算方法的原理与应用，掌握数值计算方法的相关技能，提高数值计算方法的实际应用能力。

二、实验内容1.数值微积分2.数值代数3.数值微分方程4.数值线性代数5.数值优化6.数值统计分析7.数值随机模拟8.数值傅立叶分析9.数值偏微分方程三、实验步骤1.数值微积分：通过不同的数值积分方法，计算给定函数的定积分值，并对不同数值积分方法的误差进行分析。

2.数值代数：通过使用线性代数方法，求解给定的线性方程组，并分析不同线性方程组求解方法的优劣。

3.数值微分方程：通过使用常微分方程数值解法，求解给定的微分方程，并比较不同求解方法的精度和稳定性。

4.数值线性代数：通过使用特征值分解方法，对给定的矩阵进行特征值分解，并分析不同特征值分解方法的优缺点。

5.数值优化：通过使用不同的优化方法，求解给定的优化问题，并比较不同的优化方法的效率和精度。

6.数值统计分析：通过使用不同的统计分析方法，对给定的数据进行统计分析，并分析不同的统计方法的优缺点。

7.数值随机模拟：通过使用随机模拟方法，模拟给定的概率分布，并分析不同随机模拟方法的效率和精度。

8.数值傅立叶分析：通过使用傅立叶分析方法，对给定的信号进行频谱分析，并分析不同的傅立叶分析方法的优缺点。

9.数值偏微分方程：通过使用偏微分方程数值解法，求解给定的偏微分方程，并比较不同求解方法的精度和稳定性。

四、实验结果与分析本实验中，通过对不同的数值计算方法的实验操作，我们可以更深入地理解数值计算方法的原理与应用，并掌握数值计算方法的相关技能，提高数值计算方法的实际应用能力。

同时，通过实验结果的分析，我们可以更好地比较不同数值计算方法的优缺点，为实际应用提供参考依据。

五、实验总结本实验旨在通过数值计算方法的实验操作，深入理解数值计算方法的原理与应用，掌握数值计算方法的相关技能，提高数值计算方法的实际应用能力。

数值计算的实验报告

一、实验目的1. 熟悉数值计算的基本原理和方法。

2. 掌握常用的数值计算算法及其应用。

3. 提高数值计算软件的使用能力。

4. 培养分析问题和解决问题的能力。

二、实验环境1. 操作系统：Windows 102. 编程语言：Python3. 数值计算软件：NumPy、SciPy、Matplotlib三、实验内容1. 实验一：数值积分（1）实验目的：学习数值积分方法，计算定积分的近似值。

（2）实验内容：a. 使用辛普森法则计算函数f(x) = x^2在区间[0, 1]上的定积分。

b. 使用梯形法则计算函数f(x) = e^x在区间[0, 1]上的定积分。

（3）实验步骤：a. 编写Python代码，实现辛普森法则和梯形法则。

b. 分别使用两种方法计算定积分的近似值。

c. 对比两种方法的计算结果，分析误差来源。

2. 实验二：数值微分（1）实验目的：学习数值微分方法，计算函数在某点的导数近似值。

（2）实验内容：a. 使用中心差分法计算函数f(x) = sin(x)在x = π/2处的导数近似值。

b. 使用前向差分法和后向差分法计算函数f(x) = cos(x)在x = 0处的导数近似值。

（3）实验步骤：a. 编写Python代码，实现中心差分法、前向差分法和后向差分法。

b. 分别使用三种方法计算导数的近似值。

c. 对比三种方法的计算结果，分析误差来源。

3. 实验三：线性方程组求解（1）实验目的：学习线性方程组求解方法，掌握高斯消元法和迭代法。

（2）实验内容：a. 使用高斯消元法求解线性方程组：3x + 2y - z = 72x - y + 3z = -1-x + 2y + 2z = 4b. 使用雅可比迭代法求解线性方程组：3x + 2y - z = 72x - y + 3z = -1-x + 2y + 2z = 4（3）实验步骤：a. 编写Python代码，实现高斯消元法和雅可比迭代法。

b. 分别使用两种方法求解线性方程组。