exm041323_2.m

来自「matlab mfiles实例分析」· M 代码 · 共 13 行

[t,y,S,aquad,bquad]=fzzyquad(0,2,15);	%计算0到15次谐波		<1>
AA=abs(A_sym);AA(AA<1e-10)=NaN;	%在那些实际应为0的很小数处，考虑相对误差
										%没有意义，故予以排除。
BB=abs(B_sym);BB(BB<1e-10)=NaN;
A_quad=aquad(1:7)			%数值计算所得偶函数分量系数
B_quad=bquad(1:7)			%数值计算所得奇函数分量系数
asq=abs(A_quad-A_sym)./AA	%偶函数系数与符号计算值的相对误差
bsq=abs(B_quad-B_sym)./BB	%奇函数系数与符号计算值的相对误差  

SS=[S;y]';ribbon(t,SS);%把原始波形列在图形最前方
xlabel('谐波次数\rightarrow'),ylabel('t\rightarrow')
view([46,38]),colormap(jet),shading flat,light,lighting gouraud