数值分析学号姓名2016_2上机实习报告
数值分析2016上机实习报告要求
1．应提交一份完整的实习报告。具体要求如下：
（1） 报告要排版，美观漂亮；要标明姓名、学号、专业和联系电话；
（2） 要有目录，指明题目、程序、计算结果，图标和分析等内容所在位置，作
到信息简明而完全；
（3）要有总结，全方位总结机编程计算的心得体会；
（4）尽量使报告清晰明了，一般可将计算结果、图表及对比分析放在前面，程序
清单作为附录放在后面，程序中关键部分要有中文说明或标注，指明该部分的功能和作用。
2．程序需完好保存到期末考试后的一个星期，以便老师索取用于验证、询问或质疑部分内容。
3．认真完成实验内容，可以达到既学习计算方法又提高计算能力的目的，还可以切身体会书本内容之精妙所在，期间可以得到很多乐趣。
4．拷贝或抄袭他人结果是不良行为，将视为不合格。
5．请按任课老师要求的时间和载体（电子或纸质）提交给任课老师。
6.提交电子邮件请务必按述格式命名邮件主题和文档名，以方便统计
学号姓名2016_2上机实习报告
如：2015201597张帅哥数值分析2016-2计算报告
数值分析上机试题
1． 编程计算∑∞11=n cn ， 其中c= 4.4942?10307， 给出并观察计算结果，若有问题，分析之。
2． 利用牛顿法求方程ln 2x x -=(x-ln x=2)于区间[2,4]的根，考虑不同初值下牛顿法的收敛情况。
3． 给出1251
2+=x y 在xk=0+0.25k 处的值yk, k=0,1,2,3,4. 请给出由节点 xk 确定的三次样
条插值函数S(x), 使其满足条件：S /(0)=0, S /(1)=-0.074, 分析逼近效果如何？
4. 松弛因子对SOR 法收敛速度的影响。
用SOR 法求解方程组Ax =b ,其中
3-2-.2-2-3-B ,41141......14114 要求程序中不存系数矩阵A ,分别对不同的阶数取w=1.1, 1.2, ...,1.9进行迭代，记录近似解x (k)达到||x (k)-x (k-1)||<10-6时所用的迭代次数k ，观察松弛因子对收敛速度的影响，并观察当w ≤0或w ≥2会有什么影响？
5. 用Ru n ge-Kutt a 4阶算法对初值问题y /=-20*y ，y (0)=1按不同步长求解，用于观察稳定区间的
作用，推荐两种步长h=0.1,0.2。
注：此方程的精确解为：y =e -20x
有限元分析实验报告
武汉理工大学 学生实验报告书 实验课程名称机械中的有限单元分析 开课学院机电工程学院 指导老师姓名 学生姓名 学生专业班级机电研 1502班 2015—2016 学年第2学期
实验一方形截面悬臂梁的弯曲的应力与变形分析 钢制方形悬臂梁左端固联在墙壁，另一端悬空。工作时对梁右端施加垂直向下的30KN的载荷与60kN的载荷，分析两种集中力作用下该悬臂梁的应力与应变，其中梁的尺寸为10mmX10mmX100mm的方形梁。 1.1方形截面悬臂梁模型建立 建模环境：DesignModeler 15.0。 定义计算类型：选择为结构分析。 定义材料属性：弹性模量为2.1Gpa，泊松比为0.3。 建立悬臂式连接环模型。 （1）绘制方形截面草图：在DesignModeler中定义XY平面为视图平面，并正视改平面，点击sketching下的矩形图标，在视图中绘制10mmX10mm的矩形。（2）拉伸：沿着Z方向将上一步得到的矩阵拉伸100mm，即可得到梁的三维模型，建模完毕，模型如下图1.1所示。 图1.1 方形截面梁模型 1.2 定义单元类型： 选用6面体20节点186号结构单元。 网格划分：通过选定边界和整体结构，在边界单元划分数量不变的情况下，通过分别改变节点数和载荷大小，对同一结构进行分析，划分网格如下图1.2所示：
图1.2 网格划分 1.21 定义边界条件并求解 本次实验中，讲梁的左端固定，将载荷施加在右端，施以垂直向下的集中力，集中力的大小为30kN观察变形情况，再将力改为50kN，观察变形情况，给出应力应变云图，并分析。 （1）给左端施加固定约束； （2）给悬臂梁右端施加垂直向下的集中力； 1.22定义边界条件如图1.3所示： 图1.3 定义边界条件 1.23 应力分布如下图1.4所示： 定义完边界条件之后进行求解。
计算方法上机实验报告
《计算方法》上机实验报告 班级：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 四、运行结果： 实验二：
2014级硕士研究生数值分析上机实习报告
2014级硕士研究生数值分析上机实习（第一次） 姓名：学号：学院： 实习题目：分别用二分法和Newton迭代法求方程x3■ 2x210x-20=0的根.实习目的：掌握两种解法，体会两种解法的收敛速度. 实习要求：用C程序语言编程上机进行计算，精确到8位有效数字. 报告内容： 1.确定实根的个数以及所在区间 2.将最后两次计算结果填入下表（保留8位数字）： 3.实习过程中遇到哪些问题？如何解决？有何心得体会?
4.两种解法的计算程序（此页写不下时可以加页）：
2014级硕士研究生数值分析上机实习（第二次）姓名：学号：学院： 实习题目：计算8阶三对角矩阵A=tridiag（0.235, 1.274, 0.235）的行列式.实习目的：掌握计算行列式的方法. 实习要求：首先选择一种算法，然后用C程序语言编程上机进行计算.报告内容： 1.简单描述所采用的算法: 2?计算结果： A 3.实习过程中遇到哪些问题？如何解决？有何心得体会?
4.写出C语言计算程序（此页写不下时可以加页）：
2014级硕士研究生数值分析上机实习（第三次） 姓名：学号：学院： 分别用Jacobi迭代法和Gauss-Seidel迭代法求解线性方程组实习题目： 2lx + 9.8y+ 3.4z= 6.7 <2.7x + 1.8y+ 7.2z= 2.4 8.6x + 1.5y + 3.4z = 1.9 实习目的：感受两种迭代法的收敛速度. 首先构造收敛的Jacobi迭代法和Gauss-Seidel迭代法，然后用实习要求： C程序语言编程上机进行求解，初始值均取为0,精确到4位小 数. 报告内容： 1.写出收敛的Jacobi迭代法和Gauss-Seidel迭代法:
有限元上机实验报告
有限元上机实验报告结构数值分析与程序设计 上机实验 院系: 土木工程与力学学院专业: 土木工程 班级: 姓名: 学号: 指导教师: 1、调试教材P26-30程序FEM1。 1.1、输入数据文件为: 6,4,12,6,1.0E0,0.0,1.0,0.0,1 3,1,2 5,2,4 3,2,5 6,3,5 0.0,2.0 0.0,1.0 1.0,1.0 0.0,0.0 1.0,0.0 2.0,0.0
1,3,7,8,10,12 1.2、输出数据文件为: NN NE ND NFIX E ANU T GM NTYPE 6 4 12 60.1000E+01 0.000 1.0000.0000E+00 1 NODE X-LOAD Y-LOAD 1 0.000000E+00 -0.100000E+01 2 0.000000E+00 0.000000E+00 3 0.000000E+00 0.000000E+00 4 0.000000E+00 0.000000E+00 5 0.000000E+00 0.000000E+00 6 0.000000E+00 0.000000E+00 NODE X-DISP Y-DISP 1 -0.879121E-15 -0.325275E+01 2 0.879121E-16 -0.125275E+01 3 -0.879121E-01 -0.373626E+00 4 0.117216E-1 5 -0.835165E-15 5 0.175824E+00 -0.293040E-15 6 0.175824E+00 0.263736E-15 ELEMENT X-STR Y-STR XY-STR 1 -0.879121E-01 -0.200000E+01 0.439560E+00 2 0.175824E+00 -0.125275E+01 0.256410E-15 3 -0.879121E-01 -0.373626E+00 0.307692E+00 4 0.000000E+00 -0.373626E+00 -0.131868E+00 2、修改FEM1,计算P31例2-2。
中南大学ANSYS上机实验报告
ANSYS上机实验报告 小组成员：郝梦迪、赵云、刘俊 一、实验目的和要求 本课程上机练习的目的是培养学生利用有限单元法的商业软件进行数值计算分析，重点是了解和熟悉ANSYS的操作界面和步骤，初步掌握利用ANSYS建立有限元模型，学习ANSYS分析实际工程问题的方法，并进行简单点后处理分析，识别和判断有限元分析结果的可靠性和准确性。 二、实验设备和软件 台式计算机，ANSYS10.0软件 三、基本步骤 1）建立实际工程问题的计算模型。实际的工程问题往往很复杂，需要采用适当的模型在计算精度和计算规模之间取得平衡。常用的建模方法包括：利用几何、载荷的对称性简化模型，建立等效模型。 2）选择适当的分析单元，确定材料参数。侧重考虑一下几个方面：是否多物理耦合问题，是否存在大变形，是否需要网格重划分。 3）前处理（Preprocessing）。前处理的主要工作内容如下：建立几何模型（Geometric Modeling），单元划分（Meshing）与网格控制，给定约束（Constraint）和载荷（Load）。在多数有限元软件中，不能指定参数的物理单位。用户在建模时，要确定力、长度、质量及派生量的物理单位。在建立有限元模型时，最好使用统一的物理单位，这样做不容易弄错计算结果的物理单位。建议选用kg，N，m，sec；常采用kg，N，mm，sec。 4）求解（Solution）。选择求解方法，设定相应的计算参数，如计算步长、迭代次数等。 5）后处理（Postprocessing）。后处理的目的在于确定计算模型是否合理、计算结果是否合理、提取计算结果。可视化方法（等值线、等值面、色块图）显
计算方法第二章方程求根上机报告
实验报告名称 班级：学号：姓名：成绩： 1实验目的 1）通过对二分法与牛顿迭代法作编程练习与上级运算，进一步体会二分法与牛顿迭代法的不同特点。 2）编写割线迭代法的程序，求非线性迭代法的解，并与牛顿迭代法。 2 实验内容 用牛顿法和割线法求下列方程的根 x^2-e^x=0; x*e^x-1=0; lgx+x-2=0; 3实验步骤 1)根据二分法和牛顿迭代法，割线法的算法编写相应的求根函数； 2)将题中所给参数带入二分法函数，确定大致区间； 3)用牛顿迭代法和割线法分别对方程进行求解； 3 程序设计 牛顿迭代法x0=1.0; N=100; k=0; eps=5e-6; delta=1e-6; while(1) x1=x0-fc1(x0)/fc2(x0); k=k+1; if k>N disp('Newmethod failed')
break end if(abs(x1-x0) 数值分析实验（2014,9,16~10,28） 信计1201班，人数34人 数学系机房 数值分析 计算实习报告册 专业 学号 姓名 2014~2015年第一学期 实验一 数值计算的工具 Matlab 1.解释下MATLAB 程序的输出结果 程序： t=0.1 n=1:10 e=n/10-n*t e 的结果：0 0 -5.5511e-017 0 0 -1.1102e-016 -1.1102e-016 0 0 0 2.下面MATLAB 程序的的功能是什么？ 程序： x=1;while 1+x>1,x=x/2,pause(0.02),end 用迭代法求出x=x/2,的最小值 x=1;while x+x>x,x=2*x,pause(0.02),end 用迭代法求出x=2*x,的值，使得2x>X x=1;while x+x>x,x=x/2,pause(0.02),end 用迭代法求出x=x/2,的最小值，使得2x>X 3.考虑下面二次代数方程的求解问题 02=++c bx ax 公式a ac b b x 242-+-=是熟知的，与之等价地有a c b b c x 422-+-=，对于 1,100000000,1===c b a ，应当如何选择算法。 应该用a ac b b x 242-+-=计算，因为b 做分母 4.函数)sin(x 有幂级数展开...! 7!5!3sin 7 53+-+-=x x x x x 利用幂级数计算x sin 的MATLAB 程序为 function s=powersin(x) s=0; t=x; n=1; while s+t~=s; s=s+t ； t=-x^2/（(n+1)*(n+2)）*t ； n=n+2； end t1=cputime; pause(10); t2=cputime; t0=t2-t1 (a)解释上述程序的终止准则。 当s+t=s ，终止循环。 (b)对于2/21,2/11 ,2/πππ=x 计算的进度是多少？分别计算多少项？ X=pi/2时，s =1.0000 x=11pi/2时，s=-1.0000 x=21pi/2时，s =0.9999 Cputime 分别是0.1563 0.0469 0.0156 5.考虑调和级数∑∞ =11 n n ，它是微积分中的发散级数，在计算机上计算该级数的部 分和，会得到怎么样的结果，为什么？ function s=fun(n) s=0; t=1/n; for i=1:n s=s+1/i; end 当n=100时s =5.1874 当n=80时s =4.9655 当n=50时，s =4.4992 当n=10时，s =2.9290 6.指数函数的级数展开...! 3!213 2++++=x x x e x ，如果对于0 数值分析上机实验报告 《数值分析》上机实验报告 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 a b 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)(3 2 2 5 333647>?''<-=-=''<-=-='<>+-=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 指导教师： 姓名： 学号： 专业： 联系电话： 上海交通大学 目录 序言 (3) 实验课题(一) 雅可比迭代法和高斯-塞得尔迭代法的收敛性和收敛速度 (4) 数值分析 (5) 实验课题(二) 松弛因子对SOR法收敛速度的影响 (6) 数值分析 (12) 总结 (13) 附录（程序清单） (14) 1.雅可比迭代法和高斯-塞得尔迭代法的收敛性和收敛速度 (14) 雅可比迭代法： (14) 高斯-塞得尔迭代法： (16) 2.松弛因子对SOR法收敛速度的影响 (18) 松弛法（SOR） (18) 序言 随着科学技术的发展，提出了大量复杂的数值计算问题，在实际解决这些计算问题的长期过程中，形成了计算方法这门学科，专门研究各种数学问题的数值解法（近似解法），包括方法的构造和求解过程的误差分析，是一门内容丰富，有自身理论体系的实用性很强的学科。 解决工程问题，往往需要处理很多数学模型，这就要花费大量的人力和时间，但是还有不少数学模型无法用解析法得到解。使用数值方法并利用计算机，就可以克服这些困难。事实上，科学计算已经与理论分析、科学实验成为平行的研究和解决科技问题的科学手段，经常被科技工作者所采用。作为科学计算的核心内容——数值分析（数值计算方法），已逐渐成为广大科技工作者必备的基本知识并越来越被人重视。 由于数值方法是解数值问题的系列计算公式，所以数值方法是否有效，不但与方法本身的好坏有关，而且与数值问题本身的好坏也有关，因此，研究数值方法时，不但需要研究数值方法的好坏，即数值稳定性问题，而且还需要研究数值问题本身的好坏，即数值问题的性态，以及它们的判别问题。 数值计算的绝大部分方法都具有近似性，而其理论又具有严密的科学性，方法的近似值正是建立在理论的严密性基础上，根据计算方法的这一特点。因此不仅要求掌握和使用算法，还要重视必要的误差分析，以保证计算结果的可靠性。数值计算还具有应用性强的特点，计算方法的绝大部分方法如求微分方程近似解，求积分近似值，求解超越方程，解线性方程组等都具有较强的实用性，而插值法，最小二乘法，样条函数等也都是工程技术领域中常用的，有实际应用价值的方法。 应用数值计算方法解决科学研究和工程中的实际问题，首先要建立描述具体问题的适当数学模型；其次，要选择一定的计算方法并制定相应的计算方案；并编制、设计合理的计算机程序；最后由计算机计算出数值结果。其中，计算方案的设计和计算方法的选择是上述求解过程中极为重要的环节，是程序设计和分析数值计算结果准确性的基础。 本实验所使用的是Visual C++ 6.0程序设计语言。 C++是在C语言的基础上拓展而来的。C++引入了面向对象的机制，同时又充分保留了C语言的简洁性和高效性，并且与C语言完全兼容，具有C语言高效灵活、功能强大、可移植性好等诸多优点，是面向对象程序设计的最佳语言之一。 有限元实验报告 T1013-5 20100130508 蔡孟迪 ANSYS有限元上机报告（一） 班级:T1013-5 学号:20100130508 姓名:蔡孟迪 上机题目： 图示折板上端固定，右侧受力F=1000N，该力均匀分布在边缘各节点上；板厚t=2mm 材料选用低碳钢，弹性模量E=210Gpa，μ=0.33. 一、有限元分析的目的： 1．利用ANSYS构造实体模型； 2．根据结构的特点及所受载荷的情况，确定所用单元类型；正确剖分网格并施加外界条件；3．绘制结构的应力和变形图，给出最大应力和变形的位置及大小；并确定折板角点Ａ处的应力和位移； 4．研究网格密度对Ａ处角点应力的影响； 5．若在Ａ处可用过渡圆角，研究Ａ处圆角半径对Ａ处角点应力的影响。 二、有限元模型的特点： 1．结构类型 本结构属于平面应力类型 2.单位制选择 本作业选择N（牛），mm（毫米），MPa（兆帕）。 3.建模方法 采用自左向右的实体建模方法。 4.定义单元属性及类型 1）材料属性：弹性模量：EX=2.10E5MPa, 泊松比：PRXY=0.33 2）单元类型：在Preferences选Structural，Preprocessor>ElemmentType>Add/Edit/Delete中定义单元类型为：Quad4 node 182，K3设置为：平面薄板问题（Plane strs w/thk） 3）实常数：薄板的厚度THK=2mm 5.划分网格 在MeshTool下选set,然后设置SIZE Element edge length的值，再用Mesh进行网格划分。6.加载和约束过程：在薄板的最上端施加X、Y方向的固定铰链，在薄板的最右端施加1000N 的均匀布置的载荷。 计算方法实验报告 班级： 学号： 姓名： 成绩： 1 舍入误差及稳定性 一、实验目的 （1）通过上机编程，复习巩固以前所学程序设计语言及上机操作指令； （2）通过上机计算，了解舍入误差所引起的数值不稳定性 二、实验内容 1、用两种不同的顺序计算10000 21n n -=∑，分析其误差的变化 2、已知连分数() 1 01223//(.../)n n a f b b a b a a b =+ +++，利用下面的算法计算f ： 1 1 ,i n n i i i a d b d b d ++==+ (1,2,...,0 i n n =-- 0f d = 写一程序，读入011,,,...,,,...,,n n n b b b a a 计算并打印f 3、给出一个有效的算法和一个无效的算法计算积分 1 041 n n x y dx x =+? (0,1,...,1 n = 4、设2 2 11N N j S j == -∑ ，已知其精确值为1311221N N ?? -- ?+?? （1）编制按从大到小的顺序计算N S 的程序 （2）编制按从小到大的顺序计算N S 的程序 （3）按两种顺序分别计算10001000030000,,,S S S 并指出有效位数 三、实验步骤、程序设计、实验结果及分析 1、用两种不同的顺序计算10000 2 1n n -=∑，分析其误差的变化 （1）实验步骤： 分别从1~10000和从10000~1两种顺序进行计算，应包含的头文件有stdio.h 和math.h （2）程序设计： a.顺序计算 #include 数值分析上机报告 班级：20级学隧2班 姓名：000000000 学号：00000000000 目录 1 序言 (6) 2 题目 (7) 2.1 题2 (7) 2.1.1 题目内容 (7) 2.1.2 MATLAB程序 (8) 2.1.3 计算结果 (8) 2.1.4 图形 (9) 2.1.5 分析 (14) 2.2 题3 (14) 2.2.1 题目内容 (14) 2.2.2 程序 (14) 2.2.3 计算结果 (14) 2.2.4 图形 (15) 2.2.5 分析 (16) 2.3 选做题5 (16) 2.3.1方法介绍 (17) 2.3.2计算结果及分析 (17) 3总结 (18) 4．附录 (19) 4.1 题1程序代码 (19) 4.2 题2程序代码 (22) 4.3 题3程序代码 (26) 数值分析2015上机实习报告要求 1．应提交一份完整的实习报告。具体要求如下： （1）报告要排版，美观漂亮（若是纸质要有封面，封面上）要标明姓名、学号、专业和联系电话； （2）要有序言，说明所用语言及简要优、特点，说明选用的考量； （3）要有目录，指明题目、程序、计算结果，图标和分析等内容所在位置，作到信息简明而完全； （4）要有总结，全方位总结机编程计算的心得体会； （5）尽量使报告清晰明了，一般可将计算结果、图表及对比分析放在前面，程序清单作为附录放在后面，程序中关键部分要有中文说明或标注， 指明该部分的功能和作用。 2．程序需完好保存到期末考试后的一个星期，以便老师索取用于验证、询问或质疑部分内容。 3．认真完成实验内容，可以达到既学习计算方法又提高计算能力的目的，还可以切身体会书本内容之精妙所在，期间可以得到很多乐趣。 4．拷贝或抄袭他人结果是不良行为，将视为不合格。 5．请按任课老师要求的时间和载体（电子或纸质）提交给任课老师。 一、实验目的 通过上机对有限元法的基本原理和方法有一个更加直观、深入的理解；通过对本实验所用软件平台Ansys 的初步涉及，为将来在设计和研究中利用该类大型通用CAD/CAE 软件进行工程分析奠定初步基础。 二、实验设备 机械工程软件工具包Ansys 三、实验内容及要求 1) 简支梁如图3.1.1所示，截面为矩形，高度h=200mm ，长度L=1000mm ，厚 度t=10mm 。上边承受均布载荷，集度q=1N/mm2，材料的E=206GPa ，μ=0.29。平面应力模型。 X 方向正应力的弹性力学理论解如下： 图3.1.1 ①在Ansys 软件中用有限元法探索整个梁上x σ，y σ的分布规律。 ②计算下边中点正应力x σ的最大值；对单元网格逐步加密，把x σ的计算值与理论解对比，考察有限元解的收敛性。 ③针对上述力学模型，对比三节点三角形平面单元和4节点四边形平面等参元的求解精度。 2) 一个正方形板，边长L = 1000mm ，中心有一小孔，半径R = 100mm ，左右边 受均布拉伸载荷，面力集度q = 25MPa ，如图 3.2.1所示。材料是 206E GPa =，0.3μ=，为平面应力模型。当边长L 为无限大时，x = 0截面上理论解为： ) 534()4 (6222 23-+-=h y h y q y x L h q x σ )32(2|44 220r R r R q x x ++==σ 其中R 为圆孔半径，r 为截面上一点距圆心的距离。x = 0截面上孔边（R r =）应力q x 3=σ。所以理论应力集中系数为3.0。 图3.2.1 用四边形单元分析x = 0截面上应力的分布规律和最大值，计算孔边应力集中系数，并与理论解对比。利用对称性条件，取板的四分之一进行有限元建模。 3) 如图3.3.1所示，一个外径为0.5m ，内径为0.2m ，高度为0.4m 的圆筒，圆 筒的外壁施加100MPa 的压强，圆筒的内部约束全部的自由度，材料参数是密度。 使用平面单元，依照轴对称的原理建模分析。 q 数值分析上机实习报告要求 1．应提交一份完整的实习报告。具体要求如下： （1）要有封面，封面上要标明姓名、学号、专业和联系电话； （2）要有序言，说明所用语言及简要优、特点，说明选用的考量； （3）要有目录，指明题目、程序、计算结果，图标和分析等内容所在位置，作到信息简明而完全； （4）要有总结，全方位总结机编程计算的心得体会； （5）尽量使报告清晰明了，一般可将计算结果、图表及对比分析放在前面，程序清单作为附录放在后面，程序中关键部分要有中文说明或标注，指明该部分 的功能和作用。 2．程序需完好保存到期末考试后的一个星期，以便老师索取用于验证、询问或质疑部分内容。 3．认真完成实验内容，可以达到既学习计算方法又提高计算能力的目的，还可以切身体会书本内容之精妙所在，期间可以得到很多乐趣。 4．拷贝或抄袭他人结果是不良行为，将视为不合格。 5．报告打印后按要求的时间提交给任课老师。 数值分析上机试题2 （选择其中1个题目） 5.用Ru n ge-Kutt a 4阶算法对初值问题y/=-20*y，y(0)=1按不同步长求解，用于观察稳定区间的 作用，推荐两种步长h=0.1,0.2。 注：此方程的精确解为：y=e-20x 6.实验内容 (1)实际验证梯形求积公式、Simpson求积公式、Newton-Cotes求积公式的代数精度。 (2)针对上述三个函数和积分区间[a,b]，实验观察梯形求积公式、Simpson求积公式和 Newton-Cotes求积公式的复化求积公式的实际计算效果。 y=exp(-x.^2).*sin(10*x)+4; a=1; b=3; y=sin(5*x)./x.^3;a=2*pi; b=4*pi; y=sin(5*x)./x.^3;a=2*pi;b=9.4248; 1 实验一——插值方法 实验学时：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]; 航空航天大学 《数值分析》计算实习报告 第一大题 学院：自动化科学与电气工程学院 专业：控制科学与工程 学生姓名： 学号： 教师： 电话： 完成日期： 2015年11月6日 航空航天大学 Beijing University of Aeronautics and Astronautics 实习题目： 第一题 设有501501?的实对称矩阵A ， ??? ???? ?????????=5011A a b c b c c b c b a 其中，064.0,16.0),501,,2,1(64.0)2.0sin()024.064.1(1 .0-==???=--=c b i e i i a i i 。矩阵A 的特征值为)501,,2,1(???=i i λ，并且有 ||min ||,501 150121i i s λλλλλ≤≤=≤???≤≤ 1.求1λ，501λ和s λ的值。 2.求A 的与数40 1 5011λλλμ-+=k k 最接近的特征值)39,,2,1(???=k k i λ。 3.求A 的(谱数)条件数2)A (cond 和行列式detA 。 说明： 1.在所用的算法中，凡是要给出精度水平ε的，都取12-10=ε。 2.选择算法时，应使矩阵A 的所有零元素都不储存。 3.打印以下容： （1）全部源程序； （2）特征值），，39,...,2,1(,s 5011=k k i λλλλ以及A det ,)A (cond 2的值。 4.采用e 型输出实型数，并且至少显示12位有效数字。 一、算法设计方案 1、求1λ，501λ和s λ的值。 由于||min ||,501 150121i i s λλλλλ≤≤=≤???≤≤，可知绝对值最大特征值必为1λ和501 λ其中之一，故可用幂法求出绝对值最大的特征值λ，如果λ=0，则1λ=λ，否则 501λ=λ。将矩阵A 进行一下平移： I -A A'λ= （1） 对'A 用幂法求出其绝对值最大的特征值'λ，则A 的另一端点特征值1λ或501λ为'λ+λ。 s λ为按模最小特征值，||min ||501 1i i s λλ≤≤=，可对A 使用反幂法求得。 2、求A 的与数40 1 5011λλλμ-+=k k 最接近的特征值)39,...,2,1(=k k i λ。 计算1)1,2,...,50=(i i λ-k μ，其模值最小的值对应的特征值k λ与k μ最接近。因此对A 进行平移变换： )39,,2,1k -A A k k ==（I μ （2） 对k A 用反幂法求得其模最小的特征值'k λ，则k λ='k λ+k μ。 3、求A 的(谱数)条件数2)(A cond 和行列式detA 。 由矩阵A 为非奇异对称矩阵可得： | | )(min max 2λλ=A cond （3） 其中max λ为按模最大特征值，min λ为按模最小特征值，通过第一问我们求得的λ和s λ可以很容易求得A 的条件数。 在进行反幂法求解时，要对A 进行LU 分解得到。因L 为单位下三角阵，行 列式为1，U 为上三角阵，行列式为主对角线乘积，所以A 的行列式等于U 的行列式，为U 的主对角线的乘积。 有限元上机实验报告 姓名柏小娜 学号0901510401 实验一 一 已知条件 简支梁如图所示，截面为矩形，高度h=200mm ，长度L=1000mm ，厚度t=10mm 。上边承受均布载荷，集度q=1N/mm 2，材料的E=206GPa ，μ=0.29。平面应力模型。 X 方向正应力的弹性力学理论解如下： )534()4 (6222 23-+-=h y h y q y x L h q x σ 二 实验目的和要求 （1）在Ansys 软件中用有限元法探索整个梁上x σ，y σ的分布规律。 （2）计算下边中点正应力x σ的最大值；对单元网格逐步加密，把x σ的计算值与理论解对比，考察有限元解的收敛性。 （3）针对上述力学模型，对比三节点三角形平面单元和4节点四边形平面等参元的求解精度。 三 实验过程概述 （1） 定义文件名 （2） 根据要求建立模型：建立长度为1m ，外径为0.2m ，平行四边行区域 （3） 设置单元类型、属性及厚度，选择材料属性： （4） 离散几何模型，进行网格划分 （5） 施加位移约束 （6） 施加载荷 （7） 提交计算求解及后处理 （8） 分析结果 四 实验内容分析 （1）根据计算得到应力云图，分析本简支梁模型应力分布情况和规律。主要考察x σ和y σ，并分析有限元解与理论解的差异。 由图1看出沿X 方向的应力呈带状分布，大小由中间向上下底面递增，上下底面应力方向相反。由图2看出应力大小是由两侧向中间递增的，得到X 方向 上最大应力就在下部中点，为0.1868 MPa 。根据理论公式求的的最大应力值为0.1895MPa 。由结果可知，有限元解与理论值非常接近。由图3看出Y 的方向应力基本相等，应力主要分布在两侧节点处。 图 1 以矩形单元为有限元模型时计算得出的X 方向应力云图 图 2 以矩形单元为有限元模型时计算得出的底线上各点x 方向应力图 （2）对照理论解，对最大应力点的x σ应力收敛过程进行分析。列出各次计算 应力及其误差的表格，绘制误差－计算次数曲线，并进行分析说明。 答：在下边中点位置最大应力理论值为： MPa h y h y q y x L h q x 1895.0)5 34()4(622223=-+-=σ 计算方法上机实验报告——拉格朗日插值问题 一、方法原理 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 数值分析2016上机实习报告要求 1．应提交一份完整的实习报告。具体要求如下： （1） 报告要排版，美观漂亮；要标明姓名、学号、专业和联系电话； （2） 要有目录，指明题目、程序、计算结果，图标和分析等内容所在位置，作 到信息简明而完全； （3）要有总结，全方位总结机编程计算的心得体会； （4）尽量使报告清晰明了，一般可将计算结果、图表及对比分析放在前面，程序 清单作为附录放在后面，程序中关键部分要有中文说明或标注，指明该部分的功能和作用。 2．程序需完好保存到期末考试后的一个星期，以便老师索取用于验证、询问或质疑部分内容。 3．认真完成实验内容，可以达到既学习计算方法又提高计算能力的目的，还可以切身体会书本内容之精妙所在，期间可以得到很多乐趣。 4．拷贝或抄袭他人结果是不良行为，将视为不合格。 5．请按任课老师要求的时间和载体（电子或纸质）提交给任课老师。 6.提交电子邮件请务必按述格式命名邮件主题和文档名，以方便统计 学号姓名2016_2上机实习报告 如：2015201597张帅哥数值分析2016-2计算报告 数值分析上机试题 1． 编程计算∑∞11=n cn ， 其中c= 4.4942?10307， 给出并观察计算结果，若有问题，分析之。 2． 利用牛顿法求方程ln 2x x -=(x-ln x=2)于区间[2,4]的根，考虑不同初值下牛顿法的收敛情况。 3． 给出1251 2+=x y 在xk=0+0.25k 处的值yk, k=0,1,2,3,4. 请给出由节点 xk 确定的三次样 条插值函数S(x), 使其满足条件：S /(0)=0, S /(1)=-0.074, 分析逼近效果如何？ 4. 松弛因子对SOR 法收敛速度的影响。 用SOR 法求解方程组Ax =b ,其中 3-2-.2-2-3-B ,41141......14114 要求程序中不存系数矩阵A ,分别对不同的阶数取w=1.1, 1.2, ...,1.9进行迭代，记录近似解x (k)达到||x (k)-x (k-1)||<10-6时所用的迭代次数k ，观察松弛因子对收敛速度的影响，并观察当w ≤0或w ≥2会有什么影响？ 5. 用Ru n ge-Kutt a 4阶算法对初值问题y /=-20*y ，y (0)=1按不同步长求解，用于观察稳定区间的 作用，推荐两种步长h=0.1,0.2。 注：此方程的精确解为：y =e -20x数值分析实验2014
数值分析上机实验报告
数值分析上机实习报告
有限元实验报告模板
计算方法上机实习题大作业(实验报告).
数值分析上机报告
有限元实验报告
2011上机实习报告2
计算方法实验报告册
北航数值分析计算实习报告一
ansys实验报告
计算方法上机实验报告——拉格朗日插值问题
数值分析学号姓名2016_2上机实习报告