FFT Golang 实现

最近项目要用到快速傅立叶变换，自己写了个算法，测试了下，性能和精度还可以接受

len,time= 1048576 378.186167ms

diff=-0.00000000000225974794 I0.00000000000936106748成功: 进程退出代码 0.

百万级别，变换花了378ms

逆变换误差和0.000000000009

 l :=mvm.PowerOf2(20)
   arr :=make([]complex128,l) 
   for i,_ :=range arr{
      arr[i]=complex(rand.Float64(),rand.Float64())
   }
   now :=time.Now()
   fft:= mvm.CoreFFT(arr,false)
   Println("len,time=",l,time.Now().Sub(now))
//   Println("arr=",arr)
//   Println("fft=",fft)
   reverse:=mvm.CoreFFT(fft,true)
   //Println("reverse=",reverse)   var sumDiff complex128
   for i,v :=range reverse{
      sumDiff =sumDiff+v-arr[i]
   }
   Printf("diff=%.20f I%.20f",real(sumDiff),imag(sumDiff))



				


				
时间： 2024-08-27 19:45:57 

		


		


	

	
	

		


		
FFT Golang 实现的相关文章


								

		

			

                XJTUOJ wmq的A×B Problem FFT
			

		

		

									

				wmq的A×B Problem 发布时间: 2017年4月9日 17:06   最后更新: 2017年4月9日 17:07   时间限制: 3000ms   内存限制: 512M 描述 这是一个非常简单的问题. wmq如今开始学习乘法了!他为了训练自己的乘法计算能力,写出了n个整数,并且对每两个数a,b都求出了它们的乘积a×b.现在他想知道,在求出的n(n−1)2个乘积中,除以给定的质数m余数为k(0≤k<m)的有多少个. 输入 第一行为测试数据的组数. 对于每组测试数据,第一行为2个正整数n,			

		

	

				

		

			

                对AM信号FFT的matlab仿真
			

		

		

									

				普通调幅波AM的频谱,大信号包络检波频谱分析 u(t)=Ucm(1+macos ?t)cos ?ct ma称为调幅系数 它的频谱由载波,上下边频组成 , 包络检波中二极管截去负半周再用电容低通滤波,可以得到基带信号,那么,截去负半周后的AM信号必定包含基带信号的频谱.我们可以通过matlab来验证. %已知基带信号为1hz,载波为64hz,调制系数ma=0.3,采样频率1024hz,FFT变换区间N为2048 clear; fs=1024; f=1; %1hz基带信号 fc=64; %64hz载			

		

	

				

		

			

                golang []byte转string
			

		

		

									

				golang中,字符切片[]byte转换成string最简单的方式是 package main import ( "fmt" _ "unsafe" ) func main() { bytes := []byte("I am byte array !") str := string(bytes) bytes[0] = 'i'//注意这一行,bytes在这里修改了数据,但是str打印出来的依然没变化, fmt.Println(str) } 打印信息:			

		

	

				

		

			

                多项式FFT相关模板
			

		

		

									

				自己码了一个模板...有点辛苦...常数十分大,小心使用 #include <iostream> #include <stdio.h> #include <math.h> #include <string.h> #include <time.h> #include <stdlib.h> #include <algorithm> #include <vector> using namespace std; #de			

		

	

				

		

			

                golang实现Ringbuf
			

		

		

									

				Ring buffer算法优点:高内存使用率,在缓冲buffer内存模型中,不太容易发生内存越界.悬空指针等 bug ,出了问题也容易在内存级别分析调试.做出来的系统容易保持健壮. package main import ( "bytes" "fmt" ) type Ringbuf struct { buf         []byte start, size int } func New(size int) *Ringbuf { return &Ringb			

		

	

				

		

			

                Golang Hash MD4
			

		

		

									

				//Go标准包中只有MD5的实现 //还好,github上有MD4实现. package main import (     "golang.org/x/crypto/md4"     "encoding/hex"     "fmt" ) func get_md4(buf []byte) ([] byte) { ctx := md4.New() ctx.Write(buf) return ctx.Sum(nil) } func main() {			

		

	

				

		

			

                Java程序员的Golang入门指南(上)
			

		

		

									

				Java程序员的Golang入门指南 1.序言 Golang作为一门出身名门望族的编程语言新星,像豆瓣的Redis平台Codis.类Evernote的云笔记leanote等. 1.1 为什么要学习 如果有人说X语言比Y语言好,两方的支持者经常会激烈地争吵.如果你是某种语言老手,你就是那门语言的"传道者",下意识地会保护它.无论承认与否,你都已被困在一个隧道里,你看到的完全是局限的.<肖申克的救赎>对此有很好的注脚: [Red] These walls are funny.			

		

	

				

		

			

                golang学习笔记：golang 语法篇(二)
			

		

		

									

				在语法篇(一)中学习了go中基本的数据类型.变量.常量等组成语言的基本要素,在这一节中将会学习如何将这些元素组织起来,最终写成可以执行的代码. 在这一部分包括: go中的流程控制语句: go中函数的用法: go特殊的错误处理方式: Golang中的流程控制语句 在具体编程的时候免不了需要使用一些特殊的语句实现某些功能,比如使用循环语句来进行迭代,使用选择语句控制程序的执行方式等.这些语句在任何一门程序设计语言 中都会有支持,golang中除了支持常用的循环,条件选择语句以外,还支持跳转语句,下面			

		

	

				

		

			

                Golang关键字—— if/else
			

		

		

									

				Golang中,if/else 关键字用于条件判断,如果满足条件就做某事,否则做另一件事: if age >= 18 { fmt.Println("成年人") } else { fmt.Println("未成年") } 多重判断: if score >= 90 { fmt.Println("优秀") } else if score >= 70 { fmt.Println("良好") } else if sco