用matlab进行信号滤波

《用matlab进行信号滤波》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、第二次作业姓名：陈贤波对007hz进行滤波clc;clear;fs=25000;%采样频率Ts=1/fs;%周期signal=csvread('E:matlab信号处理PPT课件作业2607.csv',0,3);%载入信号N=length(signal);%信号长度t=0:Ts:(N-1)*Ts;%时间轴ws=60*2*pi/fs;%归一化wp=100*2*pi/fs;%归一化deltaw=wp-ws;%通带长度n=ceil(1.8*pi/deltaw)+1;%计算滤波器阶数wn=(ws+wp)/pi/2

2、;%计算截止频率[b,a]=fir1(n,wn,boxcar(n+1));%一定是n+1y=filter(b,a,signal);subplot(2,1,1);plot(t,signal);xlabel('时间(s)');subplot(2,1,2);plot(t,y);xlabel('时间(s)');对80Hz进行滤波clc;clear;fs=25000;%采样频率Ts=1/fs;%周期signal=csvread('E:matlab信号处理PPT课件作业26315.csv',0,4);%载入信号N=le

3、ngth(signal);%信号长度t=0:Ts:(N-1)*Ts;%时间轴ws=400*2*pi/fs;%归一化wp=500*2*pi/fs;%归一化deltaw=wp-ws;%通带长度n=ceil(1.8*pi/deltaw)+1;%计算滤波器阶数wn=(ws+wp)/pi/2;%计算截止频率[b,a]=fir1(n,wn,boxcar(n+1));%一定是n+1y=filter(b,a,signal);subplot(2,1,1);plot(t,signal);xlabel('时间(s)');subplot(2

4、,1,2);plot(t,y);xlabel('时间(s)');对315Hz进行滤波clc;clear;fs=100000;%采样频率Ts=1/fs;%周期signal=csvread('E:matlab信号处理PPT课件作业26315.csv',0,4,[0,4,999,4]);%载入信号N=length(signal);%信号长度t=0:Ts:(N-1)*Ts;%时间轴ws=1500*2*pi/fs;%归一化wp=2000*2*pi/fs;%归一化deltaw=wp-ws;%通带长度n=ceil(1.8

5、*pi/deltaw)+1;%计算滤波器阶数wn=(ws+wp)/pi/2;%计算截止频率[b,a]=fir1(n,wn,boxcar(n+1));%一定是n+1y=filter(b,a,signal);subplot(2,1,1);plot(t,signal);xlabel('时间(s)');subplot(2,1,2);plot(t,y);xlabel('时间(s)');