b= zeros(len,1);%式系数 b(1)=1; fori=3:2:2*m+1 b(i)=1*b(i-2)/(i2)/(i-1) C= zeros(len,1);%乘积系数 c(2)=1 for i=4: 2: len c(i)=4*c(i-2)/(i-1)/(i-2) end x=fft(a,len);% Fourier变换 fft(b,len);% curler变换 z=ifft(z1);% Fourier逆变换 maxerror=O for i=l: len e=abs(z(i-c(i)) if e>maxerror maxerror-e db=zeros(len,1);%乘式系数 b(1)=1; for i=3:2:2*m+1 b(i)=-1*b(i-2)/(i-2)/(i-1); end c=zeros(len,1);%乘积系数 c(2)=1; for i=4:2:len c(i)=-4*c(i-2)/(i-1)/(i-2); end x=fft(a,len);%Fourier 变换 y=fft(b,len);%Fourier 变换 z1=x.*y; z=ifft(z1);%Fourier 逆变换 maxerror=0; for i=1:len e=abs(z(i)-c(i)); if e>maxerror maxerror=e; end end 9