STM32F103C8_DMA方式2路ADC采集雙聲道音頻FFT變換顯示頻譜
這是一個完整的項目，沒有加狗
關鍵點在于多路采集，FFT庫函數應用，這個FFT函數用匯編優化，速度比較快。
當然，本系統速度瓶頸在畫點顯示。后期我采用直接寫顯存，就比較快了。
上段忙，這也是我對同道的一點小小的貢獻
51hei圖片20210121160005.jpg (130.86 KB, 下載次數: 83)
下載附件
2021-1-21 16:00 上傳
單片機源程序如下:
- #include "stm32f10x.h"
- #include "usart.h"
- #include "adc.h"
- #include "delay.h"
- #include "Lcd_Driver.h"
- #include "LCD_Config.h"
- #include "stm32_dsp.h"
- #include "math.h"
- #define NPT 256 //FFT采樣點數
- long lBUFMAG[NPT+NPT]; //存儲求模后的數據[NPT+NPT/2]
- long lBUFOUT[NPT]; //FFT輸出序列NPT=256
- long lBUFIN0[NPT]; //FFT輸入系列
- long lBUFIN1[NPT]; //FFT輸入系列
- void dsp_column0(void);
- void dsp_column1(void);
- void powerMag(long nfill);//計算頻點幅值
- void dsp_column_1(void);
- void OLED_Fill(u8 x1,u8 y1,u8 x2,u8 y2,u16 dot);
- extern __IO uint16_t ADC_ConvertedValue[NOFCHANEL];
- // 局部變量，用于保存轉換計算后的電壓值
- float ADC_ConvertedValueLocal[NOFCHANEL];
-
- int main(void)
- {
- u16 i;
- uart_init(115200);
- ADCx_Init();
- delay_init();
- Lcd_Init();//初始化硬件SPI
- Lcd_Clear(BLUE); //清屏函數
- while(1)
- {
- ADC_SoftwareStartConvCmd(ADC_x, ENABLE);
- while(!DMA_GetFlagStatus(DMA1_FLAG_TC1));
- DMA_ClearFlag(DMA1_FLAG_TC1);
- ADC_SoftwareStartConvCmd(ADC1, DISABLE);
- ADC_ConvertedValueLocal[0] =(float) ADC_ConvertedValue[0]/4096*3.3;
- ADC_ConvertedValueLocal[1] =(float) ADC_ConvertedValue[1]/4096*3.3;
- printf("CH0 = %f V \r\n",ADC_ConvertedValueLocal[0]);
- printf("CH1 = %f V \r\n",ADC_ConvertedValueLocal[1]);
- for(i=0;i<NPT;i++)
- {// 由于沒有采用外部觸發，所以使用軟件觸發ADC轉換
- ADC_SoftwareStartConvCmd(ADC_x, ENABLE);
- while(!DMA_GetFlagStatus(DMA1_FLAG_TC1));
- DMA_ClearFlag(DMA1_FLAG_TC1);
- ADC_SoftwareStartConvCmd(ADC1, DISABLE);
- lBUFIN0[i] =(float) ADC_ConvertedValue[0];
- lBUFIN1[i] =(float) ADC_ConvertedValue[1];
- }
- cr4_fft_256_stm32(lBUFOUT,lBUFIN1,NPT);//調用STM32的DSP庫作FFT變換
- powerMag(NPT);//計算頻點幅值
- dsp_column1();//顯示x根柱條。
-
- cr4_fft_256_stm32(lBUFOUT,lBUFIN0,NPT);//調用STM32的DSP庫作FFT變換
- //(FFT輸出序列,輸入序列,NPT=256)
- powerMag(NPT);//計算頻點幅值
- dsp_column0();//顯示x根柱條。
- }
- }
- //顯示各頻點的柱條
- void dsp_column0(void)
- {
- u8 i,j=80;
- for(i=1;i<161;i+=2)
- { OLED_Fill(i,0,i,128,0x0000); //填充區域的對角坐標
- OLED_Fill(i,0,i,lBUFMAG[j],RED); //填充區域的對角坐標
- j --;
- }
- }
- //顯示各頻點的柱條
- void dsp_column1(void)
- {
- u8 i,j=80;
- for(i=0;i<160;i+=2)
- { OLED_Fill(i,0,i,128,0x0000); //填充區域的對角坐標
- OLED_Fill(i,0,i,lBUFMAG[j],GREEN); //填充區域的對角坐標
- j --;
- }
- }
- //x1,y1,x2,y2 填充區域的對角坐標
- //確保x1<=x2;y1<=y2 0<=x1<=127 0<=y1<=63
- //dot:0,清空;1,填充
- void OLED_Fill(u8 x1,u8 y1,u8 x2,u8 y2,u16 dot) //可以快
- {
- u8 x,y;
- for(x=x1;x<=x2;x++)
- {
- for(y=y1;y<=y2;y++) Gui_DrawPoint(x,y,dot);//畫一個點
- }
- }
- ////////////////////////////////////////////
- //計算各頻點的模值
- void powerMag(long nfill) //計算頻點幅值
- { int32_t lX,lY;
- uint32_t i,j;
- for (i=0; i < nfill; i++) //256
- {
- lX= (lBUFOUT[i]<<16)>>16; /* sine_cosine --> cos */
- lY= (lBUFOUT[i] >> 16); /* sine_cosine --> sin */
- {
- float X= 64*((float)lX)/32768;
- float Y = 64*((float)lY)/32768;
- float Mag = sqrt(X*X+ Y*Y)/nfill; // 先平方和,再開方sqrt
- j= (long)(Mag*65536); //存儲求模后的數據
- if(j>128) j=128;//避免顯示越界
- lBUFMAG[i] =j;
- }
- }
- }
所有程序51hei提供下載:
DMA2通道ADC_FFT.7z
(216.08 KB, 下載次數: 169)
2021-1-21 16:46 上傳
點擊文件名下載附件
下載積分: 黑幣 -5
|
[復制鏈接]
QQ好友和群
QQ空間
騰訊微博
騰訊朋友
收藏
淘帖
頂
踩
