Перенос спектра в 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 и не совпадают по фазе с исходным синусом.

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

[Самурай]

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

Угу, из-за фильтрации. Тот фильтр, что у Вас, просто давит зеркальную часть смещенного в ноль синуса. Что там происходит при обратном переносе мне не ведомо (лень формулы писать), но явно не то, что хотелось бы. Полосу пропускание фильтра надо бы увеличить, или делать комплексный фильтр.

[ataradov]

Полосу пропускание фильтра надо бы увеличить, или делать комплексный фильтр.

Да, действительно, это из-за того, что фильр с действительными коэффициентами.

В общем идею это иллюстрирует, только нужно правильно фильтр посчитать.

Но вообще так избаляться от помех - это изврат.