9实验数据及等距五点三次平滑.docVIP

901 实验数据的等距五点三次平滑 功能 用五点三次平滑公式., 对等距节点上的实验数据进行平滑,力争消除随机误差的影响, 提高实验数据的质量 . 算法简介 设 2n+1 个等距节点X,X,……X.X,X,….,X,X上的实验数据分别为Y,Y,….,Y,Y,Y,….,Y,Y.再设两节点间的等距为h,作交换t=则原节点变为t=-n ,t= -n+1,….., t= -1,t=0,t=1,….t= n-1,t=n, 我们用m次多项式来拟合得到的实验数据.设拟合多项式为 Y(t)=a+at+at+…+at 用最小二乘法来确定方程(t)中的待定系数,令, =[-Y ]=(a,a,…,a) 为使(a,a,…,a)达到最小.将它分为对a(k=0,1,….,m)求编导数,并令其为0,可得方程组 = (k=0.1,…,m) 该方程组称为正规方程组. 当n=2(5 ?个节点), m=3时, 得到一个具体得正规方程组, 由此解出, a,a,a,a 代入(1)式, 并令t=0,+1,-1+2,-2,得到五点三次平滑公式: =(69Y+4Y-6Y+4 Y-Y) (4-1) =(2 Y+27 Y+12 Y-8 Y+2 Y) (4-2) =(-3 Y+12 Y+17 Y+12 Y-3 Y) (4-3) =(2 Y-8 Y+12 Y+27 Y+2 Y) (4-4) =(- Y+4 Y-6 Y+4 Y+69 Y) (4-5) 其中为Y得改进值. 该公式要求节点个数为k5,当节点个数多于5时,为对称其见,除在两端分别用公式,(4-1).(4-2).(4-4),(4-5)外,其余都用公式(4-3)进行平滑,这就相当在每个子区间上用不同得三次最小二乘多项式进行平滑. 由上述推导知,对于等距节点,平滑公式中只用到实验数据Y,而和节点 X及节点间等距离h无关. 三.程序使用说明 输入参数 N 实验数据(X,Y)的对数 B(N)一维数组,存放N个实验数据y,i=1,2,….. ,n X(N) 一维数组,存放实验序号j,j=1,2,…..,n 2. 输出参数 X(N) 一维数组,存放实验序号 “OLD” 即原始实验数据B(N) “NEW” 经平滑后的数据,存放在一维数组A(N)中输出. ************************************************************************************************************************************************************** 第九章 数理统计 10 ******************************************* 20 * 901 实验数据的等距五点三次平滑********* 30 ******************************************* 40 NRUT “n=”; N 50 PRINT TAB(3); “例:” ; Atab(10); “n=”;N 60 IN A(N),B(N),X(N) 70 OR J=1 TO N 80 EAD B(J):X(J)=J 90 EXT J 100 GOSUB 300 110 PRINT TAB(2);”X(J)”;TAB(10);”OLD”;TAB(25):”NEW” 120 FOR J=1 TO N 130 LPRINT TAB(2);”X(J)”;TAB(8);B(J);TAB(23);A(J) 140 NEXT J 150 DATA 54,145,227,359,381,401,372,342,259,112,65 160 DATA 170 END 300 子程序 310 A(1)=(69*B(1)+4*B(2)+B(4))-6*B(3)-B(5))/70 320 A(2)=(2*(B(1)+B(5))+27*B(2)+12*B(3)-8*B(4))/35 330 FOR J=3 TO N-2 340 A(J)=(-3*(B(J-2)+B(J+2))+12*(B(J-1)+B(J+1))+17*B(J))/35 350 NEXT J 360 A(N-1)=(2*(B(N-4)+B(N))-8*B(N-3)+12*B(N-2)+27*B(N-1))/35 370 A(N)==(-B(N-4)+4*(B(N-3)+B(N-1))-6*B(N-2)+69*B(N))/70 380 RETURN 例: