[問題]是否有更快的做法?

作者: k0939620064 (skan)   2014-04-30 17:59:24
各位大大
小弟目前遇上一個問題就是用for迴圈做訊號處理時,往往都要跑個1.20分鐘以上
想請問各位大大是否有更快的方法?
以下是小弟的程式碼
for m=1:4
for n=1:500000
signal_square(n,m)=sum(sig((n-1)+1:n,m).^2);
% 將收到的訊號平方
phase_noise(n,m) = angle(signal_square(n,m))./2;
% 將平方的訊號取角度後除2 = phase noise
phase_sig((n-1)+1:n,m) = (((angle(sig((n-1)+1:n,m)))- phase_noise(n,m)));
% 原始訊號= arg(收下來的訊號)-phase noise
phase_comp((n-1)+1:n,m) = exp(j*phase_sig((n-1)+1:n,m));
% 原始訊號
end
end
作者: celestialgod (天)   2014-04-30 18:58:00
寫MEX or 用GPU
作者: sunev (Veritas)   2014-04-30 19:29:00
(n-1)+1:n ?

Links booklink

Contact Us: admin [ a t ] ucptt.com