用matlab进行信号滤波.doc

第二次作业姓名:陈贤波对007hz进行滤波clc;clear;fs=25000;%采样频率Ts=1/fs;%周期signal=csvread('E:\matlab\信号处理PPT课件\作业2\6\',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(*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)');对80Hz进行滤波clc;clear;fs=25000;%采样频率Ts=1/fs;%周期signal=csvread('E:\matlab\信号处理PPT课件\作业2\6\',0,4);%载入信号N=length(signal);%信号长度t=0:Ts:(N-1)*Ts;%时间轴ws=400*2*pi/fs;%归一化wp=500*2*pi/fs;%归一化deltaw=wp-ws;%通带长度n=ceil(*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)');对315Hz进行滤波clc;clear;fs=100000;%采样频率Ts=1/fs;%周期si