Перенос спектра в 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
Спасибо!

[ataradov]

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

Код:

CODE

clear all;

Fs = 10e3;
f1 = 1e3;

t = 0:1/Fs:100000/Fs;

s = sin(2*pi*f1*t);

s0 = s .* exp(-j*2*pi*f1*t);

h = firls(150, [0 150 1000 Fs] / Fs, [1 1 0 0]);

sf = filtfilt(h, 1, s0);

sh = sf .* exp(j*2*pi*f1*t);

figure; hold on; plot(s); plot(real(sh), 'r'); plot(imag(sh), 'g');

Тут происходит перенос помехи в 0, фильтрация ФНЧ и перенос результата назад.

Проблема - если частота помехи низкая (f1 = 20 Гц), то imag(sh) около 0, а real(sh) совпадает по фазе и амплитуде с исходным синусом. Но если частоту увеличивать, то постепенно действительная и мнимая части становятся амплитудой 0.5 и не совпадают по фазе с исходным синусом.

Похоже это из-за фильтрации, но я не могу с ходу понять почему именно.