斯大林
估计你需要向量化代码中所有循环可指一种可能的转接法 减少二环入代码很难不重写逻辑即压缩第三方循环
索引化
并
宝斯
变量大小迭代希望它帮助
函数显示yss
N1=长度%计算Xn序列大小
NFFT=2QSeil(log2N1);%计算样本数补全
x=x0s(1nFFT-N1)序列完成0-nFFT元素
N=长度(x);%计算xn样本数
b=bin2dec(fliplr(dec2bin (0:1:nFFT-1)))+1;百分比重排序样本
x=xxb信号重排序
S=log2N%计算阶段数
半=1;
For级=1:S%循环算法Log2
index=0:(2^stage):(N-1);蝴蝶算法的每一阶段
N=0:(Half-1);蝴蝶计算并保存结果
spos=n+index+1%样本索引
pow=(2^(S-stage)).*n;%复杂乘法功率
w=exp((-1i)*(2*pi).*pow/N);%复杂乘法
a=x(pos)+x(pos+Half).*w;首端蝴蝶
b=x(pos)-x(pos+Half).*w;二分
x(sp)=a%保存前端蝴蝶结果
x(sp+Half)=b% 保存蝴蝶第二部分结果
半=2*半
结束%计算下个'半'值
Y=x;
函数结果保存
结束