悬赏2000P 感恩
第一个短CODE:
%% ========
clear;clc;
%% ========
% 假设系数(done)
% 假设初始条件
C=[1
2
3
4
5];
% 其他假设
t=10;% 计算次数
%% ========
% 建立F矩阵
gg=[4 5 6]; hh=eye(9);
aa=conv2(hh,gg,'same');
F=sparse([aa]);
%% ========
% 开始计算
for k=1:t
D=F*C% 计算D矩阵
C=D;% 用D矩阵取代C矩阵
end
................................
第二个短CODE:
I=sparse(eye(10));
N=5;
y=[1
2
3
4
5];
A=[dd];
M = sparse(eye(9));
H = sparse(zeros(1,9));
Q = sparse(eye(9));
R = 1;
x0 = [5
4
3
2
1];
Pk = sparse(eye(9));
xk = x0;
for k=1:N
x_kp = A*xk; %公式 1
P_kp = A*Pk*A'+ Q; %公式 2
K = (P_kp*H') / (H*P_kp*H'+ R); %公式 3
xk = x_kp + K*(y-H*x_kp); %公式 4
Pk = (I - K*H)*P_kp; %公式 5
end
=====================================
抱歉 刚接触GPU加速运算
我自己有尝试改成GPU的形式
但是我笔电没有显卡可以测试 所以还不知道自己写的GPU是不是对
(所以明天要到网咖去测试)
为了怕自己写错= = 想先请教各位高手 这两个短程式码...要怎么改成GPU形式 才是正确