matlab - matlabバターワースフィルター

Question

バターワースフィルターを適用する必要があります。カットオフ周波数を証明するには、写真の2番目の正方形のように、「残差とカットオフ周波数」のプロットを作成する必要があります。

残差は

だから私は適用したい

-バターワースフィルター

-ローパス

• 順序n=2

• カットオフ周波数wn=[0,20] Hz

だから私の疑問は、matlabで残差を見つける方法があるかどうかです。そして、それらをカットオフ周波数に対してプロットします。

どんな助けでも大歓迎です。

Answer 1

入力信号を受け取り（それを呼び出しますx）、次のようなことを行うことができます

[b,a]=butter(2, [0 0.3]) % assuming the cutoff is 0.3 times the Nyquist frequency
y = filter(b,a,x)
n = length(x)
z = x - y
z = z.^2
z = sum(z) / n
R = sqrt(z)