clc;
clear;
close?all;

%%?參數設置
N=1024;
fs=300;
ff=input（‘input?f/Hz\n‘）
Ts=1/fs;
n=0:N-1;
fp=n/N*fs;

%%?時、頻波形
[FileNamePathName]?=?uigetfile（‘USB.dat‘‘Select?the?USB.dat?file‘）;
f?=?fullfile（PathNamefilesepFileName）;
fid?=?fopen（f‘r‘）;
data?=?fscanf（fid‘%x‘）;
fclose（fid）;
data?=?data（1:2:end）*256?+?data（2:2:end）;
figure（1）;
plot（data）;?????
data=data-mean（data）;
b=max（data）;
a=min（data）;
data=data/（b-a）;
F=fft（data）;
figure（2）;
plot（abs（F））;

%%??濾波后的時、頻波形
fp=0.05;fc=120;As=80;Ap=1;Fs=300;
wc=2*pi*fc/Fs;?wp=2*pi*fp/Fs;
wd=wc-wp;
beta=0.1102*（As-8.7）;
N=ceil（（As-7.95）/2.286/wd）;
wn=?kaiser（N+1beta）;?
ws=（wp+wc）/2/pi;
b=fir1（Nwswn）;
x=fftfilt（bdata）;
X=fft（x）;
figure（3）;
plot（x）;title（‘低通濾波后的信號‘）;
figure（4）;
plot（abs（X））;

%%?信噪比計

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2019-03-28?15:01??程序+采集數據\
?????文件????????8192??2017-08-29?15:17??程序+采集數據\0.8hz.dat
?????文件????????8192??2017-08-28?14:58??程序+采集數據\1.5hz.dat
?????文件??????448539??2019-03-28?14:59??程序+采集數據\1.docx
?????文件????????8192??2017-08-29?15:12??程序+采集數據\1hz.dat
?????文件????????8192??2017-08-29?15:04??程序+采集數據\2hz.dat
?????文件????????8192??2017-08-30?10:34??程序+采集數據\400.dat
?????文件????????8192??2017-08-30?11:00??程序+采集數據\500.dat
?????文件????????8192??2017-08-29?15:48??程序+采集數據\800m.dat
?????文件????????8192??2017-08-30?11:06??程序+采集數據\USB.dat
?????文件????????8192??2017-08-29?11:20??程序+采集數據\USB_100.dat
?????文件????????3232??2017-09-01?13:57??程序+采集數據\a1.m