午夜伦理在线-午夜伦理在线观看-午夜毛片免费看-午夜美女久久久久爽久久-午夜免费

以文本方式查看主題

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

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

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

主站蜘蛛池模板: 人与动物xxxxx| 91啦视频在线观看 | 高清毛片一区二区三区 | 色综合网亚洲精品久久久 | 啪视频在线 | 国产精品日韩欧美一区二区三区 | 欧美成人久久 | 又黄又爽一线毛片免费观看 | 国产日韩在线播放 | 国产三级毛片视频 | 日本一级毛片2021免费 | 国产资源在线看 | 一级a级国产不卡毛片 | 免费精品国产 | 看黄视频免费 | 欧美二级在线观看免费 | 青青青网 | 夜夜综合| 国产精品嫩草影院视频 | 亚洲在线看片 | 久久亚洲精品无码观看不卡 | 男女很黄很色床视频网站免 | 老司机深夜福利影院 | 国产人体视频 | 欧美aaa性bbb毛片 | 色香婷婷| 久久久久琪琪免费影院 | 国产69精品久久久久9牛牛 | 黄色一级毛片网站 | 亚洲第一页在线 | 亚洲毛片视频 | 亚洲国产欧美日韩精品一区二区三区 | 欧美午夜理伦三级在线观看 | 亚洲精品一区二区不卡 | 国产91久久久久久久免费 | 看毛片网| 真正毛片 | 麻豆精品视频网站在线观看 | 手机在线日韩高清理论片 | 成人在线观看视频免费 | 午夜国产精品无套 |