题主给出的代码，为什么会画出来一个奇怪的图形？
分析了题主提供的代码，由于运行得到的方程有奇点，所以用solve函数得到是与本题不符的复数解（±1/z1^(1/2)），那如何解决呢？可以通过下列方法来求解：
1、给出theta的范围值，即
theta0=pi/10:pi/10:2*pi
2、利用for循环语句和fsolve函数联合求解，某theta值对应于的r值。即
for i=1:length(theta0)
theta=theta0(i)
。。。。。。
fun=@(r) 表达式
r_sln=fsolve(fun,pi/10)
r0(i)=r_sln
end
3、利用polar极坐标绘图函数，绘制其函数的图形。即
polar(theta0,r0)
4、按上述要求，编程运行可以得到如下图形。
8点ditfft蝶形图的对应关系是：1. 原位运算 在DIT-FFT的蝶形图中,取第m级且两输入节点分别在第k、j行的蝶形为例,讨论DIT-FFT的原位运算规律。
2. . 倒序规律 从图3.5可以看出,按原位计算时,蝶形图的输出正好是自然顺序X(0),X(1),...,X(7),
基2FFT的蝶形图对信号进行分析和处理时最常用的工具之一。在200多年前法国数学家、物理学家傅里叶提出后来以他名字命名的傅里叶级数之后，用DFT这个工具来分析信号就已经为人们所知。历史上最伟大的数学家之一。
它是根据离散傅氏变换的奇、偶、虚、实等特性，对离散傅立叶变换的算法进行改进获得的。它对傅氏变换的理论并没有新的发现，但是对于在计算机系统或者说数字系统中应用离散傅立叶变换，可以说是进了一大步。
扩展资料：
FFT的基本思想是把原始的N点序列，依次分解成一系列的短序列。充分利用DFT计算式中指数因子 所具有的对称性质和周期性质，进而求出这些短序列相应的DFT并进行适当组合，达到删除重复计算，减少乘法运算和简化结构的目的。
此后，在这思想基础上又开发了高基和分裂基等快速算法，随着数字技术的高速发展，1976年出现建立在数论和多项式理论基础上的维诺格勒傅里叶变换算法(WFTA）和素因子傅里叶变换算法。它们的共同特点是，当N是素数时，可以将DFT算转化为求循环卷积，从而更进一步减少乘法次数，提高速度。
