[an error occurred while processing this directive]
|
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) потом сделать