16x16音乐频谱(源代码原理图) - 51单片机 - - 手机版

论坛 > 51单片机

发帖|

16x16音乐频谱(源代码原理图)

看9640|回18|收藏

楼主 ID:104503 只看他
2016-1-26 22:34

16x16音乐频谱(源代码原理图)

部分源码预览：

代码:

#include "STC12C5A.h"
#include <intrins.h>
#include <math.h>
#include "DS1302.h"
#include "Led.h"
#define uchar unsigned char
#define uint unsigned int
#define BinToHex(n) (((n>>21)&0x80)|((n>>18)&0x40)|((n>>15)&0x20)|((n>>12)&0x10)|((n>>9)& 0x08)|((n>>6)&0x04)|((n>>3)&0x02)|((n)&0x01))
#define B(n) BinToHex(0x##n##l) // 8位2进制宏
uchar code ColScan_2[16] = {0x17,0x13,0x15,0x11,0x16,0x12,0x14,0x10,
0x27,0x23,0x25,0x21,0x26,0x22,0x24,0x20}; // 74hc138 进行列扫描
float code iw[64]=
{
1.000,0,0.9952,-0.0980,0.9808,-0.1951,0.9569,-0.2903,0.9239,-0.3827,0.8819,-0.4714,0.8315,-0.5556,
0.7730,-0.6344,0.7071,-0.7071,0.6344,-0.7730,0.5556,-0.8315,0.4714,-0.8819,0.3827,-0.9239,0.2903,-0.9569,
0.1951,-0.9808,0.0980,-0.9952,0.0,-1.0000,-0.0980,-0.9952,-0.1951,-0.9808,-0.2903,0.9569,-0.3827,-0.9239,
-0.4714,-0.8819,-0.5556,-0.8315,-0.6344,-0.7730,-0.7071,-0.7071,-0.7730,-0.6344,-0.8315,-0.5556,-0.8819,-0.4714,
-0.9239,-0.3827,-0.9569,-0.2903,-0.9808,-0.1951,-0.9952,-0.0980
};
uint data LEDBuf[16];
uint idata refreshflag[16];
uchar fft_sign = 0; // 进行fft变换标志位
struct compx
{
float real;
float imag;
};
struct compx dd[65]; // FFT数据
data struct compx temp;
void delay500ms(void)
{
unsigned char a,b,c,n;
for(c=218;c>0;c--)
for(b=79;b>0;b--)
for(a=238;a>0;a--);
for(n=4;n>0;n--);
}
void Delay(uint n)
{
uint x;
while (n--)
{
x = 5000;
while (x--);
}
}
void GPIO_Init() // GPIO口的初始化
{
P1M1 = B(00000011);
P1M0 = B(00000000);
P1 = B(00000011);
P1ASF = B(00000011);
}
void Interrupt_Init(void) // 单片机中断初始化
{
TMOD = 0x01; // 高4位控制T/C1
EA = 1; // 开总中断
TH0 = 0x00; // 16位计数寄存器T1高8位
TL0 = 0x00; // 16位计数寄存器T1低8位
ET0 = 1; // T/C1中断开
TR0 = 1; // T/C1启动
}
void ADC_Init() // 集成ADC的初始化（官方函数）
{
ADC_RES = 0; //Clear previous result
ADC_CONTR = ADC_POWER | ADC_SPEEDLL;
Delay(2); //ADC power-on and delay
}
uchar GetADCResult(uchar ch) // 进行AD转换（官方函数）
{
ADC_CONTR = ADC_POWER | ADC_SPEEDLL | ch | ADC_START;
_nop_(); //Must wait before inquiry
_nop_();
_nop_();
_nop_();
while (!(ADC_CONTR & ADC_FLAG));//Wait complete flag
ADC_CONTR &= ~ADC_FLAG; //Close ADC
return ADC_RES; //Return ADC result
}
void ee(struct compx b1,uchar b2) //复数乘法
{
temp.real=b1.real*iw[2*b2]-b1.imag*iw[2*b2+1];
temp.imag=b1.real*iw[2*b2+1]+b1.imag*iw[2*b2];
}
uint mypow(uchar nbottom,uchar ntop) //乘方函数
{
uint result=1;
uchar t;
for(t=0;t<ntop;t++)result*=nbottom;
return result;
}
void fft(struct compx *xin,uchar data N) //快速傅立叶变换
{
uchar data fftnum,i,j,k,l,m,n,disbuff,dispos,dissec;
data struct compx t;
fftnum=N; //傅立叶变换点数
for(m=1;(fftnum=fftnum/2)!=1;m++);//求得M的值
for(k=0;k<=N-1;k++) //码位倒置
{
n=k;
j=0;
for(i=m;i>0;i--) //倒置
{
j=j+((n%2)<<(i-1));
n=n/2;
}
if(k<j){t=xin[1+j];xin[1+j]=xin[1+k];xin[1+k]=t;}//交换数据
}
for(l=1;l<=m;l++) //FFT运算
{
disbuff=mypow(2,l); //求得碟间距离
dispos=disbuff/2; //求得碟形两点之间的距离
for(j=1;j<=dispos;j++)
for(i=j;i<N;i=i+disbuff) //遍历M级所有的碟形
{
dissec=i+dispos; //求得第二点的位置
ee(xin[dissec],(uint)(j-1)*(uint)N/disbuff);//复数乘法
t=temp;
xin[dissec].real=xin[i].real-t.real;
xin[dissec].imag=xin[i].imag-t.imag;
xin[i].real=xin[i].real+t.real;
xin[i].imag=xin[i].imag+t.imag;
}
}
}
void processfft(uchar num)
{
uchar pt;
uint tmp;
for(pt=1;pt<65;pt++)
{
dd[pt].imag=0; //清零虚部
}
fft(dd,64); //对当前数据进行傅立叶变换
for(pt=1;pt<33;pt++)
{
dd[pt].real=sqrt(dd[pt].real*dd[pt].real+dd[pt].imag*dd[pt].imag);//取均方根
}
switch(num)
{ // 选择显示模式
case 1: {
for(pt=0;pt<16;pt++) // style 1
{
LEDBuf[pt]=0xffff;
tmp = dd[(pt+1)*2].real;
tmp = (tmp/16)+1;
LEDBuf[pt]<<=tmp;
LEDBuf[pt]=~(LEDBuf[pt]);
}
break;
}
case 2: {
for(pt=0;pt<16;pt++) // style 2
{
if(refreshflag[pt]<(dd[pt].real/16)+1)
{
LEDBuf[pt]=0xffff;
tmp = dd[(pt+1)*2].real;
tmp = (tmp/16)+1;
refreshflag[pt] = tmp;
LEDBuf[pt]<<=tmp;
}
else
{
if(refreshflag[pt]>1)
{
refreshflag[pt]--;
}
LEDBuf[pt]=0xffff;
tmp = refreshflag[pt];
LEDBuf[pt]<<=tmp;
}
LEDBuf[pt]=~(LEDBuf[pt]);
}
break;
}
case 3: { // style 3
for(pt=0;pt<16;pt++)
{
LEDBuf[pt]=0xffff;
tmp = dd[(pt+1)*2].real;
tmp = (tmp/16)+1;
if(refreshflag[pt]<tmp)
{
refreshflag[pt] = tmp;
}
else
{
if(refreshflag[pt]>1)
{
refreshflag[pt]--;
}
tmp = refreshflag[pt];
}
LEDBuf[pt]&=(0x0001<<(tmp-1));
}
break;
}
case 4: { // style 4
for(pt=0;pt<16;pt++)
{
LEDBuf[pt]=0xffff;
tmp = dd[(pt+1)*2].real;
tmp = (tmp/16)+1;
LEDBuf[pt]<<=tmp;
LEDBuf[pt]=~(LEDBuf[pt]);
if(refreshflag[pt]<tmp)
{
refreshflag[pt] = tmp;
}
else
{
if(refreshflag[pt]>1)
{
refreshflag[pt]--;
}
tmp = refreshflag[pt];
}
LEDBuf[pt]|=(0x0001<<(tmp-1));
}
}
}
}
void main()
{
uchar i,num=1;
uint n;
GPIO_Init(); // io口初始化
ADC_Init(); // 集成adc初始化
DS1302Init(); // 初始化1302
MyLove(); // 开机动画
while(GetADCResult(0)+GetADCResult(1)==0) // 没有插入耳机时
{
DisplayTime(); // 显示时间
TimeSetting(); // 时间设定
}
for(i=0;i<16;i++) // 清空数据缓存
{
refreshflag[i] = 0x0000;
}
Interrupt_Init();
while(1)
{
if(fft_sign == 1) // 读取频率为20Hz
{
fft_sign = 0;
for(i=0;i<65;i++)
{
dd[i].real=(GetADCResult(0)+GetADCResult(1))<<2; // 读取ad结果并放大4倍；
}
processfft(num); // 傅立叶变化及处理
}
for(i=0;i<16;i++)
{ // 显示
LineInput(0x0000);
P2 = ColScan_2[i];
LineInput(LEDBuf[i]);
n = 512;
while(n--);
LineInput(0x0000);
}
if(ACT_Key == 0) // 动作检测，切换效果
{
num++;
if(num>4)
{
num = 1;
}
delay500ms();
}
}
}
void Timer0() interrupt 1 //20hz
{
TL0 = 0x3c;
TH0 = 0xb0;
fft_sign = 1;
}

附件列表

16x16音乐频谱(源代码原理图).zip (2016-1-26 22:34 上传)

196.15 KB, 下载次数: 195, 下载积分: 黑币 -5

16x16音乐频谱(源代码原理图)

沙发 ID:155591 只看他
2016-12-18 15:35

能发一下实物图，或效果展示视频吗

板凳 ID:150587 只看他
2016-12-18 18:14

这音乐的代码可以修改吗

地板 ID:155775 只看他
2016-12-18 23:47

有没有效果图呀呀

5^# ID:90996 只看他
2017-6-26 20:11

想下载没有金币，看看吧

6^# ID:169527 只看他
2017-10-17 18:39

这能用吗？那个压缩文件里面的程序和原理图对应？

7^# ID:105064 只看他
2017-10-19 18:18

分享快乐

8^# ID:240944 只看他
2017-10-19 23:18

哎，没有金币。。。

9^# ID:241152 只看他
2017-10-20 08:48

非常好，谢谢分享

10^# ID:241159 只看他
2017-10-20 09:33

楼主，是用七彩灯吗？

11^# ID:201871 只看他
2017-10-21 13:44

想看下效果图

12^# ID:124995 只看他
2018-2-28 17:39

很好很不错

13^# ID:299231 只看他
2018-4-9 22:15

效果啥样儿的，还在挣币中。。。

14^# ID:314717 只看他
2018-4-24 09:11

引用:
wwtor 发表于 2018-4-9 22:15
效果啥样儿的，还在挣币中。。。

怎么挣币啊？貌似评论都要审核

15^# ID:128101 只看他
2018-8-1 21:03

感谢分享

16^# ID:214103 只看他
2020-1-22 17:22

压缩包里的原理图怎么打不开啊，

17^# ID:694488 只看他
2020-3-1 23:02

33/5000 无法打开文件“ src stc12c5a.h”

18^# ID:694488 只看他
2020-3-1 23:02

无法打开文件“ src stc12c5a.h” 怎么解决

19^# ID:214103 只看他
2020-3-14 05:16

请问我在keil3中编译的时候出现：16X16音乐频谱(源代码+原理图)\PROJECT\MAIN.C(298): 错误 C249: 'DATA': SEGMENT TOO LARGE 这个怎么解决啊