快速傅里叶离散变换FFT (更新中)

声明:参考了 y y c yyc yycblogPPT (from smwc) ,以及 w z r wzr wzrblog

目录

  • Part 1 多项式
  • Part 2 FFT概论
  • Part 3 点值与插值
  • Part 4 复数,单位根
  • Part 5

Part 1 多项式

  • 定义:对于有限数列 A 0 A_{0} A0~ n _{n} n ,形如 A ( x ) = a 0 + a 1 x + a 2 x 2 + . . . + a n x n A(x)=a_0+a_1x+a_2x^2+ ... +a_nx^n A(x)=a0+a1x+a2x2+...+anxn (即 A ( x ) = ∑ i = 0 n A i x i A(x)= \sum_{i=0}^{n} A_ix^i A(x)=i=0nAixi) 的函数称为多项式。记 A ( x ) A(x) A(x) x n x^n xn 项系数为 [ x n ] A ( x ) [x^n]A(x) [xn]A(x) ,简写为 A [ n ] A[n] A[n]
  • 运算:
    • 加法: C ( x ) = A ( x ) + B ( x ) = ∑ i = 0 n ( A i + B i ) x i C(x)=A(x)+B(x) = \sum_{i=0}^{n}(A_i+B_i)x^i C(x)=A(x)+B(x)=i=0n(Ai+Bi)xi
    • 减法: C ( x ) = A ( x ) − B ( x ) = ∑ i = 0 n ( A i − B i ) x i C(x)=A(x)-B(x) = \sum_{i=0}^{n}(A_i-B_i)x^i C(x)=A(x)B(x)=i=0n(AiBi)xi
    • 乘法: C ( x ) = A ( x ) × B ( x ) = ∑ i = 0 n ∑ j = 0 m A i B j x i + j C(x)=A(x) \times B(x) = \sum_{i=0}^{n} \sum_{j=0}^{m} A_iB_jx^{i+j} C(x)=A(x)×B(x)=i=0nj=0mAiBjxi+j
  • 卷积:某二元运算 ⊕ \oplus 的定义域和值域均为 N \mathbb{N} N ,数列 A , B A,B A,B ⊕ \oplus 卷积 为一个新数列 C C C ,有 C [ k ] = ∑ i ⊕ j = k A [ i ] B [ j ] C[k]= \sum_{i \oplus j=k} A[i]B[j] C[k]=ij=kA[i]B[j] ⊕ \oplus 可以代表任意运算符)。容易发现,多项式乘法就是“加法卷积”。

Part 2 FFT概论

对于两个多项式的乘法,有简单粗暴的 O ( n 2 ) O(n^2) O(n2) 级别方法。但这只适用于少部分题目,于是某个天才就发明的一种新方法,叫 Fast Fourier Transformation(FFT)。利用卷积思想,化乘为加,快速计算多项式乘法,只需 O ( n l o g n ) O(nlogn) O(nlogn)


Part 3 点值与插值

  • 点值:将多项式 A ( x ) A(x) A(x) 视为函数,对于常数 t 0 t_0 t0 A ( t 0 ) = ∑ i = 0 + ∞ A [ i ] t 0 i A(t_0)=\sum_{i=0}^{+\infty} A[i]t_0^i A(t0)=i=0+A[i]t0i A ( x ) A(x) A(x) t 0 t_0 t0 处的点值。
    • 作用:已知 A ( x ) , B ( x ) , C ( x ) A(x),B(x),C(x) A(x),B(x),C(x),可以快速判断 C ( x ) = A ( x ) B ( x ) C(x)=A(x)B(x) C(x)=A(x)B(x)
    • 核心:通过选择具体的常数,将复杂的多项式乘法转换成数的乘法,提升效率。
  • 插值:已知 A ( x ) A(x) A(x) 的若干点,求 A ( x ) A(x) A(x) 的系数序列。

将上面两种表示法组合起来,就有一个 “求值-插值法” ,得到多项式乘法一个新想法。即,把系数表达转换为点值表达 (DFT) -> 点值表达相乘 -> 把点值表达转换为系数表达 (IDFT)
这就是 F F T FFT FFT 的思路:

  1. 选择一个足够大的 n n n,选定 n n n 个不同的数 x 1 , x 2 , x 3 , . . . , x n x_1,x_2,x_3, ... ,x_n x1,x2,x3,...,xn
  2. 求出点值 A ( x 1 ) , A ( x 2 ) , . . . , A ( x n ) A(x_1),A(x_2),... ,A(x_n) A(x1),A(x2),...,A(xn) B ( x 1 ) , B ( x 2 ) , . . . , B ( x n ) B_(x_1),B(x_2),...,B(x_n) B(x1),B(x2),...,B(xn)
  3. 根据 C ( x i ) = A ( x i ) B ( x i ) C(x_i)=A(x_i)B(x_i) C(xi)=A(xi)B(xi) 得到 C ( x 1 ) , C ( x 2 ) , . . . , C ( x n ) C(x_1),C(x_2),...,C(x_n) C(x1),C(x2),...,C(xn) :
  4. 插值求出 C ( x ) C(x) C(x) 的系数;

Part 4 复数,单位根

考虑如何做 FFT 的第一步(即选择 n n n 个不同的数),一般直觉会选正整数等一些熟悉的数,但这没什么性质(加速不了多项式乘法)。于是单位根就登场了……

  • 复数:令虚单位 i i i 满足 i 2 = − 1 i^2=-1 i2=1,形如 a + b i a+bi a+bi ( a , b ∈ R a,b \in \mathbb{R} a,bR) 的数称为复数。
    将复数 z = a + b i z=a+bi z=a+bi 看作平面直角坐标系中的点 P ( a , b ) P(a,b) P(a,b) ,原点 O ( 0 , 0 ) O(0,0) O(0,0),称表示复数的平面直角坐标系为复平面 r = ∣ O P ∣ r=|OP| r=OP z z z模长,射线 O P OP OP x x x轴正半轴的夹角 θ \theta θ z z z辐角。有 z = r ( c o s θ + i z=r(cos \theta + i z=r(cosθ+i s i n θ ) sin \theta) sinθ)
    • 复数加法: ( a + b i ) + ( c + d i ) = ( a + c ) + ( b + d ) i (a+bi)+(c+di) = (a+c) + (b+d)i (a+bi)+(c+di)=(a+c)+(b+d)i
    • 复数减法: ( a + b i ) − ( c + d i ) = ( a − c ) + ( b − d ) i (a+bi)-(c+di) = (a-c) + (b-d)i (a+bi)(c+di)=(ac)+(bd)i
    • 复数乘法: ( a + b i ) × ( c + d i ) = a c − b d + ( a d + b c ) i (a+bi) \times (c+di) = ac-bd+(ad+bc)i (a+bi)×(c+di)=acbd+(ad+bc)i
    • 复数除法: ( a + b i ) (a+bi) (a+bi) / / / ( c + d i ) = a c + b d c 2 + d 2 + b c − a d c 2 + d 2 i (c+di)= \frac{ac+bd}{c^2+d^2} + \frac{bc-ad}{c^2+d^2}i (c+di)=c2+d2ac+bd+c2+d2bcadi
    • 复数的共轭: a + b i a+bi a+bi 的共轭为 a − b i a-bi abi
  • 单位根 n n n 次本原单位根记为 ω n \omega ^n ωn ( n ∈ N n \in \mathbb{N} nN),单位根为复数,满足 ω n = 1 \omega ^n=1 ωn=1 ,且 ω 0 , ω 1 , ω 2 , . . . , ω n − 1 \omega ^0,\omega ^1,\omega ^2 ,...,\omega ^{n-1} ω0,ω1,ω2,...,ωn1 互不相同。
    引入 单位圆(圆心在原点且半径为 1 1 1 的圆):单位圆
    记第 n n n n n n 次单位根为 ω n n − 1 。 \omega _n^{n-1}。 ωnn1 ω n k = k ( c o s 2 π n + i \omega _n^k=k(cos\frac{2\pi}{n}+i ωnk=k(cosn2π+i s i n 2 π n ) sin \frac{2\pi}{n}) sinn2π)
    ω n \omega _n ωn 是模长为 1 1 1 ,圆上的点所表示的复数模长都为 1 1 1 ,所以单位根一定在单位圆上。又有 ω n \omega _n ωn 的辐角为 2 π n \frac{2\pi}{n} n2π (即 1 n \frac{1}{n} n1 圆周),乘一次 ω n \omega _n ωn 相当于绕原点逆时针旋转 1 n \frac{1}{n} n1 圆周,乘 n n n 次恰好转完一圈回到原地,满足 “ ω n = 1 \omega ^n=1 ωn=1 ,且 ω 0 , ω 1 , ω 2 , . . . , ω n − 1 \omega ^0,\omega ^1,\omega ^2 ,...,\omega ^{n-1} ω0,ω1,ω2,...,ωn1 互不相同” 这一条件。
    w
  • 单位根的一些性质:
    1. ω n 0 = 1 \omega _n^0 =1 ωn0=1
    2. ω n k = ( ω n 1 ) k \omega _n^k = (\omega _n^1)^k ωnk=(ωn1)k
    3. ω n i × ω n j = ( ω n 1 ) i × ( ω n 1 ) j = ω n i + j \omega _n^i \times \omega_n^j =(\omega_n^1)^i \times (\omega_n^1)^j = \omega_n^{i+j} ωni×ωnj=(ωn1)i×(ωn1)j=ωni+j
    4. ω n − k = ω n n − k \omega_n^{-k} = \omega_n^{n-k} ωnk=ωnnk
    5. ω 2 n 2 k = ω n k = ω p n p k \omega_{2n}^{2k} = \omega_n^k = \omega_{pn}^{pk} ω2n2k=ωnk=ωpnpk
    6. ω n ( k + n / 2 ) = − ω n k \omega_{n}^{(k+n/2)} = - \omega_{n}^{k} ωn(k+n/2)=ωnk

Part 5

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.rhkb.cn/news/12939.html

如若内容造成侵权/违法违规/事实不符,请联系长河编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

w193基于Spring Boot的秒杀系统设计与实现

🙊作者简介:多年一线开发工作经验,原创团队,分享技术代码帮助学生学习,独立完成自己的网站项目。 代码可以查看文章末尾⬇️联系方式获取,记得注明来意哦~🌹赠送计算机毕业设计600个选题excel文…

Spark--如何理解RDD

1、概念 rdd是对数据集的逻辑表示,本身并不存储数据,只是封装了计算逻辑,并构建执行计划,通过保存血缘关系来记录rdd的执行过程和历史(当一个rdd需要重算时,系统会根据血缘关系追溯到最初的数据源&#xff…

旋钮屏设备物联网方案,ESP32-C3无线通信应用,助力设备智能化升级

在智能家居的浪潮中,旋钮屏以其独特的交互方式和便捷的操作体验,逐渐成为智能家电控制面板上的新宠儿。从智能冰箱、洗衣机到烤箱、空气炸锅等设备,旋钮屏的应用无处不在。 通过简单的旋转和按压操作,用户可以轻松调节温度、时间…

crewai框架第三方API使用官方RAG工具(pdf,csv,json)

最近在研究调用官方的工具,但官方文档的说明是在是太少了,后来在一个视频里看到了如何配置,记录一下 以PDF RAG Search工具举例,官方文档对于自定义模型的说明如下: 默认情况下,该工具使用 OpenAI 进行嵌…

嵌入式工程师必学(143):模拟信号链基础

概述: 我们每天使用的许多电子设备,以及我们赖以生存的电子设备,如果不使用电子工程师设计的实际输入信号,就无法运行。 模拟信号链由四个主要元件组成:传感器、放大器、滤波器和模数转换器 (ADC)。这些传感器用于检测、调节模拟信号并将其转换为适合由微控制器或其他数…

C++11详解(二) -- 引用折叠和完美转发

文章目录 2. 右值引用和移动语义2.6 类型分类(实践中没什么用)2.7 引用折叠2.8 完美转发2.9 引用折叠和完美转发的实例 2. 右值引用和移动语义 2.6 类型分类(实践中没什么用) C11以后,进一步对类型进行了划分&#x…

NeetCode刷题第21天(2025.2.4)

文章目录 114 Gas Station 加油站115 Hand of Straights 顺子之手116 Merge Triplets to Form Target 将 Triplelet 合并到 Form Target117 Partition Labels 分区标签118 Valid Parenthesis String 有效的括号字符串119 Insert Interval 插入间隔120 Merge Intervals 合并区间…

车载软件架构 --- 基于AUTOSAR软件架构的ECU开发流程小白篇

我是穿拖鞋的汉子,魔都中坚持长期主义的汽车电子工程师。 老规矩,分享一段喜欢的文字,避免自己成为高知识低文化的工程师: 简单,单纯,喜欢独处,独来独往,不易合同频过着接地气的生活…

Ollama本地搭建大模型

短短一夜之间,中国的AI大模型DeepSeek迅速崛起,成功引起了全球科技界的广泛关注。 deepSeek爆火时间线 DeepSeek大事记 技术突破与产品发布 2024年12月26日:DeepSeek-V3发布,知识类任务水平提升,生成吐字速度加快。…

C#结合html2canvas生成切割图片并导出到PDF

目录 需求 开发运行环境 实现 生成HTML范例片断 HTML元素转BASE64 BASE64转图片 切割长图片 生成PDF文件 小结 需求 html2canvas 是一个 JavaScript 库,它可以把任意一个网页中的元素(包括整个网页)绘制到指定的 canvas 中&#xf…

【通俗易懂说模型】线性回归(附深度学习、机器学习发展史)

🌈 个人主页:十二月的猫-CSDN博客 🔥 系列专栏: 🏀深度学习_十二月的猫的博客-CSDN博客 💪🏻 十二月的寒冬阻挡不了春天的脚步,十二点的黑夜遮蔽不住黎明的曙光 目录 1. 前言 2. …

C#面试常考随笔12:游戏开发中常用的设计模式【C#面试题(中级篇)补充】

C#面试题(中级篇),详细讲解,帮助你深刻理解,拒绝背话术!-CSDN博客 简单工厂模式 优点: 根据条件有工厂类直接创建具体的产品 客户端无需知道具体的对象名字,可以通过配置文件创建…

大模型的底层逻辑及Transformer架构

一、大模型的底层逻辑 1.数据驱动 大模型依赖海量的数据进行训练,数据的质量和数量直接影响模型的性能。通过大量的数据,模型能够学习到丰富的模式和规律,从而更好地处理各种任务。 2.深度学习架构 大模型基于深度学习技术,通常…

C++ 学习:深入理解 Linux 系统中的冯诺依曼架构

一、引言 冯诺依曼架构是现代计算机系统的基础,它的提出为计算机的发展奠定了理论基础。在学习 C 和 Linux 系统时,理解冯诺依曼架构有助于我们更好地理解程序是如何在计算机中运行的,包括程序的存储、执行和资源管理。这对于编写高效、可靠…

【C++】STL——list底层实现

目录 💕1.list的三个类介绍 💕2.list——节点类 (ListNode) 💕3.list——链表类 (List) 💕4.list——迭代器类(重点思考)(ListIterator) 💕5…

deepseek、qwen等多种模型本地化部署

想要在本地部署deepseek、qwen等模型其实很简单,快跟着小编一起部署吧 1 环境搭建 1.1下载安装环境 首先我们需要搭建一个环境ollama,下载地址如下 :Ollama 点击Download 根据自己电脑的系统选择对应版本下载即可 1.2 安装环境(window为例) 可以直接点击安装包进行安…

穷举vs暴搜vs深搜vs回溯vs剪枝系列一>黄金矿工

目录 决策树&#xff1a;代码设计代码&#xff1a; 决策树&#xff1a; 代码设计 代码&#xff1a; class Solution {boolean[][] vis;int ret,m,n;public int getMaximumGold(int[][] grid) {m grid.length;n grid[0].length;vis new boolean[m][n]; for(int i 0; i <…

基于springboot河南省旅游管理系统

基于Spring Boot的河南省旅游管理系统是一种专为河南省旅游行业设计的信息管理系统&#xff0c;旨在整合和管理河南省的旅游资源信息&#xff0c;为游客提供准确、全面的旅游攻略和服务。以下是对该系统的详细介绍&#xff1a; 一、系统背景与意义 河南省作为中国的中部省份&…

并发编程 - 线程同步(三)之原子操作Interlocked简介

上一章我们了解了3种处理多线程中共享资源安全的方法&#xff0c;今天我们将更近一步&#xff0c;学习一种针对简单线程同步场景的解决方案——Interlocked。 在此之前我们先学习一个概念——原子操作。 01、原子操作 原子操作&#xff0c;其概念源于化学领域&#xff0c;原子…

0205算法:最长连续序列、三数之和、排序链表

力扣128&#xff1a;最长连续序列 给定一个未排序的整数数组 nums &#xff0c;找出数字连续的最长序列&#xff08;不要求序列元素在原数组中连续&#xff09;的长度。 请你设计并实现时间复杂度为 O(n) 的算法解决此问题。 class Solution {public int longestConsecutive(in…