傅立叶变换

zhaoxuanyi

请问这个计算傅立叶变换的程序为什么运行不出来？？
function [A_sym,B_sym]=fzzysym(T,Nf,Nn)
%nf谐波的阶数 nn 输出数据的准确位数 A cos展开的系数 B 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)=vpa(A0,Nn);
for k=1:Nf
    A_sym(k+1)=vpa(subs(As,n,k),Nn);
    B_sym(k+1)=vpa(subs(Bs,n,k),Nn);
end
function yy=time_fun_s(ttt)  %子函数
y1=sym('Heaviside(ttt-0.5)')*(ttt-0.5);
yy=y1-sym('Heavside(ttt-1.5)')*((ttt-1.5)+1);

renshudong113

恩，不错，谢谢楼主

ChaChing

请问谢谢楼主什麽? 有灌水嫌疑!? :@)