下载文档原格式
学习中心/函授站姓名学号20232020601、大作业试题公布时间：(1)毕业班：2023年4月21日；(2)正常班：2023年5月19日；2、考试必须独立完成，如发现抄袭、雷同均按零分计；3、答案须用《西安电子科技大学网络与继续教育学院2023春期末考试答题纸》(个人专属答题纸)手写完成，要求字迹工整、卷面干净、整齐；4、在线上传时间：(1)毕业班学生于2023年4月21日至2023年5月4日在线上传大作业答卷；(2)正常班学生于2023年5月19日至2023年5月29日在线上传大作业答卷；5、拍照要求完整、清晰，一张图片对应一张个人专属答题纸(A4纸)，正确上传。
KP =、K F K P =K F(每空2分，共20分)1、在电子测量技术中，按被测量对象性质及技术处理领域的不同，电子测量包括时域测量、__________测量和数据域测量。
答案;频域2、电子测量所包括的技术有变换技术、比较技术、___________技术、处理技术及显示技术等。
答案：放大技术3、标定一量程为100mA 的电流表，经过大量数据测量以及和标准表的数据比较，发现在10mA 处的误差最大，且为0.8mA，则这块表的等级应标定为__1级_______级表。
4、有效数字体现了测量仪器的误差，现用一电压表测量某电压，测量结果的有效值为2.382V，则该测量仪器的最大误差为__0.0005V _。
5、在各类交流电压表的刻度标定中，电压表的刻度一般是按__________定标刻度的。
答案交流电压有效值6、最小二乘法是处理测量数据很重要的一种方法，一般在工程上常利用用最小二乘法的原理来__________。
答案：拟合测量曲线7、在模拟式交流电压中，常用的三种类型的交流电压表有峰值电压表、__________和有效值电压表。
答案：均值电压表8、用示波器测量某信号的周期，若选择的扫描速度越高，则屏幕显示的波形周期个数越__________。
电子测量技术大作业班级: 通信1109学号: 11211105姓名: 单赟吉专业: 通信工程指导老师: 朱云二零一三年十二月第一题： 一．研究题目：4-19：在Multisim 环境下，设计一种多斜积分式DVM ，给出原理图和仿真实验结果。
二．积分型A/D 转换电路2.1 双积分型A/D 转换电路双积分型ADC 是1种V —T 型A/D 转换器，原理电路如图12.2.2-1(a)所示，由积分器、比较器、计数器和部分控制电路组成。
工作过程如下：（1）平时（即A/D 转换之前），转换控制信号v C =0，计数器和触发器FFc 被清零，门G1、G2输出低电平，开关S 0闭合使电容C 完全放电，S 1掷下方，比较器输出v B =0，门G3关闭。
（2）v C =1时，开关S 0断开，开关S 1掷上方接输入信号V I ，积分器开始对V I 积分，输出电压为⎰-=-=tt RCV dt V RCv 0II O 1(2.1)显然v O 是1条负向积分直线，如图12.2.2-1(b)中t =0～T 1段实线所示。
与此同时，比较器输出v B =1（因v O <0），门G3开启，计数器开始计数。
（3）当积分到t =T 1=2n T cp 时（其中T cp 是时钟CP 的周期），n 位计数器计满2n 复0，FFc 置1，门G2输出高电平，开关S 1掷下方接基准电压（－V REF ），积分器开始对（－V REF ）进行积分。
设t =T 1时，v O 下降到v O =V O1，由式（3.1）1IO1T RCV V -= （2.2）)()(11REFO11REF O1O T t RCV V dt V RCV v tT -+=--=⎰（2.3）v O 波形如图3.5(b)中t =T 1～（T 1+T 2）段实线所示。
（4）当t =T 1+T 2时，v O 上升到v O =0V ，v B =0，门G3被关闭，计数器停止计数，此时计数器中保存下来的数字就是时间T 2。
电子测量技术大作业目录题目一测量数据误差处理 (1)(1)提供测试数据输入、粗大误差判别准则选择等的人机界面； (1)(2)编写程序使用说明； (1)(3)通过实例来验证程序的正确性。
(2)题目二时域反射计 (4)(1)时域反射计简介 (4)(2)时域反射计原理 (4)(3)时域反射计（TDR）组成 (5)(4)仿真与结果 (5)附录 (8)题目一测量数据误差处理2-21 参考例2-2-6的解题过程，用C语言或Matlab设计测量数据误差处理的通用程序,要求如下：(1)提供测试数据输入、粗大误差判别准则选择等的人机界面；图1 测试数据误差处理的输入(2)编写程序使用说明；本题用的是C语言编写的数据误差处理的通用程序，调试编译借助了CodeBlocks软件。
运行exe文件后，只需输入所需测试数据的数目、各数值大小并选择误差处理方式与置信概率即可得出处理结果。
在程序的子函数中已经将t a值表、肖维纳准则表及格拉布斯准则表的所有数据存入，无需人工查表填入。
其他具体程序内容可见附录。
图 2 程序运行流程图(3) 通过实例来验证程序的正确性。
例2-2-6中的原始数据如下表16=2.67()0.0858U =2.45%U V U V V σ∧=——，；异常值为；无累进性系统误差、无周期性系统误差；在95的置信概率下，置信区间为[2.66,2.72]V计算所得结果与图3显示结果近似相等，说明程序编译无误。
图3 数据处理后的结果显示题目二时域反射计6-14 在Multisim环境下，基于Tektronix TDS204虚拟示波器设计一种时域反射计，给出电路原理图和实验仿真结果。
（本题设计以时域反射计测量阻抗为例）(1)时域反射计简介时域反射计(TDR)用来测量信号在通过某类传输环境传导时引起的反射，如电路板轨迹、电缆、连接器等等。
TDR仪器通过介质发送一个脉冲，把来自“未知”传输环境的反射与标准阻抗生成的反射进行比较。
《电子测量技术》课程研究性作业姓名 :学号 :班级：指导教师 :日期 :目录一、6-14 在Multisim环境下，基于Tektronix TDS204虚拟示波器设计一种时域反射计，给出电路原理图和实验仿真结果。
(2)1.1时域反射计简介 (2)1.2时域反射计原理 (2)1.3电路原理图 (5)1.4实验仿真结果 (6)二、7-14查阅网络分析仪的技术资料，说明网络分析仪的功能和基本原理，比较网络分析仪和频谱分析仪的异同点。
(6)2.1网络分析仪的功能 (6)2.2网络分析仪的基本原理 (7)2.3频谱分析仪与网络分析仪的异同点 (7)一、6-14 在Multisim环境下，基于Tektronix TDS204虚拟示波器设计一种时域反射计，给出电路原理图和实验仿真结果。
1.1时域反射计简介时域反射计(TDR)用来测量信号在通过某类传输环境传导时引起的反射，如电路板轨迹、电缆、连接器等等。
TDR仪器通过介质发送一个脉冲，把来自“未知”传输环境的反射与标准阻抗生成的反射进行比较。
TDR 显示了在沿着一条传输线传播快速阶跃信号时返回的电压波形。
波形结果是入射阶跃和阶跃遇到阻抗偏差时产生的反射的组合。
1.2时域反射计原理时域反射计TDR是最常用的测量传输线特征阻抗的仪器，它是利用时域反射的原理进行特性阻抗的测量。
图1是传统TDR工作原理图。
TDR包括三部分组成：1) 快沿信号发生器：典型的发射信号的特征是：幅度200mv，上升时间35ps，频率250KHz方波。
2) 采样示波器：通用的采样示波器.3) 探头系统：连接被测件和TDR仪器。
测试信号的运行特征参考图2所示。
由阶跃源发出的快边沿信号注入到被测传输线上，如果传输线阻抗连续，这个快沿阶跃信号就沿着传输线向前传播。
当传输线出现阻抗变化时，阶跃信号就有一部分反射回来，一部分继续往前传播。
反射回来的信号叠加到注入的阶跃信号，示波器可采集到这个信号。
电子测量大作业【实验题目】查阅网络分析仪的技术资料,说明其功能与工作原理,比较网络分析仪与频谱分析仪的异同点。
目录:网络分析仪简介..．．.....．..．.．.．..．...．..．.．...．.．.......．.....２网络分析仪的功能.......．.．....．．.．．.．.....．.．..．．.．.....．．．...２网络分析仪的原理...．．.．.....．．..............．....．．.．..．.．...．2 频谱分析仪简介.....．．．..．...．..．.....．.........．．....．.．．.....３频谱分析仪的基本原理.....．...．．．..．．．.．....．.．．...．.......．．.．3 结论．．．....．．....．....．.....．.....．.......．...．．..........．．.．4一．网络分析仪简介矢量网络分析仪,它本身自带了一个信号发生器，可以对一个频段进行频率扫描．如果是单端口网络分析仪测量的话，将激励信号加在端口上，通过测量反射回来信号的幅度和相位,就可以判断出阻抗或者反射情况. 而对于双端口测量,则还可以测量传输参数. 由于受分布参数等影响明显,所以网络分析仪使用之前必须进行校准图１网络分析仪二．网络分析仪的功能可直接测量有源或无源、可逆或不可逆的双口和单口网络的复数散射参数，并以扫频方式给出各散射参数的幅度、相位频率特性。
自动网络分析仪能对测量结果逐点进行误差修正，并换算出其他几十种网络参数,如输入反射系数、输出反射系数、电压驻波比、阻抗（或导纳）、衰减（或增益）、相移和群延时等传输参数以及隔离度和定向度等。
三．网络分析仪的原理一个任意多端口网络的各端口终端均匹配时,由第n个端口输入的入射行波ａｎ将散射到其余一切端口并发射出去。
若第ｍ个端口的出射行波为bm,则n 口与ｍ口之间的散射参数Ｓmn=bm/an。
《电子测量技术》课程标准课程名称：电子测量技术 Electronic Measurement Technology课程性质：专业选修学分：2.5总学时：45，理论学时：36，实验(上机)学时：9适用专业：电子信息技术先修课程：模拟电子技术、数字电子技术、信号与系统、微机原理一、教学目的与要求《电子测量技术》是电子信息、自动控制、测量仪器等专业的通用技术基础课程。
包括电子测量的基本原理、测量误差分析和实际应用，主要电子仪器的工作原理，性能指标，电参数的测试方法，该领域的最新发展等。
电子测量技术综合应用了电子、计算机、通信、控制等技术。
通过本课程的学习，培养学生具有电子测量技术和仪器方面的基础知识和应用能力；通过本课程的学习，可开拓学生思路，培养综合应用知识能力和实践能力；培养学生严肃认真，求实求真的科学作风，为后续课程的学习和从事研发工作打下基础。
二、教学内容与学时分配三、各章节主要知识点与教学要求第1章序论第一节测量的基本概念一、测量的定义二、测量的意义三、测量技术第二节计量的基本概念一、计量二、单位和单位制三、计量标准四、测量标准的传递第三节电子测量技术的内容，特点和方法一、电子测量二、电子测量的内容和特点三、电子测量的一般方法第四节电子测量的基本技术一、电子测量的变换技术二、电子测量的放大技术三、电子测量的比较技术四、电子测量的处理技术五、电子测量的显示技术第五节本课程的任务重点：测量的基本概念、基本要素；单位和单位制，基准和标准，量值的传递准则。
难点：量值的传递准则教学要求：理解测量的基本概念、基本要素，测量误差的基本概念和计算方法。
理解计量的基本概念，单位和单位制，基准和标准，量值的传递准则。
理解测量的基本原理，信息获取原理和量值比较原理。
理解电子测量的实现原理：变换、比较、处理、显示技术。
第2章测量误差理论与数据处理第一节测量误差的基本概念一、有关误差的基本概念二、测量误差的基本表示方法第二节测量误差的来源与分类一、测量误差的来源二、测量误差的分类第三节测量误差的分析与处理一、随机误差的分析与处理二、系统误差的判断及消除方法三、粗大误差的分析与处理第四节测量误差的合成与分配一、测量误差的合成二、测量测量不确定度及其合成三、误差分配及最佳测量方案第五节测量数据处理一、有效数字处理二、测量结果的处理三、最小二乘法与回归分析重点：测量误差的分类估计和处理，系统误差和粗大误差的判断及处理，不确定度的评定方法。
电子测量实训报告第一篇：电子测量实训报告电子测量实训报告一.实训目的（1）熟悉常用电子仪器的功能及使用方法。
（2）掌握常用电子仪器的工作原理。
（3）掌握常用电子仪器附加功能的使用。
（4）熟练使用常用电子仪器进行数据测量。
（5）掌握常用电子元器件的测量方法，掌握电子元器件的焊接技巧和装配工艺；学会使用万用表、示波器、毫伏表、频率计、信号发生器等电子测量仪器。
掌握查找电子设备故障的一般方法。
培养学生实际动手操作能力；为学生以后参加工作打下良好的基础。
二．基本要求一、课程性质和任务陕航航空电子设备维修专业的主干专业课程。
其任务是使学生掌握从事航空电子设备维修工作所必需的电子基本工艺和基本技能，初步形成解决实际问题的能力，为学习其他专业知识和职业技能打下基础。
二、课程教学目标(一)知识教学目标1.了解电工电子仪表、仪器的基本结构及正确使用与维护；2.掌握常用电子元器件的正确识别与检测方法；3.理解常用电子电路和简单电子整机电路的分析、检测与常见故障排除方法；4.掌握电子电路安装的工艺知识。
(二)能力培养目标1.能正确使用常用电工电子仪表、仪器；2.能正确阅读分析电路原理图和设备方框图，并能根据原理图绘制简单印刷电路；3.初步学会借助工具书、设备铭牌、产品说明书及产品目录等资料，查阅电子元器件及产品有关数据、功能和使用方法；4.能按电路图要求，正确安装、调试单元电子电路、简单整机电路；5.处理电子设备的典型故障。
(三)思想教育目标1.具有热爱科学、实事求是的学风和创新意识、创新精神；2.加强爱岗敬业意识和职业道德意识。
三、教学内容和要求基础模块(一)常用电子仪器、仪表的使用与维护1.了解常用电子仪器、仪表的结构；2.理解常用电子仪器、仪表的基本功能；3.掌握常用电子仪器、仪表的使用方法和注意事项。
(二)常用电子元器件的识别与检测1.理解常用电子元器件的型号和主要参数；2.理解常用电子元器件的识别和分类方法；3.掌握用万用表检测常用电子元器件的方法。
电子测量大作业（选做题目）题目：数字频率计的设计一设计目的、意义1、设计目的掌握数字频率计的设计方法。
掌握振荡器、分频器、计数译码显示电路、单稳态等相关电路的设计。
2、设计意义在电子技术中，频率是最基本的参数之一，并且与许多电参量的测量方案、测量结果都有十分密切的关系，因此频率的测量就显得更为重要。
测量频率的方法有多种,其中电子计数器测量频率具有精度高、使用方便、测量迅速，以及便于实现测量过程自动化等优点，是频率测量的重要手段之一。
电子计数器测频有两种方式：一是直接测频法，即在一定闸门时间内测量被测信号的脉冲个数；二是间接测频法，如周期测频法。
直接测频法适用于高频信号的频率测量，间接测频法适用于低频信号的频率测量。
二设计内容1、总体方案简介本次设计的数字频率计由四部分组成：时基电路、闸门电路、逻辑控制电路以及可控制的计数、译码、显示电路。
如图2 .1所示图2.1 数字频率计原理图由555 定时器, 分级分频系统及门控制电路得到具有固定宽度T的方波脉冲做门控制信号,时间基准T称为闸门时间。
宽度为T的方波脉冲控制闸门的一个输入端B。
被测信号频率为fx,周期为Tx。
到闸门另一输入端A。
当门控制电路的信号到来后,闸门开启,周期为Tx的信号脉冲和周期为T 的门控制信号结束时过闸门,于输出端C产生脉冲信号到计数器, 计数器开始工作, 直到门控信号结束, 闸门关闭。
单稳1 的暂态送入锁存器的使能端, 锁存器将计数结果锁存, 计数器停止计数并被单稳2暂态清零。
若T=1s , 计数器显示f x=N( T 时间内的通过闸门信号脉冲个数) 。
若T=0. 1s ,通过闸门脉冲个数位N时, f x=10N( 闸门时间为0. 1s 时通过闸门的脉冲个数) 。
也就是说, 被测信号的频率计算公式是f x=N/T。
由此可见,闸门时间决定量程,可以通过闸门时基选择开关,选择T大一些,测量准确度就高一些,T小一些,则测量准确度就低。
测量误差理论和测量数据处理班级：姓名：学号：指导老师：1. 变值系统误差的判定1.1马利科夫判据马利科夫判据是常用的判别有无累进性系统误差的方法。
把 n 个等精密度测量值所对应的残差按测量顺序排列，把残差分成前后两部分求和，再求其差值。
若测量中含有累进性系统误差，则前后两部分残差和明显不同，差值应明显地异于零。
所以马利科夫判据是根据前后两部分残差和的差值来进行判断。
当前后两部分残差和的差值近似等于零，则上述测量数据中不含累进性系统误差， 若其明显地不等于零（与最大的残差值相当或更大） ，则说明上述测量数据中存在累进性系统误差。
n 为偶数时： n 为奇数时: 若 ，则存在累进性系差，否则不存在累进性系差。
1.2 阿卑-赫梅特判据通常用阿卑—赫梅特判据来检验周期性系统误差的存在。
把测量数据按测量顺序排列，将对应的残差两两相乘，然后求其和的绝对值，再与总体方差的估计相比较，若式 成立则可认为测量中存在周期性系统误差。
当我们按照随机误差的正态分布规律检查测量数据时，如果发现应该剔除的粗大误差占的比例较大时，就应该怀疑测量中含有非正态分布的系统误差。
存在变值系统误差的测量数据原则上应舍弃不用。
但是，若虽然存在变值系统误差，但残差的 最大值明显地小于测量允许的误差范围或仪器规定的系统误差范围，则测量数据可以考虑使用，在继续测量时需密切注意变值系统误差的情况。
2.粗大误差剔除的常用准则2.1莱特准则 若，则x i 为异常值剔除不用；否则不存在异常值。
莱特检验法是一种测量数据服从正态分布情况下判别异常值的方法，主要用于测量数据数量较多的情况，一般要求测量次数 n 大于 10。
2.2肖维纳准则若 ，则x i 为异常值应剔除不用。
否则不存在异常值。
肖维纳检验法也是以正态分布作为前提的，假设多次重复测量所得n 个测量值中，当残差绝对值时，则认为是粗大误差。
式中n 是系数，可通过查肖维纳准则表得到。
要注意的是肖维纳检验法是建立在测量数据服从正态分布的前提下，要求在n ＞5时使用。
电⼦测量⼤作业（基于C51单⽚机的简单数字电压表）基于单⽚机的简易数字电压表的设计⼀、概述本课题设计是⼀种基于单⽚机的简易数字电压表的设计。
该设计主要由三个模块组成：A/D转换模块，数据处理模块及显⽰模块。
A/D转换主要由芯⽚ADC0808来完成，它负责把采集到的模拟量转换为相应的数字量在传送到数据处理模块。
数据处理则由芯⽚AT89C51来完成，其负责把ADC0808传送来的数字量经过⼀定的数据处理，产⽣相应的显⽰码送到显⽰模块进⾏显⽰；此外,它还控制着ADC0808芯⽚⼯作。
⼆、主要芯⽚1、ADC0808 主要特性ADC0808是CMOS单⽚型逐次逼近式A/D转换器，带有使能控制端，与微机直接接⼝，⽚内带有锁存功能的8路模拟多路开关，可以对8路0-5V输⼊模拟电压信号分时进⾏转换，由于ADC0808设计时考虑到若⼲种模/数变换技术的长处，所以该芯⽚⾮常适应于过程控制，微控制器输⼊通道的接⼝电路，智能仪器和机床控制等领域。
ADC0808主要特性:8路8位A/D转换器，即分辨率8位；具有锁存控制的8路模拟开关；易与各种微控制器接⼝；可锁存三态输出，输出与TTL兼容；转换时间：128µs；转换精度：0.2%；单个+5V电源供电；模拟输⼊电压范围0- +5V，⽆需外部零点和满度调整；低功耗，约15mW。
ADC0808芯⽚有28条引脚，采⽤双列直插式封装，其引脚图如图3所⽰。
图1 ADC0808引脚图下⾯说明各个引脚功能:IN0-IN7（8条）：8路模拟量输⼊线，⽤于输⼊和控制被转换的模拟电压。
地址输⼊控制（4条）：ALE: 地址锁存允许输⼊线，⾼电平有效，当ALE为⾼电平时，为地址输⼊线，⽤于选择IN0-IN7上那⼀条模拟电压送给⽐较器进⾏A/D转换。
ADDA,ADDB,ADDC: 3位地址输⼊线，⽤于选择8路模拟输⼊中的⼀路，其对应关系如表1所⽰：表1 ADC0808通道选择表START：START为“启动脉冲”输⼊法，该线上正脉冲由CPU送来，宽度应⼤于100ns，上升沿清零SAR,下降沿启动ADC⼯作。
电子测量技术大作业目录题目一测量数据误差处理 (1)(1)提供测试数据输入、粗大误差判别准则选择等的人机界面； (1)(2)编写程序使用说明； (1)(3)通过实例来验证程序的正确性。
(1)题目二时域反射计 (1)(1)时域反射计简介 (1)(2)时域反射计原理 (2)(3)时域反射计（TDR）组成 (2)(4)仿真与结果 (2)附录 (2)题目一测量数据误差处理2-21 参考例2-2-6的解题过程，用C语言或Matlab设计测量数据误差处理的通用程序,要求如下：(1)提供测试数据输入、粗大误差判别准则选择等的人机界面；图 1 测试数据误差处理的输入(2)编写程序使用说明；本题用的是C语言编写的数据误差处理的通用程序，调试编译借助了CodeBlocks软件。
运行exe文件后，只需输入所需测试数据的数目、各数值大小并选择误差处理方式与置信概率即可得出处理结果。
在程序的子函数中已经将t a值表、肖维纳准则表及格拉布斯准则表的所有数据存入，无需人工查表填入。
其他具体程序内容可见附录。
图 2 程序运行流程图(3)通过实例来验证程序的正确性。
例2-2-6中的原始数据如下表1计算所得结果与图3显示结果近似相等，说明程序编译无误。
图 3 数据处理后的结果显示题目二时域反射计6-14 在Multisim环境下，基于Tektronix TDS204虚拟示波器设计一种时域反射计，给出电路原理图和实验仿真结果。
（本题设计以时域反射计测量阻抗为例）(1)时域反射计简介时域反射计(TDR)用来测量信号在通过某类传输环境传导时引起的反射，如电路板轨迹、电缆、连接器等等。
TDR仪器通过介质发送一个脉冲，把来自“未知”传输环境的反射与标准阻抗生成的反射进行比较。
TDR 显示了在沿着一条传输线传播快速阶跃信号时返回的电压波形。
波形结果是入射阶跃和阶跃遇到阻抗偏差时产生的反射的组合。
(2)时域反射计原理时域反射计TDR是最常用的测量传输线特征阻抗的仪器，它是利用时域反射的原理进行特性阻抗的测量。
图 4 TDR原理图(3)时域反射计（TDR）组成①快沿信号发生器：典型的发射信号的特征是：幅度200mv，上升时间35ps，频率250kHz方波。
②采样示波器：通用的采样示波器；③探头系统：连接被测件和TDR仪器。
(4)仿真与结果图 5 时域反射计仿真电路图 6 信号发生器设置选项图7 信号发生器上升沿时间设置图8 示波器仿真显示结果在图8中，第一条黄线为终端开路（反射系数为1）时的结果；第二条蓝线为终端有负载阻抗时的结果；第三条为终端短路（反射系数为-1）时的结果。
图9 TDR测试信号理论运行特征图图10 被测传输线特征阻抗的计算附录#include <stdio.h>#include <math.h>//求绝对值的函数，因为“abs”只适用于整数float Abs(float a){if(a>0.0)return a;return -a;}//计算平均值的子函数float Average(float a[],int n){int i;float sum=0.0;for(i=0; i<=n-1; i++){sum+=a[i];}return sum/n;}//计算方差的子函数float Variance(float a[],int n){int i;float sum=0.0;float xi=Average(a,n);for(i=0; i<=n-1; i++){sum+=(a[i]-xi)*(a[i]-xi);}return (float)sqrt(sum/(n-1));}//删除数组中指定位置的子函数void Delete(float a[],int n,int i){int j;for(j=i; j<=n-2; j++){a[j]=a[j+1];}}//粗大误差处理1//莱特检验法，出口参数为所剔除坏值的个数int Laite(float a[],int n){int t=0; //t用来装粗大误差的个数int i;int z;int flag=1;float aver;float var;while(flag){z=n;aver=Average(a,n);var=Variance(a,n);for(i=0; i<=n-1; i++){if(Abs(a[i]-aver)>3*var){printf(" %f是坏值！\n",a[i]);Delete(a,n,i);n--;t++;}}if(n==z) flag=0;}return t;}//粗大误差处理2//肖维纳准则int Xiaown(float a[],int n){float ch[38]= {0,0,0,0,0,1.65,1.73,1.79,1.86,1.92,1.96,2.00,2.04,2.07,2.10,2.13,2.16,2.18,2.20,2.22,2.24,2.26,2.28,2.30,2.32,2.33,2.34,2.35,2.37,2.38,2.39,2.45,2.50,2.58,2.64,2.74,2.81,3.02}; //肖维纳准则表int t=0; //t用来装粗大误差的个数int i;int z;int flag=1;float aver;float var;while(flag){z=n;aver=Average(a,n);var=Variance(a,n);for(i=0; i<=n-1; i++){if(Abs(a[i]-aver)>ch[n]*var){printf(" %f是坏值！\n",a[i]);Delete(a,n,i);n--;t++;}}if(z==n) flag=0;}}//粗大误差处理3//格拉布斯准则\//格拉布斯准则表1：95%float Ge95(int n){float Ge[30]= {0,0,0,1.15,1.46,1.67,1.82,1.94,2.03,2.11,2.18,2.23,2.29,2.33,2.37,2.41,2.44,2.47,2.50,2.53,2.56,2.58,2.60,2.62,2.64,2.66};if(n<=25) return Ge[n];switch (n){case 30:return 2.74;case 35:return 2.81;case 40:return 2.87;case 50:return 2.96;case 100:return 3.17;default:printf(" 输入n值错误！！！");return 0.0;}}//格拉布斯准则表2: 99%float Ge99(int n){float Ge[30]= {0,0,0,1.16,1.49,1.75,1.94,2.10,2.22,2.32,2.41,2.48,2.55,2.61,2.66,2.71,2.75,2.79,2.82,2.85,2.88,2.91,2.94,2.96,2.99,3.01};if(n<=25) return Ge[n];switch (n){case 30:return 3.10;return 3.18;case 40:return 3.24;case 50:return 3.34;case 100:return 3.58;default:printf(" 输入n值错误！！！");return 0.0;}}//格拉布斯执行函数int Glbs(float a[],int n,int m){int t=0; //t用来装粗大误差的个数int i;int z;float aver;float var;int flag=1,flag1=1,flag2=1;while(flag){if(m==1){while(flag1){z=n;aver=Average(a,n);var=Variance(a,n);for(i=0; i<=n-1; i++){if(Abs(a[i]-aver)>Ge95(n)*var){printf(" %f是坏值！\n",a[i]);Delete(a,n,i);n--;t++;}}if(z==n) flag1=0;}flag=0;return t;}else if(m==2){while(flag2){z=n;aver=Average(a,n);var=Variance(a,n);for(i=0; i<=n-1; i++){if(Abs(a[i]-aver)>Ge99(n)*var){printf(" %f是坏值！\n",a[i]);Delete(a,n,i);n--;t++;}if(z==n) flag2=0;}}flag=0;return t;}else{printf(" 输入有错！！请重新输入：\n");}}}//误差处理子函数int Wucha(float a[],int n,int x,int m){int pa=1;int t;int i;while (pa){if(x==1){pa=0;t=Laite(a,n);n=n-t;printf(" 处理后的数据如下：\n\t\t");for(i=0; i<=n-1; i++){printf("%f\t",a[i]);}return n;}else if(x==2){pa=0;t=Xiaown(a,n);n=n-t;printf(" 处理后的数据如下：\n\t\t");for(i=0; i<=n-1; i++){printf("%f\t",a[i]);}return n;}else if(x==3){pa=0;t=Glbs(a,n,m);n=n-t;printf(" 处理后的数据如下：\n\t\t");for(i=0; i<=n-1; i++){printf("%f\t",a[i]);}return n;}else printf(" 输入错误，请重新输入！\n");}}//判断有无累进性误差的子函数float Max(float a[],int n){int i;float f=Average(a,n);float x,y;x=Abs(a[0]-f);for(i=0; i<=n-1; i++){if(Abs(a[i]-f)>=x){x=Abs(a[i]-f);}}return x;}void Mlkf(float a[],int n){int i;float sum1=0.0,sum2=0.0;float f=Average(a,n);if(n%2==0) //偶数个{for(i=0; i<=n/2-1; i++)sum1=sum1+(a[i]-f);for(i=n/2; i<=n-1; i++)sum2=sum2+(a[i]-f);printf("\n");if(Abs(sum1-sum2)>=Max(a,n))printf(" 测量中存在累进性系统误差。
