zhaoxuanyi 发表于 2009-9-24 20:29

傅立叶变换

请问这个计算傅立叶变换的程序为什么运行不出来？？
function =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 发表于 2009-10-7 22:01

回复 楼主 zhaoxuanyi 的帖子

恩，不错，谢谢楼主

ChaChing 发表于 2009-10-8 00:21

回复 沙发 renshudong113 的帖子

请问谢谢楼主什麽? 有灌水嫌疑!? :@)
页: [1]
查看完整版本: 傅立叶变换