测量数据处理实验报告
测量数据处理实验报告
院系____________________ 班级____________________ 姓名____________________ 学号____________________目录（页码根据自己需要填写）一、实验项目须知.....................................二、实验一秩亏自由网上机实验........................三、实验二极大验后滤波、推估上机实验...............四、实验三最小二乘配置上机实验.....................五、实验四平差系统的统计假设检验上机实验...........六、编写实验报告要求................................实验项目须知一、课程介绍本课程是测绘工程类专业的专业选修课程之一，是在学生学习了误差理论与测量平差知识之后又开设的一门专业的数据处理课程，该课程主要讲授近现代测量数据处理理论、模型和方法，同时紧密结合计算机编程，从而提高学生的计算机应用能力、测量数据处理能力，培养学生养成严谨的科学态度，提高分析和解决问题的能力，更好的适应今后工作和学习。
二、测量数据处理实验教学要求1.掌握稳健估计方法处理测量数据；2.掌握极大验后滤波、推估理论和方法处理测量数据；3.初步掌握最小二乘配置理论和方法处理测量数据；4.初步掌握工程或科研项目中测量数据处理的理论、步骤、方法和技巧；5.掌握测量数据处理中常用的计算机语言和软件；6.服从实验指导教师的规定，认真、按时、独立完成任务每次实验结束，提交书写工整或按照规范打印实验报告。
7.在实验过程中，还应遵守纪律，爱护实验室设备、离开后将自己的随身物品带走，并记录仪器使用状况、将凳子整理好并放到电脑桌的下方。
8.上机实验后，提交实验报告，成绩评定重点考核算法、处理结果与结果分析。
三、实验选择要求可根据课堂讲课具体内容，完成4个实验，并要求写出实验报告。
声速的测量实验报告及数据处理
声速的测量实验报告及数据处理一、实验目的与原理1.1 实验目的为了研究声速的测量方法，我们进行了一次声速的测量实验。
通过实验，我们希望能够了解声速的定义、测量原理以及影响声速的因素，从而为实际应用提供理论依据。
1.2 实验原理声速是指在某种介质中，声波传播的速度。
声音是由物体振动产生的机械波，当这种振动传播到介质中时，会引起介质分子的振动，从而形成声波。
声波在介质中的传播速度与其内部分子的振动速度有关，而分子的振动速度又受到温度、压力等因素的影响。
因此，声速的测量实际上是测量介质中分子振动速度的过程。
二、实验设备与材料2.1 设备本次实验使用的设备包括：声源(用于产生声波)、麦克风(用于接收声波)、计时器(用于计算声波传播时间)、数据处理软件(用于分析实验数据)。
2.2 材料实验所使用的材料包括：水、玻璃、铝箔等。
这些材料都是常见的介质，可以用于测量声速。
三、实验步骤与数据处理3.1 实验步骤1) 将水倒入一个透明的容器中，使其充满水。
2) 将玻璃和铝箔分别放在水中。
3) 用麦克风分别对玻璃和铝箔进行录音。
4) 使用计时器记录每次录音所需的时间。
5) 重复以上步骤多次，以获得较为准确的数据。
6) 使用数据处理软件对实验数据进行分析，得出声速的测量结果。
3.2 数据处理我们需要计算每次录音所需的时间。
由于实验过程中可能会受到环境噪声的影响，因此我们需要在每次录音前先将麦克风校准，以减小误差。
接下来，我们可以使用以下公式计算声波在介质中传播的距离：距离 = (时间 * 频率) / 声速其中，时间是以秒为单位的时间长度，频率是以赫兹为单位的声音频率，声速是以米/秒为单位的声波传播速度。
通过对所有数据的分析，我们可以得到不同介质中声波传播速度的测量结果。
四、实验结果与分析根据我们的实验数据，我们得到了不同介质中声波传播速度的结果。
通过对比实验数据与理论预测值，我们发现实验结果与理论预测值基本一致，说明我们的实验方法是可行的。
测量数据处理实验报告
测量数据处理实验报告一、实验目的1.熟练运用一些数据处理软件（如Excel、Origin等）进行数据处理和绘图。
2.了解拟合方法及其在数据处理中的应用。
二、实验原理1.数据处理软件：数据处理软件是对实验数据进行处理和分析的重要工具，能够对数据进行各种运算和处理，例如均值、标准差、拟合等。
其中，Excel和Origin是常用的两种数据处理软件。
Excel是微软公司开发的电子表格软件，具有强大的计算和图表绘制功能，既方便又易于使用；Origin是常用于科学数据分析和绘图的专业软件，其功能包括数据处理、统计分析、绘图等，可满足科研工作的需求。
2.拟合方法：实验数据通常具有一定的规律性，而拟合方法就是通过某些数学模型，将实验数据“拟合”成一条直线、曲线或其他形式的函数表达式，以达到对数据的分析和评价的目的。
常用的拟合方法包括线性拟合、非线性拟合、最小二乘法等。
三、实验步骤1.分析实验数据，确定需要进行的计算和绘图类型。
3.根据需要计算实验数据的均值、标准差等统计量。
4.绘制数据的直方图、散点图等图形，以观察数据的分布规律。
5.进行拟合分析，选择适当的拟合方法（如线性拟合、非线性拟合），得到拟合曲线的表达式。
6.通过计算拟合优度等指标，评价拟合效果。
四、实验结果本实验选择使用Excel软件进行数据处理和绘图。
1. 输入数据：通过手动输入实验数据，得到如下数据表格。
2. 计算统计量：通过Excel的公式功能，可以方便地计算出各种统计量。
例如，输入“=AVERAGE(A2:A11)”可计算出数据的均值；输入“=STDEV(A2:A11)”可计算数据的标准差。
计算结果如下表所示。
3. 绘制图表：在Excel中，绘制图表的方法有多种，包括线图、柱状图、散点图等多种类型。
本实验选择绘制散点图和直方图。
散点图可以直观地反映实验数据的分布情况，而直方图则可以更清晰地展示出数据的分布规律。
通过Excel软件，绘制了如下散点图和直方图。
声速测量实验报告数据处理
声速测量实验报告数据处理实验报告：声速测量实验报告数据处理实验目的：1. 通过测量空气中声音在不同温度下的传播速度，了解声速与温度的关系；2. 通过数据处理和分析，掌握实验中常见数据处理方法。
实验原理：声速测量实验采用单向传播法，即利用一定距离内声波的扩散来测定声速。
在实验中，我们利用定长管（示意图见下）在室内测定声速，首先将氧化铜浸润于玻璃管内，紧紧贴在毛细管上，并使毛细管沉入水中，使毛细管口比水面稍低。
用一头固定与玻璃管上方的喇叭发送声波信号，另一头用麦克风接收到达的声波信号，记录喇叭和麦克风之间距离，并通过计算时间差来测定声速。
实验步骤：1. 按如上原理将实验装置搭建好，注意调整喇叭和麦克风之间的距离和位置，使其尽量接近玻璃管中心。
2. 先使用室温下测量声速，记录测量数据。
3. 然后，改变室温，测量不同温度下声速的变化。
分别记录测量数据，并且注意保持实验装置不变。
4. 完成测量后，计算和分析数据，绘制声速随温度变化的曲线。
数据计算和处理：1. 初始化在第一步中，使用测量设备记录了音波的通过时间，并将数据存储在不同的数组中。
对于空气，由于焓是一致的，所以方程式可以这样写：v = 343m/s （室温下的声速）2. 数据的转换对于数据进行简单的转换，注意峰值和峰谷之间的距离。
3. 计算根据测量数据和数据计算公式得到声速随温度变化的曲线。
我们运用了Mat lab来绘制数据图。
实验结果与分析在三种不同温度下，我们记录了空气中声音通过定长管的时间差：$\Delta t_1$ = 1.57ms $\Delta t_2$ = 1.7ms $\Delta t_3$ = 1.8ms根据上表数据可得，声波在低温下传播较快，高温时传播较慢。
此与流体密度和温度相关。
和理论值v=331.4+0.6t(m/s)相对，我们的实验结果的误差很小。
结论：本实验采用定长管单向传播法测定气体中声速的方法，通过测量来得到声速与温度的关系。
声速测量实验报告数据处理
一、实验目的1. 掌握声速测量的基本原理和方法；2. 了解声速与介质参数的关系；3. 学会使用逐差法进行数据处理。
二、实验原理声速是指声波在介质中传播的速度。
声速的测量方法有多种，本实验采用共振干涉法、相位比较法和时差法进行测量。
1. 共振干涉法：利用声波的干涉现象，通过测量相邻波腹或波节的距离，计算声速。
2. 相位比较法：通过比较声波传播过程中接收器接收到的信号与发射器激励电信号的相位差，计算声速。
3. 时差法：测量声波传播的距离和时间，根据公式计算声速。
三、实验仪器与材料1. 实验仪器：超声波发射器、超声波接收器、示波器、函数信号发生器、游标卡尺、温度计、湿度计等。
2. 实验材料：空气、实验数据表格。
四、实验步骤1. 共振干涉法：调整超声波发射器与接收器之间的距离，使接收器接收到的声波与发射器发出的声波发生干涉。
观察示波器上的波形，当出现相邻波腹或波节时，记录游标卡尺测得的距离L。
2. 相位比较法：调整超声波发射器与接收器之间的距离，使接收器接收到的信号与发射器激励电信号的相位差为0。
观察示波器上的波形，记录此时游标卡尺测得的距离L。
3. 时差法：调整超声波发射器与接收器之间的距离，记录声波传播的时间t。
根据公式v = L/t计算声速。
五、数据处理1. 共振干涉法：计算相邻波腹或波节的距离L的平均值，根据公式v = λf计算声速，其中λ为波长，f为频率。
2. 相位比较法：计算相位差为0时的距离L，根据公式v = λf计算声速。
3. 时差法：计算声波传播的距离L和时间t的平均值，根据公式v = L/t计算声速。
六、实验结果与分析1. 共振干涉法：测量得到相邻波腹或波节的距离L的平均值为L1，根据公式v = λf计算声速v1。
2. 相位比较法：测量得到相位差为0时的距离L的平均值为L2，根据公式v =λf计算声速v2。
3. 时差法：测量得到声波传播的距离L和时间t的平均值为L3和t3，根据公式v = L/t计算声速v3。
测量误差与数据处理实验报告
测量误差与数据处理实验报告实验报告格式：
标题：测量误差与数据处理实验报告
摘要：本实验旨在探究测量误差的来源及其处理方法，通过自己设计的实验进行数据采集与处理，最后得出结论并分析误差的影响。
实验结果表明，合理控制误差和精准处理数据非常重要。
1. 实验目的：
通过自己设计的实验了解测量误差的来源和处理方法，掌握精度等基本概念。
2. 实验步骤：
(1) 设计实验：以电容为例，设计了“通过变化距离来测量电容的实验”。
(2) 组装仪器：根据实验设计，组装了测量电容的仪器。
(3) 测量数据：对实验进行了多次测量，得到了电容的测量值。
(4) 数据处理：使用 Excel 等工具处理数据，计算出各项指标和
误差范围，并进行精度等级划分。
3. 实验结果：
(1) 根据数据处理结果，得到平均电容值为3.5μF，标准差为
0.2μF。
(2) 通过进行误差分析，可知测量误差来源主要包括仪器本身
误差、环境因素干扰和人为误差等多方面因素。
(3) 在误差控制和数据处理方面可采用实验平均法、精度等级
标准等方法。
4. 实验结论：
通过本实验的设计和数据处理，在实验中了解了测量误差的来源和处理方法，识别出了各方面因素影响到精度结果的准确性。
同时也提醒了我们在进行实验操作时需严格控制误差，避免产生干扰和误差现象，最终希望以此为基础，提高本人的实验操作、数据分析和综合思考能力。
声速的测量实验报告及数据处理
声速的测量实验报告及数据处理嘿伙计们，今天我们要来聊聊声速的测量实验报告及数据处理。
咱们得明白声速是什么吧？声速就是声音在空气中传播的速度，换句话说，就是我们听到的声音传到别人耳朵里需要多长时间。
好了，不多说了，让我们开始实验吧！实验目的：测量实验室内不同温度下的声音传播速度。
实验器材：麦克风、计时器、温度计、声速计、温度计。
实验步骤：1. 我们需要准备好实验器材。
把麦克风插上电源，打开开关，然后用计时器记录下从发出声音到接收到回声所需的时间。
用温度计测量实验室内的温度。
2. 接下来，我们要把声速计调整到合适的范围。
一般来说，声速计的量程是0-3499米/秒。
不过，我们这次实验的目的是测量不同温度下的声音传播速度，所以我们要把声速计调整到0-343米/秒这个范围内。
这样一来，我们就可以更准确地测量出声音在空气中传播的速度了。
3. 现在，我们可以开始实验了。
先让麦克风发出一个响亮的声音，然后用计时器记录下从发出声音到接收到回声所需的时间。
用温度计测量实验室内的温度。
重复这个过程几次，取平均值作为结果。
4. 把测得的时间和温度代入公式：声速 = (2 * 时间) / 温度，计算出声音在空气中传播的速度。
注意，这里的时间单位是秒，温度单位是摄氏度。
5. 我们可以把测得的结果整理成表格或图表的形式进行展示和分析。
通过对比不同温度下的声音传播速度，我们可以了解到什么因素会影响声音在空气中的传播速度。
好啦，实验完成啦！下面我们来分析一下实验数据。
根据我们的实验数据，我们发现随着温度的升高，声音在空气中传播的速度确实会变慢。
这是因为温度升高会导致空气分子的运动变得更加剧烈，从而使声音在空气中传播时受到更大的阻力。
所以呢，当我们感觉天气越来越热的时候，就会觉得声音变得“聒噪”了。
通过这次声速的测量实验报告及数据处理，我们不仅学到了如何测量声音在空气中传播的速度，还了解到了温度对声音传播速度的影响。
希望这些知识能帮助大家更好地理解我们周围的世界哦！。
声速测量实验报告数据处理
声速测量实验报告数据处理声速测量实验报告数据处理引言：声速是声波在介质中传播的速度，是声学研究中的重要参数之一。
本文将对声速测量实验中所得到的数据进行处理和分析，以得出准确的声速数值，并探讨实验中可能存在的误差来源和解决方法。
一、实验原理声速测量实验通常采用了经典的迈克尔逊干涉仪原理。
在实验中，将激光束分为两束，一束经过空气，另一束经过声速介质，然后再次汇聚，通过干涉现象来测量声速。
根据光程差和时间差的关系，可以计算出声速。
二、实验步骤1. 搭建迈克尔逊干涉仪实验装置，保证光路稳定和干涉现象清晰可见。
2. 在测量前，对实验装置进行校准，确保激光束的分束比例和光程差的准确性。
3. 将待测声速介质放置于其中一个光程中，注意保持介质的温度和压力稳定。
4. 调节迈克尔逊干涉仪的反射镜位置，使得干涉条纹清晰可见。
5. 记录下光程差和时间差的数值，进行多次测量以提高准确性。
三、数据处理1. 光程差与时间差的关系根据迈克尔逊干涉仪原理，光程差与时间差之间存在线性关系。
通过绘制光程差与时间差的图像，可以得到一条直线，斜率即为声速的倒数。
2. 数据拟合对实验测得的光程差和时间差数据进行拟合，可以得到最佳拟合直线。
利用拟合直线的斜率，即可计算出声速的数值。
3. 数据分析对实验测得的声速数值进行统计分析，计算平均值和标准差，以评估实验结果的准确性和可靠性。
4. 误差来源与解决方法（1）温度和压力的变化会导致声速介质的性质发生变化，进而影响测量结果。
解决方法是在实验过程中保持介质的温度和压力稳定。
（2）光程差的测量误差会直接影响声速的计算结果。
解决方法是使用高精度的测量仪器，并进行多次测量取平均值，以提高准确性。
（3）干涉条纹的清晰度会影响光程差的测量精度。
解决方法是调整反射镜的位置，使得干涉条纹清晰可见。
四、实验结果与讨论通过对实验数据的处理和分析，得到了声速的准确数值。
在实验中，我们得到了多组数据，并计算了平均值和标准差。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
院系____________________ 班级____________________ 姓名____________________ 学号____________________目录（页码根据自己需要填写）一、实验项目须知.....................................二、实验一秩亏自由网上机实验........................三、实验二极大验后滤波、推估上机实验...............四、实验三最小二乘配置上机实验.....................五、实验四平差系统的统计假设检验上机实验...........六、编写实验报告要求................................实验项目须知一、课程介绍本课程是测绘工程类专业的专业选修课程之一，是在学生学习了误差理论与测量平差知识之后又开设的一门专业的数据处理课程，该课程主要讲授近现代测量数据处理理论、模型和方法，同时紧密结合计算机编程，从而提高学生的计算机应用能力、测量数据处理能力，培养学生养成严谨的科学态度，提高分析和解决问题的能力，更好的适应今后工作和学习。
二、测量数据处理实验教学要求1.掌握稳健估计方法处理测量数据；2.掌握极大验后滤波、推估理论和方法处理测量数据；3.初步掌握最小二乘配置理论和方法处理测量数据；4.初步掌握工程或科研项目中测量数据处理的理论、步骤、方法和技巧；5.掌握测量数据处理中常用的计算机语言和软件；6.服从实验指导教师的规定，认真、按时、独立完成任务每次实验结束，提交书写工整或按照规范打印实验报告。
7.在实验过程中，还应遵守纪律，爱护实验室设备、离开后将自己的随身物品带走，并记录仪器使用状况、将凳子整理好并放到电脑桌的下方。
8.上机实验后，提交实验报告，成绩评定重点考核算法、处理结果与结果分析。
三、实验选择要求可根据课堂讲课具体内容，完成4个实验，并要求写出实验报告。
实验一 秩亏自由网上机实验一、 实验目的1. 掌握各种秩亏自由网平差的函数模型和随机模型；2. 熟悉各种秩亏自由网平差解算的若干公式；3. 结合算例能利用matlab 进行编程计算。
二、实验内容在下图水准网中，观测高差、距离和各待定点高程近似值见教材P43表2-1，分别进行下列自由网平差：（1）以6号点为固定点的经典自由网平差； （2）以重心基准的自由网平差（P χ=E ）；(3) 以1，2，5，6四个点为拟稳基准的拟稳平差.参考程序：（1）首先建立误差方程 l x B V -=∧B=[-1 1 0 0 0 ;-1 0 0 1 0;-1 0 0 0 1 ;0 -1 1 0 0 ;0 -1 0 1 0 ... ;0 -1 0 0 0 ;0 0 -1 1 0 ;0 0 -1 0 1 ;0 0 0 0 -1 ]; l=[-1;4;5;8;16;2;-13;-20;-1];S=[204.1 188.7 344.8 149.2 142.9 250 128.2 98 196.1];P=diag(10./S);N=B'*P*B;x=inv(N)*B'*P*l V=B*x-l Qx=inv(N)(2) 首先建立误差方程 l x B V -=∧B=[-1 1 0 0 0 0;-1 0 0 1 0 0;-1 0 0 0 1 0;0 -1 1 0 0 0;0 -1 0 1 0 0 ... ;0 -1 0 0 0 1;0 0 -1 1 0 0;0 0 -1 0 1 0;0 0 0 0 -1 1];l=[-1;4;5;8;16;2;-13;-20;-1]S=[204.1 188.7 344.8 149.2 142.9 250 128.2 98 196.1];P=diag(10./S); S1=[1;1;1;1;1;1]; W=B'*P*l;Xr=inv(B'*P*B+S1*S1')*WV=B*Xz-lQxr=inv(B'*P*B+S1*S1')-inv(B'*P*B+S1*S1')*S1*S1’*inv(B'*P*B+S1*S1')(3) 首先建立误差方程 l x B V -=∧B=[-1 1 0 0 0 0;-1 0 0 1 0 0;-1 0 0 0 1 0;0 -1 1 0 0 0;0 -1 0 1 0 0 ...;0 -1 0 0 0 1;0 0 -1 1 0 0;0 0 -1 0 1 0;0 0 0 0 -1 1];l=[-1;4;5;8;16;2;-13;-20;-1]S=[1;1;1;1;1;1];S1=[204.1 188.7 344.8 149.2 142.9 250 128.2 98 196.1];P=diag(10./S1); N=B'*P*B W=B ’*P*l; px=[1 1 0 0 1 1] Px=diag(px) Ss=Px*SXs=inv(B'*P*B+Ss*Ss')*WV=B*Xs-lQxs=inv(B'*P*B+Ss*Ss')*N*inv(B'*P*B+Ss*Ss')三、 实验心得和体会实验二 极大验后滤波、推估上机实验 一、 实验目的和要求1. 掌握极大验后滤波与推估的基本原理和计算公式；2. 能根据观测值和给定的数学模型计算滤波信号的协方差矩阵的估计矩阵；3. 能使用matlab 进行有关的计算；二、 实验内容:1.（基本要求题目）沿A 、B 连线在A 、1、2、3、4等五个点测定了大气温度，获得18组观测值ij x ,其结果见下表，已知各点的距离）、、、4321( j S Aj 分别为4.511,10.747,16.753,22.220（km ）.假定AB 连线上的大气温度，是一个以各点至A 的距离S 为因素的平稳随机过程，试估计此随机过程的协方差函数。
参考程序：X=[16.4 15.9 16.3 15.8 15.4;16.3 16.1 16.2 15.9 15.5;16.8 16.4 15.9 16.3 16.3; ...16.6 16.2 16.3 16.2 16.3;16.8 16.7 16.1 16.2 16.3; ...16.6 16.6 16.4 16.2 16.3;16.7 16.6 17.0 16.3 16.2; ...16.8 16.5 16.9 16.5 15.9;16.9 16.4 16.1 16.6 16.2; ...16.1 16.6 16.2 17.1 17.0;16.5 16.7 16.7 16.8 16.5; ...17.1 16.9 16.6 16.8 16.7;16.8 16.9 16.8 17.1 17.0; ...17.0 16.9 16.7 16.9 17.0;17.1 16.9 16.7 16.9 17.0; ...16.2 16.0 16.3 16.5 16.4;16.1 15.9 16.2 16.4 16.4; ...16.1 15.9 16.0 16.1 16.5;16.3 16.0 16.5 16.6 17.0;16.1 16.0 15.9 16.5 16.3]Dx1=cov(X)s1=[0;4.510;10.746;16.754;22.230]A=[ones(5,1),s1,s1.^2,s1.^3,s1.^4]s2=[0;6.236;12.242;17.709]B=[ones(4,1),s2,s2.^2,s2.^3,s2.^4]s3=[0;6.006;11.473]C=[ones(3,1),s3,s3.^2,s3.^3,s3.^4]s4=[0;5.467]D=[ones(2,1),s4,s4.^2,s4.^3,s4.^4]E=[1 0 0 0 0]Y=[A;B;C;D;E]L=[0.1168;0.0958;0.0506;0.0466;0.0394;0.1362;0.0705;0.1023;0.1047;0.1073;0.0567;0.0346; ...0.1529;0.1561;0.2203]b=inv(Y’*Y)*Y’*Lb=[b(1),b(2)*10,b(3)*100,b(4)*1000,b(5)*10000] %将km单位化为10km单位S=[0 4.511 10.747 16.753 22.220;4.511 0 6.236 12.242 17.709;10.747 6.236 0 6.006 11.473; ... 16.753 12.242 6.006 0 5.467;22.220 17.709 11.473 5.467 0]b=fliplr(b)for i=1:25S(i)=polyval(b,S(i)/10);endS %即为滤波信号协方差阵的估计矩阵Warning: Unable to create personal MATLAB work folder:C:\yonghu\????\MATLAB Warning: Userpath must be an absolute path and must exist on disk.>> X=[16.4 15.9 16.3 15.8 15.4;16.3 16.1 16.2 15.9 15.5;16.8 16.4 15.9 16.3 16.3; ...16.6 16.2 16.3 16.2 16.3;16.8 16.7 16.1 16.2 16.3; ...16.6 16.6 16.4 16.2 16.3;16.7 16.6 17.0 16.3 16.2; ...16.8 16.5 16.9 16.5 15.9;16.9 16.4 16.1 16.6 16.2; ...16.1 16.6 16.2 17.1 17.0;16.5 16.7 16.7 16.8 16.5; ...17.1 16.9 16.6 16.8 16.7;16.8 16.9 16.8 17.1 17.0; ...17.0 16.9 16.7 16.9 17.0;17.1 16.9 16.7 16.9 17.0; ...16.2 16.0 16.3 16.5 16.4;16.1 15.9 16.2 16.4 16.4; ...16.1 15.9 16.0 16.1 16.5;16.3 16.0 16.5 16.6 17.0;16.1 16.0 15.9 16.5 16.3] X =16.4000 15.9000 16.3000 15.8000 15.400016.3000 16.1000 16.2000 15.9000 15.500016.8000 16.4000 15.9000 16.3000 16.300016.6000 16.2000 16.3000 16.2000 16.300016.8000 16.7000 16.1000 16.2000 16.300016.6000 16.6000 16.4000 16.2000 16.300016.7000 16.6000 17.0000 16.3000 16.200016.8000 16.5000 16.9000 16.5000 15.900016.9000 16.4000 16.1000 16.6000 16.200016.1000 16.6000 16.2000 17.1000 17.000016.5000 16.7000 16.7000 16.8000 16.500017.1000 16.9000 16.6000 16.8000 16.700016.8000 16.9000 16.8000 17.1000 17.000017.0000 16.9000 16.7000 16.9000 17.000017.1000 16.9000 16.7000 16.9000 17.000016.2000 16.0000 16.3000 16.5000 16.400016.1000 15.9000 16.2000 16.4000 16.400016.1000 15.9000 16.0000 16.1000 16.500016.3000 16.0000 16.5000 16.6000 17.000016.1000 16.0000 15.9000 16.5000 16.3000>> Dx1=cov(X)Dx1 =0.1203 0.1023 0.0559 0.0389 0.03040.1023 0.1405 0.0716 0.0880 0.08420.0559 0.0716 0.1094 0.0504 0.03690.0389 0.0880 0.0504 0.1371 0.14280.0304 0.0842 0.0369 0.1428 0.2157>> s1=[0;4.510;10.746;16.754;22.230]s1 =4.510010.746016.754022.2300>> A=[ones(5,1),s1,s1.^2,s1.^3,s1.^4]A =1.0e+05 *0.0000 0 0 0 0 0.0000 0.0000 0.0002 0.0009 0.0041 0.0000 0.0001 0.0012 0.0124 0.1333 0.0000 0.0002 0.0028 0.0470 0.7879 0.0000 0.0002 0.0049 0.1099 2.4421>> s2=[0;6.236;12.242;17.709]s2 =6.236012.242017.7090>> B=[ones(4,1),s2,s2.^2,s2.^3,s2.^4]B =1.0e+04 *0.0001 0 0 0 0 0.0001 0.0006 0.0039 0.0243 0.1512 0.0001 0.0012 0.0150 0.1835 2.2460 0.0001 0.0018 0.0314 0.5554 9.8350>> s3=[0;6.006;11.473]s3 =6.006011.4730>> C=[ones(3,1),s3,s3.^2,s3.^3,s3.^4]C =1.0e+04 *0.0001 0 0 0 0 0.0001 0.0006 0.0036 0.0217 0.1301 0.0001 0.0011 0.0132 0.1510 1.7326>> s4=[0;5.467]s4 =5.4670>> D=[ones(2,1),s4,s4.^2,s4.^3,s4.^4]D =1.0000 0 0 0 0 1.0000 5.4670 29.8881 163.3982 893.2979>> E=[1 0 0 0 0]Y=[A;B;C;D;E]E =1 0 0 0 0Y =1.0e+05 *0.0000 0 0 0 0 0.0000 0.0000 0.0002 0.0009 0.0041 0.0000 0.0001 0.0012 0.0124 0.1333 0.0000 0.0002 0.0028 0.0470 0.7879 0.0000 0.0002 0.0049 0.1099 2.4421 0.0000 0 0 0 0 0.0000 0.0001 0.0004 0.0024 0.0151 0.0000 0.0001 0.0015 0.0183 0.22460.0000 0.0002 0.0031 0.0555 0.98350.0000 0 0 0 00.0000 0.0001 0.0004 0.0022 0.01300.0000 0.0001 0.0013 0.0151 0.17330.0000 0 0 0 00.0000 0.0001 0.0003 0.0016 0.00890.0000 0 0 0 0>>L=[0.1168;0.0958;0.0506;0.0466;0.0394;0.1362;0.0705;0.1023;0.1047;0.1073;0.0567 ;0.0346; ...0.1529;0.1561;0.2203]b=inv(Y’*Y)*Y’*LL =0.11680.09580.05060.04660.03940.13620.07050.10230.10470.10730.05670.03460.15290.15610.2203b=inv(Y’*Y)*Y’*L|Error: The input character is not valid in MATLAB statements or expressions. >> b=inv(Y'*Y)*Y'*Lb =0.1468-0.0008-0.00290.0003-0.0000>> b=[b(1),b(2)*10,b(3)*100,b(4)*1000,b(5)*10000]S=[0 4.511 10.747 16.753 22.220;4.511 0 6.236 12.242 17.709;10.747 6.236 0 6.006 11.473; ...16.753 12.242 6.006 0 5.467;22.220 17.709 11.473 5.467 0] b =0.1468 -0.0080 -0.2860 0.2774 -0.0706 S =0 4.5110 10.7470 16.7530 22.2200 4.5110 0 6.2360 12.2420 17.7090 10.7470 6.2360 0 6.0060 11.4730 16.7530 12.2420 6.0060 0 5.4670 22.2200 17.7090 11.4730 5.4670 0>> b=fliplr(b) for i=1:25S(i)=polyval(b,S(i)/10); end S b =-0.0706 0.2774 -0.2860 -0.0080 0.1468 S =0.1468 0.1076 0.0581 0.0792 0.0400 0.1076 0.1468 0.0872 0.0589 0.0823 0.0581 0.0872 0.1468 0.0898 0.0578 0.0792 0.0589 0.0898 0.1468 0.0960 0.0400 0.0823 0.0578 0.0960 0.1468 >>2.同第1题，已知,00.23km S AB =测得A 、1、2、3、4五个测站上的大气温度为T L )8.19,2.19,7.18,1.20,0.19(00000=利用推估公式试推估B点大气温度Bt。