材料分析实验报告

材料专业实验报告

题目：扫描电镜(SEM)物相分析实验学院：先进材料与纳米科技学院专业：材料物理与化学

姓名：

学号：1514122986

2016年6月30日

扫描电镜(SEM)物相分析实验

一．实验目的

1.了解扫描电镜的基本结构与原理

2.掌握扫描电镜样品的准备与制备方法

3.掌握扫描电镜的基本操作并上机操作拍摄二次电子像

4.了解扫描电镜图片的分析与描述方法

二．实验原理

1.扫描电镜的工作原理

扫描电镜（SEM）是用聚焦电子束在试样表面逐点扫描成像。试样为块状或粉末颗粒，成像信号可以是二次电子、背散射电子或吸收电子。其中二次电子是最主要的成像信号。由电子枪发射的电子，以其交叉斑作为电子源，经二级聚光镜及物镜的缩小形成具有一定能量、一定束流强度和束斑直径的微细电子束，在扫描线圈驱动下，于试样表面按一定时间、空间顺序作栅网式扫描。聚焦电子束与试样相互作用，产生二次电子发射以及背散射电子等物理信号，二次电子发射量随试样表面形貌而变化。二次电子信号被探测器收集转换成电讯号，经视频放大后输入到显像管栅极，调制与入射电子束同步扫描的显像管亮度，得到反映试样表面形貌的二次电子像。

本次实验中主要通过观察背散射电子像及二次电子像对样品进行分析表征。

1）背散射电子

背散射电子是指被固体样品原子反射回来的一部分入射电子，其中包括弹性背反射电子和非弹性背反射电子。弹性背反射电子是指被样品中原子和反弹回来的，散射角大于90度的那些入射电子，其能量基本上没有变化（能量为数千到数万电子伏）。非弹性背反射电子是入射电子和核外电子撞击后产生非弹性散射，不仅能量变化，而且方向也发生变化。非弹性背反射电子的能量范围很宽，从数十电子伏到数千电子伏。背反射电子的产生范围在100nm-1mm深度。背反射电子产额和二次电子产额与原子序数的关系背反射电子束成像分辨率一般为50-200nm（与电子束斑直径相当）。背反射电子的产额随原子序数的增加而增加，所以，利用背反射电子作为成像信号不仅能分析形貌特征，也可以用来显示原子序数衬

度，定性进行成分分析。

2）二次电子像

二次电子是指背入射电子轰击出来的核外电子。由于原子核和外层价电子间的结合能很小，当原子的核外电子从入射电子获得了大于相应的结合能的能量后，可脱离原子成为自由电子。如果这种散射过程发生在比较接近样品表层处，那些能量大于材料逸出功的自由电子可从样品表面逸出，变成真空中的自由电子，即二次电子。二次电子来自表面5-10nm的区域，能量为0-50eV。它对试样表面状态非常敏感，能有效地显示试样表面的微观形貌。由于它发自试样表层，入射电子还没有被多次反射，因此产生二次电子的面积与入射电子的照射面积没有多大区别，所以二次电子的分辨率较高，一般可达到5-10nm。扫描电镜的分辨率一般就是二次电子分辨率。二次电子产额随原子序数的变化不大，它主要取决与表面形貌。

2.扫描电镜的构造

扫描电镜由下列五部分组成，如图1（a）所示。

（a）（b）

图1 扫描电子显微镜构造示意图

各部分主要作用简介如下：

1)电子光学系统

它由电子枪、电磁透镜、光阑、样品室等部件组成，如图1（b）所示。为了获得较高

的信号强度和扫描像，由电子枪发射的扫描电子束应具有较高的亮度和尽可能小的束斑直径。常用的电子枪有三种形式：普通热阴极三极电子枪、六硼化镧阴极电子枪和场发射电子枪，其性能如表1所示。

表1 几种类型电子枪性能

普通热阴极三极电子枪104～105 20～50 ≈50 10－2

六硼化镧阴极电子枪105～106 1～10 ≈500 10－4

场发射电子枪107～108 0.01～0.1 ≈5000 10-7～10－8前两种属于热发射电子枪，后一种则属于冷发射电子枪，也叫场发射电子枪。

由表可以看出场发射电子枪的亮度最高、电子源直径最小，是高分辨本领扫描电镜的理想电子源。

电磁透镜的功能是把电子枪的束斑逐级聚焦缩小，因照射到样品上的电子束斑越小，其分辨率就越高。扫描电镜通常有三个磁透镜，前两个是强透镜，缩小束斑，第三个透镜是弱透镜，焦距长，便于在样品室和聚光镜之间装入各种信号探测器。为了降低电子束的发散程度，每级磁透镜都装有光阑；为了消除像散，装有消像散器。

样品室中有样品台和信号探测器，样品台还能使样品做平移运动。

2）扫描系统

扫描系统的作用是提供入射电子束在样品表面上以及阴极射线管电子束在荧光屏上的同步扫描信号。

3）信号检测、放大系统

样品在入射电子作用下会产生各种物理信号、有二次电子、背散射电子、特征X射线、阴极荧光和透射电子。不同的物理信号要用不同类型的检测系统。它大致可分为三大类，即电子检测器、阴极荧光检测器和X射线检测器。

4）真空系统

镜筒和样品室处于高真空下，一般不得高于1×10－2 Pa，它由机械泵和分子涡轮泵来实现。开机后先由机械泵抽低真空，约20分钟后由分子涡轮泵抽真空，约几分钟后就能达到高真空度。此时才能放试样进行测试，在放试样或更换灯丝时，阀门会将镜筒部分、电子枪室和样品室分别分隔开，这样保持镜筒部分真空不被破坏。

5）电源系统

由稳压、稳流及相应的安全保护电路所组成，提供扫描电镜各部分所需要的电源。3.扫描电镜的特点

1）有较高的放大倍数，20-20万倍之间连续可调；

2）有很大的景深，视野大，成像富有立体感，可直接观察各种试样凹凸不平表面的细微结构；

3）试样制备简单。

三．实验内容与结果分析

1.样品制备

1）固体样品：将样品用双面胶带或导电胶固定于样品台上，非导体样品需要喷镀金或铂导电层。

2 ）粉末样品：将样品均匀洒落在贴有双面胶带的样品台上，用吸耳球吹去未粘牢的颗粒，非导体样品需要喷镀金或铂导电层。

3 ）生物样品：干燥含水样品并进行固定，然后进行喷镀处理。

4 ）需作能谱分析的非导电样品应喷镀碳膜

2.实验中观察表征的样品图例

（a）x1800倍的ZnO的BES图像（b）x2300倍的ZnO的BES图像

（c）x4000倍的ZnO的SEI图像（d）x4000倍的ZnO的BES图像

如图（a）、（b）所示ZnO的BES图像：样品微区分布不同尺寸的晶粒、晶界，由于主要成分ZnO晶粒生长所需的能量较少，表现为晶粒尺寸较大；重金属元素在小晶粒中取代Zn的位置，晶粒在生长过程所需的能量较大，表现为晶粒尺寸相对较小。

对比（c）、（d）所示ZnO样品的SEI、BES图像：样品表面倾斜度较小的区域产生的二次电子较少，表现为微区图像亮度相对较暗; 样品表面倾斜度较大的区域产生的二次电子较多，表现为微区图像亮度相对明亮。

在背散射电子（BES）成像模式下，成像区域亮说明此处元素的原子序数大，成像区域暗说明此处元素的原子序数小。通过SEI图像纵向对比，亮的地方相对暗的地方更平整，微区图像（c）、（d）反映了ZnO表面缺陷的存在。

四．扫描电镜在材料分析中的应用

扫描电镜在分辨率、景深及微分析等方面具有一定优越性。利用二次电子背散射电子成像，放大倍数在10到10万倍之间连续可调，观察凹凸不平的表面时，立体感强，细节显示清晰；能对试样的亚表面实现非破坏性的剖面成像；可以实现试样表面、断面从低倍到高倍的定位分析；在样品室中试样不仅可以沿三维空间移动，还能够根据观察需要进行空间转动，以利于使用者对感兴趣的部位进行连续、系统的观察分析；扫描电子显微图像因真实、清晰，并富有立体感；且样品制备方法简单，因而发展迅速，应用广泛。目前，在半导体、金属和陶瓷材料，电子器件及生物材料的显微组织三维形态的观察研究方面获得了广泛应用。

五. 思考题

1. 通过实验你体会扫描电镜有哪些特点？

能直接观察样品的表面结构；制样方法简单，取样后可不做任何改变来观察表面形貌；样品可在样品室中作三维空间的平移和旋转；放大倍数可在大范围内连续可调；富有立体感，景深长，视野范围很大；分辨能力高，可达10nm；可对样品进行综合分析和动态观察；可通过调节图像衬度观察到清晰图像。

2. 根据实验叙述样品制备的步骤？

在样品台上粘上少量的导电胶；用棉签粘取少量干燥的固体样品后涂在导电胶上；去除多余未粘在导电胶上的粉末；若样品导电性差必须喷金。

喷金一般在真空镀膜机上进行，表面喷镀不要太厚、否则会掩盖细节，也不能太薄，不均匀，一般控制在5～10 nm为宜。需要设置适当的电流大小和真空度。

3. 说明二次电子像成像过程？

二次电子是进入样品表面的部分一次电子轰击样品原子激发的单电子。其由静电场引导，沿镜简向下加速。在镜筒中，通过一系列电磁透镜将电子束聚焦并射向样品。靠近镜简底部，在样品表面上方，扫描线圈使电子束以光栅扫描方式偏转。最后一级电磁透镜把电子束聚焦成一个尽可能小的斑点射入样品，从而激发出各种成像信号，其强弱随样品表面的形貌和组成元素不同而变化。

3. 背散射电子像的应用？

原子序数衬度是利用对样品微区原子序数或化学成分变化敏感的物理信号作为调制信号得到的一种显示微区化学成分差别的像衬度。背散射电子信号对微区原子序数或化学成分的变化敏感，都可以作为原子序数衬度或化学成分衬度。

背散射电子是被样品原子反射回来的入射电子，样品背散射系数η随元素原子序数Z 的增加而增加。即样品表面平均原子序数越高的区域，产生的背散射电子信号越强，在背散射电子像上显示的衬度越亮；反之越暗。因此可以根据背散射电子像（成分像）亮暗衬度来判断相应区域原子序数的相对高低。

通过背散射电子像检测器，就可以对样品进行shadow像、topo像和comp像的观察。

材料分析方法实验报告

篇一：材料分析方法实验报告篇二：材料分析方法课程设计报告材料分析测试方法课程设计（论文）题目：磁控溅射c/w多层膜成分及微观分析学院材料科学与工程专业材料化学班级材化082 学生王维娜学号 3080101296 指导教师陈迪春起止时间 2010.12.27-2011.1.1 年材料分析测试方法课程设计任务书课程设计内容要求：掌握高分辨透射电子显微镜样品制备方法，学习并了解真空镀膜技术-磁控溅射技术，多层膜制备过程，以及其微观结构分析，成分分析所用仪器和原理。学生（签名）月日材料分析测试方法课程设计评语指导教师（签名）年日目录材料分析测试方法 ............................................................................. .. (1) 1.1 磁控溅射 ............................................................................. (5) 1.2 x射线衍射仪 ............................................................................. . (5) 1.3 透射电子显微镜 ............................................................................. (6) 1.4 x射线光电子能谱仪（xps） ........................................................................ (7) 第二章实验方法 ............................................................................. .. (9) 2.1 tem样品的制备方法 .............................................................................

数值分析实验报告176453

实验报告插值法数学实验室数值逼近算法设计级 ____________________________ 号 ____________________________ 名 _____________________________ 实验项目名称实验室所属课程名称实验类型实验日期

实验概述：【实验目的及要求】本次实验的目的是熟练《数值分析》第二章“插值法”的相关内容，掌握三种插多项式插值，三次样条插值，拉格朗日插值，并比较三种插值方法的优劣。本次试验要求编写牛顿多项式插值，三次样条插值，拉格朗日插值的程序编码，并去实现。【实验原理】《数值分析》第二章“插值法”的相关内容，包括：牛顿多项式插值，三次样条插值，拉格朗日插值的相应算法和相关性质。【实验环境】（使用的软硬件）软件： MATLAB 2012a 硬件：电脑型号：联想 Lenovo 昭阳E46A 笔记本电脑操作系统： Win dows 8专业版处理器：In tel ( R Core ( TM i3 CPU M 350 @2.27GHz 2.27GHz 实验内容：【实验方案设计】第一步，将书上关于三种插值方法的内容转化成程序语言，用 MATLA B 现; 第二步，分别用牛顿多项式插值，三次样条插值，拉格朗日插值求解不同的问题。【实验过程】（实验步骤、记录、数据、分析）实验的主要步骤是：首先分析问题，根据分析设计 MATLA 程序，利用程序算出问题答案，分析所得答案结果，再得出最后结论。实验一：已知函数在下列各点的值为试用4次牛顿插值多项式 P 4（ x ）及三次样条函数 S （ x ）（自然边界条件）对数据进行插值。用图给出{（ X i , y i ）, X i =0.2+0.08i , i=0 , 1, 11, 10 } , P 4 （ x ）及 S （ x ）。值方法：牛顿在MATLAB 件中

数据分析实验报告

数据分析实验报告文稿归稿存档编号：[KKUY-KKIO69-OTM243-OLUI129-G00I-FDQS58-

第一次试验报告习题1.3 1建立数据集，定义变量并输入数据并保存。 2数据的描述，包括求均值、方差、中位数等统计量。分析—描述统计—频率，选择如下：输出：统计量全国居民农村居民城镇居民 N 有效 22 22 22 缺失均值 1116.82 747.86 2336.41 中值 727.50 530.50 1499.50 方差 1031026.918 399673.838 4536136.444 百分位数 25 304.25 239.75 596.25 50 727.50 530.50 1499.50 75 1893.50 1197.00 4136.75 3画直方图，茎叶图，QQ 图。（全国居民）分析—描述统计—探索，选择如下：输出：全国居民 Stem-and-Leaf Plot Frequency Stem & Leaf 5.00 0 . 56788 数据分析实验报告【最新资料，WORD 文档，可编辑修改】

2.00 1 . 03 1.00 1 . 7 1.00 2 . 3 3.00 2 . 689 1.00 3 . 1 Stem width: 1000 Each leaf: 1 case(s) 分析—描述统计—QQ图，选择如下：输出：习题1.1 4数据正态性的检验：K—S检验，W检验数据：取显着性水平为0.05 分析—描述统计—探索，选择如下：（1）K—S检验

结果：p=0.735 大于0.05 接受原假设，即数据来自正太总体。（2 ）W 检验结果：在Shapiro-Wilk 检验结果972.00 w ，p=0.174大于0.05 接受原假设，即数据来自正太总体。习题1.5 5 多维正态数据的统计量数据：

数值分析实验报告1

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 姓名：万轩实验二插值法

材料分析(SEM)实验报告

材料专业实验报告题目：扫描电镜(SEM)物相分析实验学院：先进材料与纳米科技学院专业：材料物理与化学姓名：学号：1514122986 2016年6月30日

扫描电镜(SEM)物相分析实验一．实验目的 1.了解扫描电镜的基本结构与原理 2.掌握扫描电镜样品的准备与制备方法 3.掌握扫描电镜的基本操作并上机操作拍摄二次电子像 4.了解扫描电镜图片的分析与描述方法二．实验原理 1.扫描电镜的工作原理扫描电镜（SEM）是用聚焦电子束在试样表面逐点扫描成像。试样为块状或粉末颗粒，成像信号可以是二次电子、背散射电子或吸收电子。其中二次电子是最主要的成像信号。由电子枪发射的电子，以其交叉斑作为电子源，经二级聚光镜及物镜的缩小形成具有一定能量、一定束流强度和束斑直径的微细电子束，在扫描线圈驱动下，于试样表面按一定时间、空间顺序作栅网式扫描。聚焦电子束与试样相互作用，产生二次电子发射以及背散射电子等物理信号，二次电子发射量随试样表面形貌而变化。二次电子信号被探测器收集转换成电讯号，经视频放大后输入到显像管栅极，调制与入射电子束同步扫描的显像管亮度，得到反映试样表面形貌的二次电子像。本次实验中主要通过观察背散射电子像及二次电子像对样品进行分析表征。 1）背散射电子背散射电子是指被固体样品原子反射回来的一部分入射电子，其中包括弹性背反射电子和非弹性背反射电子。弹性背反射电子是指被样品中原子和反弹回来的，散射角大于90度的那些入射电子，其能量基本上没有变化（能量为数千到数万电子伏）。非弹性背反射电子是入射电子和核外电子撞击后产生非弹性散射，不仅能量变化，而且方向也发生变化。非弹性背反射电子的能量范围很宽，从数十电子伏到数千电子伏。背反射电子的产生范围在100nm-1mm深度。背反射电子产额和二次电子产额与原子序数的关系背反射电子束成像分辨率一般为50-200nm（与电子束斑直径相当）。背反射电子的产额随原子序数的增加而增加，所以，利用背反射电子作为成像信号不仅能分析形貌特征，也可以用来显示原子序数衬

数值计算实验报告

(此文档为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);

数值分析实验报告模板

数值分析实验报告模板篇一：数值分析实验报告(一)(完整) 数值分析实验报告 1 2 3 4 5 篇二：数值分析实验报告实验报告一题目：非线性方程求解摘要：非线性方程的解析解通常很难给出，因此线性方程的数值解法就尤为重要。本实验采用两种常见的求解方法二分法和Newton法及改进的Newton法。利用二分法求解给定非线性方程的根，在给定的范围内，假设f(x,y)在[a,b]上连续，f(a)xf(b) 直接影响迭代的次数甚至迭代的收敛与发散。即若x0 偏离所求根较远，Newton法可能发散的结论。并且本实验中还利用利用改进的Newton法求解同样的方程，且将结果与Newton法的结果比较分析。前言：（目的和意义）掌握二分法与Newton法的基本原理和应用。掌握二分法的原理，验证二分法，在选对有根区间的前提下，必是收

敛，但精度不够。熟悉Matlab语言编程，学习编程要点。体会Newton使用时的优点，和局部收敛性，而在初值选取不当时，会发散。数学原理：对于一个非线性方程的数值解法很多。在此介绍两种最常见的方法：二分法和Newton法。对于二分法，其数学实质就是说对于给定的待求解的方程f(x)，其在[a,b]上连续，f(a)f(b) Newton法通常预先要给出一个猜测初值x0，然后根据其迭代公式xk?1?xk?f(xk) f'(xk) 产生逼近解x*的迭代数列{xk}，这就是Newton法的思想。当x0接近x*时收敛很快，但是当x0选择不好时，可能会发散，因此初值的选取很重要。另外，若将该迭代公式改进为 xk?1?xk?rf(xk) 'f(xk) 其中r为要求的方程的根的重数，这就是改进的Newton 法，当求解已知重数的方程的根时，在同种条件下其收敛速度要比Newton法快的多。程序设计：本实验采用Matlab的M文件编写。其中待求解的方程写成function的方式，如下 function y=f(x);

编译原理实验报告材料(预测分析报告表方法)

预测分析表方法一、实验目的理解预测分析表方法的实现原理。二、实验内容：编写一通用的预测法分析程序，要求有一定的错误处理能力，出错后能够使程序继续运行下去，直到分析过程结束。可通过不同的文法（通过数据表现）进行测试。三、实验步骤 1．算法数据构造：构造终结符数组：char Vt[10][5]={“id”,”+”……}；构造非终结符数组：char Vn[10]={ }; 构造follow集数组：char *follow[10][10]={ } (可将follow集与预测分析表合并存放) 数据构造示例（使用的预测分析表构造方法1）： /*data1.h简单算术表达式数据*/ char VN[10][5]={"E","E'","T","T'","F"}; //非终结符表 int length_vn=5; //非终结符的个数 char VT[15][5]={"id","+","*","(",")","#"}; //终结符表 int length_vt=6; //终结符的个数 char Fa[15][10]={"TE'","+TE'","","FT'","*FT'","","(E)","id"}; //产生式表:0:E->TE' 1:E'->+TE' 2:E'->空 // 3:T->FT' 4:T'->*FT' 5:T'->空 6:F->(E) 7:F->id int analysis_table[10][11]={0,-1,-1,0,-2,-2,0,0,0,0,0, -1,1,-1,-1,2,2,0,0,0,0,0, 3,-2,-1,3,-2,-2,0,0,0,0,0, -1,5, 4,-1,5, 5,0,0,0,0,0, 7,-2,-2,6,-2,-2,0,0,0,0,0}; //预测分析表,-1表示出错，-2表示该行终结符的follow集合，用于错误处理，正数表示产生式在数组Fa 中的编号，0表示多余的列。（1）预测分析表的构造方法1 给文法的正规式编号：存放在字符数组中，从0开始编号，正规式的编号即为该正规式在数组中对应的下标。如上述Fa数组表示存储产生式。构造正规式数组：char P[10][10]={“E->TE’”,”E’->+TE’”,……..}; (正规式可只存储右半部分，如E->TE’可存储为TE’，正规式中的符号可替换，如可将E’改为M ) 构造预测分析表：int analyze_table[10][10]={ } //数组元素值存放正规式的编号，-1表示出错（2）预测分析表的构造方法2 可使用三维数组 Char analyze_table[10][10][10]={ }

数据分析实验报告

《数据分析》实验报告班级：07信计0班学号：姓名：实验日期2010-3-11 实验地点：实验楼505 实验名称：样本数据的特征分析使用软件名称：MATLAB 实验目的1.熟练掌握利用Matlab软件计算均值、方差、协方差、相关系数、标准差与变异系数、偏度与峰度，中位数、分位数、三均值、四分位极差与极差； 2.熟练掌握jbtest与lillietest关于一元数据的正态性检验； 3.掌握统计作图方法； 4.掌握多元数据的数字特征与相关矩阵的处理方法；实验内容安徽省1990-2004年万元工业GDP废气排放量、废水排放量、固体废物排放量以及用于污染治理的投入经费比重见表6.1.1，解决以下问题：表6.1.1废气、废水、固体废物排放量及污染治理的投入经费占GDP比重年份万元工业GDP 废气排放量万元工业GDP 固体物排放量万元工业GDP废水排放量环境污染治理投资占GDP比重（立方米）（千克）（吨）（%）1990 104254.40 519.48 441.65 0.18 1991 94415.00 476.97 398.19 0.26 1992 89317.41 119.45 332.14 0.23 1993 63012.42 67.93 203.91 0.20 1994 45435.04 7.86 128.20 0.17 1995 46383.42 12.45 113.39 0.22 1996 39874.19 13.24 87.12 0.15 1997 38412.85 37.97 76.98 0.21 1998 35270.79 45.36 59.68 0.11 1999 35200.76 34.93 60.82 0.15 2000 35848.97 1.82 57.35 0.19 2001 40348.43 1.17 53.06 0.11 2002 40392.96 0.16 50.96 0.12 2003 37237.13 0.05 43.94 0.15 2004 34176.27 0.06 36.90 0.13 1.计算各指标的均值、方差、标准差、变异系数以及相关系数矩阵； 2.计算各指标的偏度、峰度、三均值以及极差； 3.做出各指标数据直方图并检验该数据是否服从正态分布？若不服从正态分布，利用boxcox变换以后给出该数据的密度函数； 4.上网查找1990-2004江苏省万元工业GDP废气排放量，安徽省与江苏省是否服从同样的分布？

(完整版)哈工大-数值分析上机实验报告

实验报告一题目：非线性方程求解摘要：非线性方程的解析解通常很难给出，因此线性方程的数值解法就尤为重要。本实验采用两种常见的求解方法二分法和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);

数值分析2016上机实验报告

序言数值分析是计算数学的范畴，有时也称它为计算数学、计算方法、数值方法等，其研究对象是各种数学问题的数值方法的设计、分析及其有关的数学理论和具体实现的一门学科，它是一个数学分支。是科学与工程计算（科学计算）的理论支持。许多科学与工程实际问题（核武器的研制、导弹的发射、气象预报）的解决都离不开科学计算。目前，试验、理论、计算已成为人类进行科学活动的三大方法。数值分析是计算数学的一个主要部分,计算数学是数学科学的一个分支,它研究用计算机求解各种数学问题的数值计算方法及其理论与软件实现。现在面向数值分析问题的计算机软件有：C,C++,MATLAB,Python,Fortran等。 MATLAB是matrix laboratory的英文缩写，它是由美国Mathwork公司于1967年推出的适合用于不同规格计算机和各种操纵系统的数学软件包，现已发展成为一种功能强大的计算机语言，特别适合用于科学和工程计算。目前，MATLAB应用非常广泛，主要用于算法开发、数据可视化、数值计算和数据分析等，除具备卓越的数值计算能力外，它还提供了专业水平的符号计算，文字处理，可视化建模仿真和实时控制等功能。本实验报告使用了MATLAB软件。对不动点迭代，函数逼近（lagrange插值，三次样条插值，最小二乘拟合），追赶法求解矩阵的解，4RungeKutta方法求解，欧拉法及改进欧拉法等算法做了简单的计算模拟实践。并比较了各种算法的优劣性，得到了对数值分析这们学科良好的理解，对以后的科研数值分析能力有了极大的提高。

目录序言 (1) 问题一非线性方程数值解法 (3) 1.1 计算题目 (3) 1.2 迭代法分析 (3) 1.3计算结果分析及结论 (4) 问题二追赶法解三对角矩阵 (5) 2.1 问题 (5) 2.2 问题分析（追赶法） (6) 2.3 计算结果 (7) 问题三函数拟合 (7) 3.1 计算题目 (7) 3.2 题目分析 (7) 3.3 结果比较 (12) 问题四欧拉法解微分方程 (14) 4.1 计算题目 (14) 4.2.1 方程的准确解 (14) 4.2.2 Euler方法求解 (14) 4.2.3改进欧拉方法 (16) 问题五四阶龙格-库塔计算常微分方程初值问题 (17) 5.1 计算题目 (17) 5.2 四阶龙格-库塔方法分析 (18) 5.3 程序流程图 (18) 5.4 标准四阶Runge-Kutta法Matlab实现 (19) 5.5 计算结果及比较 (20) 问题六舍入误差观察 (22) 6.1 计算题目 (22) 6.2 计算结果 (22) 6.3 结论 (23) 7 总结 (24) 附录

材料现代分析方法实验报告

力学与材料学院材料现代分析方法实验报告二 XRD图谱分析专业年级：1 姓名：1 指导老师：1 学号：1 2016年12月中国南京目录实验名称：XRD图谱分析…………………………………………… 一、实验目的……………………………………………………

二、实验要求…………………………………………………… 三、操作过程…………………………………………………… 四、结果分析与讨论……………………………………………… 实验名称：XRD图谱分析一、实验目的了解XRD基本原理及其应用，不同物相晶体结构XRD图谱的区别，熟练掌握如何来分析利用X射线测试得到的XRD图谱。二、实验要求

1、熟练掌握如何来利用软件打开、分析XRD图谱，以及输出分析结果。 2、明确不同物质的XRD图谱，掌握XRD图谱包含的晶体结构的关系，通过自己分析、数据查找和鉴别的全过程，了解如何利用软件正确分析和确定不同物相的XRD图谱，并输出分析结果。 3、实验报告的编写，要求报告能准确的反映实验目的、方法、过程及结论。三、操作过程 1、启动Jade 6.0，并打开实验数据。 2、点击图标使图谱平滑后，再连续两次点击图标扣除背景影响。 3、右击工具栏中的图标，全选左侧的项目，取消选择右侧中的Use Chemistry Filter，最后在下方选择S/M Focus on Major Phases（如图一），并点击OK。图一

4、得到物相分析，根据FOM值（越小，匹配性越高）可推断出该物相为以ZnO为主，可能含有CaF2、Al2O3、Mg(OH)2混合组成的物质（如图二），双击第一种物质可以得到主晶相的PDF卡片（如图三），点击图三版面中的Lines可以观察到不同角度处的衍射强度（如图四）。图二

数值分析实验报告

实验一、误差分析一、实验目的 1．通过上机编程，复习巩固以前所学程序设计语言及上机操作指令； 2．通过上机计算，了解误差、绝对误差、误差界、相对误差界的有关概念； 3．通过上机计算，了解舍入误差所引起的数值不稳定性。二．实验原理误差问题是数值分析的基础，又是数值分析中一个困难的课题。在实际计算中，如果选用了不同的算法，由于舍入误差的影响，将会得到截然不同的结果。因此，选取算法时注重分析舍入误差的影响，在实际计算中是十分重要的。同时，由于在数值求解过程中用有限的过程代替无限的过程会产生截断误差，因此算法的好坏会影响到数值结果的精度。三．实验内容对20,,2,1,0 =n ，计算定积分 ?+=10 5dx x x y n n . 算法1：利用递推公式 151--=n n y n y , 20,,2,1 =n , 取 ?≈-=+=1 00182322.05ln 6ln 51dx x y . 算法2：利用递推公式 n n y n y 51511-= - 1,,19,20 =n . 注意到 ???=≤+≤=10 10202010201051515611261dx x dx x x dx x , 取 008730.0)12611051(20120≈+≈y .: 四.实验程序及运行结果程序一： t=log(6)-log(5);

n=1; y(1)=t; for k=2:1:20 y(k)=1/k-5*y(k-1); n=n+1; end y y =0.0884 y =0.0581 y =0.0431 y =0.0346 y =0.0271 y =0.0313 y =-0.0134 y =0.1920 y =-0.8487 y =4.3436 y =-21.6268 y =108.2176 y =-541.0110 y =2.7051e+003 y =-1.3526e+004 y =6.7628e+004 y =-3.3814e+005 y =1.6907e+006 y =-8.4535e+006 y =4.2267e+007 程序2： y=zeros(20,1); n=1; y1=(1/105+1/126)/2;y(20)=y1; for k=20:-1:2 y(k-1)=1/(5*k)-(1/5)*y(k); n=n+1; end 运行结果：y = 0.0884 0.0580 0.0431 0.0343 0.0285 0.0212 0.0188 0.0169

数值分析实验报告

学生实验报告实验课程名称开课实验室学院年级专业班学生姓名学号开课时间至学年学期

if(A(m,k)~=0) if(m~=k) A([k m],:)=A([m k],:); %换行 end A(k+1:n, k:c)=A(k+1:n, k:c)-(A(k+1:n,k)/ A(k,k))*A(k, k:c); %消去end end x=zeros(length(b),1); %回代求解 x(n)=A(n,c)/A(n,n); for k=n-1:-1:1 x(k)=(A(k,c)-A(k,k+1:n)*x(k+1:n))/A(k,k); end y=x; format short;%设置为默认格式显示,显示5位（2）建立MATLAB界面利用MA TLAB的GUI建立如下界面求解线性方程组：详见程序。五、计算实例、数据、结果、分析下面我们对以上的结果进行测试，求解：

? ? ? ? ? ? ? ? ? ? ? ? - = ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? - - - - - - 7 2 5 10 13 9 14 4 4 3 2 1 13 12 4 3 3 10 2 4 3 2 1 x x x x 输入数据后点击和，得到如下结果：更改以上数据进行测试，求解如下方程组： 1 2 3 4 43211 34321 23431 12341 x x x x ?? ???? ?? ???? ?? ???? = ?? ???? - ?? ???? - ???? ?? 得到如下结果：

材料分析与表征方法实验报告

材料分析与表征方法实验报告热重分析实验报告一、实验目的 1.了解热重分析法的基本原理和差热分析仪的基本构造。 2.掌握热重分析仪的使用方法。二、实验原理热重分析指温度在程序控制时，测量物质质量与温度之间的关系的技术。热重分析所用的仪器是热天平，它的基本原理是，样品重量变化所引起的天平位移量转化成电磁量，这个微小的电量经过放大器放大后，送入记录仪记录；而电量的大小正比于样品的重量变化量。当被测物质在加热过程中有升华、汽化、分解出气体或失去结晶水时，被测的物质质量就会发生变化。三、实验原料一水草酸钙CaC2O4·H2O 四、实验仪器美国TA公司TGA55 升温与降温速率（K/min）0.1-100℃/min 天平灵敏度（μg）0.1μg 温度范围（°C）室温-1000℃ 五、操作条件

第一组：10℃／min空气条件下和20℃／min空气条件下，对TG和DTG 曲线进行对比。第二组：10℃／min空气条件下和10℃／min氮气条件下，对DSC进行对比。第三组：10℃／min氮气条件下，得到TG、DTG、DSC曲线。六、结果与讨论含有一个结晶水的草酸钙(242CaC.OHO)在100℃以前没有失重现象,其热重曲线呈水平状,为TG曲线的第一个平台。DTG曲线在0刻度。在100℃和200℃之间失重并出现第二个平台。DTG曲线先升后降，在108.4℃达到最大值，即失重速率的最大值。DSC曲线先降后升，在188.4℃达到最小值，即热功率的最小值。这一步的失重量占试样总质量的12.47%，相当于每mo CaC2O4·H2O失掉1mol H2O,其热分解反应为： CaC2O4·H2O CaC2O4 + H2O 在400℃和500℃之间失重并开始呈现第三个平台,DTG曲线先升后降，在

数据分析实验报告

数据分析实验报告【最新资料，WORD文档，可编辑修改】第一次试验报告习题1.3 1建立数据集，定义变量并输入数据并保存。 2数据的描述，包括求均值、方差、中位数等统计量。分析—描述统计—频率，选择如下：输出：

方差1031026.918399673.8384536136.444百分位数25304.25239.75596.25 50727.50530.501499.50 751893.501197.004136.75 3画直方图，茎叶图，QQ图。（全国居民）分析—描述统计—探索，选择如下：输出：全国居民Stem-and-Leaf Plot Frequency Stem & Leaf 9.00 0 . 122223344 5.00 0 . 56788 2.00 1 . 03 1.00 1 . 7 1.00 2 . 3 3.00 2 . 689

1.00 3 . 1 Stem width: 1000 Each leaf: 1 case(s) 分析—描述统计—QQ图，选择如下：输出：习题1.1 4数据正态性的检验：K—S检验，W检验数据：取显着性水平为0.05 分析—描述统计—探索，选择如下：（1）K—S检验单样本Kolmogorov-Smirnov 检验身高N60正态参数a,,b均值139.00

标准差7.064 最极端差别绝对值.089 正.045 负-.089 Kolmogorov-Smirnov Z.686 渐近显着性(双侧).735 a. 检验分布为正态分布。 b. 根据数据计算得到。结果：p=0.735 大于0.05 接受原假设，即数据来自正太总体。（2）W检验

数值分析实验报告资料

机电工程学院机械工程陈星星 6720150109 《数值分析》课程设计实验报告实验一函数插值方法一、问题提出对于给定的一元函数)(x f y =的n+1个节点值(),0,1,,j j y f x j n ==。试用Lagrange 公式求其插值多项式或分段二次Lagrange 插值多项式。数据如下：（1 求五次Lagrange 多项式5L ()x ，计算(0.596)f ,(0.99)f 的值。（提示：结果为(0.596)0.625732f ≈, (0.99) 1.05423f ≈）实验步骤：第一步：先在matlab 中定义lagran 的M 文件为拉格朗日函数代码为： function[c,l]=lagran(x,y) w=length(x); n=w-1; l=zeros(w,w); for k=1:n+1 v=1; for j=1:n+1 if(k~=j) v=conv(v,poly(x(j)))/(x(k)-x(j)); end end l(k,:)=v; end c=y*l; end

第二步：然后在matlab命令窗口输入： >>>> x=[0.4 0.55 0.65 0.80,0.95 1.05];y=[0.41075 0.57815 0.69675 0.90 1.00 1.25382]; >>p = lagran(x,y) 回车得到： P = 121.6264 -422.7503 572.5667 -377.2549 121.9718 -15.0845 由此得出所求拉格朗日多项式为 p（x）=121.6264x5-422.7503x4+572.5667x3-377.2549x2+121.9718x-15.0845 第三步：在编辑窗口输入如下命令： >> x=[0.4 0.55 0.65 0.80,0.95 1.05]; >> y=121.6264*x.^5-422.7503*x.^4+572.5667*x.^3-377.2549*x.^2+121.9718 *x-15.0845; >> plot(x,y) 命令执行后得到如下图所示图形，然后 >> x=0.596; >> y=121.6264*x.^5-422.7503*x.^4+572.5667*x.^3-377.2549*x.^2+121.9718 *x-15.084 y =0.6257 得到f（0.596）=0.6257 同理得到f（0.99）=1.0542

数值分析实验报告

实验报告实验项目名称估计水塔的水流量实验项目类型验证演示综合设计其他指导教师成绩一、实验目的（1）学会对实际问题的分析方法（2）学会利用所学的知识解决实际问题（3）设计出相应的算法，编制相应的应用程序二、实验内容某居民区，其自来水是有一个圆柱形水塔提供，水塔高12.2m ，塔的直径为17.4m ，水塔是由水泵根据水塔中的水位自动加水，一般水泵每天工作两次。按照设计，当水塔中的水位降低至最低水位，约8.2m 时，水泵自动启动加水。当水位升至最高水位，约10.8m 时，水泵停止工作。下表给出了某一天的测量记录，测量了28个时刻的数据，但由于水泵正向水塔供水，由3个时刻无法测量到水位（表中为—）。试建立数学模型，计算居民的用水速度和日总用水量。三、实验原理、方法（算法）、步骤时刻 0 0.921 1.843 2.949 3.871 4.978 5.900 水位 9.677 9.479 9.308 9.125 8.982 8.814. 8.686 时刻 7.006 7.928 8.967 9.981 10.925 10.945 12.032 水位 8.525 8.388 8.220 — — 10.820 10.500 时刻 12.954 13.875 14.982 15.903 1 6.826 1 7.931 19.037 水位 10.210 9.936 9.653 9.409 9.180 8.921 8.662 时刻 1 9.959 20.839 22.015 22.958 23.880 24.986 25.908 水位 8.433 8.220 — 10.820 10.591 10.354 10.180

材料研究方法与分析测试实验

本科生实验报告实验课程材料研究方法与分析测试实验学院名称材料与化学化工学院专业名称材料科学与工程(无机非金属方向) 学生姓名学生学号指导教师实验地点实验成绩二〇一四年12月15日——二〇一五年12月19日

填写说明 1、适用于本科生所有的实验报告（印制实验报告册除外）； 2、专业填写为专业全称，有专业方向的用小括号标明； 3、格式要求： ①用A4纸双面打印（封面双面打印）或在A4大小纸上用蓝黑色水笔书写。 ②打印排版：正文用宋体小四号，1.5倍行距，页边距采取默认形式（上下2.54cm，左右2.54cm，页眉1.5cm，页脚1.75cm）。字符间距为默认值（缩放100%，间距：标准）；页码用小五号字底端居中。 ③具体要求：题目（二号黑体居中）；摘要（“摘要”二字用小二号黑体居中，隔行书写摘要的文字部分，小4号宋体）；关键词（隔行顶格书写“关键词”三字，提炼3-5个关键词，用分号隔开，小4号黑体)；正文部分采用三级标题；第1章××(小二号黑体居中，段前0.5行) 1.1 ×××××小三号黑体×××××（段前、段后0.5行） 1.1.1小四号黑体（段前、段后0.5行）参考文献（黑体小二号居中，段前0.5行），参考文献用五号宋体，参照《参考文献著录规则（GB/T 7714－2005）》。

实验一扫描电镜实验（SEM）一、实验目的 1、了解扫描电子显微镜的原理、结构； 2、运用扫描电子显微镜进行样品微观形貌观察。二、实验原理扫描电镜（SEM）是用聚焦电子束在试样表面逐点扫描成像。试样为块状或粉末颗粒，成像信号可以是二次电子、背散射电子或吸收电子。其中二次电子是最主要的成像信号。由电子枪发射的电子，以其交叉斑作为电子源，经二级聚光镜及物镜的缩小形成具有一定能量、一定束流强度和束斑直径的微细电子束，在扫描线圈驱动下，于试样表面按一定时间、空间顺序作栅网式扫描。聚焦电子束与试样相互作用，产生二次电子发射以及背散射电子等物理信号，二次电子发射量随试样表面形貌而变化。二次电子信号被探测器收集转换成电讯号，经视频放大后输入到显像管栅极，调制与入射电子束同步扫描的显像管亮度，得到反映试样表面形貌的二次电子像。扫描电镜由下列五部分组成，如图1（a）所示。各部分主要作用简介如下：

数值分析实验报告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 ε 其中ε是一个非常小的数。这相当于是对（1.1）中19x 的系数作一个小的扰动。我们希望比较（1.1）和（1.2）根的差别，从而分析方程（1.1）的解对扰动的敏感性。实验内容：为了实现方便，我们先介绍两个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 =

数据分析与挖掘实验报告

《数据挖掘》实验报告目录 1.关联规则的基本概念和方法 (1) 1.1数据挖掘 (1) 1.1.1数据挖掘的概念 (1) 1.1.2数据挖掘的方法与技术 (2) 1.2关联规则 (5) 1.2.1关联规则的概念 (5) 1.2.2关联规则的实现——Apriori算法 (7) 2.用Matlab实现关联规则 (12) 2.1Matlab概述 (12) 2.2基于Matlab的Apriori算法 (13) 3.用java实现关联规则 (19) 3.1java界面描述 (19) 3.2java关键代码描述 (23) 4、实验总结 (29) 4.1实验的不足和改进 (29) 4.2实验心得 (30)

1.关联规则的基本概念和方法 1.1数据挖掘 1.1.1数据挖掘的概念计算机技术和通信技术的迅猛发展将人类社会带入到了信息时代。在最近十几年里，数据库中存储的数据急剧增大。数据挖掘就是信息技术自然进化的结果。数据挖掘可以从大量的、不完全的、有噪声的、模糊的、随机的实际应用数据中，提取隐含在其中的，人们事先不知道的但又是潜在有用的信息和知识的过程。许多人将数据挖掘视为另一个流行词汇数据中的知识发现（KDD）的同义词，而另一些人只是把数据挖掘视为知识发现过程的一个基本步骤。知识发现过程如下： ·数据清理（消除噪声和删除不一致的数据）·数据集成（多种数据源可以组合在一起）·数据转换（从数据库中提取和分析任务相关的数据） ·数据变换（从汇总或聚集操作，把数据变换和统一成适合挖掘的形式） ·数据挖掘（基本步骤，使用智能方法提取数

据模式） ·模式评估（根据某种兴趣度度量，识别代表知识的真正有趣的模式） ·知识表示（使用可视化和知识表示技术，向用户提供挖掘的知识）。 1.1.2数据挖掘的方法与技术数据挖掘吸纳了诸如数据库和数据仓库技术、统计学、机器学习、高性能计算、模式识别、神经网络、数据可视化、信息检索、图像和信号处理以及空间数据分析技术的集成等许多应用领域的大量技术。数据挖掘主要包括以下方法。神经网络方法：神经网络由于本身良好的鲁棒性、自组织自适应性、并行处理、分布存储和高度容错等特性非常适合解决数据挖掘的问题，因此近年来越来越受到人们的关注。典型的神经网络模型主要分3大类：以感知机、bp反向传播模型、函数型网络为代表的，用于分类、预测和模式识别的前馈式神经网络模型；以hopfield 的离散模型和连续模型为代表的，分别用于联想记忆和优化计算的反馈式神经网络模型；以art 模型、koholon模型为代表的，用于聚类的自组