?
fp=500;fs=400?;%注：如果fp在fs區間內，則為帶通，反之為帶阻
?ap=2;as=30;
?fsa=4096;
?wp=fp/fsa*2;ws=fs/fsa*2;
[nwn]=buttord（wpwsapas）;%濾波器階數選擇函數，確定符合要求性能濾波器的最小階數和固有頻率
[BA]=butter（nwn）;
[Hw]=freqz（BA）;%?H為頻域傳遞函數，w為所要計算相應的頻率點組成的向量
figure（1）;subplot（211）;
length（w）
plot（w*fsa/（2*pi）abs（H））;?????????????%繪制帶通頻譜圖
xlabel（‘Hz‘）;ylabel（‘頻率響應幅度‘）;
title（‘低通濾波器‘）;
grid;axis?tight;
subplot（212）;plot（w*fsa/（2*pi）angle（H））;
xlabel（‘Hz‘）;ylabel（‘angel‘）;
title（‘相位特性‘）;
grid;axis?tight;
x=load（‘D:\debugging\matlabfiles\監測診斷試驗\數據\齒輪\19點20分14秒\通道1.txt‘）;
figure（2）
subplot（211）
t=（0:length（x）-1）/fsa;
plot（tx）;
xlabel（‘t‘）;ylabel（‘x（t）‘）;
title（‘濾波前連續信號‘）;
grid;axis?tight;
y=filter（BAx）;%時域濾波
figure（2）?
subplot（212）;
plot（ty）;
xlabel（‘t‘）;ylabel（‘y（t）‘）;
title（‘濾波后連續信號‘）;
grid;axis?tight;
figure（3）
N=length（x）;
X=fft（xN）/N;w=fsa/2*linspace（01N/2+1）;
subplot（211）;plot（wabs（X（1:N/2+1）））;??%繪制頻譜圖
xlabel（‘Hz‘）;ylabel（‘頻率響應幅度‘）;
title（‘濾波前頻譜圖‘）;
grid;axis?tight;
Y=fft（yN）/N;w=fsa/2*linspace（01N/2+1）;
subplot（212）;plot（wabs（Y（1:N/2+1）））;???????%繪制頻譜圖
xlabel（‘Hz‘）;ylabel（‘頻率響應幅度‘）;
title（‘濾波后頻譜圖‘）;
grid;axis?tight;

?

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2018-08-29?16:12??低通濾波器?matlab?實現\
?????文件??????????55??2018-08-27?11:33??低通濾波器?matlab?實現\【源碼使用必讀】.url
?????文件????????1338??2012-11-19?17:37??低通濾波器?matlab?實現\lowpass.m
?????文件??????237958??2012-11-19?17:50??低通濾波器?matlab?實現\帶通特性.bmp
?????文件??????237958??2012-11-19?17:45??低通濾波器?matlab?實現\低通濾波后.bmp