《第五章线性参数的最小二乘处理.doc》由会员分享,可在线阅读,更多相关《第五章线性参数的最小二乘处理.doc(30页珍藏版)》请在金锄头文库上搜索。
1、第五章 线性参数的最小二乘处理习 题5-1研究铂-铱米尺基准器的膨胀系数时得出,在不同温度下该米尺基准器的长度的修正值可用下述公式表示:x+ty+t2z=L式中x表示在0时米尺基标准器的修正值(微米);y和z为温度系数;t为温度();L为t基准器的长度的修正值(微米)。经研究得在不同温度下米尺基准器长度的修正值如下表:i12345t17.2500.5515.36310.45914.277L160.285.7047.6141.49124.2567891017.80622.10324.63328.98634.417154.87192.64214.57252.09299.84求未知参量x,y,z的最
2、可依赖值。5-2对未知量x,y,z,组合测量的结果如下: x=0 y=0 z=0x-y=0.92,-y+x=1.35-x+z=1.00试求x,y,z的最可依赖值及其标准误差。5-3由等精度测定方程为:x+37y+1369z=36.3x+32y+1024z=41.4 x+27y+729z=47.5 x+2y+484z=54.7 x+17y+289z=63.2 x+12y+144z=72.9 x+7y+49z=83.7试用矩阵最小二乘法求x,y,z的最可依赖值及其精度。5-4交流电路的电抗x=L,在角频率1=3时,测得x为x1=0.82=2时,测得x为x2=0.23=1时,测得x为x3=-0.3试
3、求:(i) L,C及其方差;(ii) =3时(=0.1)电抗值及其方差。5-5试求下列方程给出的x,y的最大或然值及其标准误差。2x+y=5.1 x-y=1.14x-y=7.4x+4y=5.95-6测得一直线上四段长度AB、BC、CD、DE分别为24.1,35.8,30.3和33.8厘米,但已知AD准确长90厘米和BE准确长100厘米。试求AB,BC,CD,DE的最大或然值。5-7由方程组3x+y=2.9x-2y=0.92x-3y=1.9试求x,y的最大或然值及其标准误差。5-8由下面的不等精度的测定方程组,求x1,x2的最可信赖值及其标准误差。x1=0权: P1=8x2=0 P2=10x1+
4、2x2=0.25 P3=1x1-3x2=0.92 P4=55-9由下面的不等精度的测定方程组,试用矩阵最小二乘法求x,y的最大或然值及其标准误差。x-3y=-5.6权: P1=14x+y=8.1 P2=22x-y=0.5 P3=35-10由下面的测定方程组,试求x,y的最可依赖值及其标准误差。2x+y=5.1 权:P1=1x-y=1.1 P2=34x-y=7.2 P3=25-11试求满足下列方程的x,y,z及其标准误差(假设它们是等权的)。x+y+z=4.012x-y+z=1.04x+3y-2z=5.02 3x+y=4.975-12由座标点(1,0) (3,1) 和 (-1,2)到某点的距离分
5、别为3.1,2.2和3.2。试求该点座标位置的最大或然值及其标准误差。5-13对某一角度值,分两个测回进行测定,其权等于测定次数,测定值如下。试求该角度的最可信赖及其标准误差。第一测回第二测回piaipiai7345633455401345423455301345520134550234551345570134551013455505-14某平面三角形三个角被测出为A=48510,B=602524,C=70427,令假设这种测量(i)各次权相等;(ii)各次权分别为1、2、3;试求A、B、C的最大或然值。5-15数N系时间t的函数N=x1+ x2t+ x3t2测定后的N的值如下。测定是在异权情
6、况下进行的,试求x1,x2,x3的最可信赖值。i123456789ti1.51.10.70.3-0.1-0.5-1.0-1.5-2.0Ni6.203.452.001.802.404.558.8515.7024.40Pi0.7070.5005-16硝酸钠在100份水内的溶解度与温度的关系,测定为温度0410152129365168溶解度66.771.076.380.685.792.999.4113.6125.1上述关系可用直线67.5+0.87t表示(式中t为温度)。试用最小二乘法来检证。5-17由下列测定的方程组,求X、Y最可信赖及其或然误差。 X+Y=37.0 权:P1=5 2X+Y=61.
7、9 P2=43X+Y=86.7 P3=4X+2Y=49.2 P4=4X+3Y=60.6 P5=32X+3Y=86.7 P6=23X+2Y=98.4 P7=35-18由下列测定方程组,求X、Y最可信赖及其标准误差。 2X+4Y+8Z=0.16122.200X+4.840Y+10.648Z=0.19863.200X+10.240Y+32.768Z=0.50982.600X+6.760Y+17.576Z=0.28963X+9Y+27Z=0.41815-19假设有三个某种量规,其值分别为Y1、Y2、Y3。现在将它们直接地或间接地与数值已知为N的标准量规比较,比较的方案为下述三种(三种组合):(i)每一
8、个量规各与标准量规比较二次;(ii)第一个量规(Y1)与标准量规比较二次,第二个量规(Y2)与第一个量规比较二次,第三个量规(Y3)与第二量规比较二次;(iii)每一个量规各与标准量比较一次,然后它们相互按不同的组合比较一次;上述三种测量方案得到的条件方程式如下表所示:(1)(2)(3)Y1N=X1Y1N=X1Y1N=X1Y1N=X2Y1N=X2Y2N=X2Y2N=X3Y2Y1=X3Y3N=X3Y2N=X4Y2Y1=X4Y2Y1=X4Y3N=X5Y3Y2=X5Y3Y1=X5Y3N=X6Y3Y2=X6Y3Y2=X6试研究采用那一种测量方案能够获得最好的结果。(提示:可以比较不同测量方案下未知数的
9、权)。典型题解5-1 由测量方程 试求、的最小二乘法处理及其相应精度。解:方法一:列出误差方程组:分别对求偏导,并令它们的结果为0,即, 由上式可解得结果: 方法二:直接列表计算给出正规方程常数项和系数 1319132.98.72.921-214-20.90.9-1.832-349-61.93.8-5.7-1414-5-13.4-4.6可得正规方程将的结果代入分别求得:得,由题已知,得由不定乘数的方程组得得方法二:按矩阵形式计算,由误差方程 上式可以表示为即; ; 可得式中 所以即解得,将最佳估计值代入误差方程可得, 将计算得到的数据代入式中为求出估计量的标准差,首先求出不定常数 。由已知,不
10、定常数的系数与正规方程的系数相同,因而是矩阵中各元素,即则可得估计量的标准差为5-2 已知误差方程为 试给出的最小二乘法处理及其相应精度。解:根据矩阵形式,误差方程可以表示为即; ;可得式中 得即解得 将最佳估计值代入误差方程可得 得可得为求出估计量的标准差,首先求出不定乘数,不定乘数的系数与正规方程的系数相同,因而是矩阵中各元素,即则于是估计量的标准差, 5-3 测力计示值与测量时的温度的对应值独立测得如下表所示。/15182124273043.6143.6343.6843.7143.7443.78设无误差,值随的变化呈线性关系,试给出线性方程中系数和的最小二乘估计及其相应精度。解:方法一:
11、列出误差方程式,令为待估计量,则误差方程可写成为为计算方便,将数据列表如下:11522543.61654.1521832443.63785.3432144143.68917.2842457643.711049.0452772943.741180.9863090043.781313.41353195262.155900.19 根据误差方程,列出正规方程: 将表中计算出的相应系数值代入上面的正规方程得解得即将代入误差方程得:将代入上式,可得残余误差为:NNNNNN可得: 可得标准差为,N由上面所给的正规方程的系数,可列出求解不定乘数的方程组分别解得估计量的标准差为 方法二:直接利用矩阵求解,误差方程可写成即; ;可得式中所以将最佳估计值代入误差方程得可计算为求出估计量的标准差,需要求出不定乘数的系数,而不定乘数的系数与正规方程的系数相同,因而是矩阵中各元素,即则可得估计量的标准差为 5-4 研究米尺基准器的线膨胀系数,得出在不同温度时该基准器的长度修正值可用公式表示。式中为时米尺基准器的修正值(单位为):和为温度系数;为温度。在不同温度时米尺基准器的修正值如下表所示:/0.5515.36310.45914.27717.806