fzzysym.m

matlab mfiles实例分析

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);