Перенос спектра в Matlab Опции

[deneb13]

Здравствуйте! Возникла проблема. Реализую перенос спектра sin на 0-ю частоту. Для этого умножаю на экспоненту. Вопрос, почему из за переноса изменяются амплитудные соотношения?Совсем немного,но меняются!ведь по идее должны изменяться только частоты гармоник и все. Код привожу ниже
Fs=10*10^3;
f1=2*10^3;

t=0:1/Fs:1;
signal=1*sin(2*pi*f1*t);
y1=fft(signal);

figure(1)
plot(abs(y1)),grid


s1=signal.*exp(j*2*pi*f1*t);
figure(2)
y2=fft(s1);
plot(abs(y2)),grid
Спасибо!

[deneb13]

Например,в Тузове "Стат.теория приема сложных сигналов",1977 г.,стр.128.Встречала еще в некоторых книгах,к сожалению, названия не вспомню,за последнее время много информации перечитала. Для шпс сигнала в принципе метод не планировала менять. Т.е.предполагала просто к сигналу прибавить sin,а дальше действовать по тому же алгоритму. Почему fft нежелательно использовать, не могли бы объяснить?


Хочу еще задать вопрос:почему в аплитудах sin я указываю 1 и 2, а на графиках они отображаются с амплитудами 5000 и 10000.из за чего это?

[ataradov]

Например,в Тузове "Стат.теория приема сложных сигналов",1977 г.,стр.128.

Хм, ну книжки 77 года я сейчас икать не буду, конечно

Почему fft нежелательно использовать, не могли бы объяснить?

В первую очередь потому, что FFT вещь блочная, причем FFT не очень хорошо сохраняет временное представление сигнала, после его обраюотки в чатотной области. В частности - Ваш пример, где фильтрация достигается путем зануления некоторых компонент в частотной области и обратным преобразованием, так не делают в жизни. Это работает для суммы двух синусов, но не более.

Я бы рекомендовал экспериментировать сразу на более или менее реальном сигнале, так как синусы (да еще с очень красиво подогнанными частотами под окно FFT) не репрезентативны. Ваш метод в реальности исказит сигнал так, что демодулятору скорее самому проще с помехой будет справиться, чем после такой фильтрации.

Хочу еще задать вопрос:почему в аплитудах sin я указываю 1 и 2, а на графиках они отображаются с амплитудами 5000 и 10000.из за чего это?

FFT в матлабе не нормализованное, результат нужно делить на длинну выборки / 2.


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

В Вашем случае получается, что берется FFT сразу от всего сигнала, так не получится сделать на практике, размер блока скорее всего не будет превышать 16384 или около того. Вот это сразу заложите в модель. И даже на 2-х синусах с текущим подходом возникнут проблемы.