懸賞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形式 才是正確
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 36.233.109.146
※ 文章網址: https://www.ptt.cc/bbs/MATLAB/M.1559403905.A.352.html
※ 編輯: Ecampus (118.170.50.22), 06/02/2019 01:02:59
→
06/02 09:32,
5年前
, 1F
06/02 09:32, 1F
→
06/02 18:53,
5年前
, 2F
06/02 18:53, 2F
推
06/03 08:12,
5年前
, 3F
06/03 08:12, 3F
→
06/03 13:03,
5年前
, 4F
06/03 13:03, 4F