矢量化for循環，做一筆總金額

Question

問候全部

我該如何矢量化這個for循環？

t_rebuilt=linspace(0,1,length(inner_freq)); 

for ii=1:1:length(inner_freq)-1;ii=ii+1; 

    aa_sig_rebuilt=inner_freq(ii,2)*cos(2*pi*t_rebuilt*inner_freq(ii,1)+inner_freq(ii,3)); 
    aa_sig_combined=aa_sig_combined+aa_sig_rebuilt; 

end; 

我試圖與線替換它下面，看看它的工作，但我只是得到一個直線

aa_sig_rebuilt=inner_freq(ii,2).*cos(2*pi*t_rebuilt*inner_freq(ii,1)+inner_freq(ii,3));

Answer 1

看來，t_rebuilt是一個1乘N的時間點向量，您要在其中計算一組三角函數，並使用在M乘3矩陣的行中定義的參數inner_freq。然後，您想要將所有這些結果添加到單個組合信號中。

temp = 2*pi*inner_freq(:,1)*t_rebuilt; %# This matrix multiplication will 
             %# result in an M-by-N matrix 
temp = temp+repmat(inner_freq(:,3),1,numel(t_rebuilt)); %# Replicate and add 
                 %# column 3 
aa_sig_combined = inner_freq(:,2).'*cos(temp); %'# Gives you your 1-by-N 
               %# combined signal 

然後可以繪製結果如下：：可以爲循環如下做到這一點沒有一個

plot(t_rebuilt,aa_sig_combined); 

Answer 2

試試這個

sum(aa_sig_rebuilt=inner_freq(:,2).*cos(2*pi*t_rebuilt*inner_freq(:,1)+inner_freq(:,3)));