I want to use Interpole with Zero Padding

Question

Zaref Li 2024-1-5

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2066846-i-want-to-use-interpole-with-zero-padding

编辑： Hassaan 2024-1-5

Hi everyone,

I have a signal [ cos(2*pi*0.2*n) n:0:1:15]. I want Interpolate in time domain by a factor of 5 using zero padding in frequency domain. I write this code but I am not sure that it is correct.

n = 1:15;
x = cos(2*pi*0.2*n);
N = length(x);
X = fft(x, N*5);
x_interpolated = ifft(X);
figure;
subplot(2,1,1);
stem(n, x);
title('Original Signal');
subplot(2,1,2);
stem(1:length(x_interpolated), x_interpolated);
title('Interpolated Signal');

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

请先登录，再进行评论。

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

Answer 1

Image Analyst 2024-1-5

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2066846-i-want-to-use-interpole-with-zero-padding#answer_1384181

编辑：Image Analyst 2024-1-5

在 MATLAB Online 中打开

Why are you doing anything with fft? Just use interp1

xInterp = linspace(min(n), max(n), 5 * length(n));
x_interpolated =interp1(n, x, xInterp);

Or just change your step in your n:

n = 1 : 0.2 : 15;
x = cos(2*pi*0.2*n);

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

Zaref Li 2024-1-5

thank you!

请先登录，再进行评论。

Answer 2

Hassaan 2024-1-5

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2066846-i-want-to-use-interpole-with-zero-padding#answer_1384226

编辑：Hassaan 2024-1-5

在 MATLAB Online 中打开

n = 0:15; % n should start at 0

x = cos(2*pi*0.2*n);

N = length(x);

% Define the half_N variable correctly

half_N = ceil((N+1)/2);

% FFT of the original signal

X = fft(x);

% Zero-padding should be done by adding zeros in the middle of the array, not at the end

% to maintain the Hermitian symmetry of the FFT for a real signal.

X_padded = [X(1:half_N), zeros(1, (N*5)-(N)), X(half_N+1:end)];

% To maintain symmetry, ensure that the number of zeros padded is even

% This can be done by appending an extra zero if the number of points to pad is odd

num_zeros = N*5 - N;

if mod(num_zeros, 2) ~= 0

X_padded = [X(1:half_N), zeros(1, num_zeros+1), X(half_N+1:end)];

end

% IFFT to get the interpolated signal

x_interpolated = ifft(X_padded, 'symmetric'); % Use 'symmetric' to ensure the output is real

% Time vector for the interpolated signal

n_interpolated = linspace(0, max(n), length(x_interpolated));

figure;

subplot(2,1,1);

stem(n, x);

title('Original Signal');

subplot(2,1,2);

stem(n_interpolated, x_interpolated);

title('Interpolated Signal');

------------------------------------------------------------------------------------------------------------------------------------------------

If you find the solution helpful and it resolves your issue, it would be greatly appreciated if you could accept the answer. Also, leaving an upvote and a comment are also wonderful ways to provide feedback.

Professional Interests