作者:
Ecampus (7.7)
2019-06-01 23:45:03懸賞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形式 才是正確