[½»Á÷] ¡¾ÇóÖú¡¿´ËVB³ÌÐòÔõôÔËÐв»ÆðÀ´£¿ ÒÑÓÐ4È˲ÎÓë

Private Sub Command1_Click()
'PROGRAM DIR7
'Driver program for routine TOEPLZ
N = 5
N2 = 2 * N
Dim X(5), Y(5), R(10)
'ÊäÈëÒÑÖª·½³ÌϵÊý¾ØÕóµÄR
For I = 1 To 2 * N - 1
R(I) = 1 / I
Next I
'ÊäÈëÒÑÖª·½³ÌµÄÓÒ¶ËÏòÁ¿
For I = 1 To N
   Y(I) = 0.1 * I
Next I
Print Tab(5); "ÒÑÖªµÄ·½³Ì×éÓÒ¶ËÏòÁ¿"
Print Tab(12); Format$(Y(1), "##.##00")
Print Tab(12); Format$(Y(2), "##.##00")
Print Tab(12); Format$(Y(3), "##.##00")
Print Tab(12); Format$(Y(4), "##.##00")
Print Tab(12); Format$(Y(5), "##.##00")
Call TOEPLZ(R(), X(), Y(), N)
'Êä³ö·½³Ì×éµÄ½âX
Print
Print Tab(5); "¼ÆËã³ö·½³Ì×éµÄ½â"
Print Tab(9); Format$(X(1), "#####0e+00")
Print Tab(9); Format$(X(2), "#####0e+00")
Print Tab(9); Format$(X(3), "#####0e+00")
Print Tab(9); Format$(X(4), "#####0e+00")
Print Tab(9); Format$(X(5), "#####0e+00")
'½«¼ÆËã³öµÄ½âX³ËÒÔϵÊý¾ØÕó£¬ÒÔÑéÖ¤¼ÆËã½á¹ûÕýÈ·
Print
Print Tab(4); "½«¼ÆËã³öµÄ½â³ËÒÔϵÊý¾ØÕó£¬ÒÔÑéÖ¤¼ÆËã½á¹ûÕýÈ·"
Print
Print Tab(5); "½â³ËÒÔϵÊý¾ØÕó"; Tab(24); "·½³ÌµÄÓÒ¶ËÏòÁ¿"
For I = 1 To N
    SUM1 = 0
    For J = 1 To N
     SUM1 = SUM1 + R(N + 1 - J) * X(J)
    Next J
    Print Tab(10); Format$(SUM1, "#.####0"); Tab(30); Format$(Y(I), "#.###0")
   Next I
End Sub
Sub TOEPLZ(R(), X(), Y(), N)
NMAX = 100
Dim G(100), H(100)
If R(N) = 0# Then GoTo 99
X(1) = Y(1) / R(N)
If N = 1 Then Exit Sub
G(1) = R(N - 1) / R(N)
H(1) = R(N + 1) / R(N)
For M = 1 To N
    M1 = M + 1
    SXN = -Y(M1)
    SD = -R(N)
    For J = 1 To M
    SXN = SXN + R(N + M1 - J) * X(J)
    SD = SD + R(N + M1 - J) * G(M - J + 1)
   Next J
   If SD = 0# Then GoTo 99
   X(M1) = SXN / SD
   For J = 1 To M
       X(J) = X(J) - X(M1) * G(M - J + 1)
   Next J
   If M1 = N Then Exit Sub
   SGN1 = -R(N - M1)
   SHN = -R(N + M1)
   SGD = -R(N)
   For J = 1 To N
       SGN1 = SGN1 + R(N + J - M1) * G(J)
       SHN = SHN + R(N + M1 - J) * H(J)
       SGD = SGD + R(N + J - M1) * H(M - J + 1)
   Next J
   If SD = 0# Or SGD = 0# Then GoTo 99
   G(M1) = SGN1 / SGD
   H(M1) = SHN / SD
   K = M
   M2 = (M + 1) / 2
   PP = G(M1)
   QQ = H(M1)
   For J = 1 To M2
       PT1 = G(J)
       PT2 = G(K)
       QT1 = H(J)
       QT2 = H(K)
       G(J) = PT1 - PP * QT2
       G(K) = PT2 - PP * QT1
       H(J) = QT1 - QQ * PT2
       H(K) = QT2 - QQ * PT1
       K = K - 1
     Next J
    Next M
   Print "Never get here"
99    Print "Levinson method fails:singular principal minor"
End Sub
   
   

Originally posted by stephenliu89 at 2010-08-29 15:59:05:
Èý·Ö±à³Ì£¬Æß·Öµ÷ÊÔ¡£µ÷ÊÔµÄʱ¼ä»á±È±à³ÌµÄʱ¼ä³¤ºÜ¶à£¬¶øÇÒÓÐʱºò´íÎóºÜÄÑÕÒµ½£¬Â¥Ö÷²»ÒªÐ¹Æø£¬¼ÓÓÍ£¡