lisanfuliyebianhuan.rar_离散傅里叶

1. 信号傅里叶分解代码： （1） 离散傅里叶变换 N=256;dt=0.02;%数据的个数和采样间隔 n=0:N-1; t=n*dt;%序号序列和时间序列 x=sin(2*pi*t);%合成信号 m=floor(N/2)+1;%分解a，b的最大序号值为分解的N/2个参数加参数a0 a=zeros(1,m);b=zeros(1,m);%产生a，b两个为零的序列 for k=0:m-1, for i=0:N-1 a(k+1)=a(k+1)+2/N*x(i+1)*cos(2*pi*k*i/N); b(k+1)=b(k+1)+2/N*x(i+1)*sin(2*pi*k*i/N); %MATLAB中的数组序号只能从1开始 end c(k+1)=sqrt(a(k+1)^2+b(k+1)^2); end subplot(211);plot(t,x);title('原始信号');xlabel('时间/s'); f=(0:m-1)/(N*dt); subplot(212);plot(f,c); title('Fourier变换’)；xlabel('频率/Hz');ylabel('振幅')