UC知道求助关于matlab滤波功能
来源：百度知道 编辑：UC知道 时间：2024/05/25 09:48:04
现有一组拉索振动的实验数据，每隔0.002s在拉索中间点出测量位移值，得到txt文件，其中一部分数据如下：
0.0000000（时间） 6.0159998（位移值）
0.0020000 5.8280001
0.0040000 5.6399999
0.0060000 5.3920002
0.0080000 5.0939999
0.0100000 4.8660002
0.0120000 4.5679998
0.0140000 4.2519999
0.0160000 3.9200001
0.0180000 3.5139999
0.0200000 3.1500001
0.0220000 2.7460001
0.0240000 2.3239999
0.0260000 1.9040000
0.0280000 1.3900000
0.0300000 0.9140000
0.0320000 0.4040000
0.0340000 0.0000000
0.0360000 -0.4720000
0.0380000 -0.9740000
0.0400000 -1.3099999
为了消除其他频率波或者噪音的影响，现在希望通过滤波将杂质除去，小弟以前没有用过matlab滤波功能，请问高手具体如何操作啊？感激不尽
请问heaventian，能稍微详细一点吗
0.0000000（时间） 6.0159998（位移值）
0.0020000 5.8280001
0.0040000 5.6399999
0.0060000 5.3920002
0.0080000 5.0939999
0.0100000 4.8660002
0.0120000 4.5679998
0.0140000 4.2519999
0.0160000 3.9200001
0.0180000 3.5139999
0.0200000 3.1500001
0.0220000 2.7460001
0.0240000 2.3239999
0.0260000 1.9040000
0.0280000 1.3900000
0.0300000 0.9140000
0.0320000 0.4040000
0.0340000 0.0000000
0.0360000 -0.4720000
0.0380000 -0.9740000
0.0400000 -1.3099999
为了消除其他频率波或者噪音的影响，现在希望通过滤波将杂质除去，小弟以前没有用过matlab滤波功能，请问高手具体如何操作啊？感激不尽
请问heaventian，能稍微详细一点吗
低通滤波器参数：Fs=8000，fp=2500，fs=3500，Rp=1dB，As=30dB，其他滤波器可以通过与低通之间的映射关系实现。
%%模拟滤波器
%巴特沃斯——滤波器设计
wp=2*pi*2500;ws=2*pi*3500;Rp=1;As=30;
[N,wc]=buttord(wp,ws,Rp,As,'s')%计算率波器的阶数和3dB截止频率
[B,A]=butter(N,wc,'s');%计算滤波器系统函数分子分母多项式
fk=0:800/512:8000;wk=2*pi*fk;
Hk=freqs(B,A,wk);
figure
plot(fk/1000,20*log10(abs(Hk)));
grid on,xlabel('频率（kHz）'),ylabel('幅度（dB）')
title('巴特沃斯模拟滤波器')
axis([0,4,-35,5])
%%
%切比雪夫I——滤波器设计
wp=2*pi*2500;ws=2*pi*3500;Rp=1;As=30;
[N1,wp1]=cheb1ord(wp,ws,Rp,As,'s')%计算切比雪夫滤波器的阶数和通带边界频率
[B1,A1]=cheby1(N1,Rp,wp1,'s');%计算滤波器系统函数分子分母多项式
fk=0:800/512:8000;wk=2*pi*fk;
Hk=freqs(B1,A1,wk);figure,
plot(fk/1000,20*log10(abs(Hk)));
grid on,xlabel('频率（kHz）'),ylabel('幅度（dB）')
title