for 循环关联和编辑向量

Question

我需要在 2 个变量中编辑向量，使它们相互对应。

For example, 
y = 1 2 3 4 5 6
q = 8 26 1 5 1 6

Let's say I fix q and shift y by 1
y =    1 2 3 4 5 6
q = 8 26 1 5 1 6

I need vector y1 = 1 2 3 4 5 and q1 = 26 1 5 1 6
This means I need to kick out 6 in y and 8 in q respectively.

Let's say I fix q and shift y by 2
y =      1 2 3 4 5 6
q = 8 26 1 5 1 6
As before, now I need to remove 5,6 in y and 8,26 in q respectively.

我想在 for 循环中执行此操作，因为我的向量很长。现在，我正在努力为 q （这是我的 soundtwo）获取正确的向量，如下所示。有什么建议吗？

% Creating time vector, "t"
t = linspace(0,16*pi,1000);

sound1 = 5*(cos(t) + 1*(rand(size(t))-0.5));
sound2 = 8*(cos(t) + 1.5*(rand(size(t))-0.5));

% Setting the time shift "dt"
dt = 1000;

% Creating a matrix to store product later on
list = zeros(dt,1);

% For loop for different shifts
for i=1:dt

      % Now edit sound1 such that sound1 shifts while sound2 remains unchanged
      %different time shift

      sound1 = 5*(cos(t+ i ) + 1*(rand(size(t))-0.5));
      sound2 = 8*(cos(t) + 1.5*(rand(size(t))-0.5));

      % Shifting sound1
      soundone = sound1(i:numel(sound1))

      % Sound 2 unchanged, but have to assign respective vector to sound1
      soundtwo = sound2()

      multipliedsound = (soundone) .* (soundtwo);

      add = sum(multipliedsound)

      product = add  / numel(t);

      % Append product to list vector
      list(i,1) = product;

end

Answer 1

circshift 函数会旋转 y 的值，然后您可以设置旋转后的 NaN 值。将 y 的旋转值设置为 NaN 就无需更改或删除 q 中的值：具有 NaN 的乘积值也将是 NaN，然后​​被 NaN sum 忽略。下面的示例将值移位 2，但您可以轻松地替换为其他移位值。

y = [1 2 3 4 5 6];
q = [8 26 1 5 1 6];

shift_value = 2;

y_shifted = circshift(y, shift_value);
y_shifted(1:shift_value) = NaN;

product_value = y_shifted .* q;
sum_value = nansum(product_value);