clc
clear?all
%%?FFT采樣
fs=3200;
N=64;???????????????????????????????????????????????????????????????%采樣頻率和數(shù)據(jù)點(diǎn)數(shù)
n=0:N-1;
t=n/fs;?????????????????????????????????????????????????????????????%時(shí)間序列
x=0.8*sin（4.5*2*pi*50*t+pi/4）+0.3*sin（2*2*pi*50*t+pi/5）;????????????%信號(hào)
y=fft（xN）;?????????????????????????????????????????????????????????%對信號(hào)進(jìn)行快速Fourier變換
f=fs*（0:length（y）/2）/length（y）;
z=angle（y）*180/pi;
mag=abs（y）;??
fz=mag（1:length（y）/2+1）*2/N;????????????????????????????????????????%求得Fourier變換后的振幅
figure（1）
plot（ffz）;?????????????????????????????????????????????????????????%繪出Nyquist頻率之前隨頻率變化的振幅
xlabel（‘頻率/Hz‘）;
ylabel（‘振幅‘）;
title（‘FFT采樣‘）;
grid?on;

%%?BP神經(jīng)網(wǎng)絡(luò)初始化
inputNums=1;
outputNums=1;
hideNums=4;
maxcount=3000;
sampleNums=N;
precision=1e-9;
xxl1=0.02;
xxl2=250;
count=1;
error=zeros（1count）;
errorp=zeros（1sampleNums）;
w1=rand（inputNumshideNums/2）;?????????????????????
w2=rand（inputNumshideNums/2）;????????????????????
b=[1*pi*2001*pi*120];

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2013-11-07?16:08??傅立葉BP諧波分析\
?????文件????????2441??2013-11-07?16:08??傅立葉BP諧波分析\FLY_BP.m