1. lowpass 函数
注意,只有2018年之后的matlab才有lowpass, bandpass 函数。
lowpass(x,fpass,fs) % x 一维信号,fpass:截止频率,fs:采样频率。
例子:
两个频率的信号,一个是50HZ, 振幅为1, 一个是250HZ,振幅为2
fs = 1e3; %sampling rate per second
t = 0:1/fs:1; % time series
x = [1 2]*sin(2*pi*[50 250]'.*t) + randn(size(t))/10; %两个频率的信号,一个是50HZ, 一个是250HZ.
做lowpass, 以150为截止频率
lowpass(x,150,fs) %以150为频率上限,做lowpass
运行结果:
如果以60HZ 为截止频率,
若把Steepness,设置为0.999