over-lapping & adding [FFT] - Error - Index exceeds matrix dimensions

Question

MatlabGirl 2014-12-7

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/165825-over-lapping-adding-fft-error-index-exceeds-matrix-dimensions

评论： MatlabGirl 2014-12-8

Hello,

I am trying to process a signal using the overlap and add method, however I get an error, when running the loop, I am not sure what matrix is exceeding the dimensions or how to solve it, i tried adding counters to trace the issue but I still can't understand where i went wrong, here's the code:

clc;
clear all;
x = [3 9 1 2 3 4 5 6 3 4 5 6 7 8 9 8 7 8];
h = [1 2 1 1]; 
% Code to perform Convolution using Overlap Add Method
n1 = length(x); 
% M-1
n2 = length(h); 
% length of the output
N = n1+n2-1; 
% intializing the output vector
y = zeros(1,N);
% zero padded impulse response
h1 = [h zeros(1,n2-1)]
% block size 
n3 = length(h1);
% new size of output function
% now that we know the block size
y = zeros(1,N+n3-n2);
% FFT of the impulse response
%fft_size = 
H = fft(h1); %, fft_size);
count_1 = 0;
count_2 = 0;
count_3 = 0;
count_4=0;
count_5=0;
for i = 1:n2:n1
if i<=(n1+n2-1)
x1 = [x(i:i+n3-n2) zeros(1,n3-n2)];
count_1 = count_1+1;
else
x1 = [x(i:n1) zeros(1,n3-n2)];
count_2 = count_2 +1;
end
x2 = fft(x1);
x3 = x2.*H;
x4 = round(ifft(x3));
if (i==1)
    y(1:n3) = x4(1:n3);
    count_3 = count_3+1;
else
y(i:i+n3-1) = y(i:i+n3-1)+x4(1:n3);
count_4 = count_4+1;
end
count_5 = count_5+1;
end

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Geoff Hayes 2014-12-7

2
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/165825-over-lapping-adding-fft-error-index-exceeds-matrix-dimensions#answer_161630

在 MATLAB Online 中打开

Maith - if I run your code as a function (named g13.m) or as a script, I observe the following error

 Index exceeds matrix dimensions.
 Error in g13 (line 33)
 x1 = [x(i:i+n3-n2) zeros(1,n3-n2)];

So we know the line, but not which iteration. To make the debugging easier, run the following command before running your above code

dbstop if error

This will cause the debugger to pause at the line which generates the error, and you will be able to look at all variables at that point and so you can get an idea of what is happening. When I do this, I notice that this fails for when i is 17 (as an aside, you may want to use a different index variable name as i and j are also used to represent the imaginary number). The code then fails because of

x(i:i+n3-n2)

as i is 17, n3 is 7, and n2 is 4 with x only being a vector with 18 elements. With that in mind, how would you modify this code to handle the case where i+n3-n2 exceeds 18?

3 个评论
显示 1更早的评论隐藏 1更早的评论

Geoff Hayes 2014-12-8

在 MATLAB Online 中打开

Try something like this for the code that is causing the problem

 if i<=(n1+n2-1)
    % fix the index into x as the minimum of the length of x or
    % the usual equation
    endIdx   = min(length(x),i+n3-n2);
    % determine the number of zeros to pad
    numZeros = n3 - (endIdx-i+1);
    % set x1
    x1 = [x(i:endIdx) zeros(1,numZeros)]
    count_1 = count_1+1;
 else
   % do other stuff
 end

The above code guards against us trying to extract data beyond the length of x by setting the "end index" to the minimum of the length of x or a valid index according to your equation. We then determine the number of zeros to pad given the length of x1 which should be seven.

Note that you may need to try something similar for the else block which seems to create a vector of length 18 if ever we should enter this block.

MatlabGirl 2014-12-8