作业帮FFT MATLAB 代码出现一点小问题,clear; clc;N=16;%离散数组长度n=0:N-1;%时域采样xn=cos(pi*n/6);%产生离散数组M=nextpow2(N);%计算常量M,NA=[xn,zeros(1,N-length(xn))];disp('输入到各储存单元的数据：'),disp(A);J=0;%给
来源：学生作业帮助网 编辑：作业帮 时间：2024/11/28 02:21:49
![FFT MATLAB 代码出现一点小问题,clear; clc;N=16;%离散数组长度n=0:N-1;%时域采样xn=cos(pi*n/6);%产生离散数组M=nextpow2(N);%计算常量M,NA=[xn,zeros(1,N-length(xn))];disp('输入到各储存单元的数据：'),disp(A);J=0;%给](/uploads/image/z/11295530-26-0.jpg?t=FFT+MATLAB+%E4%BB%A3%E7%A0%81%E5%87%BA%E7%8E%B0%E4%B8%80%E7%82%B9%E5%B0%8F%E9%97%AE%E9%A2%98%2Cclear%3B+clc%3BN%3D16%3B%25%E7%A6%BB%E6%95%A3%E6%95%B0%E7%BB%84%E9%95%BF%E5%BA%A6n%3D0%3AN-1%3B%25%E6%97%B6%E5%9F%9F%E9%87%87%E6%A0%B7xn%3Dcos%28pi%2An%2F6%29%3B%25%E4%BA%A7%E7%94%9F%E7%A6%BB%E6%95%A3%E6%95%B0%E7%BB%84M%3Dnextpow2%28N%29%3B%25%E8%AE%A1%E7%AE%97%E5%B8%B8%E9%87%8FM%2CNA%3D%5Bxn%2Czeros%281%2CN-length%28xn%29%29%5D%3Bdisp%28%27%E8%BE%93%E5%85%A5%E5%88%B0%E5%90%84%E5%82%A8%E5%AD%98%E5%8D%95%E5%85%83%E7%9A%84%E6%95%B0%E6%8D%AE%EF%BC%9A%27%29%2Cdisp%28A%29%3BJ%3D0%3B%25%E7%BB%99)
FFT MATLAB 代码出现一点小问题,clear; clc;N=16;%离散数组长度n=0:N-1;%时域采样xn=cos(pi*n/6);%产生离散数组M=nextpow2(N);%计算常量M,NA=[xn,zeros(1,N-length(xn))];disp('输入到各储存单元的数据：'),disp(A);J=0;%给
FFT MATLAB 代码出现一点小问题,
clear; clc;
N=16;%离散数组长度
n=0:N-1;%时域采样
xn=cos(pi*n/6);%产生离散数组
M=nextpow2(N);%计算常量M,N
A=[xn,zeros(1,N-length(xn))];
disp('输入到各储存单元的数据：'),disp(A);
J=0;%给倒序数赋初值
for I=0:N-1%使原来数组次序改变
if I=K;
J=J-K; K=K/2;
end
J=J+K;
end
disp('倒序后各储存单元的数据：'),disp(A);%将变换次序以后的数组赋%值给A
WN=exp(-j*2*pi/N);
for L=1:M%%最外层循环为蝴蝶算法的总行数
disp('运算级次：'),disp(L);
B=2^(L-1);%M为第r行的组数
for R=0:B-1;
P=2^(M-L)*R;
for K=R:2^L:N-2;
T=A(K+1)+A(K+B+1)*WN^P;
A(K+B+1)=A(K+1)-A(K+B+1)*WN^P;
A(K+1)=T
end
end
disp('本级运算后储存单元的数据：'),disp(A);
end
disp('输出各储存单元的数据：'),Kx=A%输出运算结果
%telapsed = toc（tstart）；
FFT MATLAB 代码出现一点小问题,clear; clc;N=16;%离散数组长度n=0:N-1;%时域采样xn=cos(pi*n/6);%产生离散数组M=nextpow2(N);%计算常量M,NA=[xn,zeros(1,N-length(xn))];disp('输入到各储存单元的数据：'),disp(A);J=0;%给
clear; clc;
N=16; %离散数组长度
n=0:N-1; %时域采样
xn=cos(pi*n/6); %产生离散数组
M=nextpow2(N); %计算常量M,N
A=[xn,zeros(1,N-length(xn))];
disp('输入到各储存单元的数据：'),disp(A);
A=A(:,N:-1:1);
disp('倒序后各储存单元的数据：'),disp(A); %将变换次序以后的数组赋%值给A
WN=exp(-j*2*pi/N);
for L=1:M%%最外层循环为蝴蝶算法的总行数
disp('运算级次：'),disp(L);
B=2^(L-1);%M为第r行的组数
for R=0:B-1;
P=2^(M-L)*R;
for K=R:2^L:N-2;
T=A(K+1)+A(K+B+1)*WN^P;
A(K+B+1)=A(K+1)-A(K+B+1)*WN^P;
A(K+1)=T
end
end
disp('本级运算后储存单元的数据：'),disp(A);
end
disp('输出各储存单元的数据：'),Kx=A%输出运算结果
%telapsed = toc（tstart）；