- 1、本文档共37页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
计算物理学pptchp2讲述
计算物理 (Computational Physics) 2.2 平均值 方差 标准偏差 目录 计算机数和误差 实验数据的统计处理 正态分布密度函数 ?为数学期望, ?为标准误差 如?=0, ?=1 x=-5:0.01:5; y=exp(-0.5*x.^2)/sqrt(2*pi); plot(x,y) 计算物理 (Computational Physics) 2.2 平均值 方差 标准偏差 目录 计算机数和误差 实验数据的统计处理 对于物理量x的n次测量值,通常用其算术平均值作为观测量的真值的估计值。 算术平均值: 测量值对其平均值的均方偏差: 其数学期望: 故令测量值的方差为: 其数学期望为: 测量值的标准偏差: 计算物理 (Computational Physics) 2.2 平均值 方差 标准偏差 目录 计算机数和误差 实验数据的统计处理 相关MATLAB命令: 标准差:std(), 方差:var(), 最大值:max(),最小值:min(), 平均值:mean(), 求和:sum(),求积: prod(); 例: x=rand(1,100); mean(x) ans = 0.5286 std(x) ans = 0.2818 var(x) ans = 0.0794 max(x) ans = 0.9883 x=1:4 x = 1 2 3 4 sum(x) ans = 10 prod(x) ans = 24 A=[1 2 3;4 5 6] A = 1 2 3 4 5 6 sum(A) ans = 5 7 9 prod(A) ans = 4 10 18 计算物理 (Computational Physics) 2.3 错误值的剔除 目录 计算机数和误差 实验数据的统计处理 错误值:测量数据在采集、传输与记录的过程中由于很多原因可能造成数据丢失或个别数据产生不切合实际的偏差,这种数据称为错误值(或奇异项,坏值)。 拉依达方法: 当某测量值与平均值之差大于标准偏差的3倍,则剔除。 则剔除 xi 肖维勒方法: 当 其中 则剔除 xi 注:肖维勒方法和拉依达方法均属事后处理方法,不适用于实时处理过程。但他们共同的优点是简单易行。 计算物理 (Computational Physics) 2.3 错误值的剔除 目录 计算机数和误差 实验数据的统计处理 一阶差分法: 预估值: 比较判别:如果 则保留 xn否则,用 代替 xn , W 为误差窗 function y = lyd(x) % 拉依达法 sx = std(x); mx = mean(x); y = [];%空向量 for i=1:length(x) d = abs(x(i)-mx); if d = 3*sx y = [y,x(i)];%连接 end end x=randn(1,10000); hist(x,20) y=lyd(x); hist(y,20) 计算物理 (Computational Physics) 2.3 错误值的剔除 目录 计算机数和误差 实验数据的统计处理 function y = cf(x,w) % 差分法剔除错误值 % w为误差窗 n = length(x); y(1) = x(1); y(2) = x(2); for i=1:n-2 xx = 2*x(i+1)-x(i); d = abs(x(i+2)-xx); if d = w y(i+2) = xx; else y(i+2) = x(i+2); end end 例: x=rand(1,1000); hist(x,20) y=cf(x,0.1); hist(y,20) 计算物理 (Computational Physics) 思考题 目录 计算机数和误差 实验数据的统计处理 1.已知一组某物理量的测量数据, 求 x 的平均值, 方差, 标准偏差, 并用拉依达方法剔除所测数据的错误值。 2.用肖维勒方法剔除x=randn(1,10000)中的错误值,并用直方图表示。 3.用差分法剔除x=rand(1,10000)中的错误值,并用直方图表示。 肖维勒方法 function y = xwl(x,n) % 肖维勒法 sx = std
您可能关注的文档
- 计算机组成原理第二章讲述.ppt
- 计算机组成原理第5章存储系统和结构讲述.ppt
- 计算机组成原理第五章-09年讲述.ppt
- 遗传学复习重点详解.doc
- 计算机组成原理第三章存储系统[三]讲述.ppt
- 计算机组成原理题库讲述.doc
- 遗传学课后作业题目及答案详解.docx
- 遗传学试题详解.docx
- 计算机组装与维护实习报告讲述.doc
- 计算机组成原理课件_总结讲述.ppt
- DB23_T 3866-2024冰上龙舟赛事组织服务规范.docx
- DB32 2163-2012 棉纱单位可比综合电耗限额及计算方法.docx
- DB32∕T 2429-2013 棉花田间生长发育观察记载规范.docx
- DB23_T 3820-2024 工业互联网综合平台数据质量管理规范.docx
- DB23_T 3474-2023非煤智慧矿山信息系统技术规范.docx
- DB3201_T 1106-2022 工程地质层划分技术规范.docx
- DB33T 2515-2022公共机构“零碳”管理与评价规范.docx
- DB3301_T 0415-2023 装修垃圾收运处置管理规范.docx
- DB3311/T 106―2019“丽水山居”民宿服务要求与评价规范.docx
- DB3201T 1127-2022 慈善捐赠受赠管理规范.docx
最近下载
- 亿赛通电子文档安全管理系统V5.0--客户端使用手册V1113379.pdf VIP
- [标准规范]煤矿窄轨铁道维修质量标准及检查评级办法(精编).doc VIP
- 人教统编版小升初语文总复习专题七:句式变换课件(共28张PPT).ppt VIP
- 2025重庆永川区板桥镇村(社区) 后备干部选拔笔试备考试题及答案解析.docx VIP
- 深度解析《GBT 44026-2024预制舱式锂离子电池储能系统技术规范》.pptx
- 品质部各岗位职责说明书.pdf VIP
- 2022.01 深圳南山红花岭工业区城市更新项目一期建筑方案设计 华艺.pdf VIP
- 电冰箱安全技术规程培训.pptx VIP
- 07FK02 防空地下室 通风设备安装.pdf VIP
- 解决我国农民负担问题的路线图设想-中国发展.pdf
文档评论(0)