用matlab进行信号滤波

姓名：陈贤波

对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+1

y=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+1

y=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;%周期

signal=csvread('E:\\matlab\\信号处理PPT课件\\作业2\\6\\',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*pi/deltaw)+1;%计算滤波器阶数

wn=(ws+wp)/pi/2;%计算截止频率

[b,a]=fir1(n,wn,boxcar(n+1));%一定是n+1

y=filter(b,a,signal);

subplot(2,1,1);

plot(t,signal);

xlabel('时间(s)');

subplot(2,1,2);

plot(t,y);

xlabel('时间(s)');