用matlab进行信号滤波

合集下载

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

第二次作业

姓名：陈贤波

对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)');