塑料排水板-植草格-土工膜厂家-泰安东诺工程材料有限公司

Rss & SiteMap

曙海教育集團論壇 http://www.bjzhda.cn

曙海教育集團論壇
共1 條記錄, 每頁顯示 10 條, 頁簽: [1]
[瀏覽完整版]

標題:vb語言實現的卡爾曼濾波源程序

1樓
wangxinxin 發表于:2010-12-14 14:47:56
(沒有測試,輸出部分要配合相應的過程)

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
共1 條記錄, 每頁顯示 10 條, 頁簽: [1]

Copyright © 2000 - 2009 曙海教育集團
Powered By 曙海教育集團 Version 2.2
Processed in .01563 s, 2 queries.
主站蜘蛛池模板: 九九在线精品视频播放 | 色视频在线免费观看 | 日本a级视频在线播放 | 在线观看亚洲专区 | 国产亚洲精品看片在线观看 | 国产日产欧美精品 | 精品无码国产自产在线观看水浒传 | av在线亚洲男人的天堂 | 精品九九九 | 成人免费网站视频ww | 亚洲一区二区三区香蕉 | 99国产精品人妻噜啊噜 | 成人亚洲欧美日韩中文字幕 | 亚洲国产片在线观看 | 亚洲午夜精品a片久久www慈禧 | 国产精品久久久久影视青草 | 免费一级欧美片在线观免看 | 国产成人精品久久一区二区三区 | 日韩精品无码一区二区视频 | 好吊妞人成视频在线观看27du | 亚洲欧美v国产一区二区 | 亚洲一区二区三区自拍公司 | 四川丰满妇女毛片四川话 | 99国产精品免费视频 | 国精品人妻无码一区二区三区喝尿 | 国产在线小视频 | 国产成人精品.一二区 | www.国产.com| a毛片基地免费全部香蕉 | 天天做天天添婷婷我也去 | 国产精品男女 | 不卡视频在线观看 | 亚洲天堂免费在线视频 | 国产偷人妻精品一区二区在线 | 青青草国产精品人人爱99 | 一级毛片牲交大片 | 国产福利视频在线播放 | 国产精品久久久久久久久久一区 | 四虎影视国产884a精品亚洲 | 人妻哺乳奶头奶水 | 九九视频免费精品视频免费 |