[创意DIY]

基于AVR单片机的音乐频谱

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

比较简单的~算是个半成品，还要和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-7-2 19:29:57 评论举报 61^# 绘沙

0 液晶屏多少钱？

2013-7-6 08:08:03 评论举报 62^# 卢作广

0 学习中，向楼主学习了。

2013-7-6 08:18:37 评论举报 63^# 电子海豚

0 非常不错，谢谢分享

2013-7-6 08:58:22 评论举报 64^# 王洁民

0 {:12:}{:12:}

2013-7-6 16:28:58 评论举报 65^# 瓜瓜不见了

0 支持顶顶顶好哦

2013-9-7 19:24:32 评论举报 66^# 会飞的IC

0 求原理图哦。

2013-9-9 10:28:05 评论举报 67^# C60024878

0 非常非常有才{:10:}{:10:}{:10:}{:10:}

2013-9-19 13:18:48 评论举报 68^# 小杜

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

2013-9-20 22:10:36 评论举报 69^# d743145193

0 谢谢楼主分享

2013-9-21 00:23:50 评论举报 70^# hizzxx

0 借鉴一下楼主的资料，很感谢了！

2013-10-17 14:19:41 评论举报 71^# 知识阅览者

0 不错，谢谢分享

2013-10-17 15:04:11 评论举报 72^# 曹玲

0 谢谢分享

2013-10-17 15:38:14 评论举报 73^# 刘昌明

0 好东西必须学习

2013-10-22 21:56:44 评论举报 74^# 梁骅

0 感谢分享。。。

2013-10-23 18:15:32 评论举报 75^# 458088240

0 谢谢分享，，，，

2013-10-23 18:42:52 评论举报 76^# 谢涛灼

0 谢谢分享！！！

2013-11-21 16:14:33 评论举报 77^# 何广礼

0 强人一个，还厉害

2013-11-23 13:20:13 评论举报 78^# 可苦可乐LCH

0 有意思，不过没技术……慢慢学

2013-11-23 16:56:59 评论举报 79^# 烈火菜鸟

0 谢谢分享！！！！！！

2013-11-25 18:33:38 评论举报 80^# BIN六七八

声明：本文内容及配图由入驻作者撰写或者入驻合作网站授权转载。文章观点仅代表作者本人，不代表电子发烧友网立场。文章及其配图仅供工程师学习之用，如有内容图片侵权或者其他问题，请联系本站作侵删。侵权投诉

发表新帖

精选推荐

测试板RA-Eco-RA0E1-32PIN-V1.0

335 浏览 0 评论
星闪技术已经发布一年多了，应用何去何从？

565 浏览 1 评论
基于ESP32的终极智能手表，全过程！全开源！

554 浏览 0 评论
开源项目！用ESP32做一个可爱的无用机器人

494 浏览 0 评论
为了学习内核开发，大佬手搓了一个轻量级操作系统YiYiYa OS

812 浏览 0 评论

快速回复 返回顶部 返回列表

关注微信公众号

电子发烧友网

电子发烧友论坛

社区合作: 刘勇; 联系电话：15994832713; 邮箱地址：liuyong@huaqiu.com

社区管理: elecfans短短; 微信：elecfans_666; 邮箱：users@huaqiu.com

【优惠升级】华秋PCB首单最高立减100元，再返2000元优惠券！！！

返回 DIY及创意

搜索历史

基于AVR单片机的音乐频谱

评分

相关推荐

精选推荐

站长推荐 /6