Question

我一直试图在MATLAB中使用谐波产品频谱找到基本音符。我遇到了一个算法并尝试使用它。我用C大调（钢琴）和音符测试了它 C4 D4 E4 F4 G4 A4 B4 C5 B4 A4 G4 F4 E4 D4 C4

我得到除最后一个C4音符之外的所有音符的正确（几乎接近）频率值，因为此时我得到错误

矩阵尺寸必须达成一致。

这是显示错误的行

f_ym = (1*seg_fft) .* (1.0*seg_fft2) .* (1*seg_fft3)  .* (1*seg_fft4);

我不太确定这里有什么问题.. 我发现了音符开始并在每次发作时执行了FFT，并将其用于谐波产品频谱。这是执行HPS的部分

h = 1;
for i = 2:No_of_peaks

song_seg = song(max_col(i-1):max_col(i)-1);
L = length(song_seg); 
NFFT = 2^nextpow2(L); % Next power of 2 from length of y
seg_fft = fft(song_seg,NFFT);%/L;

%HPS
seg_fft = seg_fft(1 : size(seg_fft,1) / 2);
seg_fft = abs(seg_fft);

%HPS: downsampling
for i = 1:length(seg_fft)
    seg_fft2(i,1) = 1;
    seg_fft3(i,1) = 1;
    seg_fft4(i,1) = 1;
%    f_x5(i,1) = 1;
end

for i = 1:floor((length(seg_fft)-1)/2)
    seg_fft2(i,1) = (seg_fft(2*i,1) + seg_fft((2*i)+1,1))/2;
end

for i = 1:floor((length(seg_fft)-2)/3)
    seg_fft3(i,1) = (seg_fft(3*i,1) + seg_fft((3*i)+1,1) + seg_fft((3*i)+2,1))/3;    
end

for i = 1:floor((length(seg_fft)-3)/4)
    seg_fft4(i,1) = (seg_fft(4*i,1) + seg_fft((4*i)+1,1) + seg_fft((4*i)+2,1) + seg_fft((4*i)+3,1))/4;
end


%HPS, PartII: calculate product
f_ym = (1*seg_fft) .* (1.0*seg_fft2) .* (1*seg_fft3)  .* (1*seg_fft4);


%HPS, PartIII: find max
f_y1 = max(f_ym);

for c = 1 : size(f_ym)
    if(f_ym(c, 1) == f_y1)
        index = c;
    end
end

% Convert that to a frequency
f_y(h) = (index / NFFT) * FS


    h=h+1;
    f_y = abs(f_y)';

 end

嗯，我正试图在存在谐波的情况下找到基频。谐波产品频谱是实现这一目标的一种方式，也就是上述代码中实现的内容。当我进行乘法时，seg_fft的大小似乎是seg_fft2, seg_fft3, seg_fft4的一半。我不知道如何使尺寸达到相同尺寸..这就是我需要帮助的地方.. 非常感谢一些快速的帮助。 Thanx提前:)）

Answer 1

如果没有更多解释，很难准确理解你要做什么。但是在尝试对多个矩阵进行逐元素乘法时会出错。此错误消息表示并非所有矩阵都具有相同的大小，这是元素操作所必需的。要调试这种错误，最简单的方法通常是检查各种矩阵的大小。因此，在产生错误的行之前添加这些行：

disp(size(seg_fft))
disp(size(seg_fft2))
disp(size(1*seg_fft3))
disp(size(seg_fft4))

这可能会告诉您其中一个的尺寸与您的预期不同。在此之后，尝试修复你的for循环，使它们都具有相同的大小。

执行谐波产品频谱时出错

1 个答案: