铸造实验报告
一、铸造方法：湿型砂型铸造，手工造型。
二、造型材料的配方：由查《铸造技术数据手册》中的湿型砂配比表可得 1、背砂：新砂5%、旧砂94%、膨润土1%。
2、面砂：新砂40.4%、旧砂50%、膨润土4.5%、煤粉4%、重油1%、碳酸钠0.1%。
3、芯砂：新砂52%、旧砂40%、粘土6%、糖浆2%。
三、HT200化学成分的确定：由查表可得牌号
四、确定炉料配比
1、新生铁：
根据感应电炉熔炼铸铁的特性，为保证显微组织正常，炉料中生铁锭的用量不能超过20%。故选择新生铁的配比为20%，则新生铁的加入量：
150公斤x20%=30公斤
2、废钢
为了使炉料含碳量足够，废钢的配比为23%，则废钢的加入量为：
150公斤x23%=34.5公斤
3、回炉料
回炉料的加入量为：150公斤-30公斤-34.5公斤=85.5公斤五、计算炉料中各元素的应有含量 1、炉料应含碳量
铁水所需的平均含碳量（铁水C）应等于毛坯所需的含碳量（1C），即铁水C=1C，碳熔炼烧损为1%，则
炉料C=铁水C/（1-0.01）=1C/0.99
因为毛坯所需的含碳量1C已知为3.3%，所以
炉料C=铁水C/1-0.01=1C/0.99=3.3/0.99=3.33%
验算炉料实际含碳量：
新生铁带进的碳量：2C=4.0x20%=0.8% 回炉料带进的碳量：3C=4.15x57%=2.37% 废钢带进的碳量：4C=1.0x23%=0.23%
所以炉料实际含碳量C=2C+3C+4C=0.8%+2.37%+0.23%=3.4%
2、炉料应含硅量：
铁水所需的平均含硅量（铁水Si）应等于毛坯所需的含硅量（1Si），即铁水Si=1Si，硅熔炼烧损为6%，则
炉料Si=铁水Si/（1-0.06）=1Si/0.94
因为毛坯所需的含硅量1Si已知为1.95%，所以
炉料Si=铁水Si/（1-0.06）=1Si/0.94=1.95/0.94=2.07%
验算炉料实际含硅量：
新生铁带进的含硅量：2Si=1.85%x20%=0.37% 回炉料带进的含硅量：3Si=1.5%x57%=0.855% 废钢带进的含硅量：4Si=0.26%x23%=0.05%
所以炉料实际含硅量Si=2Si+3Si+4Si=0.37%+0.855%+0.06%=1.285% 炉料中尚缺硅量=炉料应含硅量-炉料中实际含硅量=2.07%-1.285%=0.785% 3、炉料应含锰量
铁水所需的平均含锰量（铁水Mn）应等于毛坯所需的含锰量（1Mn），即铁水Mn=1Mn，锰熔炼烧损为8%，则
因为毛坯所需的含锰量1Mn已知为0.8%，所以
炉料Mn=铁水Mn/(1-0.08)=0.8//0.92=0.87%
验算炉料实际含锰量：
新生铁带进的锰量：2Mn=0.03%x20%=0.006% 回炉料带进的锰量：3Mn=0.7%x57%=0.399% 废钢带进的锰量：4Mn=0.25%x23%=0.058%
所以炉料中实际含锰量Mn=2Mn+3Mn+4Mn=0.006%+0.399%+0.058%=0.463%
炉料中尚缺锰量=炉料应含锰量-炉料实际含锰量=0.87%-0.463%=0.407% 六、计算铁合金加入量 1、硅铁加入量
补加含硅75%的硅铁的百分数=炉料中尚缺的硅量（%）/硅铁的含硅量（%）=0.785/75%=1.05% 即每100公斤炉料补加1.05公斤含硅75%的硅铁则150公斤炉料补加含硅75%的硅铁量为： 1.5x1.05公斤=1.57公斤
2、锰铁加入量
补加含锰65%的锰铁的百分数=炉料中尚缺的锰铁（%）/锰铁的含量（%）=0.407%/65%=0.626% 即每100公斤炉料补加0.626公斤含锰65%锰铁则则150公斤炉料补加含锰65的锰铁量为： 1.5x0.626公斤=0.939公斤
《计算方法》实验报告 姓名： 班级： 学号： 实验日期: 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 四、运行结果： 实验二：
第1章铸造工艺设计基础 § 1-1零件结构的铸造工艺性分析 § 1-2铸造工艺方案的确定 § 1-3铸造工艺参数的确定 § 1-4砂芯设计 铸造生产周期较长，工艺复杂繁多。为了保证铸件质量，铸造工作者应根据铸件特点，技术条件和生产批量等制订正确的工艺方案，编制合理的铸造工艺流程，在确保铸件质量的 前提下，尽可能地降低生产成本和改善生产劳动条件。本章主要介绍铸造工艺设计的基础知 识，使学生掌握设计方法，学会查阅资料，培养分析问题和解决问题的能力。 § 1-1零件结构的铸造工艺性分析 铸造工艺性，是指零件结构既有利于铸造工艺过程的顺利进行，又有利于保证铸件质量。 还可定义为：铸造零件的结构除了应符合机器设备本身的使用性能和机械加工的要求外，还应符合铸造工艺的要求。这种对铸造工艺过程来说的铸件结构的合理性称为铸件的铸造工艺性。 另定义：铸造工艺性是指零件的结构应符合铸造生产的要求，易于保证铸件品质，简化 铸造工艺过程和降低成本。 铸造工艺性不好，不仅给铸造生产带来麻烦，不便于操作，还会造成铸件缺陷。因此，为了简化铸造工艺，确保铸件质量，要求铸件必须具有合理的结构。 一、铸件质量对铸件结构的要求 1 .铸件应有合理的壁厚 某些铸件缺陷的产生，往往是由于铸件结构设计不合理而造成的。采用合理的铸件结构，可防止许多缺陷。 每一种铸造合金，都有一个合适的壁厚范围，选择得当，既可保证铸件性能（机械性能）要求，又便于铸造生产。在确定铸件壁厚时一般应综合考虑以下三个方面：保证铸件达到所需要的强度和刚度；尽可能节约金属；铸造时没有多大困难。 （1 ）壁厚应不小于最小壁厚 在一定的铸造条件下，铸造合金能充满铸型的最小壁厚称为该铸造合金的最小壁厚。为了避免铸件的浇不足和冷隔等缺陷，应使铸件的设计壁厚不小于最小壁厚。各种铸造工艺条件下，铸件最小允许壁厚见表7-1?表7-5 表1-1砂型铸造时铸件最小允许壁厚（单位：mm） 合金种类铸件最大轮廓尺寸为下列值时/ mm
实验一误差分析 实验1.1（病态问题） 实验目的：算法有“优”与“劣”之分，问题也有“好”与“坏”之别。对数值方法的研究而言，所谓坏问题就是问题本身对扰动敏感者，反之属于好问题。通过本实验可获得一个初步体会。 数值分析的大部分研究课题中，如线性代数方程组、矩阵特征值问题、非线性方程及方程组等都存在病态的问题。病态问题要通过研究和构造特殊的算法来解决，当然一般要付出一些代价（如耗用更多的机器时间、占用更多的存储空间等）。 问题提出：考虑一个高次的代数多项式 显然该多项式的全部根为1,2,…,20共计20个，且每个根都是单重的。现考虑该多项式的一个扰动 其中ε（1.1）和（1.221,,,a a 的输出b ”和“poly ε。 （1（2 （3）写成展 关于α solve 来提高解的精确度，这需要用到将多项式转换为符号多项式的函数poly2sym,函数的具体使用方法可参考Matlab 的帮助。 实验过程： 程序： a=poly(1:20); rr=roots(a); forn=2:21 n form=1:9 ess=10^(-6-m);
ve=zeros(1,21); ve(n)=ess; r=roots(a+ve); -6-m s=max(abs(r-rr)) end end 利用符号函数：（思考题一）a=poly(1:20); y=poly2sym(a); rr=solve(y) n
很容易的得出对一个多次的代数多项式的其中某一项进行很小的扰动，对其多项式的根会有一定的扰动的，所以对于这类病态问题可以借助于MATLAB来进行问题的分析。 学号：06450210 姓名：万轩 实验二插值法
本科实验报告 课程名称：计算机数值方法 实验项目：方程求根、线性方程组的直接解 法、线性方程组的迭代解法、代数插值和最 小二乘拟合多项式 实验地点：行勉楼 专业班级： ******** 学号： ********* 学生姓名： ******** 指导教师：李誌，崔冬华 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 探究结论：
简单计算器 姓名: 周吉祥 实验目的：模仿日常生活中所用的计算器，自行设计一个简单的计算器程序，实现简单的计算功能。 实验内容： (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; }
重庆工商大学 《统计学》实验报告 实验课程：统计学 _ 指导教师：陈正伟 _ 专业班级： 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 同理可得
砂型铸造工艺流程 砂型铸造工艺流程图 制作木模-造型-熔化-浇注-落砂-冒口拆除-检验入库 熔模铸造工艺 失蜡铸造现在称为熔模铸造。这是一种很少切割或不切割的铸造工艺，是铸造行业的一项优秀技术。它被广泛使用。它不仅适用于各种类型和合金的铸造，而且可以生产出比其他铸造方法具有更高尺寸精度和表面质量的铸件，甚至复杂的、耐高温的、难以加工的、其他铸造方法难以铸造的铸件也可以通过熔模精密铸造来铸造。 熔模铸造是在古代蜡模铸造的基础上发展起来的。作为一个古老的文明，中国是最早使用这项技术的国家之一。早在公元前几百年，中国古代劳动人民就创造了这种失传的铸蜡技术，用来铸造钟鼎和具有各种精美图案和文字的器皿，如春秋时期曾侯乙墓的青铜板。曾侯乙墓雕像板的底座是多条龙缠绕在一起，首尾相连，上下交错，形成一个中间镂空的多层云纹图案。这些图案很难用普通的铸造工艺来制作，而失蜡法的铸造工艺可以利用石蜡无强度、易雕刻的特点，用普通的工具雕刻出与曾侯乙墓的雕像板相同的石蜡工艺品，然后加入浇注系统，经过上漆、脱蜡、浇注，得到精美的曾侯乙雕像板 现代熔模铸造法在20世纪40年代实际应用于工业生产当时，航空喷气发动机的发展要求制造具有复杂形状、精确尺寸和光滑表面的耐热合金部件，如叶片、叶轮和喷嘴。由于耐热合金材料难以加工，零件形状复杂，因此不可能或难以用其他方法制造。因此，需要找到一
种新的精确的成型工艺。因此，现代熔模铸造法借鉴了古代传下来的失蜡铸造法，通过对 材料和工艺的改进，在古代工艺的基础上取得了重要的发展。因此，航空工业的发展促进了熔模铸造的应用，熔模铸造的不断改进也为航空工业进一步提高性能创造了有利条件。 中国在20世纪50年代和60年代开始将熔模铸造应用于工业生产此后，这种先入为主的铸造技术得到了极大的发展，并已广泛应用于航空、汽车、机床、船舶、内燃机、燃气轮机、电信仪器、武器、医疗器械、切割工具等制造业，以及工艺品的制造。所谓的 熔模铸造工艺简单地指用易熔材料(如蜡或塑料)制作易熔模型(称为熔模或模型)，在其上涂覆几层特殊的耐火涂层，干燥并硬化形成整体外壳，然后用蒸汽或温水将外壳上的模型熔化，然后将外壳放入砂箱中，在其周围填充干砂，最后将模具放入穿透式烘烤器中进行高温烘烤(例如，当使用高强度外壳时，脱模后的外壳可以不造型直接烘烤)、模具或外壳 熔模铸件尺寸精度高，一般可达CT4-6(砂型铸造CT10~13，压铸CT5~7)。当然，由于熔模铸造工艺过程复杂，影响铸件尺寸精度的因素很多，如模具材料的收缩、熔模的变形、加热和冷却过程中模壳的线性变化、合金的收缩率以及铸件在凝固过程中的变形等。因此，普通熔模铸件的尺寸精度相对较高，但其一致性仍有待提高(使用中高温蜡材料的铸件的尺寸一致性有待提高)用 压制熔体模具时，采用型腔表面光洁度高的型材，因此熔体模具的
(此文档为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语言实现算法的程序代码 #include
计算方法实验报告格式 小组名称: 组长姓名(班号)： 小组成员姓名(班号): 按贡献排序情况: 指导教师评语: 小组所得分数: 一个完整的实验，应包括数据准备、理论基础、实验内容及方法，最终对实验结果进行分析，以达到对理论知识的感性认识，进一步加深对相关算法的理解，数值实验以实验报告形式完成，实验报告格式如下： 一、实验名称 实验者可根据报告形式需要适当写出． 二、实验目的及要求 首先要求做实验者明确，为什么要做某个实验，实验目的是什么，做完该实验应达到什么结果，在实验过程中的注意事项，实验方法对结果的影响也可以以实验目的的形式列出． 三、算法描述（实验原理与基础理论） 数值实验本身就是为了加深对基础理论及方法的理解而设置的，所以要求将实验涉及到的理论基础，算法原理详尽列出． 四、实验内容 实验内容主要包括实验的实施方案、步骤、实验数据准备、实验的算法以及可能用到的仪器设备． 五、程序流程图 画出程序实现过程的流程图，以便更好的对程序执行的过程有清楚的认识，在程序调试过程中更容易发现问题． 六、实验结果 实验结果应包括实验的原始数据、中间结果及实验的最终结果，复杂的结果可以用表格
形式列出，较为简单的结果可以与实验结果分析合并出现． 七、实验结果分析 实验结果分析包括对对算法的理解与分析、改进与建议． 数值实验报告范例 为了更好地做好数值实验并写出规范的数值实验报告，下面给出一简单范例供读者参考． 数值实验报告 小组名称: 小组成员(班号): 按贡献排序情况: 指导教师评语: 小组所得分数: 一、实验名称 误差传播与算法稳定性． 二、实验目的 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
?砂型铸造的基本过程https://www.doczj.com/doc/8513421299.html,/20111213/62031.html ?砂型铸造有六个基本步骤: 1) 把模样放入砂中制成一个模具。 2) 在浇注系统中把原型和砂子接合起来。 3) 把模样去掉。 4) 把模具的空隙用熔化了的金属填充起来。 5) 让金属冷却。 6) 把砂型模具敲掉取出铸件。 砂型铸造案例 项目导入：轴承座铸件的造型工艺方案。 铸件简图：轴承座如图2-1所示。 铸件材料：HT150。 体积参数：轮廓尺寸240mm′65mm′75mm，铸件重量约5kg。 生产性质：单件生产。 项目要求：确定铸件的造型工艺方案并完成造型操作。
图2-1 轴承座 将液体金属浇入用型砂捣实成的铸型中，待凝固冷却后，将铸型破坏，取出铸件的铸造方法称为砂型铸造。砂型铸造是传统的铸造方法，它适用于各种形状、大小及各种常用合金铸件的生产。套筒的砂型铸造过程如图2-2所示，主要工序包括制造模样型芯盒、制备造型材料、造型、制芯、合型、熔炼、浇注、落砂、清理与检验等。 图2-2 套筒的砂型铸造过程 铸件生产前需根据零件图绘制出铸造工艺图，铸造工艺图是在零件图上用各种工艺符号及参数表示出铸造工艺方案的图形。其中包括：浇注位置，铸型分型面，型芯的数量、形状、尺寸及其固定方法，加工余量，收缩率，浇注系统，起模斜度，冒口和冷铁的尺寸和布置等。铸造工艺图是指导模样(型芯盒)设计、生产准备、铸型制造和铸件检验的基本工艺文件。砂型铸造主要工序包括： (1) 根据零件图制造模样和型芯盒； (2) 配制性能符合要求的型(芯)砂； (3) 用模样和型芯盒进行造型和造芯； (4) 烘干型芯(或砂型)并合型； (5) 熔炼金属并进行浇注； (6) 落砂、清理和检验。 2.1.1 常用造型工模具 1. 砂箱
实验一 误差分析 实验（病态问题） 实验目的：算法有“优”与“劣”之分，问题也有“好”与“坏”之别。对数值方法的研究而言，所谓坏问题就是问题本身对扰动敏感者，反之属于好问题。通过本实验可获得一个初步体会。 数值分析的大部分研究课题中，如线性代数方程组、矩阵特征值问题、非线性方程及方程组等都存在病态的问题。病态问题要通过研究和构造特殊的算法来解决，当然一般要付出一些代价（如耗用更多的机器时间、占用更多的存储空间等）。 问题提出：考虑一个高次的代数多项式 )1.1() ()20()2)(1()(20 1∏=-=---=k k x x x x x p 显然该多项式的全部根为1,2,…,20共计20个，且每个根都是单重的。现考虑该多项式的一个扰动 )2.1(0 )(19=+x x p ε 其中ε是一个非常小的数。这相当于是对（）中19x 的系数作一个小的扰动。我们希望比较（）和（）根的差别，从而分析方程（）的解对扰动的敏感性。 实验内容：为了实现方便，我们先介绍两个Matlab 函数：“roots ”和“poly ”。 roots(a)u = 其中若变量a 存储n+1维的向量，则该函数的输出u 为一个n 维的向量。设a 的元素依次为121,,,+n a a a ，则输出u 的各分量是多项式方程 01121=+++++-n n n n a x a x a x a 的全部根；而函数 poly(v)b =
的输出b 是一个n+1维变量，它是以n 维变量v 的各分量为根的多项式的系数。可见“roots ”和“poly ”是两个互逆的运算函数。 ;000000001.0=ess );21,1(zeros ve = ;)2(ess ve = ))20:1((ve poly roots + 上述简单的Matlab 程序便得到（）的全部根，程序中的“ess ”即是（）中的ε。 实验要求： （1）选择充分小的ess ，反复进行上述实验，记录结果的变化并分析它们。 如果扰动项的系数ε很小，我们自然感觉（）和（）的解应当相差很小。计算中你有什么出乎意料的发现表明有些解关于如此的扰动敏感性如何 （2）将方程（）中的扰动项改成18x ε或其它形式，实验中又有怎样的现象 出现 （3）（选作部分）请从理论上分析产生这一问题的根源。注意我们可以将 方程（）写成展开的形式， ) 3.1(0 ),(1920=+-= x x x p αα 同时将方程的解x 看成是系数α的函数，考察方程的某个解关于α的扰动是否敏感，与研究它关于α的导数的大小有何关系为什么你发现了什么现象，哪些根关于α的变化更敏感 思考题一：（上述实验的改进） 在上述实验中我们会发现用roots 函数求解多项式方程的精度不高，为此你可以考虑用符号函数solve 来提高解的精确度，这需要用到将多项式转换为符号多项式的函数poly2sym,函数的具体使用方法可参考Matlab 的帮助。
南京信息工程大学实验（实习）报告 一、实验目的: 用最小二乘法将给定的十个点拟合成三次多项式。 二、实验步骤： 用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．分析与论证：通过上面的实验，收集到一些数据，观察到一些现象，对其分析，得出事实与假设的关系，通过归纳、概括等方法，得到结论。
实验报告 一、求方程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-1 零件结构的铸造工艺性分析 铸造工艺性，是指零件结构既有利于铸造工艺过程的顺利进行，又有利于保证铸件质量。 还可定义为：铸造零件的结构除了应符合机器设备本身的使用性能和机械加工的要求外，还应符合铸造工艺的要求。这种对铸造工艺过程来说的铸件结构的合理性称为铸件的铸造工艺性。 另定义：铸造工艺性是指零件的结构应符合铸造生产的要求，易于保证铸件品质，简化铸造工艺过程和降低成本。 铸造工艺性不好，不仅给铸造生产带来麻烦，不便于操作，还会造成铸件缺陷。因此，为了简化铸造工艺，确保铸件质量，要求铸件必须具有合理的结构。 一、铸件质量对铸件结构的要求 1．铸件应有合理的壁厚 某些铸件缺陷的产生，往往是由于铸件结构设计不合理而造成的。采用合理的铸件结构，可防止许多缺陷。 每一种铸造合金，都有一个合适的壁厚范围，选择得当，既可保证铸件性能（机械性能）要求，又便于铸造生产。在确定铸件壁厚时一般应综合考虑以下三个方面：保证铸件达到所需要的强度和刚度；尽可能节约金属；铸造时没有多大困难。 （1）壁厚应不小于最小壁厚 在一定的铸造条件下，铸造合金能充满铸型的最小壁厚称为该铸造合金的最小壁厚。为了避免铸件的浇不足和冷隔等缺陷，应使铸件的设计壁厚不小于最小壁厚。各种铸造工艺条件下，铸件最小允许壁厚见表7-1～表7-5 合金种类铸件最大轮廓尺寸为下列值时/㎜ ﹤200200-400400-800800-12501250-2000﹥ 2000 碳素铸钢 低合金钢 高锰钢 不锈钢、耐热钢灰铸铁 孕育铸铁 （HT300以上）球墨铸铁8 8-9 8-9 8-11 3-4 5-6 3-4 9 9-10 10 10-12 4-5 6-8 4-8 11 12 12 12-16 5-6 8-10 8-10 14 16 16 16-20 6-8 10-12 10-12 16～18 20 20 20-25 8-10 12-16 12-14 20 25 25 - 10-12 16-20 14-16铸件最大轮廓为下列值时mm 实验报告一 题目：非线性方程求解 摘要：非线性方程的解析解通常很难给出，因此线性方程的数值解法就尤为重要。本实验采用两种常见的求解方法二分法和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; else b=c; end R=b-a;%求出误差 k=k+1; end x=c%给出解 Newton法及改进的Newton法源程序：clear %%%% 输入函数 f=input('请输入需要求解函数>>','s') %%%求解f(x)的导数 df=diff(f);铸造工艺设计基础
(完整版)哈工大-数值分析上机实验报告
相关主题
文本预览