在线中文-在线在线ccc66-在线永久免费观看的毛片-在线永久免费观看的a站视频-久久精品最新免费国产成人-久久精品综合一区二区三区

以文本方式查看主題

-  曙海教育集團論壇  (http://www.bgl88.cn/bbs/index.asp)
--  VB語言  (http://www.bgl88.cn/bbs/list.asp?boardid=77)
----  vb語言實現的卡爾曼濾波源程序  (http://www.bgl88.cn/bbs/dispbbs.asp?boardid=77&id=2608)

--  作者:wangxinxin
--  發布時間:2010-12-14 14:47:56
--  vb語言實現的卡爾曼濾波源程序
(沒有測試,輸出部分要配合相應的過程)

Dim num1, num2, num3
Dim a(), a1(), e() As Single
Dim alf, xgm(), p(), k1(), k2, k(), p1(), p2(), c(), c1, k3(), k4(), cc() As Single
Private Sub Command1_Click()
CommonDialog1.ShowOpen
Open CommonDialog1.FileName For Input As #1
num3 = Val(Text1.Text) + 2
num2 = Val(Text2.Text)
ReDim a(num3, num2)
For j = 1 To num3
For i = 1 To num2
Input #1, a(j, i)
Next i
Next j
Close #1

End Sub

Private Sub Command2_Click()
CommonDialog1.ShowSave
Open CommonDialog1.FileName For Append As #2

For i = 1 To num2
Print #2, cc(i, 1), cc(i, 2), cc(i, 3)
Next i

Close #2
Shell "c:\\windows\\notepad " & CommonDialog1.FileName
End Sub

Private Sub Command3_Click()

num1 = Val(Text1.Text) + 1
num2 = Val(Text2.Text)
num3 = Val(Text1.Text) + 2
ReDim a1(num2), e(num1, num2) As Single
ReDim xgm(num1), p(num1, num1), k1(num1), k(num1), p1(num1, num1), p2(num1, num1), c(num1), k3(num1), k4(num1), cc(num2, num1) As Single





c(1) = c(2) = c(3) = 0
r = 0.000001
alf = 10
For j = 1 To num1
For i = 1 To num2
e(j, i) = a(j, i)
a1(i) = a(num3, i)
Next i
Next j
For x = 1 To num1
xgm(x) = alf * Sqr(r / e(x, 1))
p(x, x) = xgm(x) ^ 2
Next x

For q = 1 To num2
For i = 1 To num1
k3(i) = k4(i) = 0
Next i
k2 = 0
For i = 1 To num1
For j = 1 To num1
p1(i, j) = p2(i, j) = 0
Next j
Next i
c1 = 0

For m = 1 To num1
k3(m) = 0
For n = 1 To num1
k3(m) = k3(m) + e(n, q) * p(n, m)
Next n
Next m

For m = 1 To num1
k4(m) = 0
k4(m) = k4(m) + 1000 * k3(m) * e(m, q)
Next m

For i = 1 To num1
k2 = k2 + k4(i)
Next i

For v = 1 To num1
k(v) = 0
k(v) = k3(v) * (k2 / 1000 + r) ^ -1
Next v

For i = 1 To num1
For j = 1 To 3
p1(i, j) = k(i) * e(j, q)
Next j
Next i

For l = 1 To num1
For i = 1 To num1
p2(l, i) = 0
For j = 1 To num1

p2(l, i) = 1000 * p1(l, j) * p(j, i) + p2(l, i)
Next j
Next i
Next l

For i = 1 To num1
For j = 1 To num1
p(i, j) = p(i, j) - p2(i, j) / 1000
Next j
Next i

For i = 1 To num1
c1 = c1 + e(i, q) * c(i)
Next i
For i = 1 To num1
c(i) = c(i) + k(i) * (a1(q) - c1)

Next i


For i = 1 To num1

cc(q, i) = c(i)
Next i
Next q

End Sub

Private Sub Command4_Click()
End
End Sub

主站蜘蛛池模板: 亚洲激情五月 | 国产精品视频你懂的网址 | 亚洲三区视频 | 亚洲福利二区 | 亚洲国产欧美日韩精品一区二区三区 | 久久一区二区三区四区 | 久久精品亚洲视频 | 日韩在线观看视频网站 | 丁香婷婷色综合 | 制服丝袜天堂 | 在线亚洲观看 | 九色97| 国产美女福利在线 | 国产精品免费久久 | 亚洲视频免费在线 | 日本一区二区三区免费看 | 久久免费99精品久久久久久 | 亚洲a成人| 麻豆精品久久精品色综合 | 在线播放一区 | 在线视频这里只有精品 | 综合久久99久久99播放 | 男女男精品视频在线播放 | 国产成人综合久久亚洲精品 | 在线播放你懂 | 亚洲毛片免费视频 | gogo大胆高清免费视频 | 精品久久久久久亚洲精品 | 中文字幕精品视频在线观看 | a级毛片在线免费观看 | 一区二区国产精品 | 欧美第一福利 | 在线视频国产网址你懂的在线视频 | 精品视频自拍 | 亚洲国产成人久久99精品 | 国产在线观看一区 | 免费国产最新进精品视频 | 不卡国产视频 | 五月天激情婷婷大综合 | 欧美日在线观看 | 欧美综合图区亚欧综合图区 |