代凤英 16 通信 16020023004clear;%信号Fs=500;T=1/Fs;N=512;t=(0:N-1)*T;x=cos(2*pi*50*t)+1.5*sin(2*pi*200*t);subplot(2,1,1)plot(t,x)title('输入信号');%系数矩阵M=log2(N);xishu=zeros(N,M+1);%计算旋转因子k=0:1:N-1;W=exp(-1i*2*pi/N);Wk=W.^k;%先码位倒置得到序号oder=zeros(N,2);oder(:,1)=0:1:N-1;for i=1:1:N oder(i,2)=bin2dec(fliplr(dec2bin(oder(i,1),M)));end%得到序号后变换位置x_=zeros(1,N);for i=1:1:N x_(i)=x(oder(i,2)+1);endxishu(:,1)=x_;%循环 M 次，计算得到结果for i=0:1:M-1%i 为级数 step=2^i; number=N/2^(i+1); %每次计算之前对系数...
