Перенос спектра в 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]

Помеха в дальнейшем так и будет sin-узкополосная.Сигнал будет широкополосный. Извините,запутала вас. До этого я коэффициенты фильтра рассчитывала с помощью fdatool, а сейчас сделала с помощью firls() . Фильтрация осуществлена с помощью filter(). Привожу код
Fs=10*10^3;
f1=2*10^3;
f2=4*10^3;
t=0:1/Fs:1-1/Fs;

signal=1*sin(2*pi*f1*t)+2*sin(2*pi*f2*t);
y1=fft(signal);

figure(1)
plot(abs(y1)),grid
title('Спектр исходного сигнала');


s1=signal.*exp(-j*2*pi*f1*t);
figure(2)
y2=fft(s1);
plot(abs(y2)),grid
title('Спектр перенесенного сигнала');


f=[0 0.1 0.15 1];
a=[1 1 0 0];
b=firls(80,f,a);


B=b';
z=filter(B,1,s1);


y2(10:(length(y2)-10))=0;
figure(3)
plot(abs(y2)),grid
s_f = ifft(y2);
title('Спектр фильтрованного сигнала ОБНУЛЕНИЕМ');

s_obr=z.*exp(j*2*pi*f1*t);
spektr=fft(s_obr);
f_spektr=angle(s_obr);
figure(4)
plot(abs(spektr)),grid
title('АХ обратного переноса сигнала');


rez=signal-s_obr;
spektr_rez=fft(rez);


z0=fft(z);
figure(5)
plot(abs(z0)),grid
title('Спектр фильтрованного сигнала');

figure(6)
plot(abs(spektr_rez)),grid
title('КОМПЕНСАЦИЯ');