[创意DIY]

基于AVR单片机的音乐频谱

2013-5-1 16:53:20 18427 单片机

比较简单的~算是个半成品，还要和TDA2030功放机结合的~
源程序：
/*
* MusicSpectrum.c
*
* Created: 2013/4/30 19:34:05
*  Author: 詹磊
*/
#include
#include
#include "../Include/main.h"
#include "../Include/ADC.h"
#include "../Include/TFT.h"
#include "../Include/FFT.h"
compx Data[16];
void DrawingColumnar(uchar x,uchar y,uchar Wide,uchar High,int Data,int FullData,uint ColumnarColor)
{
uint i,tempData1,tempData2;
uchar tempData;
tempData=High-(uint)High*Data/FullData;
tempData1=tempData*Wide;
tempData2=((uint)High)*Wide;
y=159-y;
TFT_RamAddSet(x,y-High,x+Wide-1,y);
for (i=0;i {
TFT_WriteByteDAT(BackgroundColor>>8); //高八位
TFT_WriteByteDAT(BackgroundColor); //低八位
}
for (i=tempData1;i {
TFT_WriteByteDAT(ColumnarColor>>8); //高八位
TFT_WriteByteDAT(ColumnarColor); //低八位
}
}
int main(void)
{
uint temp;
uchar i;
ADCInit();
TFT_Init();
FFT_Init();
TFT_putstr(0,0,"--ZHAN LEI--n--MusicSpectrum--n--2013/04/30--",DataRed);
while(1)
{
if (time==0)
{
for (i=0;i<16;i++)
{

///*
Data[i].real=255+sin(PI*2*i/FFT_N)*20
+sin(PI*2*i/FFT_N*2)*40
+sin(PI*2*i/FFT_N*3)*80
+sin(PI*2*i/FFT_N*4)*120
+sin(PI*2*i/FFT_N*5)*110
+sin(PI*2*i/FFT_N*6)*100
+sin(PI*2*i/FFT_N*7)*90
+sin(PI*2*i/FFT_N*7.9)*20;
//*/
Data[i].imag=0;
/*
Data[i].real=(ADC_Buffer[i]-448)<<3;
*/
}
FFT(Data);
for (i=0;i<16;i++)
{
DrawingColumnar(8*i,0,2,63,Data[i].real,1023,DataGreen);
}
ADCSRA |=(1< }
}
}

/*
* FFT.h
*
* Created: 2013/5/1 0:15:37
*  Author: zhanddkk
*/
#ifndef FFT_H_
#define FFT_H_
#include "../Include/main.h"

#define PI 3.1415926535897932384626433832795028841971  //定义圆周率值
#define FFT_N 16 //频谱分数
#define FFT_m 4 //最大级数(2^FFT_m=FFT_N)

typedef struct COMPX
{
float real;
float imag;
}compx;

extern void FFT_Init();
extern void FFT(compx *xin);
#endif /* FFT_H_ */

/*
* FFT.c
*
* Created: 2013/5/1 0:16:05
*  Author: zhanddkk
*/
#include
#include "../Include/FFT.h"
float sin_tab[FFT_N],cos_tab[FFT_N];
//-------------------------------------------//
//函数名：FFT初始化函数
//入口：void
//出口：void
//功能：建立sin()、cos()表
//-------------------------------------------//
void FFT_Init()
{
uchar i;
for (i=0;i {
sin_tab[i]=sin(PI*2*i/FFT_N);
cos_tab[i]=cos(PI*2*i/FFT_N);
}
}
//-------------------------------------------//
//函数名：FFT运算函数
//入口：compx *xin:需要计算的FFT数据表结构体指针
//出口：void
//功能：FFT运算，最后得到峰值功率频谱
//-------------------------------------------//
void FFT(compx *xin)
{
uchar i,j=0,k;
uchar L,b,p;
uchar Nv2,Nm1;
float TR,TI;
//float tempData;
float tempData_rc,tempData_is,tempData_rs,tempData_ic;
float t;
Nv2=FFT_N/2;                //变址运算，即把自然顺序变成倒位序，采用雷德算法
Nm1=FFT_N-1;
//--倒序运算--//
for(i=0;i {
if(i {
t=xin[j].real;
xin[j].real=xin[i].real;
xin[i].real=t;
}
k=Nv2;
while(k<=j)
{
j=j-k;
k=k/2;
}
j=j+k;
}
//--基2的FFT蝶形运算--//
for (L=1;L<=FFT_m;L++) /* Loop_1 L是L级蝶形运算 (2^FFT_m=FFT_N)*/
{
b=1<<(L-1); /* b=2^(L-1) b是进行蝶形运算的两个数据的距离 */
for (j=0;j {
p=1<<(FFT_m-L); /* p=2^(FFT_m-L) 旋转因子计算 */
p*=j;
for (k=j;k {
TR=xin[k].real;
TI=xin[k].imag;

tempData_rc=xin[k+b].real*cos_tab[p];
tempData_is=xin[k+b].imag*sin_tab[p];
tempData_rs=xin[k+b].real*sin_tab[p];
tempData_ic=xin[k+b].imag*cos_tab[p];
/*
tempData=xin[k+b].real;
xin[k].real=xin[k].real+xin[k+b].real*cos_tab[p]+xin[k+b].imag*sin_tab[p];
xin[k].imag=xin[k].imag-xin[k+b].real*sin_tab[p]+xin[k+b].imag*cos_tab[p];
xin[k+b].real=TR-xin[k+b].real*cos_tab[p]-xin[k+b].imag*sin_tab[p];
xin[k+b].imag=TI+tempData*sin_tab[p]-xin[k+b].imag*cos_tab[p];
*/
///*
xin[k].real=xin[k].real+tempData_rc+tempData_is;
xin[k].imag=xin[k].imag-tempData_rs+tempData_ic;
xin[k+b].real=TR-tempData_rc-tempData_is;
xin[k+b].imag=TI+tempData_rs-tempData_ic;
//*/
}
}
}
//--功率峰值计算--//
for (i=0;i {
xin[i].real=sqrt(xin[i].real*xin[i].real+xin[i].imag*xin[i].imag);
}
}

IMAG0510.jpg (260.75 KB, 下载次数: 6)

IMAG0512.jpg (265.32 KB, 下载次数: 6)

评分

参与人数 1	威望 +10	+10	积分 +20	收起理由
电子眼	+ 10	+ 10	+ 20	原创技术先锋,加分鼓励!

查看全部评分

2013-5-1 16:53:20　　评论 淘帖0 举报

相关推荐

97 个讨论

0

2013-5-10 21:07:46 评论举报 41^# 韦珍

0 厉害滴童鞋……学习啦，谢谢分享。

2013-5-10 21:53:20 评论举报 42^# 日心尧亘

0 支持一下，感谢分享

2013-5-10 22:15:02 评论举报 43^# chenengood

0 进入看看。

2013-5-11 00:04:27 评论举报 44^# shcnashn

0 厉害

2013-5-11 10:51:17 评论举报 45^# 罗雨明

0 {:1:}{:1:}

2013-5-12 22:14:45 评论举报 46^# 网博

0

2013-5-13 21:49:09 评论举报 47^# 桑旺兴

0 不错啊！支持支持！

2013-5-13 22:45:06 评论举报 48^# ljs_0437

0 ，强大楼主，果断收藏学习

2013-5-14 16:05:26 评论举报 49^# jinniuxing

0 看看，好好学习

2013-5-15 13:10:52 评论举报 50^# hqyuan

0 惭愧还没写过FFT 。。学习呢 c语言也有复数是吧？

2013-5-16 19:45:55 评论举报 51^# Coo_o

0 顶一个

2013-5-18 20:21:43 评论举报 52^# 可口可

0 赞一个。。。。。

2013-5-22 08:46:04 评论举报 53^# 杨润靖

0 好，非常好。

2013-5-22 13:46:09 评论举报 54^# 郭义斌

0 顶一个！！！！

2013-5-31 17:56:10 评论举报 55^# 天一处

0 人才就是这么低调，但是做的东西却一点也不低调

2013-6-1 16:15:12 评论举报 56^# 张文强

0 支持。。。。。。。。。

2013-6-1 20:42:01 评论举报 57^# 天凡幽星

0 {:1:}{:1:}{:1:}{:1:}{:9:}

2013-6-2 11:58:16 评论举报 58^# 李步国

0 关于转换后彩屏驱动显示的函数，可以加点注释吗？学的不好看不懂。谢谢了！

2013-6-9 18:20:07 评论举报 59^# uoudio