出炉了:基于ATMEGA16的128点FFT变换，声音频谱分析。

fan313220 · 发表于 2009-8-9 13:17:08

mark

lichangan · 发表于 2009-8-9 15:13:42

记号

mcu_M3 · 发表于 2009-8-9 20:49:23

mark

roderickbao · 发表于 2009-8-10 15:19:50

好东西标记了

kafei523 · 发表于 2009-8-22 20:05:44

楼主：是否漏掉了
#include "misc.h"
#include "real_fft.h"
我没有找到

shusheng · 发表于 2009-8-22 20:17:36

mark

googse · 发表于 2009-8-23 23:10:59

MARK

xuquanrui · 发表于 2009-8-31 20:07:59

记号

leeber · 发表于 2009-8-31 20:33:52

mark！

yqhpaladin123 · 发表于 2009-8-31 21:28:37

mark,谢谢楼主分享！

xhero · 发表于 2009-9-7 12:43:27

也想用AVR实现FFT

zhonghua_li · 发表于 2009-9-11 22:43:29

real_fft.h 文件内容
#define SEQ_NUM 128
#define FFT_L 7 // log2(SEQ_NUM)-1

============================================================

fmuls16 = mulcd * muler; // 快速算法

wen007 · 发表于 2009-9-15 01:51:22

楼主，电路图有没有？

gxy508 · 发表于 2009-9-15 08:54:09

mark

denky · 发表于 2009-9-15 09:38:46

mark

51hubao · 发表于 2009-9-20 15:08:25

yurifeng · 发表于 2009-9-20 16:45:42

MARK

betwins · 发表于 2009-11-6 09:31:51

请教下是否可以用数字滤波器的方法来做这个音频频谱显示，这是我找到的数字带通滤波器的公式:
Bandpass

Transfer Function

H(s) = s/(s2+s/Q+1)

Coefficient
alpha = sin(ω0) * sinh( ln(2)/2 * bandwidth * ω0/sin(ω0)

a0 = 1 + alpha

a1 = -2*cos(ω0)

a2 = 1 - alpha

b0 = alpha * gainLinear

b1 = 0

b2 = -alpha * gainLinear
比如说我呀想显示 300～500 Hz 的低频信号，只需要用这个数字带通滤波器滤出来 300～500Hz 的信号，算出强度显示.
是否可以这样做呢？

vipcff · 发表于 2009-12-14 16:53:08

mark

iloveh · 发表于 2009-12-14 17:07:42

mark

gaoxvli · 发表于 2009-12-31 00:34:17

还会回来

flyingcys · 发表于 2010-1-4 15:57:14

关注！

ottife · 发表于 2010-1-4 22:06:30

MARK

yunqian09 · 发表于 2010-1-5 00:01:40

mark

scat-01 · 发表于 2010-1-5 11:42:40

mark,需要的时候一定研究下~~~

joni · 发表于 2010-1-5 20:02:33

mark

xml2028 · 发表于 2010-1-5 22:48:03

xjf20072608 · 发表于 2010-1-5 23:01:45

mark

wangff2531 · 发表于 2010-3-30 14:55:40

mark

yusufu · 发表于 2010-3-30 15:36:28

好资料，mark~~~

thinki · 发表于 2010-4-1 16:09:51

太强大了！！！

flywater · 发表于 2010-4-1 19:13:15

mark 了

smartzou · 发表于 2010-4-1 19:17:30

mark too

zzm24 · 发表于 2010-4-3 19:06:01

mark

zhonghua_li · 发表于 2010-4-29 22:54:36

看来帖子的人气不错，对本人是士气有极大的鼓舞。你们的回帖就是我前进的动力。
改天把这个东西step by step的方式讲解一下。

huike · 发表于 2010-4-29 23:42:22

楼主好人！

7c399 · 发表于 2010-4-30 00:02:14

哇~~~高手啊！用AVR居然可以实现FFT变换

lv998127 · 发表于 2010-4-30 00:17:21

FFT mark

yusufu · 发表于 2010-4-30 00:31:52

支持楼主，等你的step by step讲解~~

wcm_e · 发表于 2010-4-30 00:39:29

学习

thtlj · 发表于 2010-4-30 06:56:25

很少见的顶一次贴，谢谢。

hdxet · 发表于 2010-4-30 16:13:33

回复【137楼】zhonghua_li 蓝色天空
看来帖子的人气不错，对本人是士气有极大的鼓舞。你们的回帖就是我前进的动力。
改天把这个东西step by step的方式讲解一下。
-----------------------------------------------------------------------

值得称赞！

guangan854214 · 发表于 2010-4-30 17:21:01

关注！

longquan · 发表于 2010-4-30 17:21:10

看见算法就头疼

bg3mdo · 发表于 2010-4-30 18:22:31

我用飞思卡尔的pc9s12xs128编过fft 做电能分析呵呵

kxhui · 发表于 2010-4-30 19:11:30

这一堆数组咋样的出来

span42 · 发表于 2010-5-6 23:07:47

mark

papa123 · 发表于 2010-5-7 00:41:44

关注楼主中。。。

dmhs31 · 发表于 2010-5-12 23:29:17

期待楼主step by step讲解，学习FFT中

sqmm · 发表于 2010-5-13 15:17:57

原帖在这里,大家别到处找了.
http://elm-chan.org/works/akilcd/report_e.html
还有人设计出了下变频电路,使得这个设计能够支持高频信号的频谱分析.

里面还有PWM播放wave文件的设计,也很不错,我做了一个效果很好.

flyingcys · 发表于 2010-5-14 16:20:00

关注！

bjj9217 · 发表于 2010-5-20 09:38:26

速度怎样？

kingsabbit · 发表于 2010-5-20 11:41:24

借鉴

hiferr · 发表于 2010-5-20 12:10:12

做个记号

liumaojun_cn · 发表于 2010-5-20 12:31:37

mark

qinzutaim · 发表于 2010-5-20 17:26:13

不错的网站,不错的帖子!

lghtjpu · 发表于 2010-5-20 18:26:42

mask

lysn19 · 发表于 2010-5-21 15:42:48

回复【楼主位】zhonghua_li 蓝色天空
-----------------------------------------------------------------------

请问怎么在12864显示输入信号的波形，采样间隔怎么考虑，楼主是将输入信号序列的虚部全部置零？实部就是需转换的幅度？经过FFT变换之后频谱横轴的频率怎么确定？
刚开始研究这个，请楼主指条明路

myhonour · 发表于 2010-5-30 08:48:55

MARK

p2pahz · 发表于 2010-6-4 21:09:06

make,呵呵

60333tim · 发表于 2010-6-4 21:17:47

mark

number007cool · 发表于 2010-6-5 10:47:38

marK!

shipingjing · 发表于 2010-6-5 11:04:36

记号

getix · 发表于 2010-6-5 12:46:09

mark

yanrz · 发表于 2010-6-5 23:31:48

太牛逼了！！！
佩服！！

411412 · 发表于 2010-6-7 20:21:27

记号

fedora10 · 发表于 2010-6-8 11:02:26

mark

xdco · 发表于 2010-6-8 11:27:45

mark

lixingling · 发表于 2010-6-25 21:44:17

汉宁窗口里的数据是如何得来的

xiubing65 · 发表于 2010-6-27 22:32:03

哦？还加窗了啊？能告诉我你加窗的目的吗？

xiubing65 · 发表于 2010-6-27 22:38:05

回复【172楼】bx88w
-----------------------------------------------------------------------

回复【172楼】bx88w
mark
-----------------------------------------------------------------------

Malab里面抽样出来的，可以使用公式也可以使用wvtool。

xiubing65 · 发表于 2010-6-27 22:38:37

楼主，什么编译器啊？

googse · 发表于 2010-8-16 14:34:26

顶，LZ是个像神一样的人

esdart · 发表于 2010-8-16 15:02:11

收藏，FFT

lsany · 发表于 2010-8-28 10:39:29

非常好~~mark！

hhl871113 · 发表于 2010-10-18 20:43:56

mark

xiejob · 发表于 2010-10-22 20:11:23

LZ这个能否改成用LED灯显示的那种？

fshunj · 发表于 2010-11-1 10:25:38

mark

jielove2003 · 发表于 2010-11-1 12:39:42

mark

n104391455 · 发表于 2010-12-6 11:54:03

mark

fshunj · 发表于 2010-12-6 12:07:45

关注楼主讲解

GNMXD · 发表于 2010-12-18 22:08:29

关注楼主的 step by step

lxd_youxi · 发表于 2011-1-29 10:21:41

脚印，等楼主详解，最好是在春节前

yl604922959 · 发表于 2011-3-1 16:10:48

标记一下...

liguan1024 · 发表于 2011-3-3 20:48:19

mark

lylatyou · 发表于 2011-3-4 09:46:51

窗函数好像是布莱克曼吧！
FFT还是用DSP吧，C2000都很快，全汇编实现，1024点32位FFT，1秒钟可以做上千次
用AVR配合DSP做通讯IO，LCD显示这些，系统就很酷了。

IGO_AVR · 发表于 2011-3-4 12:22:37

MARK

aaaa4444 · 发表于 2011-3-4 13:49:26

想问一下那些数组是旋转因子吗？是怎么得来的直接用就行吗
看了圈圈的程序数组如下不一样呀不明白
#ifdef USE_TABLE
#if LENGTH==32
//计算余弦用的表，放大了1024倍（32点FFT时用）
const short int COS_TABLE[LENGTH/2]=
{
  1024,1004,946,851,724,569,392,200,0,-200,
  -392,-569,-724,-851,-946
};

//计算正弦用的表，放大了1024倍（32点FFT时用）
const short int SIN_TABLE[LENGTH/2]=
{
  0,200,392,569,724,851,946,1004,1024,1004,
  946,851,724,569,392
};
#endif

#if LENGTH==64
//计算余弦用的表，放大了1024倍（64点FFT时用）
const short int COS_TABLE[LENGTH/2]=
{
  1024,1019,1004,980,946,903,851,792,724,650,

  569,483,392,297,200,100,0,-100,-200,-297,

  -392,-483,-569,-650,-724,-792,-851,-903,-946,-980,

  -1004,-1019
};

//计算正弦用的表，放大了1024倍（64点FFT时用）
const short int SIN_TABLE[LENGTH/2]=
{
  0,100,200,297,392,483,569,650,724,792,

  851,903,946,980,1004,1019,1024,1019,1004,980,

  946,903,851,792,724,650,569,483,392,297,

  200,100
};
#endif

#if LENGTH==128
//计算余弦用的表，放大了1024倍（128点FFT时用）
const short int COS_TABLE[LENGTH/2]=
{
  1024,1023,1019,1013,1004,993,980,964,946,926,

  903,878,851,822,792,759,724,688,650,610,

  569,526,483,438,392,345,297,249,200,150,

  100,50,0,-50,-100,-150,-200,-249,-297,-345,

  -392,-438,-483,-526,-569,-610,-650,-688,-724,-759,

  -792,-822,-851,-878,-903,-926,-946,-964,-980,-993,

  -1004,-1013,-1019,-1023
};

//计算正弦用的表，放大了1024倍（128点FFT时用）
const short int SIN_TABLE[LENGTH/2]=
{
  0,50,100,150,200,249,297,345,392,438,

  483,526,569,610,650,688,724,759,792,822,

  851,878,903,926,946,964,980,993,1004,1013,

  1019,1023,1024,1023,1019,1013,1004,993,980,964,

  946,926,903,878,851,822,792,759,724,688,

  650,610,569,526,483,438,392,345,297,249,

  200,150,100,50
};
#endif

#if LENGTH==256
//计算余弦用的表，放大了1024倍（256点FFT时用）
const short int COS_TABLE[LENGTH/2]=
{
  1024,1024,1023,1021,1019,1016,1013,1009,1004,999,

  993,987,980,972,964,955,946,936,926,915,

  903,891,878,865,851,837,822,807,792,775,

  759,742,724,706,688,669,650,630,610,590,

  569,548,526,505,483,460,438,415,392,369,

  345,321,297,273,249,224,200,175,150,125,

  100,75,50,25,0,-25,-50,-75,-100,-125,

  -150,-175,-200,-224,-249,-273,-297,-321,-345,-369,

  -392,-415,-438,-460,-483,-505,-526,-548,-569,-590,

  -610,-630,-650,-669,-688,-706,-724,-742,-759,-775,

  -792,-807,-822,-837,-851,-865,-878,-891,-903,-915,

  -926,-936,-946,-955,-964,-972,-980,-987,-993,-999,

  -1004,-1009,-1013,-1016,-1019,-1021,-1023,-1024
};

//计算正弦用的表，放大了1024倍（256点FFT时用）
const short int SIN_TABLE[LENGTH/2]=
{
  0,25,50,75,100,125,150,175,200,224,

  249,273,297,321,345,369,392,415,438,460,

  483,505,526,548,569,590,610,630,650,669,

  688,706,724,742,759,775,792,807,822,837,

  851,865,878,891,903,915,926,936,946,955,

  964,972,980,987,993,999,1004,1009,1013,1016,

  1019,1021,1023,1024,1024,1024,1023,1021,1019,1016,

  1013,1009,1004,999,993,987,980,972,964,955,

  946,936,926,915,903,891,878,865,851,837,

  822,807,792,775,759,742,724,706,688,669,

  650,630,610,590,569,548,526,505,483,460,

  438,415,392,369,345,321,297,273,249,224,

  200,175,150,125,100,75,50,25
};
#endif

#if LENGTH==512
//计算余弦用的表，放大了1024倍（512点FFT时用）
const short int COS_TABLE[LENGTH/2]=
{
  1024,1024,1024,1023,1023,1022,1021,1020,1019,1018,

  1016,1015,1013,1011,1009,1007,1004,1002,999,996,

  993,990,987,983,980,976,972,968,964,960,

  955,951,946,941,936,931,926,920,915,909,

  903,897,891,885,878,872,865,858,851,844,

  837,830,822,815,807,799,792,784,775,767,

  759,750,742,733,724,715,706,697,688,678,

  669,659,650,640,630,620,610,600,590,579,

  569,558,548,537,526,516,505,494,483,472,

  460,449,438,426,415,403,392,380,369,357,

  345,333,321,309,297,285,273,261,249,237,

  224,212,200,187,175,163,150,138,125,113,

  100,88,75,63,50,38,25,13,0,-13,

  -25,-38,-50,-63,-75,-88,-100,-113,-125,-138,

  -150,-163,-175,-187,-200,-212,-224,-237,-249,-261,

  -273,-285,-297,-309,-321,-333,-345,-357,-369,-380,

  -392,-403,-415,-426,-438,-449,-460,-472,-483,-494,

  -505,-516,-526,-537,-548,-558,-569,-579,-590,-600,

  -610,-620,-630,-640,-650,-659,-669,-678,-688,-697,

  -706,-715,-724,-733,-742,-750,-759,-767,-775,-784,

  -792,-799,-807,-815,-822,-830,-837,-844,-851,-858,

  -865,-872,-878,-885,-891,-897,-903,-909,-915,-920,

  -926,-931,-936,-941,-946,-951,-955,-960,-964,-968,

  -972,-976,-980,-983,-987,-990,-993,-996,-999,-1002,

  -1004,-1007,-1009,-1011,-1013,-1015,-1016,-1018,-1019,-1020,

  -1021,-1022,-1023,-1023,-1024,-1024
};

//计算正弦用的表，放大了1024倍（512点FFT时用）
const short int SIN_TABLE[LENGTH/2]=
{
  0,13,25,38,50,63,75,88,100,113,

  125,138,150,163,175,187,200,212,224,237,

  249,261,273,285,297,309,321,333,345,357,

  369,380,392,403,415,426,438,449,460,472,

  483,494,505,516,526,537,548,558,569,579,

  590,600,610,620,630,640,650,659,669,678,

  688,697,706,715,724,733,742,750,759,767,

  775,784,792,799,807,815,822,830,837,844,

  851,858,865,872,878,885,891,897,903,909,

  915,920,926,931,936,941,946,951,955,960,

  964,968,972,976,980,983,987,990,993,996,

  999,1002,1004,1007,1009,1011,1013,1015,1016,1018,

  1019,1020,1021,1022,1023,1023,1024,1024,1024,1024,

  1024,1023,1023,1022,1021,1020,1019,1018,1016,1015,

  1013,1011,1009,1007,1004,1002,999,996,993,990,

  987,983,980,976,972,968,964,960,955,951,

  946,941,936,931,926,920,915,909,903,897,

  891,885,878,872,865,858,851,844,837,830,

  822,815,807,799,792,784,775,767,759,750,

  742,733,724,715,706,697,688,678,669,659,

  650,640,630,620,610,600,590,579,569,558,

  548,537,526,516,505,494,483,472,460,449,

  438,426,415,403,392,380,369,357,345,333,

  321,309,297,285,273,261,249,237,224,212,

  200,187,175,163,150,138,125,113,100,88,

  75,63,50,38,25,13
};
#endif

#if LENGTH==1024
//计算余弦用的表，放大了1024倍（1024点FFT时用）
const short int COS_TABLE[LENGTH/2]=
{
  1024,1024,1024,1024,1024,1024,1023,1023,1023,1022,

  1022,1022,1021,1021,1020,1020,1019,1018,1018,1017,

  1016,1016,1015,1014,1013,1012,1011,1010,1009,1008,

  1007,1006,1004,1003,1002,1000,999,998,996,995,

  993,992,990,989,987,985,983,982,980,978,

  976,974,972,970,968,966,964,962,960,958,

  955,953,951,948,946,944,941,939,936,934,

  931,928,926,923,920,917,915,912,909,906,

  903,900,897,894,891,888,885,882,878,875,

  872,868,865,862,858,855,851,848,844,841,

  837,834,830,826,822,819,815,811,807,803,

  799,796,792,788,784,779,775,771,767,763,

  759,755,750,746,742,737,733,729,724,720,

  715,711,706,702,697,692,688,683,678,674,

  669,664,659,654,650,645,640,635,630,625,

  620,615,610,605,600,595,590,584,579,574,

  569,564,558,553,548,543,537,532,526,521,

  516,510,505,499,494,488,483,477,472,466,

  460,455,449,443,438,432,426,421,415,409,

  403,398,392,386,380,374,369,363,357,351,

  345,339,333,327,321,315,309,303,297,291,

  285,279,273,267,261,255,249,243,237,230,

  224,218,212,206,200,194,187,181,175,169,

  163,156,150,144,138,132,125,119,113,107,

  100,94,88,82,75,69,63,57,50,44,

  38,31,25,19,13,6,0,-6,-13,-19,

  -25,-31,-38,-44,-50,-57,-63,-69,-75,-82,

  -88,-94,-100,-107,-113,-119,-125,-132,-138,-144,

  -150,-156,-163,-169,-175,-181,-187,-194,-200,-206,

  -212,-218,-224,-230,-237,-243,-249,-255,-261,-267,

  -273,-279,-285,-291,-297,-303,-309,-315,-321,-327,

  -333,-339,-345,-351,-357,-363,-369,-374,-380,-386,

  -392,-398,-403,-409,-415,-421,-426,-432,-438,-443,

  -449,-455,-460,-466,-472,-477,-483,-488,-494,-499,

  -505,-510,-516,-521,-526,-532,-537,-543,-548,-553,

  -558,-564,-569,-574,-579,-584,-590,-595,-600,-605,

  -610,-615,-620,-625,-630,-635,-640,-645,-650,-654,

  -659,-664,-669,-674,-678,-683,-688,-692,-697,-702,

  -706,-711,-715,-720,-724,-729,-733,-737,-742,-746,

  -750,-755,-759,-763,-767,-771,-775,-779,-784,-788,

  -792,-796,-799,-803,-807,-811,-815,-819,-822,-826,

  -830,-834,-837,-841,-844,-848,-851,-855,-858,-862,

  -865,-868,-872,-875,-878,-882,-885,-888,-891,-894,

  -897,-900,-903,-906,-909,-912,-915,-917,-920,-923,

  -926,-928,-931,-934,-936,-939,-941,-944,-946,-948,

  -951,-953,-955,-958,-960,-962,-964,-966,-968,-970,

  -972,-974,-976,-978,-980,-982,-983,-985,-987,-989,

  -990,-992,-993,-995,-996,-998,-999,-1000,-1002,-1003,

  -1004,-1006,-1007,-1008,-1009,-1010,-1011,-1012,-1013,-1014,

  -1015,-1016,-1016,-1017,-1018,-1018,-1019,-1020,-1020,-1021,

  -1021,-1022,-1022,-1022,-1023,-1023,-1023,-1024,-1024,-1024,

  -1024,-1024
};

//计算正弦用的表，放大了1024倍（1024点FFT时用）
const short int SIN_TABLE[LENGTH/2]=
{
  0,6,13,19,25,31,38,44,50,57,

  63,69,75,82,88,94,100,107,113,119,

  125,132,138,144,150,156,163,169,175,181,

  187,194,200,206,212,218,224,230,237,243,

  249,255,261,267,273,279,285,291,297,303,

  309,315,321,327,333,339,345,351,357,363,

  369,374,380,386,392,398,403,409,415,421,

  426,432,438,443,449,455,460,466,472,477,

  483,488,494,499,505,510,516,521,526,532,

  537,543,548,553,558,564,569,574,579,584,

  590,595,600,605,610,615,620,625,630,635,

  640,645,650,654,659,664,669,674,678,683,

  688,692,697,702,706,711,715,720,724,729,

  733,737,742,746,750,755,759,763,767,771,

  775,779,784,788,792,796,799,803,807,811,

  815,819,822,826,830,834,837,841,844,848,

  851,855,858,862,865,868,872,875,878,882,

  885,888,891,894,897,900,903,906,909,912,

  915,917,920,923,926,928,931,934,936,939,

  941,944,946,948,951,953,955,958,960,962,

  964,966,968,970,972,974,976,978,980,982,

  983,985,987,989,990,992,993,995,996,998,

  999,1000,1002,1003,1004,1006,1007,1008,1009,1010,

  1011,1012,1013,1014,1015,1016,1016,1017,1018,1018,

  1019,1020,1020,1021,1021,1022,1022,1022,1023,1023,

  1023,1024,1024,1024,1024,1024,1024,1024,1024,1024,

  1024,1024,1023,1023,1023,1022,1022,1022,1021,1021,

  1020,1020,1019,1018,1018,1017,1016,1016,1015,1014,

  1013,1012,1011,1010,1009,1008,1007,1006,1004,1003,

  1002,1000,999,998,996,995,993,992,990,989,

  987,985,983,982,980,978,976,974,972,970,

  968,966,964,962,960,958,955,953,951,948,

  946,944,941,939,936,934,931,928,926,923,

  920,917,915,912,909,906,903,900,897,894,

  891,888,885,882,878,875,872,868,865,862,

  858,855,851,848,844,841,837,834,830,826,

  822,819,815,811,807,803,799,796,792,788,

  784,779,775,771,767,763,759,755,750,746,

  742,737,733,729,724,720,715,711,706,702,

  697,692,688,683,678,674,669,664,659,654,

  650,645,640,635,630,625,620,615,610,605,

  600,595,590,584,579,574,569,564,558,553,

  548,543,537,532,526,521,516,510,505,499,

  494,488,483,477,472,466,460,455,449,443,

  438,432,426,421,415,409,403,398,392,386,

  380,374,369,363,357,351,345,339,333,327,

  321,315,309,303,297,291,285,279,273,267,

  261,255,249,243,237,230,224,218,212,206,

  200,194,187,181,175,169,163,156,150,144,

  138,132,125,119,113,107,100,94,88,82,

  75,69,63,57,50,44,38,31,25,19,

  13,6
};
#endif
#endif

jacky1982512 · 发表于 2011-3-4 16:01:11

LZ厉害啊。

geniusjia · 发表于 2011-3-4 17:20:31

mark

lylatyou · 发表于 2011-3-8 09:25:16

回转因子和倒胃表好像都是查表来的吧

xinzhi1986 · 发表于 2011-3-8 09:57:13

好东西，mark！~

strongtang · 发表于 2011-3-8 15:37:01

mark

wpnx · 发表于 2011-3-8 17:53:14

mark

raohao101 · 发表于 2011-4-11 22:08:55

mark

pangfen · 发表于 2011-4-11 22:18:32

mark!

xingong · 发表于 2011-4-12 14:33:48

看不懂

yfgww · 发表于 2011-4-12 16:25:15

mark

xuhaikun · 发表于 2011-4-13 00:11:33

呵呵，我都不知道是不是已经留名了。。。

zhuhanliang · 发表于 2011-4-13 08:10:18

mark...