向量化求和的 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 矩阵 的行中定义内部频率。然后您希望将所有这些结果添加到单个组合信号中。您可以在不使用 for 循环的情况下执行此操作，如下所示：

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)));