Obtain filter x,y points (Designed with notchfilt)

Question

Jesus Sanchez 2021-9-21

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1457969-obtain-filter-x-y-points-designed-with-notchfilt

回答： Pratyush Roy 2021-9-28

Hello everyone,

I am trying to do some Mean Square Error between a signal that I measured and the "ideal" shape of it, which corresponds to a stopband filter. I am able to design the filter by using notchspec, design and notchfilt. However, I would like to get the data points plotted on the image by fvtool, to perform MSE. I have not been able to find how to do it at all.

I would really appreciate some advice!

The code that I use to generate the filter:

F0 = 105e9;   % Notch
Fs = 250e9; % Sampling freq
Q = 10;
BW = 10e9;
orderFilter = 2;
notchspec = fdesign.notch('N,F0,Q',orderFilter,F0,Q,Fs);
notchfilt = design(notchspec,'IIR','SystemObject',true);
fvtool(notchfilt,'Color','white');
%ylim([-30 0])
%xlim([75 110])

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

请先登录，再进行评论。

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

Answer 1

Pratyush Roy 2021-9-28

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1457969-obtain-filter-x-y-points-designed-with-notchfilt#answer_796894

在 MATLAB Online 中打开

Hi Jesus,

As per my understanding, you want to get the data ponts of the ideal response plotted on the image generated by fvtool.

This can be done once we get the fvtool handle. You can refer to the code snippet given below for details:

F0 = 105e9;   % Notch
Fs = 250e9; % Sampling freq
Q = 10;
BW = 10e9;
orderFilter = 2;
notchspec = fdesign.notch('N,F0,Q',orderFilter,F0,Q,Fs);
notchfilt = design(notchspec,'IIR','SystemObject',true);
h = fvtool(notchfilt,'Color','white');
%% h is the fvtool handle. We will get the axis objects from the fvtool and 
% use that to superimpose the ideal response
s = get(h);
hchildren = s.Children;
haxes = hchildren(strcmpi(get(hchildren,'type'),'axes'));
hline = get(haxes,'children');
x = get(hline,'XData');
y = get(hline,'YData');
y1 = randn(1,8192); % Using random noise for example. You can replace that with your own ideal response
                    % data points here
line(haxes,x,y1)