线胀系数实验报告模板

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

《计算方法》实验报告 姓名： 班级： 学号： 实验日期: 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;

固体热膨胀系数的测量实验报告图文稿

固体热膨胀系数的测量 实验报告 集团文件发布号：（9816-UATWW-MWUB-WUNN-INNUL-DQQTY-

固体热膨胀系数的测量 班级： 姓名： 学号： 实验日期： 一、实验目的 测定金属棒的线胀系数，并学习一种测量微小长度的方法。 二、仪器及用具 热膨胀系数测定仪(尺读望远镜、米尺、固体线膨胀系数测定仪、铜棒、光杠杆、温度计等) 三、实验原理 1．材料的热膨胀系数 线膨胀是材料在受热膨胀时，在一维方向上的伸长。在一定的温度范围内，固体受热后，其长度都会增加，设物体原长为L ，由初温t1加热至末温t2，物体伸长了 △L,则有 ()12t t L L -=?α （1） （2） 此式表明，物体受热后其伸长量与温度的增加量成正比，和原长也成正比。比例系数称为固体的线胀系数。一般情况下，固体的体胀系数为其线胀系数的3倍。 2．线胀系数的测量 在式（1）中△L 是个极小的量，这样微小的长度变化，普通米尺、游标卡尺的精度是不够的，可采用千分尺、读数显微镜、光杠杆放大法、光学干涉法等。考虑到测量方便和测量精度，我们采用光杠杆法测量。光杠杆系统是由平面镜及底座，望远镜和米尺组成的。光杠杆放大原理如下图所示： 当金属杆伸长△L 时，从望远镜中叉丝所对标尺刻度前后为b1、b2，这时()12t t L L -?= α

有： 带入（2）式得固体线膨胀系数为： 四、实验步骤及操作 1.单击登陆进入实验大厅 2.选择热力学试验单击 3.双击固体热膨胀系数的测量进入实验界面 4.在实验界面单击右键选择“开始实验” 5.调节平面镜至竖直状态 6.进行望远镜调节，调节方位、聚焦、目镜是的标尺刻线清晰，调节 中丝读数为0.0mm,并打开望远镜视野 7.单击铜棒测量长度，单击温度计显示铜棒温度，打开电源加热，记 录每升高10度时标尺读数直至温度升高到90度止 8.单击卷尺，分别测量l、D， 9．以t为横轴，b为纵轴作b-t关系曲线，求直线斜率。 10.代入公式计算线膨胀系数值。 由图得k=0.3724 五、实验数据记录与处理 六、思考题 1．对于一种材料来说，线胀系数是否一定是一个常数为什么 答：不是。因为同一材料在不同的温度区域，其线性系数是不同的，有实验结果的事实可证明。 2．你还能想出一种测微小长度的方法，从而测出线胀系数吗？ 答：目前想不到更好地方法。 3. 引起测量误差的主要因素是什么？ 答：仪器的精准度，操作过程中的不可避免性的失误，温度变化的控制，铜棒受热不均匀等。

计算方法上机实验报告

. / 《计算方法》上机实验报告 班级：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 四、运行结果：

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

本科实验报告 课程名称：计算机数值方法 实验项目：方程求根、线性方程组的直接解 法、线性方程组的迭代解法、代数插值和最 小二乘拟合多项式 实验地点：行勉楼 专业班级： ******** 学号： ********* 学生姓名： ******** 指导教师：李誌，崔冬华 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．材料的热膨胀系数 线膨胀是材料在受热膨胀时，在一维方向上的伸长。在一定的温度范围内，固体受 热后，其长度都会增加，设物体原长为L，由初温t1加热至末温t2，物体伸长了 △L,则有 () 1 2 t t L L- = ?α（1）（2） 此式表明，物体受热后其伸长量与温度的增加量成正比，和原长也成正比。比例系 数称为固体的线胀系数。一般情况下，固体的体胀系数为其线胀系数的3倍。 2．线胀系数的测量 在式（1）中△L是个极小的量，这样微小的长度变化，普通米尺、游标卡尺的精度是不够的，可采用千分尺、读数显微镜、光杠杆放大法、光学干涉法等。考虑到测 量方便和测量精度，我们采用光杠杆法测量。光杠杆系统是由平面镜及底座，望远 镜和米尺组成的。光杠杆放大原理如下图所示： () 1 2 t t L L - ? = α

当金属杆伸长△L时，从望远镜中叉丝所对标尺刻度前后为b1、b2，这时有：带入（2）式得固体线膨胀系数为： 四、实验步骤及操作 1.单击登陆进入实验大厅 2.选择热力学试验单击 3.双击固体热膨胀系数的测量进入实验界面 4.在实验界面单击右键选择“开始实验” 5.调节平面镜至竖直状态 6.进行望远镜调节，调节方位、聚焦、目镜是的标尺刻线清晰，调节中丝读 数为0.0mm,并打开望远镜视野 7.单击铜棒测量长度，单击温度计显示铜棒温度，打开电源加热，记录每升 高10度时标尺读数直至温度升高到90度止 l L D b b? = - 2 1 2 () D l b b L 2 1 2 - = ? () ()k DL l t t DL b b l 2 2 1 2 1 2= - - = α

数值计算实验报告

(此文档为word格式，下载后您可任意编辑修改！) 2012级6班###(学号)计算机数值方法 实验报告成绩册 姓名:宋元台 学号： 成绩：

数值计算方法与算法实验报告 学期： 2014 至 2015 第 1 学期 2014年 12月1日课程名称: 数值计算方法与算法专业:信息与计算科学班级 12级5班 实验编号： 1实验项目Neton插值多项式指导教师：孙峪怀 姓名：宋元台学号：实验成绩： 一、实验目的及要求 实验目的： 掌握Newton插值多项式的算法，理解Newton插值多项式构造过程中基函数的继承特点，掌握差商表的计算特点。 实验要求： 1. 给出Newton插值算法 2. 用C语言实现算法 二、实验内容 三、实验步骤(该部分不够填写.请填写附页)

1.算法分析： 下面用伪码描述Newton插值多项式的算法： Step1 输入插值节点数n，插值点序列{x(i),f(i)},i=1,2,……,n,要计算的插值点x. Step2 形成差商表 for i=0 to n for j=n to i f(j)=((f(j)-f(j-1)(x(j)-x(j-1-i)); Step3 置初始值temp=1,newton=f(0) Step4 for i=1 to n temp=(x-x(i-1))*temp*由temp(k)=(x-x(k-1))*temp(k-1)形成 (x-x(0).....(x-x(i-1)* Newton=newton+temp*f(i); Step5 输出f（x）的近似数值newton(x)=newton. 2.用C语言实现算法的程序代码 #includeMAX_N) { printf("the input n is larger than MAX_N,please redefine the MAX_N.\n"); return 1; } if(n<=0) { printf("please input a number between 1 and %d.\n",MAX_N); return 1; } printf("now input the (x_i,y_i)i=0,...%d\n",n); for(i=0;i<=n;i++) { printf("please input x(%d) y(%d)\n",i,i);

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; result=0; sort=0; append=0; }

金属线胀系数的测定实验报告

实验5 金属线胀系数的测定 测量固体的线胀系数，实验上归结为测量在某一问题范围内固体的相对伸长量。此相对伸长量的测量与杨氏弹性模量的测定一样，有光杠杆、测微螺旋和千分表等方法。而加热固体办法，也有通入蒸气法和电热法。一般认为，用电热丝同电加热，用千分表测量相对伸长量，是比较经济又准确可靠的方法。 一、实验目的 1.学会用千分表法测量金属杆长度的微小变化。 2.测量金属杆的线膨胀系数。 二、实验原理 一般固体的体积或长度，随温度的升高而膨胀，这就是固体的热膨胀。设物体的温度改变t ?时，其长度改变量为L ?,如果t ?足够小，则t ?与L ?成正比，并且也与物体原长L 成正比，因此有 t L L ?=?α （1） 式（1）中比例系数α称为固体的线膨胀系数，其物理意义是温度每升高1℃时物体的伸长量与它在0℃时长度之比。设在温度为0℃时，固体的长度为0L ,当温度升高为t ℃时，其长度为t L ，则有 t L L L t α=-00/)( 即 )1(0t L L t α+= （2） 如果金属杆在温度为1t ，2t 时，其长度分别为1L ，2L ，则可写出 )1(101t L L α+= （3） )1(202t L L α+= (4) 将式（3）代入式（4），又因1L 与2L 非常接近，所以，1/12=L L ,于是可得到如下

结果： )(12112t t L L L --=α （5） 由式（5），测得1L ，2L ，1t 和2t ，就可求得α值。 三、仪器介绍 （一）加热箱的结构和使用要求 1.结构如图5-1。 2.使用要求 （1）被测物体控制于mm 4008?φ尺寸； （2）整体要求平稳，因伸长量极小，故仪器不应有振动； （3）千分表安装须适当固定（以表头无转动为准）且与被测物体有良好的接触（读数在0.2~0.3mm 处较为适宜，然后再转动表壳校零）； （4）被测物体与千分表探头需保持在同一直线。 （二）恒温控制仪使用说明

固体线胀系数测定

SUES大学物理选择性实验讲义Typeset by L A T E X2ε 固体线胀系数测定? 一实验目的 本实验通过固体线胀系数测定仪测定不同金属的线胀系数，要求达到： 1．掌握使用千分表和温度控制仪的操作方法； 2．分析影响测量精度的诸因素； 3．观察合金材料在金相组织发生变化温度附近，出现线膨胀量的突变现象。二实验原理 绝大多数物质具有“热胀冷缩”的特性，这是由于物体内部分子热运动加剧或减弱造成的。这个性质在工程结构的设计中，在机械和仪表的制造中，在材料的加工(如焊接)中都应考虑到。否则，将影响结构的稳定性和仪表的精度，考虑失当，甚至会造成工程结构的毁损，仪表的失灵以及加工焊接中的缺陷和失败等等。 固体材料的线膨胀是材料受热膨胀时，在一维方向上的伸长。线胀系数是选用材料的一项重要指标，在研制新材料中，测量其线胀系数更是必不可少的。SLE-1固体线胀系数测定仪通过加热温度控制仪，精确地控制实验样品在一定的温度下，由千分表直接读出实验样品的伸长量，实现对固体线胀系数测定。 SLE-1固体线胀系数测定仪的恒温控制由高精度数字温度传感器与HTC-1加热温度控制仪组成，可加热温度控制在室温至80.0?C之间。HTC-1加热温度控制?修订于2009年2月4日 1

仪自动检测实测温度与目标温度的差距，确定加热策略，并以一定的加热输出电压维持实测温度的稳度，分别由四位数码管显示设定温度和实验样品实测温度,读数精度为±0.1?C。专用加热部件的加热电压为12V。 物质在一定温度范围内，原长为l的物体受热后伸长量?l与其温度的增加量?t近似成正比，与原长l也成正比，即：?l=α·l·?t。式中α为固体的线胀系数。实验证明：不同材料的线膨胀系数是不同的。本实验配备的实验样品为铁棒、铜棒、铝棒(加工成6×400mm的圆棒)。 三仪器技术指标 1、温度读数精度：±0.1?C。 2、温度控制稳定度：±0.1?C/10分钟。 3、温度设定范围：?5.0?C～+85?C，四位数码管显示。 4、实验样品实测温度：室温至82.0?C，四位数码管显示。 5、伸长量测量精度：0.001mm，量程：0～1mm。 6、HTC-1加热温度控制仪使用条件 1)输入电源：220V±10%50Hz～60Hz 2)湿度：<85% 3)温度：0～40?C 4)功耗：<70W 四仪器组成 由SLE-1固体线胀系数测定仪实验装置和HTC-1加热温度控制仪组成。 1、实验仪器如图1： 2、实验条件 2

计算方法实验报告

计算方法实验报告（四） 方程和方程组的迭代解法 一、实验问题 利用简单迭代法，两种加速技术，牛顿法，改进牛顿法，弦割法求解习题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;

计算方法实验报告格式

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

形式列出，较为简单的结果可以与实验结果分析合并出现． 七、实验结果分析 实验结果分析包括对对算法的理解与分析、改进与建议． 数值实验报告范例 为了更好地做好数值实验并写出规范的数值实验报告，下面给出一简单范例供读者参考． 数值实验报告 小组名称: 小组成员(班号): 按贡献排序情况: 指导教师评语: 小组所得分数: 一、实验名称 误差传播与算法稳定性． 二、实验目的 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，则

金属线胀系数

金属线胀系数的测定 实验目的：1）学会用千分表法测量金属杆长度的微小变化 2）测量金属杆的线胀系数，并判断此金属为何种金属 实验仪器： 实验原理：大家都知道热胀冷缩的现象，一般固体的长度或体积会随着温度的升高而膨胀，这就是固体的热膨胀。 设物体的温度改变Δt 时，其长度改变量为ΔL,如果Δt 足够小，则Δt 与ΔL 成正比，并且也与物体的原长有关系。因此它们三个量之间有： ΔL=αL Δt 式中的比例系数α称为固体的线胀系数，其物理意义是温度每升高1℃时其伸长量与它在0℃时长度的比。设金属在0℃时的长度是L0,当温度升高为t ℃时其长度为Lt,则有： （Lt-L0)/L0=αt 即Lt=L0(1+αt) 如果金属杆在温度为t1,t2时的长度分别为L1，L2，则可加热箱 恒温控制仪

以得到： L1=L0(1+αt1),L2=L0(1+αt2) 因为L1，L2非常接近，所以得到下式： α=（L2-L1)/L0(t2-t1) 由上式测得L1，L2，t1,t2就可以测得α值了。 实验过程： 1）接好电源和各个接口。 2）打开恒温控制仪，记录室温t1,再设定温度最大值，再记录此时千分表读数n1，最后按下确定键开始加热。 （实验所用金属杆0℃时长度为400mm) 3）每隔5℃读一次数tn ，同时记录千分表读数n n 。 4) 将数据整理填入设计好的表格中，待处理。 实验数据记录与处理： t1=21℃ L0=400mm n1=0.4012mm tn/℃ 26 31 36 41 46 51 tn-t1/℃ 5 10 15 20 25 30 n n /mm 0.4630 0.5119 0.553 0.591 0.624 0.658 n n -n1/mm 0.062 0.111 0.152 0.19 0.223 0.26

数学计算方法实验报告

数学计算方法实验报告 习题二 2.估计用二分法求方程f(x)=x3+4x2-10=0在区间[1,2]内根的近似值，为使方程不超过10时所需的二分次数。f(x k) 程序过程： function two (tolerance) a=1;b=2;counter=0; while (abs(b-a)>tolerance) c=(a+b)/2; fa=a^3+4*a^2-10;

fb=b^3+4*b^2-10; fc=c^3+4*c^2-10; if ((fa==0|fb==0)) disp(counter); elseif (fa*fc<0) b=c;counter=counter+1; elseif (fb*fc<0) a=c;counter=counter+1; elseif (fb==0) disp(counter); end end solution=(a+b)/2; disp(solution); disp(counter); 实验结果： 6.取x0=1.5,用牛顿迭代法求第三中的方程根.f(x)=x3+4x2-10=0的近似值（精确到||x k+1-x k|≦10-5,并将迭代次数与3题比较。 程序过程： function six (g) a=1.5; fa=a^3+4*a^2-10;

ga=3*a^2+8*a; b=a-fa/ga; k=1; while(abs(b-a)>g) a=b; fa=a^3+4*a^2-10; ga=3*a^2+8*a; b=a-fa/ga; k=k+1; end format long; disp(a); disp(k); 实验结果：程序结果计算结果 8.用弦割法求方程f(x)=x3-3x2-x+9=0在区间[-2,-1]内的一个实根近似值x k,|f(x k)|≦10-5. 程序过程： function eight (t) a=-2; b=-1; fa=a^3-3*a^2-a+9; fb=b^3-3*b^2-b+9; c=b-fb*(b-a)/(fb-fa); k=1; while(abs(c-b)>t) a=b; b=c; fa=a^3-3*a^2-a+9; fb=b^3-3*b^2-b+9; c=b-fb*(b-a)/(fb-fa); k=k+1; end

计算方法实验报告 拟合

南京信息工程大学实验（实习）报告 一、实验目的: 用最小二乘法将给定的十个点拟合成三次多项式。 二、实验步骤： 用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 结论： 一般情况下，拟合函数使得所有的残差为零是不可能的。由图形可以看出最小二乘解决了残差的正负相互抵消的问题，使得拟合函数更加密合实验数据。 优点：曲线拟合是使拟合函数和一系列的离散点与观测值的偏差平方和达到最小。 缺点：由于计算方法简单，若要保证数据的精确度，需要大量的数据代入计算。

计算方法实验报告

实验报告 一、求方程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)

线膨胀系数实验报告参考

线胀系数测量实验报告参考稿 【实验目的】 1．学习并掌握测量金属线膨胀系数的一种方法。 2．学会用千分表测量长度的微小增量。 【实验仪器】 FB712型金属线膨胀系数测量仪一台，千分表（1-0-0.001mm ）一个，待测铜管一根。 【实验原理】 材料的线膨胀是材料受热膨胀时，在一维方向的伸长。线胀系数是选用材料的一项重要指标。特别是研制新材料，少不了要对材料线胀系数做测定。 如图所示，待测铜管的线胀系数为： () t L L ???= α 式中L 为温度为1t 摄氏度时的管长，L ?为管受热后温度从1t 升高到2t 时的伸长量，t ?为管受热前后的温度升高量 (12t t t -=?) 。 该式所定义的线胀系数的物理意义是固体材料在()21t , t 温度区域内，温度每升高一度时材料的相对伸长量，其单位为()1 C -?。 【实验内容和步骤】 1．把样品铜管安装在测试架上。连接好加热皮管，打开电源开关，以便从仪器面板水位显示器上观察水位情况。水箱容积大约为ml 750。 3.加水步骤:先打开机箱顶部的加水口和后面的溢水管口塑料盖，用漏斗从加水口往系统内加水，管路中的气体将从溢水管口跑出，直到系统的水位计仅有上方一个红灯亮，其余都转变为绿灯时，可以先关闭溢水管口塑料盖。接着可以按下强制冷却按钮，让循环水泵试运行，由于系统内可能存在大量气泡，造成水位计显示虚假水位，只有利用循环水泵试运行过程，把系统内气体排出，这时候水位下降，仪器自动保护停机。 4．设置好温度控制器加热温度:金属管加热温度设定值可根据金属管所需要的实际温度值设置。 5．将铜管（或铝管）对应的测温传感器信号输出插座与测试仪的介质温度传感器插座相连接。将千分尺装在被测介质铜管（或铝管）的自由伸缩端固定位置上，使千分表测试端与被测介质接触，为了保证接触良好，一般可使千分表初读数为mm 2.0左右，只要把该数值作为初读数对待，不必调零。(如认为有必要，可以通过转动表面，把千分尺主指针读数基本调零，而副指针无调零装置。) 6．正常测量时，按下加热按钮（高速或低速均可，但低速档由于功率小，一般最多只能加热到C 50?左右)，观察被测金属管温度的变化，直至金属管温度等于所需温度值（例如C 35?）。.

数据分析方法与技术- 实验报告模板

《数据分析方法与技术》上机实验——实验1描述性统计方法 学号： 姓名： 日期：

实验项目（一）：描述性统计方法 一、实验内容 1．实验目的 掌握常用的描述性图表展示方法的原理及操作，包括：频数分布表、分组频数表、列联表、茎叶图、箱线图、误差图、散点图等； 掌握常用的描述性统计方法的原理及操作，包括：算术平均值、中位数、众数、四分位数、极差、平均差、方差、标准差、标准分数、离散系数等。 2. 实验内容和要求 实验内容：基于标准数据集，属性描述性图表展示方法（数分布表、分组频数表、列联表、茎叶图、箱线图、误差图、散点图等），对统计指标（算术平均值、中位数、众数、极差、平均差、方差、标准差、标准分数、离散系数、偏态峰态）进行计算。 实验要求：掌握各种描述性统计指标的计算思路及其在SPSS或EXCEL环境下的操作方法，掌握输出结果的解释。 二、实验过程 1、数据集介绍 1.数据库标题:鲍鱼数据 2.该数据库共计4177行数据 3.该数据有八个属性（包含性别共有九项） 4.以下是关于属性的描述，包括属性的名称，数据类型，测量单元和一个简短的描述： Name Data TypeMeas.Description ---- --------- ----- ----------- Sex nominal M, F, and I (infant)鲍鱼宝宝 Length continuousmm Longest shell measurement最长壳 Diameter continuousmm perpendicular to length垂直长度 Height continuousmm with meat in shell有肉的壳高度 Whole weightcontinuousgramswhole abalone整个鲍鱼 Shucked weightcontinuousgramsweight of meat肉的重量 Viscera weightcontinuousgramsgut weight (after bleeding)放血后内脏重 Shell weightcontinuousgramsafter being dried弄干后重量 Rings integer +1.5 gives the age in years +1.5=年龄 5.数据的值域

计算方法实验报告册

实验一——插值方法 实验学时：4 实验类型：设计 实验要求：必修 一 实验目的 通过本次上机实习，能够进一步加深对各种插值算法的理解；学会使用用三种类型的插值函数的数学模型、基本算法，结合相应软件（如VC/VB/Delphi/Matlab/JAVA/Turbo C ）编程实现数值方法的求解。并用该软件的绘图功能来显示插值函数，使其计算结果更加直观和形象化。 二 实验内容 通过程序求出插值函数的表达式是比较麻烦的，常用的方法是描出插值曲线上尽量密集的有限个采样点，并用这有限个采样点的连线，即折线，近似插值曲线。取点越密集，所得折线就越逼近理论上的插值曲线。本实验中将所取的点的横坐标存放于动态数组[]X n 中，通过插值方法计算得到的对应纵坐标存放 于动态数组[]Y n 中。 以Visual C++.Net 2005为例。 本实验将Lagrange 插值、Newton 插值和三次样条插值实现为一个C++类CInterpolation ，并在Button 单击事件中调用该类相应函数，得出插值结果并画出图像。CInterpolation 类为 class CInterpolation { public : CInterpolation();//构造函数 CInterpolation(float *x1, float *y1, int n1);//结点横坐标、纵坐标、下标上限 ~ CInterpolation();//析构函数 ………… ………… int n, N;//结点下标上限，采样点下标上限 float *x, *y, *X;//分别存放结点横坐标、结点纵坐标、采样点横坐标 float *p_H,*p_Alpha,*p_Beta,*p_a,*p_b,*p_c,*p_d,*p_m;//样条插值用到的公有指针，分别存放 i h ，i α，i β，i a ，i b ，i c ，i d 和i m }; 其中，有参数的构造函数为 CInterpolation(float *x1, float *y1, int n1) { //动态数组x1，y1中存放结点的横、纵坐标，n1是结点下标上限（即n1+1个结点） n=n1; N=x1[n]-x1[0]; X=new float [N+1]; x=new float [n+1]; y=new float [n+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;

金属线膨胀系数测量实验报告

梧州学院学生实验报告 成绩： 指导教师： 专业： 班别： 实验时间： 实验人： 学号： 同组实验人： 实验名称：金属线膨胀系数测量 实验目的：1、学习测量金属线膨胀系数的一种方法。 2、学会使用千分表。 实验仪器: 型号规格 单位 数量 备注 FB7 1 2型金属线膨 胀系数测定仪 台 1 被测件测试架 台 1 千分表 只 1 传感器连接线 根 2 L=80c m 红黑各一根 小漏斗 只 1 电源线 根 1 实验讲义(说明书)] 本 1 注意事项：1、做实验前必须精读FB712型金属线膨胀系数测定仪的使用说明书，正规操作 2 、注意千分表的使 用规范。 FB712型金属线膨胀系数测量仪实验装置示意图 【实验原理】 材料的线膨胀是材料受热膨胀时，在一维方向的伸长。线胀系数是选用材料的一项重要指标。 特别是研制新材料，少不了要对材料线胀系数做测定。 固体受热后其长度的增加称为线膨胀。经验表明，在一定的温度范围内，原长为 L 的物体, 受热后其伸长量厶L 与其温度的增加量△ t 近似成正比，与原长L 亦成正比，即： △ L=a ? L ?△ t (1) 式中的比例系数a 称为固体的线膨胀系数（简称线胀系数）。大量实验表明，不同材料的线胀系数 不同，塑料的 47 -J?V 叱-■: <■："負号 ■'a ^_A s'.Vi Pf jW 丹 >￥ -i~ ■ "I irtf I - *■ 4 !■":■_! 牡二盂:J 豪迂二辽山输咤或典: &：：?,、性％世*巴电冷忙即卜亠：.豆凳;其 応宓云I 恣心加[文 图&匹丁型金属线勝胀無数测定仪实物黑片 强制风冷 低速如撰 高速&]壇 盥控设齧 放水阀 H 水fr 匕 千分表 铝骨 FT1碱度传感黯 循环水管 削* 口 金廉管温度扬示 甥管 爲虔倩号践 S 度 指