用Matlab計算jacobian矩陣解析解

2022-08-22 01:45:09 字數 439 閱讀 7169

做擴充套件卡爾曼濾波(ekf)的時候需要用到jacobian矩陣。有時手工求解難度較大這時可以用matlab自動求出jacobian矩陣的解析解。

以雷達觀測矩陣為例為例

syms x y vx vy %定義符號變數

jacobian([sqrt(x2+y2), atan(y/x),(xvx + yvy)/(sqrt(x2+y2))], [x; y; vx;vy])

輸出的結果比較複雜, 進行多項式簡化

ss = vx/(x^2 + y2)(1/2) - (x(vxx + vyy))/(x^2 + y2)(3/2)

simplify(ss)

輸出結果為:

(vxy^2 - vyxy)/(x^2 + y2)(3/2)
對比文獻上的公式是一樣

用MATLAB計算函式的積分

很基礎但是很實用,自己寫的時候知道一些指令,但是搭配在一起就用的不太好 a x 36 x 11 n quadrat x c.2 omega0.4 x.4.1 2.d.omega0 x 2 1 1.omega0 x 2 2 2.d.omega0 x 2 n quadrat 是x的函式 這種都是常見的語...

用matlab程式設計實現計算器

這學期開了軟體工程工具及其應用這麼課,任老師給我們教了matlab程式語言,做了課程實驗,自己用matlab實現了計算器,發現matlab確實是一門很有用得計算機程式語言,它把數值計算和視覺化環境整合到一起,非常直觀,而且提供了大量的函式,功能強大,工應用範圍廣泛。calculate.m 函式檔案,...

matlab 時間計算

matlab 時間函式 2008 11 13 22 51 time and dates.時間個日期 current date and time.當前時間和日期 now 當前日期時間數字 current date and time as date number.date 當前日期字串 current ...