Как создать фильтр в матлабе? Опции

[Stefan1]

Доброго времени суток!
Подскажите код в матлабе для полосового фильтра с полосой пропускания 0.99 - 1.01 ГГц. Сигнал синусоидальный амплитудно модулированный несколькими частотами (низкие частоты+высшие гармоники). Несущая частота сигнала - 1 ГГц.

Сообщение отредактировал Stefan1 - Nov 28 2012, 06:41

[V_G]

Помещаете рассчитанные коэффициенты КИХ-фильтра в массив Num, после чего вводите:
y=filter(Num,1,x), где y - выходной сигнал, x - входной сигнал. Если фильтр - БИХ, вместо единицы пишете коэффициенты а

ЗЫ. А аппаратно-то чем 1 ГГц фильтровать собрались?

[Stefan1]

Помещаете рассчитанные коэффициенты КИХ-фильтра в массив Num, после чего вводите:
y=filter(Num,1,x), где y - выходной сигнал, x - входной сигнал. Если фильтр - БИХ, вместо единицы пишете коэффициенты а

Что-то матлаб ругается, пишет:
"Dimension argument must be a positive integer scalar in the range 1 to 2^31" в этой записи: y=filter(Hd,1,x)

Коэффициенты Num таким образом считаются?:

[N,Wn,BETA,TYPE] = kaiserord([Fstop1 Fpass1 Fpass2 Fstop2]/(Fs/2), [0 1 0], [Dstop1 Dpass Dstop2]);
b = fir1(N, Wn, TYPE, kaiser(N+1, BETA), flag);
Hd = dfilt.dffir( В );

y=filter(Hd,1,x)
plot(t2,yidi,'r')

ЗЫ. А аппаратно-то чем 1 ГГц фильтровать собрались?

Мне это вообще нужно для выделения первой гармоники из спектра сигнала, т.к. простыми методами - не получается из-за сложности сигнала. Хочу посчитать мощность усилителя по первой гармонике.

Сообщение отредактировал Stefan1 - Nov 28 2012, 11:49

[V_G]

y=filter(Hd,1,x)

y=filter(b,1,x) , если массив коэффициентов - b

[Stefan1]

y=filter(b,1,x) , если массив коэффициентов - b

Благодарю за помощь!