matlab信号分析滤波求包络

%原信号时域图
clear all;
clc;
load data.mat;
fs=1000;
n=length（ydata）;
t=tdata;
x=ydata;
figure（1）;
subplot（211）;
plot（tx）;
grid on ;
axis（[0 0.1 -2 2]）;
xlabel（‘时间‘）;
ylabel（‘加速度‘）;
title（‘Y信号时域图‘）;

%求Y信号相关参数
rms=sqrt（sum（x.^2））/n
bf=max（x）;
cf=bf/rms
k=kurtosis（x）

%Y信号FFT后的波形图
fmin=0.2;
fmax=10;
nfft=2^nextpow2（n）;
ni=round（fmin*nfft/fs+1）;
na=round（fmax*nfft/fs+1）;
y=fft（xnfft）;
a=zeros（1nfft）;
a（ni:na）=y（ni:na）;
a（nfft-na+1:nfft-ni+1）=y（nfft-na+1:nfft-ni+1）;
Y=ifft（anfft）;
y2=（real（Y（1:n）））‘;
subplot（212）;
plot（tdatay2）;
grid on
xlabel（‘时间‘）;
ylabel（‘加速度‘）;
title（‘Y信号时域图（滤波后）‘）;
axis（[0 0.1 -2 2]）
 
figure（2）
f=fs*（0:n/2-1）/n;
plot（fabs（y（1:n/2）））;
grid on 
xlabel（‘频率‘）;
ylabel（‘振幅‘）;
title（‘Y信号截断频谱图‘）

%进行IIR带通滤波器设计
rp=1;
rs=30;
wp=[50400];
ws=[40410];
ft=‘bandpass‘;
[Nwn]=cheb2ord（wp/（fs/2）ws/（fs/2）rprs）;
[bbab]=cheby2（Nrswnf