How to plot the frequency Response of a FIR filter, whose transfer function depends on sampling frequency?

Question

Alberto Carboni 2015-7-23

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/230956-how-to-plot-the-frequency-response-of-a-fir-filter-whose-transfer-function-depends-on-sampling-freq

回答： Rakshith Sharma Srinivasa 2015-8-4

采纳的回答： Rakshith Sharma Srinivasa

Hello,

I am here because I want to plot the frequency response of a FIR filter, that is

H(e^jw) = 1/3 * (1 + e^-jwk + e^-j2wk)

the problem is that k = M/6, where M = Fs / f0 (being Fs: sampling frequency, f0 = fundamental frequency).

I am using freqz, but I cannot figure out how to include the information about k.

I cannot simply change k, since the inputs I want to pick must be shifted by pi/3 and that depends on Fs.

Thanks,

Alberto

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

请先登录，再进行评论。

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

Answer 1

Rakshith Sharma Srinivasa 2015-8-4

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/230956-how-to-plot-the-frequency-response-of-a-fir-filter-whose-transfer-function-depends-on-sampling-freq#answer_188263

在 MATLAB Online 中打开

Hi Alberto,

You can include information regarding ‘K’ while defining the numerator and denominator coefficients in the ‘freqz’ command. In this case, the numerator coefficients would look like this: B = [1/3...0….0.. 1/3 (at index K+1) ..0…0…1/3(at index 2K+1)]. For example, for K = 3,

B = [1/3 0 0 1/3 0 0 1/3];

The code snippet below delays a sine wave of frequency 100Hz sampled at 6KHz by 10 samples. 10 samples in this case corresponds to a delay of pi/3.

Fs = 6000;
f = 100;
t = 1/Fs:1/Fs:.1;
x = sin(2*pi*f*t);
plot(x);
M = Fs/f;
K = (M/6);
B = zeros(1,2*K);
A = 1;
indices = [1,K+1,2*K+1];
B(indices) = 1/3;
y = filter(B,A,x);
hold on, plot(y)
hold off
figure, freqz(B,A)