合肥工业大学 计算机专业 计算方法实验报告

合肥工业大学计算机与信息学院

实验报告

课程：计算方法

专业班级：

学号：

姓名：

Java界面

其实都不难按照程序流程图就可以完成了

实验一 插值与拟合

一、实验目的

(1) 明确插值多项式和分段插值多项式各自的优缺点；

(2) 编程实现三次样条插值算法，分析实验结果体会高次插值产生的龙格现象；

(3) 理解最小二乘拟合，并编程实现线性拟合，掌握非线性拟合转化为线性拟合的方法 (4) 运用常用的插值和拟合方法解决实际问题。

二、实验内容

(1)对于f(x)=1/(1+x*x)实现三次样条插值 (2)实现最小二乘法的直线拟合

数据如下：

j x 165 123 150 123 141 j y

187

126

172

125

148

三、基本原理（计算公式）

(1)三次样条插值在每个内节点上具有2阶导数。

(2)最小二乘法拟合直线为y=a+bx ，而a ，b 有如下等式（N 为给出的数据点的总个数）

i i y x b N =+∑a ； ∑∑∑=+i i i i y x x b x 2a

四、算法设计与实现（流程图，关键点）

最小二乘法直线拟合：输入数据后，按照公式计算a，b。用得到的拟合直线计算预测点的近似函数值。

五、输入与输出

(1)三次样条插值

输入：区间长度，n+1个数据点，预测点

输出：预测点的近似函数值，精确值，及误差

(2)最小二乘法直线拟合

输入：n个数据点，预测点

输出：预测点的近似函数值

六、结果讨论和分析

代码

三次样条插值

#include

#include

#define N 10

using namespace std;

double u0(double x){

return (x-1)*(x-1)*(2*x+1);

}

double u1(double x){

return x*x*(3-2*x);

}

double v0(double x){

return x*(x-1)*(x-1);

}

double v1(double x){

return x*x*(x-1);

}

double s3(double x,double y,double y1,double m,double m1,double h){ return u0(x)*y+u1(x)*y1+h*v0(x)*m+h*v1(x)*m1;

}

double f(double x){

return 1/(1+x*x);

}

int main(){

ifstream fin;

fin.open("E:\\t.txt");

if(!fin){

cout<<"error opening input stream"<

system("pause");

return 0;

}

double x[N+1],y[N+1],m[N+1],A[N],B[N],C[N];

double h[N];

double a[N],b[N];

double f0,fn;

double temp;

int i;

for(i=0;i<=N;i++){

fin>>x[i]>>y[i];

}

fin>>f0>>fn;

h[0]=x[1]-x[0];

for(i=1;i

h[i]=x[i+1]-x[i];

a[i]=h[i-1]/(h[i-1]+h[i]);

b[i]=3*((1-a[i])*(y[i]-y[i-1])/h[i-1]+a[i]*(y[i+1]-y[i])/h[i]);

}

m[1]=b[1]-(1-a[1])*f0;

m[N-1]=b[N-1]-a[N-1]*fn;

for(i=2;i

m[i]=b[i];

}

for(i=1;i

B[i]=2;

C[i]=a[i];

}

for(i=2;i

A[i]=1-a[i];

}

C[1]=C[1]/B[1];

m[1]=m[1]/B[1];

double t;

for(i=2;i!=N-2;i++){

t=B[i]-C[i-1]*A[i];

C[i]=C[i]/t;

m[i]=(m[i]-m[i-1]*A[i])/t;

}

m[N-1]=(m[N-1]-m[N-2]*A[N-1])/(B[N-1]-C[N-2]*A[N-1]);

for(i=N-2;i>0;i--){

m[i]=m[i]-C[i]*m[i+1];

}

cout<<"please:(输入插值节点在"<

while(cin>>temp){

double tt=temp;

if(tempx[N]){

cout<<"插值节点为"<

continue;

}

for(i=1;i<=N;i++){

if(temp

}

temp=(temp-x[i-1])/h[i-1];

temp=s3(temp,y[i-1],y[i],m[i-1],m[i],h[i-1]);

cout<<"插值节点为"<

}

system("pause");

fin.close();

return 0;

}

最小二乘法的直线拟合

#include

#include

#define n 5

using namespace std;

double sum(double x[],int k){

int i;

double sum=0;

for(i=0;i

sum=sum+x[i];

return sum;

}

double sum2(double x[],int k){

int i;

double sum=0;

for(i=0;i

sum=sum+x[i]*x[i];

return sum;

}

double sumxy(double x[],double y[],int k){

int i;

double sum=0;

for(i=0;i

sum=sum+x[i]*y[i];

return sum;

}

int main(){

ifstream fin;

fin.open("E:\\t.txt");

if(!fin){

cout<<"error opening input stream"<

system("pause");

return 0;

}

double x[n],y[n],a,b;

double x0,y0;

int i;

for(i=0;i

fin>>x[i]>>y[i];

}

b=(n*sumxy(x,y,n)-sum(x,n)*sum(y,n))/(n*sum2(x,n)-sum(x,n)*sum(x,n));

a=(sum(y,n)-b*sum(x,n))/n;

cout<<"最小二乘法直线拟合得到a: "<

cout<<"请输入插值节点x：";

while(cin>>x0){

y0=a+b*x0;

cout<<"当x="<

cout<<"请输入插值节点x：";

}

system("pause");

fin.close();

return 0;

}

实验二数值积分

一、实验目的

(1)熟悉复化梯形方法、复化Simpson方法、梯形递推算法、龙贝格算法；

(2)能编程实现龙贝格算法和中点加速；

(3)理解并掌握自适应算法和收敛加速算法的基本思想；

(4)分析实验结果体会各种方法的精确度，建立计算机求解定积分问题的感性认识

二、实验内容 (1)用龙贝格算法计算

dx x

x

?

1

sin (2)用中点加速方法计算x

e 的一阶导数

三、基本原理（计算公式）

(1)龙贝格算法

梯形递推公式∑-=++=1

02

12)(22n k k n n x f h T T

加权平均公式：

n n n S T T =--1442 n n n C S S =--144222 n n

n R C C =--1

44323 (2)中点加速

中点公式: G(h)=(f(a+h)-f(a-h))/2/h 加权平均:G1(h)=4*G(h/2)/3-G(h)/3 G2(h)=16*G1(h/2)/15-G1(h)/15 G3(h)=64*G2(h/2)/63-G2(h)/63

四、算法设计与实现（流程图，关键点）

中点加速：输入数据后根据公式计算导数值

五、输入与输出

图2.2梯形递推算法流程图

图2.3龙贝格算法流程图

(1) 用龙贝格算法计算

dx x

x

1

sin 输入：积分区间，误差限

输出：序列Tn ，Sn,Cn,Rn 及积分结果 (2)用中点加速方法计算x

e 的一阶导数 输入：求导节点，步长

输出：求得的导数值，精确值

六、结果讨论和分析

代码

龙贝格算法

#include

#include

#include

using namespace std;

double f(double x){

if(x==0)return 1;

return sin(x)/x;

}

int main(){

ifstream fin;

fin.open("E:\\t.txt");

if(!fin){

cout<<"error opening input stream"<

system("pause");

return 0;

}

double a,b,e,t1,t2,s1,s2,c1,c2,r1,r2;

double x,h,s;

fin>>a>>b>>e;

cout<<"积分区间为["<

cout<<"k T2 S2 C2 R2"<

h=b-a;

t1=(f(a)+f(b))*h/2;

cout<<0<<" "<

int k;

for(k=1;k<=10;k++,h=h/2,t1=t2,s1=s2){

s=0;

x=a+h/2;

do{

s=s+f(x);

x=x+h;

}while(x

t2=t1/2+h*s/2;

s2=t2+(t2-t1)/3;

if(k==1){

cout<

continue;

}

c2=s2+(s2-s1)/15;

if(k==2){

cout<

c1=c2;

continue;

}

r2=c2+(c2-c1)/63;

cout<

if(k==3){

r1=r2;

c1=c2;

continue;

}

if(fabs(r2-r1)

cout<<"数值积分结果为"<

break;

}

r1=r2;

c1=c2;

}

system("pause");

return 0;

}

中点加速算法

#include

#include

#include

using namespace std;

double f(double x){

return exp(x);

}

double f1(double x){

return exp(x);

}

double g(double x,double h){

return (f(x+h)-f(x-h))/2/h;

}

double g1(double x,double h){

return 4*g(x,h/2)/3-g(x,h)/3;

}

double g2(double x,double h){

return 16*g1(x,h/2)/15-g1(x,h)/15;

}

double g3(double x,double h){

return 64*g2(x,h/2)/63-g2(x,h)/63;

}

int main(){

ifstream fin;

fin.open("E:\\t.txt");

if(!fin){

cout<<"error opening input stream"<

system("pause");

return 0;

}

double a,h;

while(fin>>a>>h)

cout<<"当x="<

system("pause");

fin.close();

return 0;

}

实验三非线性方程求根迭代法

一、实验目的

(1)熟悉非线性方程求根简单迭代法，牛顿迭代及牛顿下山法

(2)能编程实现牛顿下山法

(3)认识选择迭代格式的重要性

(4)对迭代速度建立感性的认识；分析实验结果体会初值对迭代的影响

二、实验内容

用牛顿下山法解方程013

=--x x （初值为0.6）

三、基本原理（计算公式）

求非线性方程组的解是科学计算常遇到的问题，有很多实际背景．各种算法层出不穷，其中迭代是主流算法。只有建立有效的迭代格式，迭代数列才可以收敛于所求的根。因此设计算法之前，对于一般迭代进行收敛性的判断是至关重要的。牛顿法也叫切线法，是迭代算法中典型方法，只要初值选取适当，在单根附近，牛顿法收敛速度很快，初值对于牛顿迭代 至关重要。当初值选取不当可以采用牛顿下山算法进行纠正。 一般迭代：)(1k k x x φ=+ 0)()(=?=x f x x φ

牛顿公式：)

()

(1k k k k x f x f x x '-

=+ 牛顿下山公式：)

()

(1k k k k x f x f x x '-=+λ

下山因子 ，

，，，322

1

21211=λ 下山条件

|)(||)(|1k k x f x f <+

四、算法设计与实现（流程图，关键点）

图3.2牛顿下山算法流程

???

?

五、输入与输出

输入：初值，误差限，迭代最大次数，下山最大次数输出：近似根各步下山因子

六、结果讨论和分析

代码

牛顿下山法

#include

#include

#include

using namespace std;

double f(double x){//函数

return x*x*x-x-1;

}

double f1(double x){//一阶导数

return 3*x*x-1;

}

int main(){

ifstream fin;

fin.open("E:\\t.txt" );

if(!fin){

cout << "Error opening input stream" << endl;

system("pause");

return 0;

}

double x0,x1,e,j,temp;

int N,M,k,i;

fin>>x0>>e>>M>>N;

cout<<"迭代初值为"<

x1=0;

for(k=0;k

i=0;

j=1;

temp=f(x0)/f1(x0);

do{

x1=x0-temp/j;

cout<<"x0:"<

i++;

j=j*2;

if(i>=M)break;

}while(fabs(f(x1))>=fabs(f(x0)));

if(i>=M){

cout<<"重新选择x0"<

break;

}

if(fabs(x1-x0)

cout<<"迭代成功,求得结果为"<

break;

}

x0=x1;

}

if(k>=N)cout<<"迭代失败"<

if(f1(x0)==0)cout<<"一阶导数为0"<

fin.close();

system("pause");

return 0;

}

实验四求解线性方程组

一、实验目的

(1)熟悉求解线性方程组的有关理论和方法；

(2)能编程实现高斯-塞德尔迭代法、列主元高斯消去法、LU分解法

(3)通过测试，进一步了解各种方法的优缺点

(4) 根据不同类型的方程组，选择合适的数值方法

二、实验内容

(1)用高斯-塞德尔迭代法求???

??=+--=-+-=--2

.453.82102.7210321321321x x x x x x x x x

(2)列主元高斯消去法求??

?

??=+--=-+-=--2

.453.82102

.7210321321321x x x x x x x x x

(3)LU 分解法求解方程组

三、基本原理（计算公式）

线性方程组大致分迭代法和直接法。只有收敛条件满足时，才可以进行迭代。雅可比及高斯-塞德尔是最基本的两类迭代方法，最大区别是迭代过程中是否引用新值进行剩下的计算。消元是最简单的直接法，并且也十分有效的，列主元高斯消去法对求解一般的线性方程组都适用，同时可以用来求矩阵对应的行列式。约当消去实质是经过初等行变换将系数矩阵化为单位阵，主要用来求矩阵的逆。在使用直接法，要注意从空间、时间两方面对算法进行优化。

高斯-塞德尔迭代：???

? ??--=∑∑+=-=++n i j k j ij i j k j ij i ii

k i

x a x a b a x 1)(11)1()

1(1 n i ,...,2,1= 列主元高斯消去法：列主元;0||max ||≠=≤≤ik n

i k lk a a

消元)

...,,1,()()()1()()()1(n k j i b

m b b a

m a a k k

ik k i k i k kj

ik k ij k ij +=???-=-=++ 回代)1...,,()(1)()(n i a

x a

b x i ii

n

i j j

i ij

i i

i =-

=

∑+=

图4.3列主元的高斯消去流程图

图4.1G-S 迭代算法流程图

四、算法设计与实现（流程图，关键点）

LU分解法：依次求得L、U、y和x

五、输入与输出

(1)用高斯-塞德尔迭代法

输入：系数矩阵A，最大迭代次数N，初始向量，误差限e 输出：解向量

(2)列主元高斯消去法

输入：系数矩阵A

输出：解向量

(3)LU分解法

输入：系数矩阵A

输出：解向量

六、结果讨论和分析

《计算方法》课内实验报告

《计算方法》实验报告 姓名： 班级： 学号： 实验日期: 2011年10月26日

一、实验题目： 数值积分 二、实验目的： 1．熟悉matlab 编写及运行数值计算程序的方法。 2．进一步理解数值积分的基础理论。 3．进一步掌握应用不同的数值积分方法求解给定的积分并给出数据结果及误差分析。 三、实验内容： 1．分别用复合梯形求积公式及复合辛普森求积公式计算积分xdx x ln 10 ? , 要求计算精度达到410-,给出计算结果并比较两种方法的计算节点数. 2．用龙贝格求积方法计算积分dx x x ?+3 021,使误差不超过510-. 3．用3=n 的高斯-勒让德公式计算积分?3 1 sin x e x ,给出计算结果. 4．用辛普森公式(取2==M N ) 计算二重积分.5 .00 5 .00 dydx e x y ? ? - 四、实验结果： 1.（1）复合梯形法: 将区间[a,b]划分为n 等份，分点n k n a b h kh a x k ,2,1,0,,=-=+=在每个区间[1,+k k x x ](k=0,1,2,···n-1)上采用梯形公式，则得 )()]()([2)()(1 11 1 f R x f x f h dx x f dx x f I n n k k k b a n k x x k k ++===∑?∑? -=+-=+ 故)]()(2)([21 1 b f x f a f h T n k k n ++=∑-=称为复合梯形公式 计算步长和划分的区间 Eps=1E-4 h1=sqrt(Eps/abs(-(1-0)/12*1/(2+1))) h1 =0.0600 N1=ceil(1/h1) N1 =17 用复合梯形需要计算17个结点。 复合梯形： function T=trap(f,a,b,n) h=(b-a)/n;

计算方法上机实验报告

《计算方法》上机实验报告 班级：XXXXXX 小组成员：XXXXXXX XXXXXXX XXXXXXX XXXXXXX 任课教师：XXX 二〇一八年五月二十五日

前言 通过进行多次的上机实验，我们结合课本上的内容以及老师对我们的指导，能够较为熟练地掌握Newton 迭代法、Jacobi 迭代法、Gauss-Seidel 迭代法、Newton 插值法、Lagrange 插值法和Gauss 求积公式等六种算法的原理和使用方法，并参考课本例题进行了MATLAB 程序的编写。 以下为本次上机实验报告，按照实验内容共分为六部分。 实验一： 一、实验名称及题目： Newton 迭代法 例2.7(P38):应用Newton 迭代法求 在 附近的数值解 ，并使其满足 . 二、解题思路： 设'x 是0)(=x f 的根,选取0x 作为'x 初始近似值,过点())(,00x f x 做曲线)(x f y =的切线L ,L 的方程为))((')(000x x x f x f y -+=,求出L 与x 轴交点的横坐标) (') (0001x f x f x x - =,称1x 为'x 的一次近似值,过点))(,(11x f x 做曲线)(x f y =的切线,求该切线与x 轴的横坐标) (') (1112x f x f x x - =称2x 为'x

的二次近似值,重复以上过程,得'x 的近似值序列{}n x ,把 ) (') (1n n n n x f x f x x - =+称为'x 的1+n 次近似值，这种求解方法就是牛顿迭代法。 三、Matlab 程序代码： function newton_iteration(x0,tol) syms z %定义自变量 format long %定义精度 f=z*z*z-z-1; f1=diff(f);%求导 y=subs(f,z,x0); y1=subs(f1,z,x0);%向函数中代值 x1=x0-y/y1; k=1; while abs(x1-x0)>=tol x0=x1; y=subs(f,z,x0); y1=subs(f1,z,x0); x1=x0-y/y1;k=k+1; end x=double(x1) K 四、运行结果： 实验二：

c计算器实验报告

简单计算器 姓名: 周吉祥 实验目的：模仿日常生活中所用的计算器，自行设计一个简单的计算器程序，实现简单的计算功能。 实验内容： (1)体系设计： 程序是一个简单的计算器，能正确输入数据，能实现加、减、乘、除等算术运算，运算结果能正确显示，可以清楚数据等。 (2)设计思路： 1)先在Visual C++ 6.0中建立一个MFC工程文件，名为calculator. 2)在对话框中添加适当的编辑框、按钮、静态文件、复选框和单选框 3)设计按钮，并修改其相应的ID与Caption. 4)选择和设置各控件的单击鼠标事件。 5)为编辑框添加double类型的关联变量m_edit1. 6)在calculatorDlg.h中添加math.h头文件，然后添加public成员。 7)打开calculatorDlg.cpp文件，在构造函数中，进行成员初始化和完善各控件 的响应函数代码。

(3)程序清单： 添加的public成员： double tempvalue; //存储中间变量 double result; //存储显示结果的值 int sort; //判断后面是何种运算：1.加法 2.减法 3.乘法 4.除法 int append; //判断后面是否添加数字 成员初始化： CCalculatorDlg::CCalculatorDlg(CWnd* pParent /*=NULL*/) : CDialog(CCalculatorDlg::IDD, pParent) { //{{AFX_DATA_INIT(CCalculatorDlg) m_edit1 = 0.0; //}}AFX_DATA_INIT // Note that LoadIcon does not require a subsequent DestroyIcon in Win32 m_hIcon = AfxGetApp()->LoadIcon(IDR_MAINFRAME); tempvalue=0;

太原理工大学数值计算方法实验报告

本科实验报告 课程名称：计算机数值方法 实验项目：方程求根、线性方程组的直接解 法、线性方程组的迭代解法、代数插值和最 小二乘拟合多项式 实验地点：行勉楼 专业班级： ******** 学号： ********* 学生姓名： ******** 指导教师：李誌，崔冬华 2016年 4 月 8 日

y = x*x*x + 4 * x*x - 10; return y; } float Calculate(float a,float b) { c = (a + b) / 2; n++; if (GetY(c) == 0 || ((b - a) / 2) < 0.000005) { cout << c <<"为方程的解"<< endl; return 0; } if (GetY(a)*GetY(c) < 0) { return Calculate(a,c); } if (GetY(c)*GetY(b)< 0) { return Calculate(c,b); } } }; int main() { cout << "方程组为：f(x)=x^3+4x^2-10=0" << endl; float a, b; Text text; text.Getab(); a = text.a; b = text.b; text.Calculate(a, b); return 0; } 2.割线法： // 方程求根（割线法）.cpp : 定义控制台应用程序的入口点。// #include "stdafx.h" #include"iostream"

心得体会 使用不同的方法，可以不同程度的求得方程的解，通过二分法计算的程序实现更加了解二分法的特点，二分法过程简单，程序容易实现，但该方法收敛比较慢一般用于求根的初始近似值，不同的方法速度不同。面对一个复杂的问题，要学会简化处理步骤，分步骤一点一点的循序处理，只有这样，才能高效的解决一个复杂问题。

物理实验报告格式范文

物理实验报告格式范文 一、实验目的 二、实验仪器和器材(要求标明各仪器的规格型号) 三、实验原理：简明扼要地阐述实验的理论依据、计算公式、画出电路图或光路图 四、实验步骤或内容：要求步骤或内容简单明了 五、数据记录：实验中测得的原始数据和一些简单的结果尽可能用表格形式列出,并要求正确表示有效数字和单位 六、数据处理：根据实验目的对测量结果进行计算或作图表示,并对测量结果进行评定,计算误差或不确定度. 七、实验结果：扼要地写出实验结论 八、误差分析：当实验数据的误差达到一定程度后,要求对误差进行分析,找出产生误差的原因. 九、问题讨论：讨论实验中观察到的异常现象及可能的解释,分析实验误差的主要来源,对实验仪器的选择和实验方法的改进提出建议,简述自己做实验的心得体会,回答实验思考题. 物理探究实验：影响摩擦力大小的因素 技能准备：弹簧测力计，长木板，棉布，毛巾，带钩长方体木块，砝码，刻度尺，秒表。 知识准备： 1. 二力平衡的条件：作用在同一个物体上的两个力，如果大小相等，方向相反，并且在同一直线上，这两个力就平衡。 2. 在平衡力的作用下，静止的物体保持静止状态，运动的物体保持匀速直线运动状态。 3. 两个相互接触的物体，当它们做相对运动时或有相对运动的趋势时，在接触面上会产生一种阻碍相对运动的力，这种力就叫摩擦力。 4. 弹簧测力计拉着木块在水平面上做匀速直线运动时，拉力的大小就等于摩擦力的大小，拉力的数值可从弹簧测力计上读出，这样就测出了木块与水平面之间的摩擦力。

探究导引 探究指导： 关闭发动机的列车会停下来，自由摆动的秋千会停下来，踢出去的足球会停下来，运动的物体之所以会停下来，是因为受到了摩擦力。 运动物体产生摩擦力必须具备以下三个条件：1.物体间要相互接触，且挤压;2.接触面要粗糙;3.两物体间要发生相对运动或有相对运动的趋势。三个条件缺一不可。 摩擦力的作用点在接触面上，方向与物体相对运动的方向相反。由力的三要素可知：摩擦力除了有作用点、方向外，还有大小。 提出问题：摩擦力大小与什么因素有关? 猜想1：摩擦力的大小可能与接触面所受的压力有关。 猜想2：摩擦力的大小可能与接触面的粗糙程度有关。 猜想3：摩擦力的大小可能与产生摩擦力的两种物体间接触面积的大小有关。 探究方案： 用弹簧测力计匀速拉动木块，使它沿长木板滑动，从而测出木块与长木板之间的摩擦力;改变放在木块上的砝码，从而改变木块与长木板之间的压力;把棉布铺在长木板上，从而改变接触面的粗糙程度;改变木块与长木板的接触面，从而改变接触面积。 物理实验报告 .化学实验报告 .生物实验报告 .实验报告格式 .实验报告模板 探究过程： 1. 用弹簧测力计匀速拉动木块，测出此时木块与长木板之间的摩擦力：0.7N 2. 在木块上加50g的砝码，测出此时木块与长木板之间的摩擦力：0.8N 3. 在木块上加200g的砝码，测出此时木块与长木板之间的摩擦力：1.2N 4. 在木板上铺上棉布，测出此时木块与长木板之间的摩擦力：1.1N 5. 加快匀速拉动木块的速度，测出此时木块与长木板之间的摩擦力：0.7N 6. 将木块翻转，使另一个面积更小的面与长木板接触，测出此时木块与长木板之间的摩擦力：0.7N 探究结论：

实验报告格式

重庆工商大学 《统计学》实验报告 实验课程：统计学 _ 指导教师：陈正伟 _ 专业班级： 08 经济学 学生姓名：程剑波 学生学号： 2008011133 __

实验项目 实验日期实验地点80608 实验目的掌握统计学的基本计算方法和分析方法。 实验内容一、统计图绘制；二、动差、偏度系数、峰度系数的计算；三、趋势性的绘制； 四、相关分析与回归分析；五、时间数列的动态指标分析；六、循环变动的测 算分析。 通过统计学（2009.9.10-2009.12.15）实验报告如下： 一、统计图绘制； （一）过程： （二）结果： （三）分析： 二、动差、偏度系数、峰度系数的计算； （一）过程： （二）结果： （三）分析： 三、趋势性的绘制； （一）过程： （二）结果： （三）分析： 四、相关分析与回归分析； （一）过程： （二）结果： （三）分析：

五、时间数列的动态指标分析 （一）过程： （二）结果： （三）分析： 六、循环变动的测算分析。 （一）过程： （二）结果： （三）分析： 体会： 参考实验报告： 重庆工商大学数学与统计学院 综合评价方法及应用 实验报告

实验课程：非参数统计 _ 指导教师：陈正伟 _ 专业班级： 06市调2班 学生姓名：何春 学生学号： 2006004151 _

实验报告一 实验项目变异系数法相关系数法熵值发坎蒂雷法 实验日期2009-4-30 实验地点80608 实验目的 通过本实验本要求掌握综合评价指标体系中各个指标重要性权数的重要意义；掌握权数确定的定性和定量技术和技能；解决实际综合评价中重要性权数确定的处理技能。 实验内容 根据资料使用变异系数法、相关系数法、熵值法和坎蒂雷方法分别确定各个指标的权数。并进行权数比较分析。 检验方法的选择及实验步骤及结果： 1用变异系数求各个指标的权数： 基本步骤：（1）先求各个指标的均值Xi 和标准差 Si (2)接着求各个指标的变异系数Vi=Si/Xi (3)对Vi作作归一化处理，及得各个指标的权数 结果如下： 从这个表中可以看到最后一列的权数最大，即人均创造总收入这个指标在这项评价上的分辨信息丰富，这个指标的数值能明确区分开各个评价被评价对象差异。同理，第四列的权数最小，也就是说各个被评价对象在某项指标上的数值差异较小，那么这项指标区分开各评价对象的能力较弱。 2 用相关系数法求各个指标的权数： 基本步骤：（1）计算各个指标之间的相关系数矩阵 （2）构造分块矩阵 R1(去掉相关系数矩阵的第一行和第一列)R2 R3 R4 R5 R6 同理可得

数电实验二：简易计算器(实验报告)

数电实验2实验报告 1、设计修改方案 （1）加入编码器连接4选一数据选择器，控制进行运算的种类 （2）修改了输出端数据选择器的程序，使得当计算器没有任何输入时，结果显示保持为0，并且利用芯片自身的灭零管脚，让显示结果中，当十位为零时，十 位的零不显示。

2、实验数据及分析 （1）修改后电路图（附后） （2）仿真波形 设置输入2个4位二进制数为0110（十进制6）和0010（十进制2），计算方式控制SW[3:0]设为0111，即模拟除法操作，加入时钟信号。 ①模拟除法波形： 可以看到十位（商）的数码管显示中，1、2、3、4、7段亮，显示为数字3，而个位（余数）显示1、2、3、4、5、6段亮，显示数字0，相当于计算出6除2商3余0。满足计算要求。 ②模拟乘法波形：（SW[3:0]设为1011，其他输入同上）

可以看到个位的数码管显示中，1、4、5、6段亮，显示为C(化为十进制为12)，而十位一直显示1、2、3、4、5、6段亮，显示数字0，相当于计算出6乘2等于0C，即等于12。当改变输入4和2是，显示结果为8,。满足计算要求。 ③模拟加法波形：（SW[3:0]设为1101，其他输入同上） 可以看到个位的数码管显示中，1、2、3、4、5、6、7段全亮，显示为数字8，而十位一直显示1、2、3、4、5、6段亮，显示数字0，相当于计算出6加2等于08，即等于8。满足计算要求。 ④模拟减法波形：（SW[3:0]设为1110，其他输入同上） 可以看到个位的数码管显示中，2、3、6、7段亮，显示为数字3，而十位一直显示1、2、3、4、5、6段亮，显示数字0，相当于计算出6减2等于03，即等于3。满足计算要求。 从上面加减乘除四种功能运算的波形仿真可以看出，本实验设计能够正确完成对输入数字的上述四种运算。满足题目要求。

合肥工业大学信息隐藏实验报告 完整代码版.

计算机与信息学院 信息隐藏实验报告 专业班级 信息安全13-1班 学生姓名及学号 马骏 2013211869 课程教学班号 任课教师 郑淑丽 实验指导教师 郑淑丽 实验地点 20 ~20 学年第学期

实验1 BMP位图信息隐藏 一、实验目的 学习BMP格式文件，并编程实现对位图文件信息隐藏 二、实验要求 将TXT文件嵌入BMP 文件中 三、问题描述 1、BMP位图文件的格式？ 2、有哪几种方法隐藏信息，分别采用什么样的数据结构 3、随机选取如何避免“碰撞”的出现 四、算法思想 1、BMP位图文件格式 0000h~0001h 2字节-------------------------bm的ASC码 0002h~0005h 4字节-------------------------文件大小102718字节 0006h~0009h 4字节-------------------------全为0 000Ah~000Dh 4字节-------------------------偏移量118字节 000Eh~0011h 4字节-------------------------位图信息块大小40字节 0012h~0015h 4字节-------------------------宽450 0016h~0019h 4字节-------------------------高450 001Ah~001Bh 2字节-------------------------恒为01h 00h 001Ch~001Dh 2字节-------------------------颜色所占二进制位数值04h 00h=4 16色位图 001Eh~0021h 4字节-------------------------压缩方式=0无压缩 0022h~0025h 4字节-------------------------图像数据区大小102600字节 0026h~0029h 4字节-------------------------水平每米多少像素39个 002Ah~002Dh 4字节-------------------------垂直每米多少像素39个 002Eh~0031h 4字节-------------------------图像所用颜色数=0 0032h~0035h 4字节-------------------------重要颜色数=0 0036h~0076h 64字节-------------------------颜色表

计算方法实验报告格式

计算方法实验报告格式 小组名称: 组长姓名(班号)： 小组成员姓名(班号): 按贡献排序情况: 指导教师评语: 小组所得分数: 一个完整的实验，应包括数据准备、理论基础、实验内容及方法，最终对实验结果进行分析，以达到对理论知识的感性认识，进一步加深对相关算法的理解，数值实验以实验报告形式完成，实验报告格式如下： 一、实验名称 实验者可根据报告形式需要适当写出． 二、实验目的及要求 首先要求做实验者明确，为什么要做某个实验，实验目的是什么，做完该实验应达到什么结果，在实验过程中的注意事项，实验方法对结果的影响也可以以实验目的的形式列出． 三、算法描述（实验原理与基础理论） 数值实验本身就是为了加深对基础理论及方法的理解而设置的，所以要求将实验涉及到的理论基础，算法原理详尽列出． 四、实验内容 实验内容主要包括实验的实施方案、步骤、实验数据准备、实验的算法以及可能用到的仪器设备． 五、程序流程图 画出程序实现过程的流程图，以便更好的对程序执行的过程有清楚的认识，在程序调试过程中更容易发现问题． 六、实验结果 实验结果应包括实验的原始数据、中间结果及实验的最终结果，复杂的结果可以用表格

形式列出，较为简单的结果可以与实验结果分析合并出现． 七、实验结果分析 实验结果分析包括对对算法的理解与分析、改进与建议． 数值实验报告范例 为了更好地做好数值实验并写出规范的数值实验报告，下面给出一简单范例供读者参考． 数值实验报告 小组名称: 小组成员(班号): 按贡献排序情况: 指导教师评语: 小组所得分数: 一、实验名称 误差传播与算法稳定性． 二、实验目的 1．理解数值计算稳定性的概念． 2．了解数值计算方法的必要性． 3．体会数值计算的收敛性与收敛速度． 三、实验内容 计算dx x x I n n ? += 1 10 ，1,2,,10n = ． 四、算法描述 由 dx x x I n n ? += 1 10 ，知 dx x x I n n ?+=--101110，则

计算方法实验报告

计算方法实验报告（四） 方程和方程组的迭代解法 一、实验问题 利用简单迭代法，两种加速技术，牛顿法，改进牛顿法，弦割法求解习题5-1，5-2，5-3中的一题，并尽可能准确。 选取5-3：求在x=1.5附近的根。 二、问题的分析（描述算法的步骤等） （1）简单迭代法算法： 给定初始近似值,求的解。 Step 1 令i=0； Step 2 令(计算)； Step 3 如果,则迭代终止，否则重复Step 2。 （2）Aitken加速法算法 Step 1 令k=0，利用简单迭代算法得到迭代序列； Step 2 令-(计算得到一个新的序列，其中k=0,1,2…)；Step 3 如果,则迭代终止，否则重复Step 2。 (3)插值加速法算法 Step 1 令k=0，利用简单迭代算法得到迭代序列； Step 2 令+(计算得到一个新的序列，其中k=1,2,3…)； Step 3 如果,则迭代终止，否则重复Step 2。 (4)牛顿法算法

Step 1给定初始近似值； Step 2令,其中k计算得到的序列； Step 3如果,则迭代终止，否则重复Step 2。 （5）改进牛顿法的算法 Step 1给定初始近似值； Step 2令,其中k迭代计算得到的序列； Step 3如果,则迭代终止，否则重复Step 2。 （6）弦割法算法（双点弦割法） Step 1给定初始近似值，； Step 2令其中k计算得到的序列； Step 3如果,则迭代终止，否则重复Step 2。 三、程序设计 (1)简单迭代法 利用迭代公式进行迭代运算。 #include #include #include double fun(double x) { double c=1+x*x; return pow(c,1/3.0); } void main() { double x=1.5; double y=0; double D=1;

单片机实验报告 计算器

单片机原理及其应用实验报告基于51单片机的简易计算器的设计 班级：12电子1班 姓名：金腾达 学号：1200401123 2015年1月6日

摘要 一个学期的51单片机的课程已经随着期末的到来落下了帷幕。“学以致用”不仅仅是一句口号更应该是践行。本设计秉承精简实用的原则，采用AT89C51单片机为控制核心，4X4矩阵键盘作为输入，LCD1602液晶作为输出组成实现了基于51单片机的简易计算器。计算器操作方式尽量模拟现实计算器的操作方式，带有基本的运算功能和连续运算能力。并提供了良好的显示方式，与传统的计算器相比，它能够实时显示当前运算过程和上一次的结果，更加方便用户记忆使用。本系统制作简单，经测试能达到题目要求。 关键词：简易计算器、单片机、AT89C51、LCD1602、矩阵键盘

目录 一、系统模块设计......................................................................................... 错误！未定义书签。 1.1 单片机最小系统 (1) 1.2 LCD1602液晶显示模块 (1) 1.3 矩阵按键模块 (2) 1.4 串口连接模块 (1) 二、C51程序设计 (2) 2.1 程序功能描述及设计思路 (2) 2.1.1按键服务函数 (2) 2.1.2 LCD驱动函数 (2) 2.1.3 结果显示函数 (2) 2.1.4状态机控制函数 (2) 2.1.5串口服务函数 (2) 2.2 程序流程图 (3) 2.2.1系统总框图 (3) 2.2.2计算器状态机流程转换图 (3) 三、测试方案与测试结果 (4) 3.1测试方案 (4) 3.3 测试结果及分析 (7) 4.3.1测试结果（仿真截图） (7) 4.3.2测试分析与结论 (7) 四、总结心得 (7) 五、思考题 (8) 附录1：整体电路原理图 (9) 附录2：部分程序源代码 (10)

c计算器实验报告

c计算器实验报告集团标准化工作小组 #Q8QGGQT-GX8G08Q8-GNQGJ8-MHHGN#

简单计算器 姓名: 周吉祥 实验目的：模仿日常生活中所用的计算器，自行设计一个简单的计算器程序，实现简单的计算功能。 实验内容： (1)体系设计： 程序是一个简单的计算器，能正确输入数据，能实现加、减、 乘、除等算术运算，运算结果能正确显示，可以清楚数据等。 (2)设计思路： 1)先在Visual C++ 中建立一个MFC工程文件，名为calculator. 2)在对话框中添加适当的编辑框、按钮、静态文件、复选框和 单选框 3)设计按钮，并修改其相应的ID与Caption. 4)选择和设置各控件的单击鼠标事件。 5)为编辑框添加double类型的关联变量m_edit1. 6)在中添加头文件，然后添加public成员。 7)打开文件，在构造函数中，进行成员初始化和完善各控件的 响应函数代码。 (3)程序清单： 添加的public成员： double tempvalue; 法 2.减法 3.乘法 4.除法

int append; //判断后面是否添加数字 ●成员初始化： CCalculatorDlg::CCalculatorDlg(CWnd* pParent /*=NULL*/) : CDialog(CCalculatorDlg::IDD, pParent) { //{{AFX_DATA_INIT(CCalculatorDlg) m_edit1 = ; //}}AFX_DATA_INIT // Note that LoadIcon does not require a subsequent DestroyIcon in Win32 m_hIcon = AfxGetApp()->LoadIcon(IDR_MAINFRAME); tempvalue=0; result=0; sort=0; append=0; } ●各控件响应函数代码： void CCalculatorDlg::OnButton1() //按钮“1” { // TODO: Add your control notification handler code here if(append==1)result=0;

计算方法实验报告 拟合

南京信息工程大学实验（实习）报告 一、实验目的: 用最小二乘法将给定的十个点拟合成三次多项式。 二、实验步骤： 用matlab编制以函数为基的多项式最小二乘拟合程序，并用于对下列数据作三次多项式最小二乘拟合(取权函数wi=1) x -2.5 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2 y -2.30 -1 -0.14 -0.25 0.61 1.03 1.75 2.75 4.42 6.94 给定直线方程为：y=1/4*x3+1/2*x2+x+1 三、实验结论： 最小二乘法：通过最小化误差的平方和寻找数据的最佳函数匹配。利用最小二乘法可以简便地求得未知的数据，并使得这些求得的数据与实际数据之间误差的平方和为最小。最小二乘法还可用于曲线拟合。 一般地。当测量数据的散布图无明显的规律时，习惯上取n次代数多项式。 程序运行结果为： a = 0.9731 1.1023 0.4862 0.2238 即拟合的三次方程为：y=0.9731+1.1023x+0.4862*x2+0.2238*x3

-2.5 -2-1.5-1-0.5 00.51 1.52 2.5 -4-20246 81012 x 轴 y 轴 拟合图 离散点 y=a(1)+a(2)*x+a(3)*x.2+a(4)*x.3 结论： 一般情况下，拟合函数使得所有的残差为零是不可能的。由图形可以看出最小二乘解决了残差的正负相互抵消的问题，使得拟合函数更加密合实验数据。 优点：曲线拟合是使拟合函数和一系列的离散点与观测值的偏差平方和达到最小。 缺点：由于计算方法简单，若要保证数据的精确度，需要大量的数据代入计算。

实验报告的书写格式模版

实验报告的书写格式模版 有关实验报告的书写格式 一、完整实验报告的书写 完整的一份实验报告一般包括以下项目：实验名称： 实验目的： 实验器材： 实验原理： 实验步骤： 实验数据记录（表格）及处理： 实验结论（结果推导）： 实验讨论或分析等。 二、实验报告书写方法  1、实验名称：就是这个实验是做什么的。 2、实验目的：一般都写掌握什么方法啊；了解什么啊；知道什么啊；会什么啊；…… 等。 3、实验器材：就是做这个实验需要的所有器材（仪器）。 4、实验原理：就是这个实验是根据什么来做的，一般书上会写，抄一下也就可以啦。 5、实验步骤：就是你做实验的过程，开始操作时，（1）做什么； （2）做什么；（3）做什么；……

6、实验数据记录（表格）及处理：根据实验中涉及以及实验得到的数据，设计表格，将有关数据填在表格相应的位置；数据处理，就是该计算的，按要求计算后填入表格对应位置。 7、实验结论（结果推导）：就是做这个实验要得到的结果。 8、分析于讨论：写你的实验结果是否适合真实值？如果有误差要分析产生误差的原因，还有实验的一些比较关键的步骤的注意事项等。 对于初中生或小学生来说，书写的实验报告也可简单一点，有时也可不要分析于讨论，也可不写实验原理等。 三、探究实验书写一般有七个环节 1．提出问题：就是在生活中发现、提出问题。 2．猜想与假设：发现问题，就要弄清楚问题，在没有搞清楚之前总有基本的猜测和设想，这就是猜想与假设。 3．制定计划与设计实验：有了猜想，就有了实验的目的，再根据实验的目的设计实验方案，制定实验计划，包括取得证据的途径和方法，确定收集证据的范围。包括实验的理论依据（实验原理）、实验器材、实验步骤等。 4．进行实验与收集证据：上一步是动脑、思维活动，这一步是手脑并用的实验过程。 5．分析与论证：通过上面的实验，收集到一些数据，观察到一些现象，对其分析，得出事实与假设的关系，通过归纳、概括等方法，得到结论。

简易计算器设计实验报告

简易计算器设计实验报告 一．设计任务及要求 1.1实验任务： 根据计算器的原理设计一个具有加减乘除功能的简易计算器。如：5+3*4/8=4。 1.2 实验基本要求： （1）实现最大输入两位十进制数字的四则运算（加减乘除）。 （2）能够实现多次连算（无优先级，从左到右计算结果）。 如：12+34*56-78/90+9=36 （3）最大长度以数码管最大个数为限，溢出报警。 二．实验设计方案 （1）用QuartusII的原理图输入来完成系统的顶层设计。 （2）用VHDL编写以及直接拖模块来各功能模块。 （3）通过2个脉冲分别实现个位数和十位数的输入。 （4）通过选择每次的输出数值，将输出值反馈到运算输入端 （4）通过除法运算实现十六进制到十进制的转换输出。 其具体实现流程图如下：

三系统硬件设计 FPGA： EP2C5T144C8目标板及相应外围硬件电路。（从略） 四系统软件设计 1.数据输入模块 原理：用VHDL创建模块，通过两个脉冲分别对两个数码管进行输入控制，再通过相应运算模块将两个独立数据转化成两位十进制数字。 2.运算模块 原理：用VHDL创建模块，四种运算同步运行，通过按键加、减、乘、除选择输出对应的计算结果，当按键等号来时，将所得结果反馈给运算模块输入端。具体实现代码见附录二。 3.输出模块 原理：用VHDL创建模块，通过按键等号来控制显示运算对象还是运算结果，当等号按下时，输出计算结果，否则显示当前输入的数据，并且通过除法模块将十六进制转化为十进制。当输出结果溢出是LED0亮，同时数码管显示都为零。部分实现见附录二。 五实验调试 输入数据12，再按加法键，输入第二个数字25，按等号键，数码管显示37；按灭加法、等号键，输入第二个数据2，依次按等号键，减法键，数码管显示35；同上，按灭减法键、等号键，输入第三个数据7，依次按等号键，除法键，数码管显示5；按灭除法键、等号键，输入第四个数据99，依次按等号键，乘法键，数码管显示495，按灭乘法键、等号键，当前显示为99，依次按等号键、乘法键，数码管显示49005，同上进行若干次之后，结果溢出，LED0亮，同时数码管显示都为零。当输出为负数时，LED0灯变亮，同时数码管显示都为零。六实验结论 本实验基本实现了计算器的加减乘法运算功能，但是存在一个突出的缺陷，就是当输出结果时，必须先按等号键导通数据反馈，再按运算键选择输出结果。这与实际应用的计算器存在很大的差距。但是，本设计可以通过等号键实现运算对象和运算结果之间的切换。

合工大汇编语言程序设计实验报告

合肥工业大学计算机与信息学院 实验报告 课程：汇编语言程序设计专业班级：**************** 学号：********** 姓名：***** 目录

实验一 (3) 实验二 (7) 实验三 (12) 实验四 (22)

实验一Debug程序的使用 一.实验目的 1、熟悉DEBUG程序中的命令，学会在DEBUG下调试运行汇编语言源程序。 2、掌握8086/8088的寻址方式及多字节数据的处理方法。 二.实验内容 1、利用DEBUG程序中的“E”命令，将两个多字节数“003F1AE7H”和“006BE5C4H”分别送入起始地址为DS:0200H和DS:0204H两个单元中。 2、分别用直接寻址方式和寄存器间接寻址方式编写程序段，实现将DS:0200H 单元和DS:0204H单元中的数据相加，并将运算结果存放在DS:0208H单元中。要求： 本次实验的内容均在DEBUG下完成，实现数据的装入、修改、显示；汇编语言程序段的编辑、汇编和反汇编；程序的运行和结果检查。 三.实验过程和程序 实验内容一： e ds:0200 E7 1A 3F 00 e ds:0204 C4 E5 6B 00 实验内容二： （1）直接寻址方式 MOV AX,[0200] MOV BX,[0202] ADD AX,[0204] ADC BX,[0206] MOV [0208],AX MOV [020A],BX (2)寄存器间接寻址方式 MOV SI,0200H MOV DI,0204H MOV BX,0208H MOV AX,[SI] MOV DX,[SI+2] ADD AX,[DI] ADC DX,[DI+2]

计算方法实验报告

实验报告 一、求方程f(x)=x^3-sinx-12x+1的全部根, ε=1e -6 1、 用一般迭代法; 2、 用牛顿迭代法; 并比较两种迭代的收敛速度。 一、首先，由题可求得：12cos 3)(2 ' --=x x x f . 其次，分析得到其根所在的区间。 ① 令()0=x f ，可得到x x x sin 1123 =+-. ② 用一阶导数分析得到1123 +-x x 和x sin 两个函数的增减区间；再用二阶导数分析得到 两个函数的拐点以及凹凸区间. ③ 在直角坐标轴上描摹出01123 =+-x x 和0sin =x 的图，在图上可以看到他们的交点，然后估计交点所在的区间，即是所要求的根的区间。经过估计，得到根所在的区间为 []3,4--，[]1,0和[]4,3. 1、 一般迭代法 （1）算法步骤： 设ε为给定的允许精度，迭代法的计算步骤为： ① 选定初值0x .由()0=x f 确定函数()x g ，得等价形式()x g x =. ② 计算()0x g .由迭代公式得()01x g x =. ③ 如果ε≤-01x x ，则迭代结束，取1x 为解的近似值；否则，用1x 代替0x ，重复步骤②和步骤③. （2）程序代码： ① 在区间[]3,4--内， 代码： clc

x0=-3.5; %初值0x iter_max=100; %迭代的最大次数 ep=1e-6; %允许精度 ε k=0; while k<=iter_max %k 从0开始到iter_max 循环 x1=(sin(x0)+12*x0-1).^(1/3); %代入0x ，算出1x 的值 if abs(x1-x0)

实验报告 编写一个简易实用的计算器。

开课实验室：现代信息交流中心40*开课时间:2012年**月**日实验报告:年月日 一、实验目的：熟悉Java图形用户界面的开发原理与技术，包括各种组件和容器技术。 二、实验内容： 编写一个简易实用的计算器。 三、源代码清单： package exp3; import Result extends Panel { TextField tf; Result() { tf = new TextField("0", 25); add(tf); } } class Key extends Panel { Button backSpace, clearEnter, cancel,lBracket,rBracket; Button one, two, three, four, five, six, seven, eight, nine, zero; Button add, subtract, multiplay, divide; Button sqrt, reciprocal, mod; Button positiveAndNegative, point, equal; Key() { setLayout(new GridLayout(5, 5)); backSpace = new Button("←"); (new Color(255, 0, 0)); add(backSpace); clearEnter = new Button("CE"); (new Color(255, 0, 0)); add(clearEnter); cancel = new Button("C"); (new Color(255, 0, 0));

合工大电路实验报告2016.

电路分析基础 实验指导 2016．3

实验一常用电子仪器使用 1.1万用表 万用表是一种多用途的电工仪表，最常用的万用表，具有测量直流电压、交流电压、直流电流、交流电流、电阻值等功能。万用表可分为指针式和数字式万用表，图1为指针式万用表，图2为数字式万用表，万用表的型号很多，但它们的结构基本相似，使用方法也基本相同。 一、使用方法 1、电阻的测量方法 步骤一：连接表棒。将表棒插入万用表的正负接线柱上； 步骤二：调零。将表棒的两端相接触，调节欧姆条正电位器，使指针准确地指在欧姆刻度的零位上； 步骤三：档位选择。将转换开关旋至欧姆档的范围内；开始尽量选大的范围，测量后根据阻值再进行选择适当的范围。 步骤四：测量电阻。将表棒分开去测量未知电阻的阻值。 2、直流、交流电压和直流电流的测量方法 步骤一：连接表棒。 步骤二：档位选择。如果测直流电压将转换开关旋至直流电压档的范围内，测交流电压就旋至交流电压档范围内，测直流电流就将旋至直流电流档的范围内； 步骤三：量程选择。估算后选择档位的适当量程； 步骤四：测量并读数。将表棒接入电路，根据指针偏转，读出待测电压的大小。 二、注意事项 1、如果用数字表测量时，被测的量小，档位放的太大，就会出现小数点的位数太多。如果被测量的量大，档位放的太小，这时测量结果数字就为“1”，这说明档位太小，被测的量太大，数值溢出了，需换大档位。 2、测量电阻时，若用指针表，首先要将两只表笔短路，用调零旋钮将表针调到零，然后再测量。测量时，两手不应同时接触电阻两端，否则相当于在被测电阻两端并联一个人体电阻，而产生误差。 3、测量出的电阻值是档位值乘上指针的读数。如果用数字表测量电阻时，电阻值可以直接读出。 4、每次使用前应将转换开关调节在正确位置上再开始测量。测量电路中的电压时，万用表要并联在被测支路上。测量电路中的电流时，万用表要串联在被测支路中。 5、应该养成良好的使用习惯，即每当万用表使用完毕，应将转换开关放在最高交流电压档位上。

计算方法实验报告 插值

实验名称：插值计算 1引言 在生产和科研中出现的函数是多种多样的。常常会遇到这样的情况：在某个实际问题中，虽然可以断定所考虑的函数f（x）在区间[a,b]上存在且连续，但却难以找到它的解析表达式，只能通过实验和观测得到在有限个点上的函数值。用这张函数表来直接求出其他点的函数值是非常困难的，在有些情况下，虽然可以写出f(x)的解析表达式，但由于结构十分复杂，使用起来很不方便。面对这些情况，构造函数P(x)作为f(x)的近似，插值法是解决此类问题比较古老却目前常用的方法，不仅直接广泛地应用与生产实际和科学研究中，而且是进一步学习数值计算方法的基础。 设函数y=f(x)在区间[a,b]上连续，且在n+1个不同的点a≤x0,x1……,xn≤b上分别取值y0,y1……,yn. 插值的目的就是要在一个性质优良、便于计算的函数φ中，求一简单函数P(x),使P(xi)=yi(i=0,1…,n)而在其他点x≠xi上，作为f(x)的近似。 通常，称区间[a,b]为插值区间，称点x0,x1,…,xn为插值节点，上式为插值条件，称函数类φ为插值函数类，称P(x)为函数f(x)在节点x0,x1,…,xn处的插值函数，求插值函数P(x)的方法称为插值法。 2实验目的和要求 用matlab定义分段线性插值函数、分段二次插值函数、拉格朗日插值函数，输入所给函 数表，并利用计算机选择在插值计算中所需的节点，计算f(0.15),f(0.31),f(0.47)的近似值。

3算法描述 1.分段线性插值流程图

2.分段二次插值流程图

3.拉格朗日插值流程图

4程序代码及注释 1.分段线性插值

计算方法上机实验报告——拉格朗日插值问题

计算方法上机实验报告——拉格朗日插值问题 一、方法原理 n次拉格朗日插值多项式为：Ln(x)=y0l0(x)+y1l1(x)+y2l2(x)+…+ynln(x) n=1时，称为线性插值，L1(x)=y0(x-x1)/(x0-x1)+y1(x-x0)/(x1-x0)=y0+(y1-x0)(x-x0)/(x1-x0) n=2时，称为二次插值或抛物线插值，精度相对高些 L2(x)=y0(x-x1)(x-x2)/(x0-x1)/(x0-x2)+y1(x-x0)(x-x2)/(x1-x0)/(x1-x 2)+y2(x-x0)(x-x1)/(x2-x0)/(x2-x1) 二、主要思路 使用线性方程组求系数构造插值公式相对复杂，可改用构造方法来插值。 对节点xi(i=0,1,…,n)中任一点xk(0<=k<=n)作一n次多项式lk(xk)，使它在该点上取值为1，而在其余点xi(i=0,1,…,k-1,k+1,…,n)上为0，则插值多项式为Ln(x)=y0l0(x)+y1l1(x)+y2l2(x)+…+ynln(x) 上式表明：n个点xi(i=0,1,…,k-1,k+1,…,n)都是lk(x)的零点。可求得lk 三．计算方法及过程：1.输入节点的个数n 2.输入各个节点的横纵坐标 3.输入插值点 4.调用函数，返回z 函数语句与形参说明 程序源代码如下： 形参与函数类型 参数意义 intn 节点的个数 doublex[n]（double*x） 存放n个节点的值 doubley[n]（double*y） 存放n个节点相对应的函数值 doublep 指定插值点的值 doublefun() 函数返回一个双精度实型函数值，即插值点p处的近似函数值 #include #include usingnamespacestd; #defineN100 doublefun(double*x,double*y,intn,doublep); voidmain() {inti,n; cout<<"输入节点的个数n："; cin>>n;