function?[Toriginal_spectrumprocessed_spectrumfrequency_spectrumselected_index]=find_periods（xnumber_of_harmonicspeak_diff_samples）

T=[];
selected_index=[];
data=x;

Data=fftshift（fft（datamax（length（data）1023）））;
Fmax=0.5;
M=（length（Data）-1）/2;
f_interval=linspace（0FmaxM+1）;
Selected_Data=Data（M+1:end）;


Temp=Selected_Data;
for?i=1:number_of_harmonics
????[~n1]=max（abs（Temp））;
????try?
????????Temp（n1-peak_diff_samples:n1+peak_diff_samples）=0;?????%?ignore?the?maximum?close?to?the?one?that?you?just?found
????catch
????????Temp（n1）=0;
????end
????T（i）=1/f_interval（n1）;
????selected_index（i）=n1;
end


original_spectrum=Selected_Data;
processed_spectrum=Temp（1:length（original_spectrum））;
frequency_spectrum=f_interval;

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件???????10960??2018-08-12?14:47??TSAF\correlation_image.mat
?????文件???????11343??2018-08-12?14:47??TSAF\data_image.mat
?????文件?????????777??2018-08-12?14:47??TSAF\find_periods.m
?????文件?????????373??2018-08-12?14:47??TSAF\find_period_candidate.m
?????文件????????1304??2018-08-12?14:47??TSAF\license.txt
?????文件???????14067??2018-08-12?14:47??TSAF\model_image.mat
?????文件????????9146??2018-08-12?14:47??TSAF\test?folder\accidents.xlsx
?????文件???????12917??2018-08-12?14:47??TSAF\test?folder\energy.xlsx
?????文件???????17866??2018-08-12?14:47??TSAF\test?folder\test.xlsx
?????文件???????95725??2018-08-12?14:47??TSAF\TSAF.fig
?????文件???????48539??2018-08-12?14:47??TSAF\TSAF.m
?????文件????????1480??2018-08-12?14:47??license.txt