[问题]是否有更快的做法?

楼主: 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