24小时热门版块排行榜

>论坛更新日志 (3351)
>虫友互识 (795)
>休闲灌水 (117)
>导师招生 (116)
>文献求助 (99)
>博后之家 (96)
>基金申请 (95)
>论文投稿 (89)
>教师之家 (48)
>硕博家园 (44)
>药学 (29)
>考博 (28)
>招聘信息布告栏 (15)
>公派出国 (12)
>考研 (12)
>海外博后 (11)

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

zyj8119

木虫 (著名写手)

应助: 65 (初中生)
贵宾: 0.003
金币: 915.1
散金: 1440
红花: 35
帖子: 2936
在线: 1329.4小时
虫号: 664177
注册: 2008-11-29
性别: GG
专业: 理论和计算化学

[交流] 【求助】此VB程序怎么运行不起来？已有4人参与

CODE:

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