バタワースフィルター (LPF) の周波数スペクトルの書き方

Question

N/A 2019-10-30

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/488249-lpf

评论： N/A 2019-10-31

采纳的回答： Yoshio

Trace_1.mat

以下は、バタワースフィルター(LPF)のプログラムです。Trace_1は10MHzのsin波のデータ(横軸：時間、縦軸：電圧)です。

x1=Trace_1;

fc=1*10^6;%フィルターのカットオフ周波数

fs=100*10^6;%サンプリング周波数

[b,a]=butter(50,fc/(fs/2));

%freqz(b,a)

x2=filter(b,a,x1);

plot(x2);

バタワースフィルターとx2 (バタワースフィルターにx1を通した後の信号) の周波数スペクトル ( 横軸：周波数(Hz)、縦軸：利得(dB) )の出力方法を教えていただけませんか？ご回答お待ちしております。

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

Yoshio 2019-10-31

t = Trace_1(:,1);
x = Trace_1(:,2);
plot(t,x)
fs = 1/(t(2)-t(1))

fs =

2.0000e+09

となっていますので、サンプリング周波数の設定は2GHzだと思います。

またfilter(b,a,x1)としていまうと、時系列ベクトルtまでフィルタすることになります。

请先登录，再进行评论。

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

Answer 1

Yoshio 2019-10-31

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/488249-lpf#answer_399148

先のコメントにも書きましたが、まずご自身のデータTrace_1について理解してください。その上で以下のプログラムが参考になればと思います。

t = Trace_1(:,1);
x = Trace_1(:,2);
plot(t,x)
fs = 1/(t(2)-t(1))
fc = 10^6;
[b,a]=butter(3,fc/(fs/2));
freqz(b,a)
y=filter(b,a,x);
figure
plot(t,[x y]);
grid on;
legend('x','y')
figure
[p,f] = pspectrum([x y],fs);
plot(f,pow2db(p))
grid on
xlabel('Frequency (Hz)')
ylabel('Power Spectrum (dB)')