clear?all;
clc;

load?data.txt;???%原信號

subplot（311）
plot（data）;
title（‘原始信號‘）;
grid?on;
Fs=999;?????%采樣頻率
fc=50;??????%50Hz
x=0:1/Fs:1;??
y1=0.2*sin（2*pi*fc*x）‘+data;??
subplot（312）
plot（xy1）;
title（‘加入50Hz正弦噪聲‘）;
grid?on;

fp=[45?55];fs=[49?51];?????????%fp?通帶上下截止頻率??%fs?阻帶上下截止頻率
wp=fp*2/Fs;ws=fs*2/Fs;??%wp?通帶邊界頻率?ws?阻帶邊界頻率
rp=3;???%通帶最大衰減
as=14;??%阻帶最小衰減
[nwc]=buttord（wpwsrpas）;???%?n?濾波器階數?wc?3dB截止頻率?
[ba]=butter（nwc‘stop‘）;????%b?系統函數分子
y2=filter（bay1）;????????????%一維數字濾波器??輸入X為濾波前序列，Y為濾波結果序列，B/A?提供濾波器系數，B為分子，?A為分母?
subplot（313）
plot（xy2）;
title（‘巴特沃斯‘）;
grid?on?

%濾波波形頻譜
figure（2）
L=1000;
NFFT?=?2^nextpow2（L）;
Y?=?fft（y2NFFT）/L;
f?=?Fs/2*linspace（01NFFT/2+1）;
subplot（311）?
plot（f2*abs（Y（1:NFFT/2+1）））?
title（‘頻譜‘）


[h1w1]=freqz（ba1000）;?????%?w是頻率，n是FFT點數??計算系統的頻率響應???????????????
subplot（312）?????
plot（w1/piabs（h1））????????%幅度特性曲線
title（‘離散系統幅頻特性曲線‘）
subplot（313）?????
plot（w1/piangle（h1））;???????%頻率特性曲線
title（‘離散系統相頻特性曲線‘）

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????9002??2013-03-04?21:24??data.txt

?????文件???????1199??2013-03-04?21:27??test1.m

-----------?---------??----------?-----??----

????????????????10201????????????????????2