function [A_sym,B_sym]=fzzysym(T,Nf,Nn)
% 采用符號計算求[0,T]內時間函數的三角級數展開系數。
%				函數的輸入輸出都是數值量
% 	Nf		諧波的階數
%	Nn		輸出數據的準確位數
%	A_sym	第1元素是直流項，其后元素依次是1,2,3...次諧波cos項展開系數
%	B_sym	第2,3,4,...元素依次是1,2,3...次諧波sin項展開系數
syms ttt n
if nargin<2;Nf=6;end
if nargin<3;Nn=32;end
yy=time_fun_s(ttt);
A0=int(yy,ttt,0,T)/T;
As=int(yy*cos(2*pi*n*ttt/T),ttt,0,T);
Bs=int(yy*sin(2*pi*n*ttt/T),ttt,0,T);
A_sym(1)=double(vpa(A0,Nn));
for k=1:Nf
   A_sym(k+1)=double(vpa(subs(As,n,k),Nn));
   B_sym(k+1)=double(vpa(subs(Bs,n,k),Nn));
end
%-------------------------------------------
function yy=time_fun_s(ttt)
% 該函數是fzzysym.m的子函數。它由符號變量和表達式寫成。
y1=sym('Heaviside(ttt-0.5)')*(ttt-0.5);
yy=y1-sym('Heaviside(ttt-1.5)')*((ttt-1.5)+1);