中国矿业大学（北京）理学院
数值分析实验报告
实验名称 不动点迭代法求方程的近似根 实验时间 2012.3.20
组长签名
龙纯鹏 班级 信息与计算
科学（1）班
学号
11107200110 成绩
组员签名
11107200101
11107200102
11107200103 11107200119
11107200120
一、实验目的,内容 二、相关背景知识介绍 三、代码
四、数值结果 五、计算结果的分析 六、计算中出现的问题，解决方法及体会
一、实验目的、内容
实验目的：熟悉掌握不动点迭代方法的思想方法，并熟悉运用MATLAB 编写相关代码求解方程的近似根； 内容：先确定方程x
e x 5
1=
的一个收敛的有根区间[a,b], 然后用不动点迭代法求x
e x 5
1=
在此有根区间的近似根,初值0x 自己确定,要求根精确到510-,并求迭代次数。 二、相关背景知识介绍 （1）算法原理或计算公式 :
不动点：将方程0)(=x f 写成等价的形式)(x x ?=.若要求*x 满足*)(x f =0,则小*x =
)(x x ?=;反之，若*)(*x x ?=，则满足0*)(=x f ，则称*x 为函数*)(x ?)的一个不动点。
不动点迭代法：求满足)(x f 的零点就等价于求)(x ?的不动点，选择一个初始近似值x0，将其代入x =)(x ?的右端，即可求得:
)(01x x ?=
)(12x x ?=
…….
可以如此反复迭代计算: )(1k k x x ?=+),2,1,0( =k
.如图:
)
(x y ?=
则x =)(x ?称为迭代函数。如果对任何0x 属于[a,b],由)(1k k x x ?=+),2,1,0( =k 得到的序列{k x }有极限，则称迭代方程)(1k k x x ?=+),2,1,0( =k 收敛，且)(*x x ?=为φ(x)的不动点，故称 )(1k k x x ?=+),2,1,0( =k 为不动点迭代法。
迭代法的基本思路是一种主次逼近的方法，其基本思想是将隐式方程0)(=x f 归结为一组显式的计算公式 )(1k k x x ?=+),2,1,0( =k ，也就是说，迭代过程的实质上是一个逐步显式化的过程。
（2）程序设计思路： 先确定方程x e x 51=
的一个收敛的有根区间[a,b],则x e x 5
1)(=?’ ),()(x x x f ?-=设则)(x f 在[a,b]上连续，
由以上条件，若 )()(a a a f ?-= 0≤ )()(b b b f ?-= 0≤
则方程x
e x 5
1=
的一个收敛的根在区间[a,b]上，迭代过程如下：取初值1=x , *
012x x x x O
)(01x x ?= =
051x e 得1151e x = )(12x x ?== )5
1(
e ? ……………… )(1k k x x ?=+),2,1,0( =k . 三、代码（Matlab ）
clear x0=1
e=10^(-5) k=1
x1=sqrt(0.2*e^x0) while (abs(x0-x1)>e) k=k+1 x0=x1
x1=sqrt(0.2*e^x0) end x0
k = 261
x0 = 0.1691 x1 = 0.1690 k = 262
x0 = 0.1690 x1 = 0.1691 四、数值结果
在matlab 运算结果如下表： k 0x
k 0x
k 0x
100 101 102 103 104 105
0.1753 0.1630 0.1750 0.1633 0.1747 0.1636
143 145 154 158 160 178
0.1709 0.1673 0.1676 0.1702 0.1678 0.1689
250 252 253 257 260 266
0.1689 0.1691 0.1690 0.1691 0.1690 0.1691
=*x 0.1691
K=266
五、计算结果的分析 方程x
e x 5
1=
在区间[a,b]的近似根为=*x 0.1691 六、计算中出现的问题，解决方法及体会
1.着重理解不动点迭代法德思想方法，体会其用法，慢慢改进，有着深刻的映像；
2.逐渐熟悉Matlab 的用法，致使解决问题效率上升，掌握学好相关基础知识;
3.对于同样的方程, 迭代函数的构造有关,不同的迭代格式,有不同的结果; 当迭代函数不收敛时,即发散的时候,迭代没有意义。
教 师 评 语
指导教师：
2013 年 3 月 20 日
结构力学（二）上机试验结构力学求解器的使用 上机报告 班级： 姓名： 学号： 日期：
实验三、计算结构的影响线 1．实验任务 （1）作以下图示梁中截面D 的内力D M 、QD F 的影响线。 观览器：D M 的影响线 观览器：QD F 的影响线 D |F=1 3 365
编辑器： 结点,1,0,0 结点,2,3,0 结点,3,6,0 结点,4,12,0 结点,6,6,1 结点,5,17,1 单元,1,2,1,1,0,1,1,1 单元,2,3,1,1,1,1,1,1 单元,3,4,1,1,1,1,1,0 单元,3,6,1,1,0,1,1,0 单元,6,5,1,1,0,1,1,0 结点支承,1,3,0,0,0 结点支承,4,1,0,0 结点支承,5,3,0,0,0 影响线参数,-2,1,1,3 影响线参数,-2,1,1,2 End
作以下图示梁中截面D 的内力D M 、QD F 的影响线。 观览器： D M 的影响线 QD F 的影响线
编辑器： 结点,1,0,0 结点,2,2,0 结点,3,4,0 结点,4,6,0 结点,5,8,0 结点,6,0,1 结点,7,8,1 结点,8,2,1 结点,9,4,1 结点,10,6,1 单元,1,2,1,1,0,1,1,1 单元,2,3,1,1,1,1,1,1 单元,3,4,1,1,1,1,1,1 单元,4,5,1,1,1,1,1,0 单元,1,6,1,1,1,1,1,0 单元,6,8,1,1,0,1,1,0 单元,8,9,1,1,0,1,1,0 单元,9,10,1,1,0,1,1,0 单元,10,7,1,1,0,1,1,0 单元,7,5,1,1,0,1,1,0
插值法数值上机实验报告 实验题目： 利用下列条件做插值逼近，并与R (x) 的图像比较 考虑函数：R x y=1 1+x2 （1）用等距节点X i=?5+i,i=0,1,...,10.给出它的10次Newton插值多项式的图像； π),i=0,1,...,20.给出它的20次Lagrange插值多项式（2）用节点X i=5cos(2i+1 42 的图像； （3）用等距节点X i=?5+i,i=0,1,...,10.给出它的分段线性插值函数的图像；（4）用等距节点X i=?5+i,i=0,1,...,10.给出它的三次自然样条插值函数的图像； （5）用等距节点X i=?5+i,i=0,1,...,10.给出它的分段三次Hermite插值函数的图像； 实验图像结果：
实验结果分析： 1.为了验证Range现象，我还特意做了10次牛顿插值多项式和20次牛顿插值多项式的对比图像，结果如下图（图对称，只截取一半） 可以看出，Range现象在高次时变得更加明显。这也是由于高次多项式在端点处的最值随次数的变大很明显。可以料定高次多项式在两侧端点处剧烈震荡，在更小的间距内急剧上升然后下降，Range现象非常明显。
2.分析实验（2）的结果，我们会惊讶地发现，由于取21个点逼近，原本预料的Range现象会很明显，但这里却和f(x)拟合的很好。（即下图中Lagrange p(x)的图像）。可是上图中取均匀节点的20次牛顿多项式逼近的效果在端点处却很差。料想是由于节点X i=5cos2i+1 42 π ,i=0,1,...,20 取得很好。由书上第五章的 知识，对于函数y=1 1+x ，y 1 2对应的cherbyshev多项式的根恰好为X i= 5cos2i+1 42 π ,i=0,1,...,20 。由于所学限制，未能深入分析。 （3）比较三次样条插值图像和Hermit插值图像对原函数图像的逼近情形。见下图：
结构力学实验报告 班级12土木2班 姓名 学号
实验报告一 实验名称 在求解器中输入平面结构体系 一实验目的 1、了解如何在求解器中输入结构体系 2、学习并掌握计算模型的交互式输入方法； 3、建立任意体系的计算模型并做几何组成分析； 4、计算平面静定结构的内力。 二实验仪器 计算机，软件：结构力学求解器 三实验步骤 图2-4-3 是刚结点的连接示例，其中图2-4-3a 中定义了一个虚拟刚结点和杆端的连接码；各个杆端与虚拟刚结点连接后成为图2-4-3b 的形式，去除虚拟刚结点后的效果为图2-4-3c 所示的刚结点；求解器中显示的是最后的图2-4-3c。图2-4-4 是组合结点的连接示例，同理，无需重复。铰结点是最常见的结点之一，其连接示例在图2-4-5 中给出。这里，共有四种连接方式，都等效于图2-4-5e 中的铰结点，通常采用图2-4-5a 所示方式即可。值得一提的是，如果将三个杆件固定住，图2-4-5b~d 中的虚拟刚结点也随之被固定不动，而图2-4-5a 中的虚拟刚结点仍然存在一个转动自由度，可以绕结点自由转动。这是一种结点转动机构，在求解器中会自动将其排除不计①。结点机构实际上也潜存于经典的结构力学之中，如将一个集中力矩加在铰结点上，便可以理解为加在了结点机构上（犹如加在可自由转动的销钉上），是无意义的。 综上所述，求解器中单元对话框中的“连接方式”是指各杆端与虚拟刚结点的连接方式，而不是杆件之间的连接方式。这样，各杆件通过虚拟刚结点这一中介再和其他杆件间接地连接。这种处理的好处是可以避免结点的重复编码（如本书中矩阵位移法中所介绍的），同时可以方便地构造各种
姓名 实验报告成绩 评语： 指导教师（签名） 年 月 日 说明：指导教师评分后，实验报告交院（系）办公室保存。 实验一 方程求根 一、 实验目的 用各种方法求任意实函数方程0)(=x f 在自变量区间[a ，b]上，或某一点附近的实根。并比较方法的优劣。 二、 实验原理 (1)、二分法 对方程0)(=x f 在[a ，b]内求根。将所给区间二分，在分点 2a b x -=判断是否0)(=x f ；若是，则有根2a b x -=。否则，继续判断是否0)()(
+)(0x f 0))(('0=-x x x f 设0)('0≠x f ,则=x -0x )(') (00x f x f 。取x 作为原方程新的近似根1x ，然后将1x 作为0x 代入上式。迭代公式为：=+1 k x -0x )(')(k k x f x f 。 三、 实验设备：MATLAB 软件 四、 结果预测 （1）11x = （2）5x = （3）2x =0,09052 五、 实验内容 （1）、在区间[0,1]上用二分法求方程0210=-+x e x 的近似根，要求误差不超 过3105.0-?。 （2）、取初值00=x ，用迭代公式=+1 k x -0x )(') (k k x f x f ，求方程0210=-+x e x 的近似根。要求误差不超过3105.0-?。 （3）、取初值00=x ，用牛顿迭代法求方程0210=-+x e x 的近似根。要求误差 不超过3105.0-?。 六、 实验步骤与实验程序 （1） 二分法 第一步：在MATLAB 软件，建立一个实现二分法的MATLAB 函数文件如下： function x=agui_bisect(fname,a,b,e) %fname 为函数名，a,b 为区间端点，e 为精度 fa=feval(fname,a); %把a 端点代入函数，求fa fb=feval(fname,b); %把b 端点代入函数，求fb if fa*fb>0 error('两端函数值为同号'); end
结构力学求解器学习报告 一、实习目的 结构力学上机实习使训练学生使用计算机进行结构计算的重要环节。通过实习，学生可以掌握如何使用计算机程序进行杆系结构的分析计算，进一步掌握结构力学课程的基本理论和基本概念。在此基础上，通过阅读有关程序设计框图，编写、调试结构力学程序，学生进一步提高运用计算机进行计算的能力，为后续课程的学习、毕业设计及今后工作中使用计算机进行计算打下良好的基础。 二、实习时间 大三上学期第19周星期一至星期五。 三、实习内容 本次实习以自学为主，学习如何使用结构力学求解器进行结构力学问题的求解，包括：二维平面结构（体系）的几何组成、静定、超静定、位移、内力、影响线、自由振动、弹性稳定、极限荷载等。对所有这些问题，求解器全部采用精确算法给出精确解答。 四、心得体会 第一天上机时，张老师对结构力学求解器的使用方法进行了简单的介绍，然后就是学生自己自学的时间了。每个学生都有自己对应的题目要完成，在完成这些题目的同时，我也逐渐对结构力学求解器的运用更加自如。 从刚开始的生疏到最后的熟练运用，我遇到了不少问题：①第一次使用在有些问题上拿不定注意，例如，在材料性质那一栏，我不知
道是EA和EI的取值②第一次接触这个软件，在使用过程中不知道该如何下手，题目条件的输入顺序也很模糊。③经常会忘记添加荷载的单位，导致计算结果出现问题。④对于有些命令不能很明确的知道其用法，致使在使用时经常出错。在面对这些问题时，我一般都会向同学和老师寻求帮助，直到最终将问题解决。 通过这几天的上机实习，不仅让我进一步掌握了结构力学的知识，同时，还使我对结构力学求解器有了更深入的了解： 1. 结构力学求解器首先是一个计算求解的强有效的工具。对于任意平面的结构，只要将参数输进求解器，就可以得到变形图和内力图，甚至还可以求得临界荷载等问题。 2.即便是结构力学的初学者，只要会用求解器，也可以用求解器来方便地求解许多结构的各类问题，以增强对结构受力特性的直观感受和切实体验。 3.书本中的方法并非所有类型的问题都可以解决，例如，不规则分布的荷载以及超静定结构用传统方法比较困难，但用求解器就较为简单。而且，用求解器求解问题时可以不忽略轴向变形等书本中忽略的条件，与实际更加相符。 4.求解器可以用静态图形显示结构简图、变形图、内力图，还可以用动画显示机构模态、振型等动态图形。利用复制到剪贴板的功能，可以将结构简图、变形图、内力图以点阵图或矢量图的形式粘贴到word文档中，并可以方便地进行再编辑。
实验报告一 平面刚架内力计算程序APF 实验目的：(1)分析构件刚度与外界温度对结构位移的影响，如各杆刚度改变对内力分布的影响、温度因数对内力分布的影响。 (2)观察并分析刚架在静力荷载及温度作用下的内力和变形规律，包括刚度的变化，结构形式的改变，荷载的作用位置变化等因素对内力及变形的影响。对结构静力分析的矩阵位移法的计算机应用有直观的了解 (3)掌握杆系结构计算的《结构力学求解器》的使用方法。通过实验加深对静定、超静定结构特性的认识。 实验设计1： 计算图示刚架当梁柱刚度12I I 分别为15、11、15、1 10时结构的内力和位移，由此分析当刚架在水平荷 载作用下横梁的水平位移与刚架梁柱 比（1 2I I ）之间的关系。(计算时忽略轴向变形)。 数据文件： (1)变量定义,EI1=1,EI2=0.2(1,5,10) 结点,1,0,0 结点,2,0,4 结点,3,6,4 结点,4,6,0 单元,1,2,1,1,1,1,1,1 单元,2,3,1,1,1,1,1,1 单元,3,4,1,1,1,1,1,1 结点支承,1,6,0,0,0,0 结点支承,4,6,0,0,0,0 结点荷载,2,1,100,0 单元材料性质,1,1,-1,EI1,0,0,-1 单元材料性质,2,2,-1,EI2,0,0,-1 单元材料性质,3,3,-1,EI1,0,0,-1 (2)变量定义,EI1=5(1,0.2,0.1),EI2=1 结点,1,0,0 结点,2,0,4 结点,3,6,4 结点,4,6,0 单元,1,2,1,1,1,1,1,1 单元,2,3,1,1,1,1,1,1 单元,3,4,1,1,1,1,1,1 结点支承,1,6,0,0,0,0 结点支承,4,6,0,0,0,0 结点荷载,2,1,100,0 单元材料性质,1,1,-1,EI1,0,0,-1 单元材料性质,2,2,-1,EI2,0,0,-1 单元材料性质,3,3,-1,EI1,0,0,-1 主要计算结果： 位移：
数值分析实验报告记录
————————————————————————————————作者：————————————————————————————————日期：
数值分析实验报告 （第二章） 实验题目： 分别用二分法、牛顿迭代法、割线法、史蒂芬森迭代法求方程 的根，观察不同初始值下的收敛性，并给出结论。 问题分析： 题目有以下几点要求： 1.不同的迭代法计算根，并比较收敛性。 2.选定不同的初始值，比较收敛性。 实验原理： 各个迭代法简述 二分法：取有根区间的重点，确定新的有根区间的区间长度仅为区间长度的一版。对压缩了的有根区间重复以上过程，又得到新的有根区间，其区间长度为的一半，如此反复，……，可得一系列有根区间，区间收敛到一个点即为根。 牛顿迭代法：不动点迭代法的一种特例，具有局部二次收敛的特性。迭代格式为 割线法：是牛顿法的改进，具有超线性收敛的特性，收敛阶为1.618. 迭代格式为 史蒂芬森迭代法：采用不动点迭代进行预估校正。至少是平方收敛的。迭代格式为 这里可采用牛顿迭代法的迭代函数。 实验内容：
1.写出该问题的函数代码如下： function py= f(x) syms k; y=(k^2+1)*(k-1)^5; yy=diff(y,k); py(1)=subs(y,k,x); py(2)=subs(yy,k,x); end 2.分别写出各个迭代法的迭代函数代码如下： 二分法： function y=dichotomie(a,b,e) i=2; m(1)=a; while abs(a-b)>e t=(a+b)/2; s1=f(a); s2=f(b); s3=f(t); if s1(1)*s3(1)<=0 b=t; else a=t; end m(i)=t; i=i+1; end y=[t,i+1,m]; end 牛顿迭代法： function y=NewtonIterative(x,e) i=2; en=2*e;m(1)=x; while abs(en)>=e s=f(x); t=x-s(1)/s(2); en=t-x; x=t; m(i)=t; i=i+1; end y=[x,i+1,m]; end 牛顿割线法： function y=Secant(x1,x2,e) i=3; m(1)=x1,m(2)=x2; while abs(x2-x1)>=e s1=f(x1); s2=f(x2); t=x2-(x2-x1)*s2(1)/(s2(1)-s1( 1)); x1=x2; x2=t; m(i)=t; i=i+1; end
结构力学 桁架结构受力性能实验报告 学号：1153377 姓名：周璇 专业：土木工程 实验时间：2016年05月04日周三，中午12:30-13:30 实验指导教师：陈涛 理论课任课教师：陈涛
一、实验目的 （1）参加并完成规定的实验项目内容，理解和掌握结构的实验方法和实验结果，通过 实践掌握试件的设计、实验结果整理的方法。 （2）进行静定、超静定结构受力的测定和影响线的绘制。 二、结构实验 （一）空间桁架受力性能概述 桁架在受结点荷载时，两边支座处产生反力，桁架中各杆件产生轴力，如图1.1为在抛物线桁架结点分别加载时结构示意图。用Q235钢材，桁架跨度6?260=1560mm ，最大高度260mm 。杆件之间为铰接相连。杆件直径为8mm 。 图1.1 （二）实验装置 图1.2为框架结构侧向受力实验采用的加载装置，25kg 挂钩和25kg 砝码。采用单结点集中力加载，由砝码、挂钩施加拉力，应变片测算待测杆件应变。结构尺寸如图1.2所示。 图1.2 （三）加载方式 简单多次加载，将挂钩和砝码依次施加在各个结点，待应变片返回数据稳定后，进行采集。采集结束后卸下重物，等待应变片数值降回初始值后再向下一节点施加荷载，重复采集操作。 （四）量测内容 需要量测桁架待测杆件的应变值在前后四对桁架杆布置单向应变片，具体布置位置如图 1.2 所示，即加粗杆件上黏贴应变片。 三、实验原理 对桁架上的5个位置分别施加相同荷载，记录不同条件下各杆件的应变值。 由公式 2 4 F A E d A σσεπ? ?=? =???=?
可以得到 24 d E F πε = 其中： F ——杆件轴力 E ——Q235钢弹性模量 d ——杆件直径 ε ——杆件应变值 σ ——杆件应力 A ——杆件横截面积 因而可以求得各杆件轴力，进而得到不同杆件的轴力影响线。 四、实验步骤 （1）将载荷挂在加载位置1，待应变片返回数据稳定后，采集相应应变数据。 （2）待应变片数值降回初始值后，重复（1）中操作，将荷载分别挂在加载位置2,3,4,5，分别采集记录各自对应的各杆件应变数据。 五、实验结果与整理 将对应位置杆件应变值取平均值，得到所示一榀桁架四根杆件的应变值如表2.2所示。
姓名实验报告成绩 评语： 指导教师（签名） 年月日
说明：指导教师评分后，实验报告交院（系）办公室保存。 实验一 方程求根 一、 实验目的 用各种方法求任意实函数方程0)(=x f 在自变量区间[a ，b]上，或某一点附近的实根。并比较方法的优劣。 二、 实验原理 (1)、二分法 对方程0)(=x f 在[a ，b]内求根。将所给区间二分，在分点2a b x -= 判 断是否0)(=x f ；若是，则有根 2a b x -= 。否则，继续判断是否0)()(
（1）11x =0.09033 （2）5x =0.09052 （3）2x =0,09052 五、 实验内容 （1）、在区间[0,1]上用二分法求方程0210=-+x e x 的近似根，要求误差不 超过 3 105.0-?。 （2）、取初值00=x ，用迭代公式=+1k x -0x ) (') (k k x f x f ，求方程0210=-+x e x 的 近似根。要求误差不超过 3 105.0-?。 （3）、取初值00=x ，用牛顿迭代法求方程0210=-+x e x 的近似根。要求误 差不超过 3 105.0-?。 六、 实验步骤与实验程序 （1） 二分法 第一步：在MATLAB 7.0软件，建立一个实现二分法的MATLAB 函数文件agui_bisect.m 如下： function x=agui_bisect(fname,a,b,e) %fname 为函数名，a,b 为区间端点，e 为精度 fa=feval(fname,a); %把a 端点代入函数，求fa fb=feval(fname,b); %把b 端点代入函数，求fb if fa*fb>0 error('两端函数值为同号'); end %如果fa*fb>0，则输出两端函数值为同号 k=0 x=(a+b)/2 while(b-a)>(2*e) %循环条件的限制
中国矿业大学力学与建筑工程学院 2013~2014学年度第二学期 《结构力学A1》上机实验报告 学号 班级 姓名 2014年5月26日
一、单跨超静定梁计算(50分) 1. 计算并绘制下面单跨超静定梁的弯矩图和剪力图。（20分） q =12N/m q =8N/m q =8N/m q=?8m 1 2 3 2. 如果按照梁跨中弯矩相等的原则，将梁上的荷载换算成均布荷载，则均布荷载应为多少？（10分） 2m 1m 1m 1m 1m 1m q=? 8m 3. 如果按照梁端部弯矩相等的原则，将梁上的荷载换算成均布荷载，则均布荷载应为多少？（10分） 4. 如果按照梁端部剪力相等的原则，将梁上的荷载换算成均布荷载，则均布荷载应为多少？（10分） 二、超静定刚架计算（50分） 1.刚架各杆EI 如图所示,计算刚架的弯矩图，剪力图和轴力图。（30分）
2. 若EI=106 （Nm 2 ），计算刚架一层梁和二层梁的水平位移。（20分）
弯矩图： y x 12345678 ( 1 )( 2 )( 3 )( 4 )( 5 )( 6 )( 7 ) -40.96 -16.29 3.04 19.04 25.04 19.04 3.04 -16.29 -40.96 剪力图： y x 12345678 ( 1 )( 2 )( 3 )( 4 )( 5 )( 6 )( 7 ) 26.00 22.00 18.00 12.00 -12.00 -18.00 -22.00 -26.00
解：跨中弯矩M1=25.04Nm（下部受拉）均布荷载q作用在梁上时，跨中弯矩为 M2=1/24*q*(l^2)(下部受拉) ∵M1=M2, ∴q=9.39N/m 如图所示： y x 12 ( 1 ) -50.08-50.08
. 牛顿插值法一、实验目的：学会牛顿插值法，并应用算法于实际问题。 x?x)f(二、实验内容：给定函数，已知： 4832401.2)?.?1449138f(2.f.f(20)?1.414214(2.1) 549193.)?1f(2.4516575(f2.3)?1. 三、实验要求：以此作为函数2.15插值多项式在处的值，用牛顿插值法求4 次Newton（ 1）2.15?N(2.15)。在MATLAB中用内部函数ezplot绘制出的近似值4次Newton插值多项式的函数图形。 （2）在MATLAB中用内部函数ezplot可直接绘制出以上函数的图形，并与作出的4次Newton插值多项式的图形进行比较。 四、实验过程： 1、编写主函数。打开Editor编辑器，输入Newton插值法主程序语句： function [y,L]=newdscg(X,Y,x) n=length(X); z=x; A=zeros(n,n);A(:,1)=Y';s=0.0; p=1.0; for j=2:n for i=j:n A(i,j)=(A(i,j-1)- A(i-1,j-1))/(X(i)-X(i-j+1)); end end C=A(n,n); for k=(n-1):-1:1 C=conv(C,poly(X(k))); d=length(C);C(d)=C(d)+A(k,k); end y(k)= polyval(C, z); L(k,:)=poly2sym(C); 0 / 3 . %%%%%%%%%%%%%%%%%% t=[2,2.1,2.2,2.3,2.4]; fx=sqrt(t); wucha=fx-Y; 以文件名newdscg.m保存。 2、运行程序。 （1）在MATLAB命令窗口输入： >> X=[2,2.1,2.2,2.3,2.4]; Y =[1.414214,1.449138,1.483240,1.516575,1.549193]; x=2.15;[y,P]=newdscg(X,Y,x) 回车得到：
实验报告一 平面刚架内力计算程序APF 日期： 2013.4.19 实验地点： 综合楼503 实验目的： 1、通过实验加深对静定、超静定结构特性的认识。如各杆刚度改变对内力分布的影响、温度和沉陷变形因数的影响等。 2、观察并分析刚架在静力荷载及温度作用下的内力和变形规律，包括刚度的变化，结构形式的改变，荷载的作用位置变化等因素对内力及变形的影响。对结构静力分析的矩阵位移法的计算机应用有直观的了解。 3、掌握杆系结构计算的《求解器》的使用方法。 实验设计1： 别为15 、11、15、110 时结构的内力和位移，由此 分析当刚架在水平荷载作用下横梁的水平位移与刚架梁柱比（1 2I I ）之间的关系。(计算时忽略轴 向变形)。 一、 数据文件： （1）TITLE, 实验一 变量定义,EI1=1 变量定义,EI2=0.2（1, 5, 10） 结点,1,0,0 结点,2,0,4 结点,3,6,0 结点,4,6,4 单元,1,2,1,1,1,1,1,1 单元,3,4,1,1,1,1,1,1 单元,2,4,1,1,1,1,1,1 结点支承,1,6,0,0,0,0 结点支承,3,6,0,0,0,0 结点荷载,2,1,100,0 单元材料性质,1,2,-1,EI1,0,0,-1 单元材料性质,3,3,-1,EI2,0,0,-1 END
二、主要计算结果： 位移： (2)令I2=1时，I1=5,1,0.2,0.1 弯矩： (1) 令I1=1时，I2=0.2,1,5,10 ①梁柱刚度比I2：I1为1：5时的刚架弯矩图如下②梁柱刚度比I2：I1为1：1时的刚架弯矩图如下
③梁柱刚度比I2：I1为5：1时的刚架弯矩图如下④梁柱刚度比I2：I1为10：1时的刚架弯矩图如下
迭代法实验报告 一． 实验目的：掌握迭代方法的用处 二． 实验环境：Cfree5.0 三． 实验时间：2013年6月20日 四． 实验地点：电子信息楼1201教室 五． 实验内容：运用编程实现迭代方法可以更好的解线性方程组，得到线性方程的解。 六． 实验理论依据： 高斯-赛德尔（Gauss-Seidel ）迭代公式 我们注意到在雅可比迭代法中并没有对新算出的分量11k x +，12k x +， ， 11k i x +-进行充分利用．不妨设想，在迭代收敛的条件下，我们把 (1)()()()11211331111(1)()()()22112332222(1)()()()1122,111()1(1(k k k k n n k k k k n n k k k k n n n n n n nn x a x a x a x b a x a x a x a x b a x a x a x a x b a +++--?=---+???=---+?????=---+?? 式中第一个方程算出的11k x +立即投入到第二个方程中，代替()1k x 进行计算，当12 k x +算出后代替()2k x 马上投入到第三个方程中计算，依次进行下去，这样也许会得到 更好的收敛效果．根据这种思路建立的一种新的迭代格式，我们称为高斯-赛德尔（Gauss-Seidel ）迭代公式， 高斯=赛德尔迭代法的分量形式：
(1)()()()11211331111(1)(1)()()22112332222(1)(1)(1)(1)1122,111()1(1(k k k k n n k k k k n n k k k k n n n n n n nn x a x a x a x b a x a x a x a x b a x a x a x a x b a +++++++--?=---+???=---+?????=---+?? 高斯-赛德尔迭代法的矩阵形式： (1)(),(0,1,2,)k k x Bx f k +=+= 其中 1()B D L U -=- ，1()f D L b -=- B 称为高斯-赛德尔迭代矩阵，f 称为高斯-赛德尔迭代常量.. 七． 运行代码如下： #include"stdio.h" #include"math.h" int main() { bool pan1=true; int n,n1,n2=0,k=0; double num[100][100],L[100][100],U[100][100],x[100],y[100],num1=0,b[100],D[100][100],x1[200][200],x2[200][200]; printf("\n"); printf("*******************************高斯迭代法解如下********************************"); printf("输入要输入矩阵的阶数为(按Enter 输入矩阵数字)：");//
实验报告：牛顿差值多项式&三次样条 问题：在区间[-1,1]上分别取n=10、20用两组等距节点对龙格函数2 1()25f x x 作多项式插值及三次样条插值，对每个n 值，分别画出插值函数及()f x 的图形。 实验目的：通过编程实现牛顿插值方法和三次样条方法，加深对多项式插值的理解。应用所编程序解决实际算例。 实验要求： 1． 认真分析问题，深刻理解相关理论知识并能熟练应用； 2． 编写相关程序并进行实验； 3． 调试程序，得到最终结果； 4． 分析解释实验结果； 5． 按照要求完成实验报告。 实验原理： 详见《数值分析 第5版》第二章相关内容。 实验内容： （1）牛顿插值多项式 1.1 当n=10时： 在Matlab 下编写代码完成计算和画图。结果如下： 代码： clear all clc x1=-1:0.2:1; y1=1./(1+25.*x1.^2); n=length(x1); f=y1(:); for j=2:n for i=n:-1:j f(i)=(f(i)-f(i-1))/(x1(i)-x1(i-j+1)); end end syms F x p ; F(1)=1;p(1)=y1(1); for i=2:n F(i)=F(i-1)*(x-x1(i-1)); p(i)=f(i)*F(i);
end syms P P=sum(p); P10=vpa(expand(P),5); x0=-1:0.001:1; y0=subs(P,x,x0); y2=subs(1/(1+25*x^2),x,x0); plot(x0,y0,x0,y2) grid on xlabel('x') ylabel('y') P10即我们所求的牛顿插值多项式，其结果为：P10（x）=-220.94*x^10+494.91*x^8-9.5065e-14*x^7-381.43*x^6-8.504e-14*x^5+123.36*x^4+2.0202e-1 4*x^3-16.855*x^2-6.6594e-16*x+1.0 并且这里也能得到该牛顿插值多项式的在[-1，1]上的图形，并和原函数进行对比（见Fig.1）。 Fig.1 牛顿插值多项式（n=10）函数和原函数图形 从图形中我们可以明显的观察出插值函数在两端点处发生了剧烈的波动，产生了极大的误差，即龙格现象，当n=20时，这一现象将更加明显。 1.2 当n=20时： 对n=10的代码进行修改就可以得到n=20时的代码。将“x1=-1:0.2:1;”改为“x1=-1:0.1:1;”即可。运行程序，我们得到n=20时的牛顿插值多项式，结果为：P20(x)= 260188.0*x^20 - 1.0121e6*x^18 + 2.6193e-12*x^17 + 1.6392e6*x^16 + 2.248e-11*x^15 - 1.4429e6*x^14 - 4.6331e-11*x^13 + 757299.0*x^12 + 1.7687e-11*x^11 - 245255.0*x^10 + 2.1019e-11*x^9 + 49318.0*x^8 + 3.5903e-12*x^7 - 6119.2*x^6 - 1.5935e-12*x^5 + 470.85*x^4 + 1.3597e-14*x^3 - 24.143*x^2 - 1.738e-14*x + 1.0 同样的，这里得到了该牛顿插值多项式的在[-1，1]上的图形，并和原函数进行对比（见Fig.2）。
用牛顿迭代法求非线性方程的根 一、 实验题目 求方程()013=--=x x x f 在5.1附近的根。 二、 实验引言 （1）实验目的 1. 用牛顿迭代法求解方程的根 2. 了解迭代法的原理 3. 改进和修缮迭代法 （2）实验意义 牛顿迭代法就是众多解非线性方程迭代法中比较普遍的一种，求解方便实用。 三、 算法设计 （1）基本原理 给定初始值0x ，ε为根的容许误差，η为()x f 的容许误差，N 为迭代次数的容许值。 1．如果()0='x f 或迭带次数大于N ，则算法失败，结束；否则执行2. 2．计算()() 0001x f x f x x '-=. 3．若ε<-21x x 或()η<1x f ，则输出1x ，程序结束；否则执行4. 4．令10x x =，转向1. （2）流程图
四、程序设计program nndd01 implicit none real,parameter::e=0.005 real,parameter::n=9 real::x1 real::x0=1.5 integer::k real,external::f,y do k=1,9 if (y(x0)==0) then write(*,*)"失败" else x1=x0-f(x0)/y(x0) if (abs(x1-x0) else x0=x1 end if end if end do end function f(x) implicit none real::f real::x f=x*x*x-x-1 return end function function y(x) implicit none real::y real::x y=3*x*x-1 return end function 五、求解结果 3 1.324718 4 1.324718 5 1.324718 6 1.324718 7 1.324718 8 1.324718 9 1.324718 六、算法评价及讨论 1.在求解在1.5处附近的根，不难发现在输入区间左端值为1时 需要迭代6次，而输入区间左端值为1.5时，却只要4次。初 实验报告内容： 一：不动点迭代法解方程 二：牛顿插值法的MATLAB实现 完成日期：2012年6月21日星期四 数学实验报告一 日期：2012-6-21 所以，确定初值为x0=1 二：不断迭代 算法： 第一步：将f(x0)赋值给x1 第二步：确定x1-x0的绝对值大小，若小于给定的误差值，则将x1当做方程的解，否则回到第一步 编写计算机程序： clear f=inline('0.5*sin(x)+0.4'); x0=1; x1=f(x0); k=1; while abs(x1-x0)>=1.0e-6 x0=x1; x1=f(x0); k=k+1; fprintf('k=%.0f,x0=%.9f,x1=%.9f\n',k,x0,x1) end 显示结果如下： k=2,x0=0.820735492,x1=0.765823700 k=3,x0=0.765823700,x1=0.746565483 k=4,x0=0.746565483,x1=0.739560873 k=6,x0=0.736981783,x1=0.736027993 k=7,x0=0.736027993,x1=0.735674699 k=8,x0=0.735674699,x1=0.735543758 k=9,x0=0.735543758,x1=0.735495216 k=10,x0=0.735495216,x1=0.735477220 k=11,x0=0.735477220,x1=0.735470548 k=12,x0=0.735470548,x1=0.735468074 k=13,x0=0.735468074,x1=0.735467157 >>。。。 以下是程序运行截图： 结构力学上机实验报告 专业建筑工程 班级一班 学号xxx 姓名xx 20 年月日 一、用求解器进行平面体系几何构造分析 （桁架或组合结构） 报告中应包括以下内容： 求解过程 命令文档 分析结果 命令文档： 结点,1,0,0 结点,2,0,2.4 结点,3,1,1 结点,4,2,2.4 结点,5,2.8,0 结点,6,2.8,1 结点,7,3.6,2.4 结点,8,4.6,1 结点,9,5.6,2.4 单元,1,2,1,1,0,1,1,0 单元,2,3,1,1,0,1,1,0 单元,3,4,1,1,0,1,1,0 单元,2,4,1,1,0,1,1,0 单元,1,5,1,1,0,1,1,0 单元,5,6,1,1,0,1,1,0 单元,6,4,1,1,0,1,1,0 单元,1,3,1,1,0,1,1,0 单元,6,7,1,1,0,1,1,0 单元,7,9,1,1,0,1,1,0 单元,4,7,1,1,0,1,1,0 单元,7,8,1,1,0,1,1,0 单元,8,9,1,1,0,1,1,0 结点,10,5.6,0 单元,9,10,1,1,0,1,1,0 单元,10,8,1,1,0,1,1,0 单元,5,10,1,1,0,1,1,0 结点支承,10,1,0,0 结点支承,2,1,-90,0 结点支承,1,2,-90,0,0 分析结果： 二、用求解器确定截面单杆 插图 报告中应包括以下内容： 求解过程 命令文档： 结点,1,0,0 结点,2,0,6 结点,3,5,6 结点,4,10,6 结点,5,10,0 单元,1,2,1,1,1,1,1,1 单元,2,3,1,1,1,1,1,0 单元,3,4,1,1,0,1,1,1 单元,4,5,1,1,1,1,1,1 结点支承,1,2,-90,0,0 结点支承,5,2,0,0,0 单元荷载,3,3,1,0,1,90 牛顿迭代法实验报告 1．功能 本程序采用牛顿法，求实系数高次代数方程 f(x)=a0x n+a1x n-1+…+a n-1x+a n=０（a n≠０）（1） 的在初始值x0附近的一个根。 2.使用说明 （1）函数语句 Y=NEWTON_1(A,N,X0,NN,EPS1) 调用M文件newton_1.m。 （2）参数说明 A n+1元素的一维实数组，输入参数，按升幂存放方程系数。 N整变量，输入参数，方程阶数。 X0 实变量，输入参数，初始迭代值。 NN整变量，输入参数，允许的最大迭代次数。 EPS1实变量，输入参数，控制根的精度。 3．方法简介 解非线性议程f(x)=0的牛顿法是把非线性方程线性化的一种近似方法。把f(x)在x0点附近展开成泰勒级数 f(x)=f(x0)+(x－x0)fˊ(x0)+(x－x0)2 !2) (0x f'' +… 取其线性部分，作为非线性方程f(x)=0的近似方程，则有 f(x0)+fˊ(x0)(x－x0)=0 设fˊ(x0)≠0则其解为 x1=x0－f(x0)/fˊ(x0) 再把f(x)在x1附近展开成泰勒级数，也取其线性部分作f(x)=0的近似方程。若f(x1)≠0，则得 x2=x1－f(x1)/fˊ(x1) 这样，得到牛顿法的一个迭代序列 x n+1=x n－f(x n)/fˊ(x n) 4．newton_1.m程序 function y=newton_1(a,n,x0,nn,eps1) x(1)=x0； b=1； i=1； while(abs(b)>eps1*x(i)) i=i+1； x(i)=x(i-1)-n_f(a,n,x(i-1))/n_df(a,n,x(i-1))； b=x(i)-x(i-1)； if(i>nn)error(ˊnn is fullˊ)； return； end end y=x(i)； i 5．程序附注 （1）程序中调用n_f.m和n_df.m文件。n_f.m是待求根的实数代数方程的函数,n_df.m 是方程一阶导数的函数。由使用者自己编写。 （2）牛顿迭代法的收敛速度：如果f(x)在零点附近存在连续的二阶微商，ξ是f(x)的一个重零点，且初始值x0充分接近于ξ，那么牛顿迭代是收敛的，其收敛速度是二阶的，即平方收敛速度。 6．例题 用牛顿法求下面方程的根 f(x)=x3+2x2+10x-20 7．运行结果 >>a=[1,2,10,-20] ； >>n=3； >>x0=1； >>nn=1000； >>eps1=1e-8； >>y=newton_1(a,n,x0,nn,eps1) 《数值分析》课程实验报告 用拉格朗日和牛顿插值法求解函数值 算法名称用拉格朗日和牛顿插值法求函数值 学科专业xxxxx 作者姓名xxxx 作者学号xxxxx 作者班级xxxxxx xxx大学 二〇一五年十二月 《数值分析》课程实验报告 得到的近似值为。 拉格朗日插值模型简单，结构紧凑，是经典的插值法。但是由于拉格朗日的插值多项式和每个节点都有关，当改变节点个数时，需要重新计算。且当增大插值阶数时容易出现龙格现象。 2.牛顿插值法 在命令窗口输入： x=[ ]; y=[ ]; xt=; [yt,N]=NewtInterp(x,y,xt) z=::2; yz=subs(N,'t',z); figure; plot(z,sqrt(z),'--r',z,yz,'-b') hold on plot(x,y,'marker','+') hold on plot(xt,yt,'marker','o') h=legend('$\sqrt{x}$','牛顿','$(x_k,y_k)$','$x=$'); set(h,'Interpreter','latex') xlabel('x') ylabel('y') 得到结果及图像如下： yt = N = - *t^4 + *t^3 - *t^2 + *t + 得到√的近似值为，插值函数为 N =- *t^4 + *t^3 - *t^2 + *t + ， 其计算精度是相当高的。 Lagrange插值法和Newton插值法解决实际问题中关于只提供复杂的离散数据的函数求值问题，通过将所考察的函数简单化，构造关于离散数据实际函数f(x)的近似函数P(x)，从而可以计算未知点出的函数值，是插值法的基本思路。 实际上Lagrange插值法和Newton插值法是同一种方法的两种变形，其构造拟合函数的思路是相同的，而实验中两个实际问题用两种算法计算出结果是相同的。matlab(迭代法-牛顿插值)Word版
结构力学求解器实验报告
牛顿迭代法的实验报告
数值分析课程实验报告-拉格朗日和牛顿插值法
相关主题
文本预览