Отправлено homekvn 19 июля 2006 г. 17:44 В ответ на: Не-е, фазе пофигу, рвется она или нет, поскольку экспонента от нее, умноженную на j все равно... Блин! Вот на j я и забыл домножить! Балда! Все получилось! отправлено <font color=gray>homekvn</font> 19 июля 2006 г. 17:25

clear;

f_src = fir1(128, 0.3, 'high');

F_SRC = fft(f_src);

A_SRC = abs(F_SRC);

PH_SRC_new = imag(hilbert( log(A_SRC(1: ceil (length(A_SRC)) ) ) ));

F_SRC_new = A_SRC.*exp(j*PH_SRC_new);

f_src_new = real(ifft(F_SRC_new));

f_src_new(2:end) = f_src_new(end:-1:2);

figure(10);
plot(f_src_new);
grid on;
hold on;
plot(f_src, 'g');
zoom on;
hold off;

figure(12);
plot(angle(fft(f_src)), 'g');
hold on;
grid on;
zoom on;
plot(angle(fft(f_src_new)), 'b');
hold off;

% Можно здесь еще и freqz(f_src_new) потом сделать