永发信息网

在线跪求,想用matlab给一大组数据进行滤波的,网上找到了滤波函数,我的数据在一个文本文档里,是一堆采

答案:2  悬赏:70  手机版
解决时间 2021-03-20 18:04
样得到的离散数据。我要怎么把我的数据带进以下的程序啊?求大侠了。呜呜呜 程序:
利用fir1函数及hamming窗设计低通滤波器:
f1=100;f2=200;%待滤波正弦信号频率
fs=2000;%采样频率
m=(0.3*f1)/(fs/2);%定义过度带宽
M=round(8/m);%定义窗函数的长度
N=M-1;%定义滤波器的阶数
b=fir1(N,0.5*f2/(fs/2));%使用fir1函数设计滤波器
%输入的参数分别是滤波器的阶数和截止频率
figure(1)
[h,f]=freqz(b,1,512);%滤波器的幅频特性图
%[H,W]=freqz(B,A,N)当N是一个整数时函数返回N点的频率向量和幅频响应向量
plot(f*fs/(2*pi),20*log10(abs(h)))%参数分别是频率与幅值
xlabel('频率/赫兹');ylabel('增益/分贝');title('滤波器的增益响应');
figure(2)
subplot(211)
t=0:1/fs:0.5;%定义时间范围和步长
s=sin(2*pi*f1*t)+sin(2*pi*f2*t);%滤波前信号
plot(t,s);%滤波前的信号图像
xlabel('时间/秒');ylabel('幅度');title('信号滤波前时域图');
subplot(212)
Fs=fft(s,512);%将信号变换到频域
AFs=abs(Fs);%信号频域图的幅值
f=(0:255)*fs/512;%频率采样
plot(f,AFs(1:256));%滤波前的信号频域图
xlabel('频率/赫兹');ylabel('幅度');title('信号滤波前频域图');
figure(3)
sf=filter(b,1,s);%使用filter函数对信号进行滤波
%参数分别为滤波器系统函数的分子和分母多项式系数向量和待滤波信号输入
subplot(211)
plot(t,sf)%滤波后的信号图像
xlabel('时间/秒');ylabel('幅度');title('信号滤波后时域图');
axis([0.2 0.5 -2 2]);%限定图像坐标范围
subplot(212)
Fsf=fft(sf,512);%滤波后的信号频域图
AFsf=abs(Fsf);%信号频域图的幅值
f=(0:255)*fs/512;%频率采样
plot(f,AFsf(1:256))%滤波后的信号频域图
xlabel('频率/赫兹');ylabel('幅度');title('信号滤波后频域图');
最佳答案
可以直接读取你的文本文档的,如[A]=textread('d:\My Documents\a.txt','%f');
全部回答
建议你试一下巴特沃斯(butterworth)滤波器吧,先二阶,再到4阶,应该能满足你的要求。
我要举报
如以上问答信息为低俗、色情、不良、暴力、侵权、涉及违法等信息,可以点下面链接进行举报!
大家都在看
有什么测试能知道自己的女朋友爱你有多深?
东莞寄衣服到西安
新光盘无法完成格式化 新光盘如何格式化吗
滨州高新恒晟窗业有限公司地址在什么地方,想
各国著名广播电台有哪些?
血小板分布宽度偏高:20.9fL.正常吗
海南省东方市什么时候建的市
以前学校电信宽带用校园极速卡可以开无线网,
双眉紧蹙的意思
阳光广场王府井购物中心南2门地址在什么地方
50元人民币尾数四个1的直钱吗
三月份描写春天汾河景色的作文
发动机型号是什么意思
优品惠购在什么地方啊,我要过去处理事情
银行卡用来扣驾照分之后会绑定吗
推荐资讯
S206/Y118(路口)怎么去啊,有知道地址的么
用油煎鱼油中起泡是怎么回事?
2300左右最好的笔记本电脑
海澜之家东明服装店我想知道这个在什么地方
在专卖店买的匡威帆布鞋为什么鞋头大小不一,
红萝卜可和红枣连子白果构杞一起食用吗
剁椒凉皮(辣味)的做法,剁椒凉皮(辣味)怎
静怡爱情诗
03款中华俊杰冷车油门加不上去是怎么回事
一笔画出EXO
孩子的梳子怎么选,哪种合适啊?
每日工作总结是什么
正方形一边上任一点到这个正方形两条对角线的
阴历怎么看 ?