[请益] 悬赏5000P 请教我把MATLAB改成C语言= =

楼主: Ecampus (7.7)   2019-06-07 14:29:47
我知道有matlab coder可以写成C
但是第一次碰 很不熟悉= = 苦手ING
以下matlab程式码:
% 17M for the example of free body falling
% Using previous data to modified
clear all
clc
I=sparse(eye(6,6));
N=10;
O=[ 这里10个数];
y=O';
gg=[1 2 3]; hh=eye(6,6);
aa=sparse(conv2(hh,gg,'same'));
aa=[aa];
M = sparse(eye(6));
H = sparse(zeros(1,6));
H(1,3) = 1;
Q = sparse(160*eye(6));
R = 1;
x0 = [这里6个数];
%初始值设定成x0
X_befored = zeros(6,N);
X_saved = zeros(6,N);
Pk = sparse(eye(6)); %初始矩阵
xk = x0; %卡尔曼滤波初始值
for k=2:N
x_kp = A*xk; %公式 1
%将公式1 计算得到的在k时间卡尔曼滤波估计值 xk ,存在X_saved的第 k-1 行
X_befored(:, k-1) = x_kp;
P_kp = A*Pk*A'+ Q; %公式 2
K = (P_kp*H') / (H*P_kp*H'+ R); %公式 3
xk = x_kp + K*(y(k-1)-H*x_kp); %公式 4
%将公式4 计算得到的在k时间卡尔曼滤波估计值 xk ,存在X_saved的第 k-1 行
X_saved(:, k-1) = xk;
Pk = (I - K*H)*P_kp; %公式 5
end
%提取卡尔曼估计值中的,第3列
xhat_3 = X_saved(3,:);
....请问上面这组程式码 怎么用matlab coder转成C
楼主: Ecampus (7.7)   2019-06-07 15:07:00
我还在想说 是不是每个参数都要用CODER定义成函数= =?
作者: sunev (Veritas)   2019-06-07 18:07:00
vector用sparse没什么意思,公式3不是直接取值吗?
楼主: Ecampus (7.7)   2019-06-07 23:17:00
vector因为里面很多零...所以我设成sparse公式3是直接取值 式的
作者: sunev (Veritas)   2019-06-08 11:44:00
只针对matlab的部份,有不少效能上可能可以改进的地方
楼主: Ecampus (7.7)   2019-06-08 12:40:00
!?是吗 可以请问一下大概的方向吗? 我已经改了几次@@
作者: sunev (Veritas)   2019-06-08 13:00:00
就如我上面所说的,取值直接取,向量不要用sparse
楼主: Ecampus (7.7)   2019-06-08 23:23:00
啊 我这个向量其实数量很多 而且超过99.9%都是零所以用sparse比较快经过尝试了XD取值直接取会比较快吗 我公式1~5都是要迭代的 迭代十万次有时间前后关系的
作者: sunev (Veritas)   2019-06-08 23:43:00
我觉得直接取值会比较快,所有索涉到H的公式都可以换掉

Links booklink

Contact Us: admin [ a t ] ucptt.com