»論壇 › 嵌入式/單片機(jī)論壇 › MCU綜合技術(shù)區(qū) › PIC單片機(jī) › PIC單片機(jī)FFT音樂頻譜Proteus仿真程序

PIC單片機(jī)FFT音樂頻譜Proteus仿真程序

查看數(shù): 7836 | 評論數(shù): 8 | 收藏 6

關(guān)燈 | 提示：支持鍵盤翻頁<-左右->

帖子模式

feeling1791

發(fā)布時間: 2020-1-3 11:50

正文摘要:

因為貼子附件大小問題，wav音樂文件沒有附上，請自行添加到目錄，打開FFT.pdsprj文件，然后雙擊Audio圖標(biāo)指定音樂wav位置，開始運行即可，稍后再發(fā)下電腦里以前做的東東

回復(fù)

ID:1136495 發(fā)表于 2024-11-27 10:32

audio無法使用怎么辦呀，proteus8.17

ID:1066888 發(fā)表于 2023-3-16 21:44

FFT自己編寫的C，不錯的內(nèi)容，學(xué)習(xí)一下

ID:651494 發(fā)表于 2022-12-20 15:14

王tt 發(fā)表于 2020-7-15 11:10
采樣數(shù)據(jù)存儲列表內(nèi)數(shù)據(jù)是怎么計算的呢？

把所有數(shù)據(jù)弄到EXECL，做個圖形就明白了

ID:512651 發(fā)表于 2020-11-16 17:58

感謝樓主共享學(xué)習(xí)資料！�。�

ID:801012 發(fā)表于 2020-7-15 11:10

hsl32 發(fā)表于 2020-7-9 21:37

采樣數(shù)據(jù)存儲列表內(nèi)數(shù)據(jù)是怎么計算的呢？

ID:767423 發(fā)表于 2020-7-9 21:37

單片機(jī)真好玩呀發(fā)表于 2020-6-13 00:49
樓主能做一個用89c51實現(xiàn)音樂頻譜的嗎

#include<STC12C5A60S2.H>
typedef unsigned char uint8;
typedef unsigned int uint16;
#define channel 0x01 //設(shè)置AD通道為 P1.0
#define ADC_POWER (1<<7)
#define ADC_SPEEDHH (0x03<<5)
#define ADC_START (1<<3)
#define PADCH (1<<5)
#define ADC_FLAG (1<<4)
#define P_3 P3 //32列掃描控制
sbit sw=P1^2;//74HC154組切換控制位
//放大128倍后的sin整數(shù)表（128）
code char SIN_TAB[128]={0,6,12,18,24,30,36,42,48,54,59,65,70,75,80,85,89,94,98,102,
105,108,112,114,117,119,121,123,124,125,126,126,126,126,126,
125,124,123,121,119,117,114,112,108,105,102,98,94,89,85,80,75,
70,65,59,54,48,42,36,30,24,18,12,6,0,-6,-12,-18,-24,-30,-36,-42,
-48,-54,-59,-65,-70,-75,-80,-85,-89,-94,-98,-102,-105,-108,-112,
-114,-117,-119,-121,-123,-124,-125,-126,-126,-126,-126,-126,-125,
-124,-123,-121,-119,-117,-114,-112,-108,-105,-102,-98,-94,-89,-85,
-80,-75,-70,-65,-59,-54,-48,-42,-36,-30,-24,-18,-12,-6};
//放大128倍后的cos整數(shù)表（128）
code char COS_TAB[128]={127,126,126,125,124,123,121,119,117,114,112,108,105,102,98,94,
89,85,80,75,70,65,59,54,48,42,36,30,24,18,12,6,0,-6,-12,-18,-24,
-30,-36,-42,-48,-54,-59,-65,-70,-75,-80,-85,-89,-94,-98,-102,
-105,-108,-112,-114,-117,-119,-121,-123,-124,-125,-126,-126,-126,
-126,-126,-125,-124,-123,-121,-119,-117,-114,-112,-108,-105,-102,
-98,-94,-89,-85,-80,-75,-70,-65,-59,-54,-48,-42,-36,-30,-24,-18,
-12,-6,0,6,12,18,24,30,36,42,48,54,59,65,70,75,80,85,89,94,98,102,
105,108,112,114,117,119,121,123,124,125,126,126};
//采樣存儲序列表
code char LIST_TAB[128] = {0,64,32,96,16,80,48,112,8,72,40,104,24,88,56,120,4,68,36,100,20,84,52,116,
12,76,44,108,28,92,60,124,2,66,34,98,18,82,50,114,10,74,42,106,26,90,58,
122,6,70,38,102,22,86,54,118,14,78,46,110,30,94,62,126,1,65,33,97,17,81,49,
113,9,73,41,105,25,89,57,121,5,69,37,101,21,85,53,117,13,77,45,109,29,93,61,
125,3,67,35,99,19,83,51,115,11,75,43,107,27,91,59,123,7,71,39,103,23,87,55,
119,15,79,47,111,31,95,63,127};
uint8 LED[32]={0x80,0x00,0xC0,0x00,0xE0,0x00,0xF0,0x00,
0xF8,0x00,0xFC,0x00,0xFE,0x00,0xFF,0x00,
0xff,0x80,0xff,0xC0,0xff,0xE0,0xff,0xF0,
0xff,0xF8,0xff,0xFC,0xff,0xFE,0xff,0xFF,};
uint8 COUNT=15,COUNT1=0,ADC_Count=0,LINE=15,G=0,T;
uint8 i,j,k,b,p,anum,led=0;
int Temp_Real,Temp_Imag,temp; // 中間臨時變量
uint16 TEMP1,max=0;
int xdata Fft_Real[128];
int xdata Fft_Image[128]; // fft的虛部
uint8 xdata LED_TAB2[16]; //記錄漂浮物是否需要停頓一下
uint8 xdata LED_TAB[16]; //記錄紅色柱狀
uint8 xdata LED_TAB1[16]; //記錄漂浮點
/////////////shop60610364.taobao.com/////////
void FFT()//基2fft
{
for( i=1; i<=7; i++) /* for(1) */
{
b=1;
b <<=(i-1); //碟式運算，用于計算隔多少行計算例如第一極 1和2行計算
for( j=0; j<=b-1; j++) /* for (2) */
{
p=1;
p <<= (7-i);
p = p*j;
for( k=j; k<128; k=k+2*b) /* for (3) */
{
Temp_Real=Fft_Real[k];
Temp_Imag=Fft_Image[k];
temp=Fft_Real[k+b];
Fft_Real[k]=Fft_Real[k]+((Fft_Real[k+b]*COS_TAB[p])>>7)+((Fft_Image[k+b]*SIN_TAB[p])>>7);
Fft_Image[k]=Fft_Image[k]-((Fft_Real[k+b]*SIN_TAB[p])>>7)+((Fft_Image[k+b]*COS_TAB[p])>>7);
Fft_Real[k+b]=Temp_Real-((Fft_Real[k+b]*COS_TAB[p])>>7)-((Fft_Image[k+b]*SIN_TAB[p])>>7);
Fft_Image[k+b]=Temp_Imag+((temp*SIN_TAB[p])>>7)-((Fft_Image[k+b]*COS_TAB[p])>>7); // 移位.防止溢出. 結(jié)果已經(jīng)是本值的 1/64
Fft_Real[k] >>= 1;
Fft_Image[k] >>= 1;
Fft_Real[k+b] >>= 1;
Fft_Image[k+b] >>= 1;
}
}
}
for(j=0;j<16;j++)//16分頻
{
TEMP1=((((Fft_Real[j+1]* Fft_Real[j+1]))+((Fft_Image[j+1]*Fft_Image[j+1])))>>1);//求各頻段幅值
if(TEMP1<1)TEMP1=0;
LED_TAB[j]=TEMP1;
if(LED_TAB[j]>max)max=LED_TAB[j];
}
if(max>16)//分級量化
{
max/=16;
for(j=0;j<16;j++)LED_TAB[j]/=max;
}
}
void delay2us(void) //2us延時
{
unsigned char a,b;
for(b=1;b>0;b--)
for(a=3;a>0;a--);
}
void Init()
{
//-----------------------------------------------------------------------------------
P1ASF = 0x02; //0000,0010, 將 P1.1 置成模擬口
AUXR1 &=0xFB; //1111,1011, 令 ADRJ=0
EADC=1; //AD中斷打開
ADC_CONTR = ADC_POWER | ADC_SPEEDHH | ADC_START | channel;//1110 1001 1打開 A/D （ADC_POWER）轉(zhuǎn)換電源；11速度為90周期一次；
//0中斷標(biāo)志清零；1啟動adc(ADC_START);001AD通道打開（這里為P1.1）;
//-----------------------------------------------------------------------------------
P2M0=1;
TMOD=0X12;
TH0=0xb0;
TL0=0xb0;
TH1=0xfe;
TL1=0Xd4;
ET0=1; //定時器0 打開
TR0=0; //關(guān)閉定時器
ET1=1;
TR1=1;
PT1=0;
PT0=1;
IPH=PADCH;
IP=PADC; //中斷優(yōu)先級
EA=1; //總中斷打開
}
//----------------------------------------------------------------------------------------
void ADC_Finish() interrupt 5
{ ADC_CONTR &= !ADC_FLAG;
Fft_Real[LIST_TAB[ADC_Count]]=(int)((ADC_RES)<<1)+(ADC_RESL>>1)*2;//按LIST_TAB表里的順序，進(jìn)行存儲采樣值
if(ADC_Count<=127)ADC_Count++;
else {EADC=0;TR0=0;}
}
//---------------------------------------------------------------------------------------
void LED_Display() interrupt 3 //中斷一次顯示一列
{
TH1=0xfe;
TL1=0Xd4;
G++;
if(G==33)G=1;
switch(G) //往點陣屏填充一列的數(shù)據(jù)
{
case 1:sw=1;anum=LED_TAB[0];P_3=0x00;break;
case 2:anum=LED_TAB[0];P_3=0x01;break;
case 3:anum=LED_TAB[1];P_3=0x02;break;
case 4:anum=LED_TAB[1];P_3=0x03;break;
case 5:anum=LED_TAB[2];P_3=0x04;break;
case 6:anum=LED_TAB[2];P_3=0x05;break;
case 7:anum=LED_TAB[3];P_3=0x06;break;
case 8:anum=LED_TAB[3];P_3=0x07;break;
case 9:anum=LED_TAB[4];P_3=0x08;break;
case 10:anum=LED_TAB[4];P_3=0x09;break;
case 11:anum=LED_TAB[5];P_3=0x0a;break;
case 12:anum=LED_TAB[5];P_3=0x0b;break;
case 13:anum=LED_TAB[6];P_3=0x0c;break;
case 14:anum=LED_TAB[6];P_3=0x0d;break;
case 15:anum=LED_TAB[7];P_3=0x0e;break;
case 16:anum=LED_TAB[7];P_3=0x0f;break;
case 17:sw=0;anum=LED_TAB[8];P_3=0x00;break;
case 18:anum=LED_TAB[8];P_3=0x10;break;
case 19:anum=LED_TAB[9];P_3=0x20;break;
case 20:anum=LED_TAB[9];P_3=0x30;break;
case 21:anum=LED_TAB[10];P_3=0x40;break;
case 22:anum=LED_TAB[10];P_3=0x50;break;
case 23:anum=LED_TAB[11];P_3=0x60;break;
case 24:anum=LED_TAB[11];P_3=0x70;break;
case 25:anum=LED_TAB[12];P_3=0x80;break;
case 26:anum=LED_TAB[12];P_3=0x90;break;
case 27:anum=LED_TAB[13];P_3=0xa0;break;
case 28:anum=LED_TAB[13];P_3=0xb0;break;
case 29:anum=LED_TAB[14];P_3=0xc0;break;
case 30:anum=LED_TAB[14];P_3=0xd0;break;
case 31:anum=LED_TAB[15];P_3=0xe0;break;
case 32:anum=LED_TAB[15];P_3=0xf0;break;
}
for(led=0;led<=anum;led++) //遞增
{ P0=LED[2*led+1];
P2=LED[2*led];
delay2us();
} //遞減
for(;led>=0;led--)
{ P0=LED[2*led+1];
P2=LED[2*led];
delay2us();
}
}
//------------------------------------------------------------------------------------------------------------
void Ad_Control() interrupt 1 //控制采樣率
{
ADC_CONTR = ADC_POWER | ADC_SPEEDHH| ADC_START | channel; //開始AD采集
}
//===============================================================================================================
void main()
{
uint8 i;
Init();
while(1)
{
ADC_Count=0;
TR0=1;
EADC=1; //開啟定時器中斷0，開啟ADC
while(ADC_Count<128);
for(i=0;i<128;i++)
{
Fft_Image[i]=0;
}
FFT();//FFT運算并轉(zhuǎn)換為各頻段幅值
}
}

復(fù)制代碼

ID:745514 發(fā)表于 2020-6-13 00:49

樓主能做一個用89c51實現(xiàn)音樂頻譜的嗎

ID:679307 發(fā)表于 2020-1-3 17:01

真是高手，學(xué)習(xí)了

帳號		自動登錄	找回密碼
密碼			立即注冊