КИХ с произвольной АЧХ

Все что касается фильтрации
dagl
Сообщения: 2
Зарегистрирован: 25 сен 2015, 16:28

КИХ с произвольной АЧХ

Сообщение dagl »

Доброго времени суток, столкнулся с задачей реализации выравнивающего фильтра для ЦАП, с амплитудной характеристикой 1/sinc(x), руководствуясь статьей при вычисление ИХ получаются комплексные значения, предполагаю, что я неверно задал ФЧХ. прикладываю кусок программы. Спасибо за внимание.

F_cap = 720e6; %ширина полосы пропускания ЦАП
f = 0:Fd/360:Fd-Fd/360;%вектор частот для ЦАП
K1=((sinc((f)/F_cap)));%АЧХ ЦАП
N=16
for k = 0:N-1
if (k<(N/2))
FI(k+1) = (-(N-1)*k/N)*pi;
else
FI(k+1) = (-(N-1)*(k-N)/N)*pi;
end
end
H1 = 1./K1([1:16]*22);
komp = H1.*exp(j*FI);
h = ifft(komp)

Аватара пользователя
Бахурин Сергей
Администратор
Сообщения: 1114
Зарегистрирован: 05 окт 2010, 19:55
Контактная информация:

Re: КИХ с произвольной АЧХ

Сообщение Бахурин Сергей »

чему равно Fd в вашем примере?

ВОт пример расчета фильтра

Код: Выделить всё

close all; clear all; clc;

Fd = 720e6; %ширина полосы пропускания ЦАП

N=16;
f = 0:N/2-1;

K=sinc(f*0.1);%АЧХ ЦАП

H = [K, 0, fliplr(K(2:N/2))];


k = 0:N/2-1;
PHI =-(N-1)*pi*k / N;

k = N/2 : N-1;

PHI = [PHI,  -(N-1)*(k-N)*pi/N];


figure; plot(H, 'o')

h = ifft(H.*exp(1i*PHI))

figure; plot(1:N, real(h), 1:N, imag(h))


dagl
Сообщения: 2
Зарегистрирован: 25 сен 2015, 16:28

Re: КИХ с произвольной АЧХ

Сообщение dagl »

720Мгц

Аватара пользователя
Бахурин Сергей
Администратор
Сообщения: 1114
Зарегистрирован: 05 окт 2010, 19:55
Контактная информация:

Re: КИХ с произвольной АЧХ

Сообщение Бахурин Сергей »

в вашем примере не выполняется условие симметричности АЧХ вашего фильтра. Как результат получаете комплексный фильтр.

Ответить