基于傅立葉級數(shù)牛頓外插法的經(jīng)濟參數(shù)預(yù)測-資料下載頁
2025-06-18 18:00本頁面
【正文】 uspP=1DO WHILE .T. amp。amp。 選擇權(quán)重系數(shù) amp。amp。 加一控件，生成權(quán)重系數(shù) amp。amp。 K=0 FOR I=1 TO N K=K+P(I) NEXT FOR J=2 TO M2 FRI(N+1,J)=0 FOR I=1 TO N FRI(N+1,J)=FRI(N+1,J)+FRI(I,J)*P(I) amp。amp。 各年份傅立葉系數(shù)乘權(quán)重系數(shù) NEXT FRI(N+1,J)=FRI(N+1,J)/K NEXT FRI(N+1,1)=1 amp。amp。平均值a0置為1 Q=1 S=0 S(1)=N FOR I=2 TO 7 amp。amp。 正則方程系數(shù) FOR J=1 TO N Q(J)=Q(J)*J S(I)=S(I)+Q(J) NEXT NEXT T=0 FOR J=1 TO N amp。amp。 正則方程常數(shù)項 T(1)=T(1)+FRI(J,1) T(2)=T(2)+FRI(J,1)*J T(3)=T(3)+FRI(J,1)*J*J T(4)=T(4)+FRI(J,1)*J*J*J NEXT FOR L=2 TO 4 amp。amp。 外插值計算，T(5)=線性插值結(jié)果, T(6)=拋物線插值結(jié)果, T(7)=三次插值結(jié)果 FOR I=1 to L amp。amp。 生成正則方程，將Si,Ti寫入數(shù)組A矩陣 FOR J=1 TO L A(I,J)=S(I+J1) NEXT A(I,L+1)=T(I) NEXT LineEqua(@A,L) amp。amp。 用最小二乘法求曲線擬合結(jié)果 T(L+3)=0 amp。amp。 計算外插值,分別用1,2,3次函數(shù)預(yù)測第 N+1 年的平均值 K=1 FOR J=1 TO L T(5)=T(5)+K*A(J,L+1) K=K*(N+1) NEXT NEXT FOR J=1 TO M2 FRI(N+1,J)=FRI(N+1,J)*T(5)amp。amp。 傅立葉系數(shù)乘預(yù)測的平均值 NEXT FOR J=1 TO M1 amp。amp。 做傅立葉逆變換，生成各月值 C(N+1,J)=FRI(N+1,1) FOR I=1 TO MC amp。amp。 計算COS 項 C(N+1,J)=C(N+1,J)+FRI(N+1,I+1)*COS(I*X(J)/M1*PI2) NEXT FOR I=1 TO MS C(N+1,J)=C(N+1,J)+FRI(N+1,I+MC+1)*SIN(I*X(J)/M1*PI2) NEXT NEXT FOR J=1 TO M1 C(N+2,J)=C(N+1,J)*T(6)/T(5) amp。amp。二次預(yù)測值按一次預(yù)測值做比例計算 C(N+3,J)=C(N+1,J)*T(7)/T(5) amp。amp。三次.......... NEXT W= WAIT 擬合完成? Y/N TO W IF UPPER(W)=Y EXIT ENDIFENDDOUSE TEMPFOR I=1 TO 3 F(1)=YEAR1(N+I) FOR J=1 TO M1 F(J+1)=C(I+N,J) NEXT APPE BLAN GATH FROM FNEXTF=0FOR J=2 TO M1 F(J)=P(J1)NEXTAPPE BLANGATH FROM FRETURNPROCEDURE PR1FUNC LineEqua(C,N) amp。amp。 線性方程組求解函數(shù)，高斯主元素法,N 為未知數(shù)個數(shù)，解在第 N+1 列返回PRIVATE I,J,K,M,TFOR M=1 TO N T=ABS(C(M,M)) K=M FOR I=M+1 TO N amp。amp。 尋找主元素 IF TABS(C(I,M)) T=ABS(C(I,M)) K=I ENDIF NEXT IF MK amp。amp。 M列發(fā)現(xiàn)絕對值更大的元素，應(yīng)交換到主元素行 FOR J=M TO N+1 T=C(M,J) C(M,J)=C(K,J) C(K,J)=T NEXT ENDIF IF ABS(C(M,M))1E7 amp。amp。 方程無解，返回，N=0 為標(biāo)志 N=0 C(1,1)=1E20 RETURN ENDIF FOR I=M+1 TO N amp。amp。 消元 T=C(I,M)/C(M,M) FOR J=M+1 TO N+1 C(I,J)=C(I,J)+T*C(M,J) NEXT NEXTNEXTC(N,N+1)=C(N,N+1)/C(N,N) amp。amp。 得到第 N 個 未知數(shù)的值FOR I=N1 TO 1 STEP 1 amp。amp。 回代 FOR J=I+1 TO N C(I,N+1)=C(I,N+1)C(I,J)*C(J,N+1) NEXT C(I,N+1)=C(I,N+1)/C(I,I)NEXTRETURN第一章 引言 1167。 數(shù)學(xué)趨勢預(yù)測的目的及意義 1第二章 時間序列波動性特征和持續(xù)增長趨勢的述 2167。 波動特征及傅立葉級數(shù) 2167。 時間序列的持續(xù)增長趨勢與冪函數(shù) 3第三章 數(shù)學(xué)模型的建立 6167。 不等距節(jié)點的離散傅立葉變換 6167。 最小二乘法意義下的曲線擬合擬合求解 7167。 傅立葉系數(shù)的加權(quán)處理 9167。 外插值計算及擬合求解 11 Newton基函數(shù) 11 多項式擬合求解 12167。 線性方程組求解 14 高斯主元素消去法 14第四章 模型應(yīng)用 18第五章 總結(jié) 20致謝 21參考文獻(xiàn) 22附論文主程序 23
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1