バターワースのバンドバスフィルタ　

2018 10 25
1 个回答
回答已采纳
更新时间：2018 10 28
6 次查看（30 天）

0 个投票

Fc1=0.5;
Fc2=50;
N = 30;
バターワースでローカット周波数0.5Hz ハイカット周波数60Hzの バンドパスフィルタですが、フィルタの遅延を補正すう方法を教えていただきたいです。
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Hd = butter(Fc1,Fc2,N);
freqz(Hd);
filt_data = filter(Hd, buttLoop); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
filt_data = filter(Hd, buttLoop);
filterの関数を
filt_data = filtfilt(Hd, buttLoop);
に変更すればよいのでしょうか。
ノッチフィルタではfiltfiltで上手く補正できたのですが。
このようなエラーが出ます。
エラー: filtfilt (line 78)
入力引数が不足しています。
エラー: tougou3 (line 84)
filt_data = filtfilt(Hd, buttLoop);
もし、ご存知の方がいたら教えてください。

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

 采纳的回答

Hirokazu Tanaka
Hirokazu Tanaka 2018-10-25

0 个投票

Fc1=0.5;
Fc2=50;
N = 30;
Hd = butter(Fc1,Fc2,N);
は、実行すると
エラー: butter (line 62)
カットオフ周波数は (0,1) の区間内になければなりません。
とのエラーになりますが、
Hd = butter(Fc1,Fc2,N);
は自作の関数でしょうか？それとも https://jp.mathworks.com/help/signal/ref/butter.html ?

3 个评论
显示 1更早的评论 隐藏 1更早的评论

neshin
neshin 2018-10-26
ありがとうございます。
回答になっているかわかりませんが。
Hd = butter(Fc1,Fc2,N); は、tdatoolで作った関数です。 function Hd = butte と作成された関数に引数(Fc1,Fc2,N) を与えて自作関数にしています。
以下のコードがfdatoolで作成された関数で、 引数を(Fc1,Fc2,N)を入れてカットオフ周波数を呼び出しもとで、 変更できるようにしています。
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function Hd = butter(Fc1,Fc2,N)
%BUTTER 離散時間フィルター オブジェクトを返します。
% MATLAB Code
% Generated by MATLAB® 9.4 and Signal Processing Toolbox 8.0.
% Generated on: 23-Oct-2018 14:50:10
% Butterworth Bandpass filter designed using FDESIGN.BANDPASS.
% All frequency values are in Hz.
Fs = 200; % Sampling Frequency
% N = 10; % Order
% Fc1 = 0.5; % First Cutoff Frequency
%Fc2 = 60; % Second Cutoff Frequency
% Construct an FDESIGN object and call its BUTTER method. h = fdesign.bandpass('N,F3dB1,F3dB2', N, Fc1, Fc2, Fs); Hd = design(h, 'butter');
% [EOF]
Hirokazu Tanaka
Hirokazu Tanaka 2018-10-28
filtfilt 関数のページ https://jp.mathworks.com/help/signal/ref/filtfilt.html を見ると、いくつか呼び出し方がありますが、
y = filtfilt(d,x)
を使う場合は d は digitalFilter オブジェクト である必要があるみたいです。エラーの原因はここかと思います。
見せて頂いた butter 関数を使うと
>> whos Hd
Name Size Bytes Class Attributes
Hd 1x1 dfilt.df2sos
と、digitalFilterオブジェクトではありません。
ちょうど filtfilt関数ページにあるサンプルコードでもバターワースフィルターの例がありますし、
https://jp.mathworks.com/help/signal/ref/designfilt.html
も参考になると思います。
neshin
neshin 2018-10-28
ありがとうございます。 参考にさせて頂きます。

请先登录，再进行评论。

更多回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 フィルターの設計 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!