Почему, когда я добавляю сигнал нулями, fft исходного сигнала и сигнала, дополненного нулями, больше не объединяются вокруг одной и той же точки?

Я думаю, что неправильно реализую код. Я думаю, что исходный и дополненный нулями сигнал будут объединены вокруг одной и той же точки с одинаковыми пиками. Я неправильно понимаю это или проблема в моем коде?

clc;clear;
N=257; %number of points in the signal
f=330.5; %frequency of signal
fs=1024; %sampling frequency 
Ts=1/fs; %sampling period
ts=0:Ts:(N-1)/fs; %duration of signal
x=sin(f*ts);%generation of sampled signal
X=fftshift(fft(x)); %shifted FFT of signal
figure(5)
stem(abs(X))
M=2048; %number of points desired in the new signal that will be zero padded
zerovec=zeros(1,(M-N)); %creating enough 0's to add to the end of the original signal to achieve the desired length
x1=[x zerovec]; %concatenating original signal and 0's to get zero padded signal 
X1=fftshift(fft(x1)); %fft of zero padded signal

figure()
stem(abs(X)) %discrete plot of original signal
hold on 
stem(abs(X1)) %discrete plot of zero padded signal

matlab signal-processing fft dft

10.04.2019 21:25

Стоит ли изучать PHP в 2023-2024 годах?

Привет всем, сегодня я хочу высказать свои соображения по поводу вопроса, который я уже много раз получал в своем сообществе: "Стоит ли изучать PHP в...

Поведение ключевого слова "this" в стрелочной функции в сравнении с нормальной функцией

В JavaScript одним из самых запутанных понятий является поведение ключевого слова "this" в стрелочной и обычной функциях.

Приемы CSS-макетирования - floats и Flexbox

Здравствуйте, друзья-студенты! Готовы совершенствовать свои навыки веб-дизайна? Сегодня в нашем путешествии мы рассмотрим приемы CSS-верстки - в...

Тестирование функциональных ngrx-эффектов в Angular 16 с помощью Jest

В системе управления состояниями ngrx, совместимой с Angular 16, появились функциональные эффекты. Это здорово и делает код определенно легче для...

Концепция локализации и ее применение в приложениях React ⚡️

Локализация - это процесс адаптации приложения к различным языкам и культурным требованиям. Это позволяет пользователям получить опыт, соответствующий...

Пользовательский скаляр GraphQL

Листовые узлы системы типов GraphQL называются скалярами. Достигнув скалярного типа, невозможно спуститься дальше по иерархии типов. Скалярный тип...

242

Перейти к ответу Данный вопрос помечен как решенный

Ответы 1

Ответ принят как подходящий

При заполнении сигнала нулями его частотный спектр становится более плотным. В некотором смысле вы интерполируете в частотной области при заполнении нулями пространственной области.

Если вы начертите два частотных спектра с правильными частотами вдоль оси x, вы увидите, что они перекрываются:

N=257;
f=330.5;
fs=1024;
Ts=1/fs;
ts=0:Ts:(N-1)/fs;
x=sin(f*ts);

X=fftshift(fft(x));
F=0:fs/N:fs-fs/N;     % <<< NEW!

M=2048;
zerovec=zeros(1,(M-N));
x1=[x zerovec];
X1=fftshift(fft(x1));
F1=0:fs/M:fs-fs/M;    % <<< NEW!

figure()
stem(F,abs(X))        % <<< NEW! using F
hold on 
stem(F1,abs(X1))      % <<< NEW! using F1

plot