采用winograd fft 算法写的radix5的代码，用MATLAB验证不对

Question

请问哪里出了问题啊

clc;
close all;
clear all;
%x=rand(1,5);
x=[1 2 3 4 5];
y=zeros(1,5);
% y=radix_5(x(1),x(2),x(3),x(4),x(5));
%第一级
t1=x(2)+x(5);
t2=x(3)+x(4);
t3=x(2)-x(5);
t4=x(3)-x(4);
%第二级
t5=t1+t2;
t6=t1-t2;
t7=t3+t4;
%第三级
c0=1;
c1=(1/2*(cos(2*pi/5)+cos(4*pi/5))-1);
c2=1/2*(cos(2*pi/5)-cos(4*pi/5));
c3=(i)*sin(4*pi/5);
c4=(i)*(sin(4*pi/5)+sin(2*pi/5));
c5=(i)*(sin(2*pi/5)-sin(4*pi/5));
m0=x(1)+t5;%乘不乘c0无所谓
m1=t5*c1;
m2=t6*c2;
m3=t7*c3;
m4=t4*c4;
m5=t3*c5;
%第四级
s1=m0+m1;
s2=s1+m2;
s4=s1-m2;
s3=m3-m4;
s5=m3+m5;
%第五级
y=zeros(1,5);
y(1)=m0;
y(2)=s2+s3;
y(5)=s2-s3;
y(3)=s4+s5;
y(4)=s4-s5;

Answer 1

上午好☀️☀️☀️️
本答案参考ChatGPT-3.5

你的代码主要使用了Winograd FFT算法来计算根号五的序列变换，但是并没有直接检查它的正确性。这是因为这个方法主要用于特定的计算任务，如大整数除法和平方根运算等，而不是用于一般的硬件问题。

然而，你可以通过添加一些条件来检测这段代码是否有效：

% 添加条件，检查输入是否为实数
if ~isreal(x) || ~isreal(y)
    disp('Invalid input. x and y must be real numbers.');
    return;
end

% 添加条件，检查是否满足递归的终止条件
if max(abs(t1)) > 1 || max(abs(t2)) > 1 || max(abs(t3)) > 1 || max(abs(t4)) > 1 || max(abs(t5)) > 1 || ...
       max(abs(c0)) > 1 || max(abs(c1)) > 1 || max(abs(c2)) > 1 || max(abs(c3)) > 1 || max(abs(c4)) > 1 || ...
       max(abs(c5)) > 1 || max(abs(s1)) > 1 || max(abs(s2)) > 1 || max(abs(s3)) > 1 || max(abs(s4)) > 1 || ...
       max(abs(s5)) > 1 || max(abs(y(1))) > 1 || max(abs(y(2))) > 1 || max(abs(y(5))) > 1)
    disp('Invalid input. The root of five is not valid.');
    return;
end

以上这些条件可以根据实际情况进行修改。如果代码运行正常，并且符合上述要求，那么就可以确认这段代码是有效的。否则，就需要重新考虑你的硬件环境和应用需求，可能需要对代码进行相应的调整或者重构。