各位大大
小弟目前遇上一个问题就是用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