function?[a0amnphn]?=?dfseries（xx?yy?tol）
%?????[a0amnphn]?=?dfseries（xx?yy?tol）
????%tol:total?terms?of?（amnphn）
????%a0:mean?term?a0/2
????%amnphn:amplitude?and?phase?of?an*cos（wt）+bn*sin（wt）
????if?nargin?==?2
????????tol=5;
????end
????xxn=length（xx）;
????ll=（xx（xxn）-xx（1））/2;
????a0=1/ll*trapz（xxyy）;
????an=zeros（tol1）;bn=an;
????for?ii=1:1:tol
????????an（ii）=1/ll*trapz（xxyy.*cos（ii*pi*xx/ll））;
????????bn（ii）=1/ll*trapz（xxyy.*sin（ii*pi*xx/ll））;
????end
????amn=（an.^2+bn.^2）.^0.5.*sign（bn）;phn=atan（an./bn）;
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件?????????551??2015-03-25?12:44??dfseries.m
?????文件????????2454??2015-03-25?12:21??dynmeshk.m