Comparing FIR1 bultiin function for generating coefficients with my Function.

Question

Ahsen Noor Subhopoto 2019-7-11

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/471239-comparing-fir1-bultiin-function-for-generating-coefficients-with-my-function

评论： Ahsen Noor Subhopoto 2019-8-3

Hi, I am trying to create my own FIR1 by using Window functions but somehow I am not getting the exact results. I have shared the codes with corresponding outputs here;

%Window function is Rectangl, W(n) = 1;%
order = 10;
L = order + 1; %length of a Window
Tau = ceil((order-1)/2); % Tau = 5
Fc = 0.99 %for simplicity that is close to 1
%Builtin Function%
Rectangl = fir1(order,Fc,'low',rectwin(L));
%My function%
n = 0:1:order;
Hd_n = (sin((n-Tau)*Fc))./(pi*(n-Tau));
Hd_n(Tau+1) = Fc/pi;
[Hd_n' Rectangl']
%end
Output:
Both outputs are different. Like, Hd_n(1) = -0.0619, whereas Reactangl(1) = 0.0099.

Am I doing something wrong ? Your suggestions and debugging will be helpful for me.

Thank you,

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Kaashyap Pappu 2019-7-26

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/471239-comparing-fir1-bultiin-function-for-generating-coefficients-with-my-function#answer_384929

The difference between the 2 filter coefficients arises for the following reasons:

“fir1” normalizes the coefficients as described here . Adding another argument to function changes this: fir1(order,Fc,'low',rectwin(L),'noscale').
The Hd_n vector has coefficients can be corrected by changing (sin((n-Tau)*Fc)) to (sin(pi*(n-Tau)*Fc)). The lack of pi expands the area under curve.
Hd_n(Tau+1) = Fc corrects the magnitude at center frequency.