MATLAB學(xué)習(xí)筆記之傅里葉變換2
1.3.3 源代碼分析
（1）創(chuàng)建74HC595驅(qū)動
<1>創(chuàng)建74hc595.h并輸入以下代碼。
#ifndef _74HC595_H_
#define _74HC595_H_
#include "sys.h"
sbit HC595_CLK_X = P0^3 ; //X軸595時鐘
sbit HC595_DAT_X = P0^4 ; //X軸595數(shù)據(jù)
sbit HC595_RCK_X = P0^5 ; //X軸595鎖存
sbit DB16_X = P0^6 ; //增補(bǔ)X軸腳
sbit HC595_CLK_Y = P0^0 ; //Y軸595時鐘
sbit HC595_DAT_Y = P0^1 ; //Y軸595數(shù)據(jù)
sbit HC595_RCK_Y = P0^2 ; //Y軸595鎖存
void HC595_X_Dat( u16 dat ) ; //74HC595橫坐標(biāo)數(shù)據(jù)
void HC595_Y_Dat( u16 dat ) ; //74HC595縱坐標(biāo)數(shù)據(jù)
#endif
<2>創(chuàng)建74hc595.c并輸入以下代碼。
#include " 74hc595.h "
/********************************************************
Name :HC595_X_Dat
Function :74HC595橫坐標(biāo)數(shù)據(jù)
Paramater :
dat:寫入的數(shù)據(jù)
Return :None
********************************************************/
void HC595_X_Dat( u16 dat )
{
u8 i ;
for( i=0 ; i<16 ; i++ )
{
if( dat&0x8000 )
HC595_DAT_X = 1 ;
else
HC595_DAT_X = 0 ;
dat <<= 1 ;
HC595_CLK_X = 0 ;
HC595_CLK_X = 1 ;
}
HC595_RCK_X = 0 ;
HC595_RCK_X = 1 ;
}
/********************************************************
Name :HC595_Y_Dat
Function :74HC595縱坐標(biāo)數(shù)據(jù)
Paramater :
dat:寫入的數(shù)據(jù)
Return :None
********************************************************/
void HC595_Y_Dat( u16 dat )
{
u8 i ;
for( i=0 ; i<16 ; i++ )
{
if( dat&0x8000 )
HC595_DAT_Y = 1 ;
else
HC595_DAT_Y = 0 ;
dat <<= 1 ;
HC595_CLK_Y = 0 ;
HC595_CLK_Y = 1 ;
}
HC595_RCK_Y = 0 ;
HC595_RCK_Y = 1 ;
}
（2）移植FFT算法
<1>創(chuàng)建fft.h并輸入以下代碼。
#ifndef _FFT_H_
#define _FFT_H_
#include "sys.h"
/*********************************************************************************************************
函 數(shù) 列 表
*********************************************************************************************************/
void FFT_Init( void ) ; //FFT初始化
void FFT_Sample( void ) ; //FFT采樣計算(128位)
#endif
<2>創(chuàng)建fft.c并輸入以下代碼。
/*********************************************************************************************************
FFT 算 法 驅(qū) 動 代 碼
*********************************************************************************************************/
#include "fft.h"
#include "intrins.h"
#include "math.h"
#include "74hc595.h"
//采樣存儲序列表
unsigned char code BRTable[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};
char code sin_tabb[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};
char code cos_tabb[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};
xdata unsigned char result[128];
xdata unsigned char temp[128];
xdata unsigned char num[128];
unsigned char timernum;//用于分離
int xdata FftReal[128];
int xdata FftImage[128];
u16 code HC_FFT_TAB_X[] = {
0x0001, 0x0002, 0x0004, 0x0008, 0x0010, 0x0020, 0x0040, 0x0080,
0x0100, 0x0200, 0x0400, 0x0800, 0x1000, 0x2000, 0x4000, 0x8000,
} ; //FFT列選
u16 code HC_FFT_TAB_Y[] = {
0xFFFE, 0xFFFC, 0xFFF8, 0xFFF0, 0xFFE0, 0xFFC0, 0xFF80, 0xFF00,
0xFE00, 0xFC00, 0xF800, 0xF000, 0xE000, 0xC000, 0x8000, 0x0000,
} ; //FFT行選模式
/********************************************************
Name :ADC_Init
Function :ADC模塊初始化
Paramater :None
Return :None
********************************************************/
void ADC_Init()
{
P1ASF = 0x01 ;
ADC_CONTR = 0xE0 ;
}
/********************************************************
Name :ADC_Read_Byte
Function :ADC讀取10位字節(jié)
Paramater :None
Return :轉(zhuǎn)換結(jié)果的10位數(shù)據(jù)
********************************************************/
float ADC_Read_Byte()
{
ADC_CONTR = 0xE8 ;
_nop_() ;
_nop_() ;
_nop_() ;
_nop_() ;
while ( !( ADC_CONTR&0x10 ) ) ;
return ( ADC_RES*4+ADC_RESL ) ;
}
/********************************************************
Name :FFT_Init
Function :FFT初始化
Paramater :None
Return :None
********************************************************/
void FFT_Init()
{
ADC_Init();
IE |= 0x82 ;
TMOD |= 0x01 ;
TH0 = ( 65535-2048 ) / 256 ;
TL0 = ( 65535-2048 ) % 256 ;
TR0 = 1 ;
P0M0 = 0x40 ;
}
/********************************************************
Name :FFT_process
Function :下落遲滯
Paramater :None
Return :None
********************************************************/
void FFT_process()
{
u8 i ;
for( i=0; i<17; i++ )
{
if( result[ i ]>7 )+( ( FftImage[ k+bb ]*sin_tabb[ p ] )>>7 ) ;
FftImage[ k ] = FftImage[ k ] - ( ( FftReal[ k+bb ]*sin_tabb[ p ] )>>7 )+( ( FftImage[ k+bb ]*cos_tabb[ p ] )>>7 ) ;
//對稱幅值計算
FftReal[ k+bb ] = TR -( ( FftReal[ k+bb ]*cos_tabb[ p ] )>>7 )-( ( FftImage[ k+bb ]*sin_tabb[ p ] )>>7 ) ;
FftImage[ k+bb ] = TI + ( ( temp*sin_tabb[ p ] )>>7 )-( ( FftImage[ k+bb ]*cos_tabb[ p ] )>>7 ) ;
//除以2
FftReal[ k ] >>= 1 ;
FftImage[ k ] >>= 1 ;
FftReal[ k+bb ] >>= 1 ;
FftImage[ k+bb ] >>= 1 ;
}
}
}
//轉(zhuǎn)換為余弦形式
for( i=0; i<17; i++ )
{
ulReal = FftReal[ i+1 ] ; //獲取實(shí)部
ulReal *= ulReal ; //計算實(shí)部平方
ulImage = FftImage[ i+1 ] ; //獲取虛部
ulImage *= ulImage ; //計算虛部平方
result[ i ] = sqrt( ulReal+ulImage )*4 ; //計算模值
}
}
/********************************************************
Name :FFT_Disp
Function :FFT數(shù)據(jù)顯示
Paramater :None
Return :None
********************************************************/
void FFT_Disp()
{
FFT_process() ;
if( result[ timernum ]>=16 )
result[ timernum ] = 15 ;
if( timernum<16 )
{
HC595_X_Dat( HC_FFT_TAB_X[ timernum ] ) ;
DB16_X = 0 ;
}
else
{
HC595_X_Dat( 0x0000 ) ;
DB16_X = 1 ;
}
HC595_Y_Dat( HC_FFT_TAB_Y[ result[ timernum ] ] ) ;
timernum ++ ;
if( timernum==17 )
timernum = 0 ;
}
/********************************************************
Name :TIM0
Function :定時器0中斷服務(wù)函數(shù)
Paramater :None
Return :None
********************************************************/
void TIM0() interrupt 1
{
u8 i , t ;
TH0 = ( 65535-2048 ) / 256 ;
TL0 = ( 65535-2048 ) % 256 ;
t ++ ;
if( t==17 )
{
for( i=0 ; i<17 ; i++ )
temp[ i ] = result[ i ] ;
t = 0;
}
FFT_Disp() ;
}
聲明：本文內(nèi)容及配圖由入駐作者撰寫或者入駐合作網(wǎng)站授權(quán)轉(zhuǎn)載。文章觀點(diǎn)僅代表作者本人，不代表電子發(fā)燒友網(wǎng)立場。文章及其配圖僅供工程師學(xué)習(xí)之用，如有內(nèi)容侵權(quán)或者其他違規(guī)問題，請聯(lián)系本站處理。
舉報投訴
-
matlab
+關(guān)注
關(guān)注
188文章
2998瀏覽量
233464 -
線性
+關(guān)注
關(guān)注
0文章
200瀏覽量
25547 -
傅里葉變換
+關(guān)注
關(guān)注
6文章
442瀏覽量
43051
發(fā)布評論請先 登錄
相關(guān)推薦
熱點(diǎn)推薦
進(jìn)群免費(fèi)領(lǐng)FPGA學(xué)習(xí)資料！數(shù)字信號處理、傅里葉變換與FPGA開發(fā)等
進(jìn)群免費(fèi)領(lǐng)FPGA學(xué)習(xí)資料啦！小編整理了數(shù)字信號處理、傅里葉變換與FPGA開發(fā)等FPGA必看資料，需要的小伙伴可以加小助手（微信：elecfans123）或進(jìn) QQ 群：913501156 群免費(fèi)領(lǐng)
發(fā)表于 04-07 16:41
圖像頻率域分析之傅里葉變換
文章目錄傅里葉變換基礎(chǔ)傅里葉級數(shù)傅里葉積分傅里葉變換一維連續(xù)傅里葉變換一維離散傅里葉變換二維離散傅里葉變換正
發(fā)表于 05-22 07:41
學(xué)習(xí)傅里葉變換意義和方法
學(xué)習(xí)傅里葉變換需要面對大量的數(shù)學(xué)公式，數(shù)學(xué)功底較差的同學(xué)聽到傅里葉變換就頭疼。事實(shí)上，許多數(shù)學(xué)功底好的數(shù)字信號處理專業(yè)的同學(xué)也不一定理解傅里葉變換的真實(shí)含義，不能做到學(xué)以致用！事實(shí)上，
發(fā)表于 06-28 07:31
DSP變換運(yùn)算-傅里葉變換
第24章 DSP變換運(yùn)算-傅里葉變換本章節(jié)開始進(jìn)入此教程最重要的知識點(diǎn)之一傅里葉變換。關(guān)于傅里葉變換，本章主要是把傅里葉相關(guān)的基礎(chǔ)知識進(jìn)行必要的介紹，沒有這些基礎(chǔ)知識的話，后面
發(fā)表于 08-03 06:14
小波變換比傅里葉變換好在哪里_小波變換與傅里葉變換詳解
小波變換與傅里葉變換有什么區(qū)別嗎？小波變換與傅里葉變換哪個好？我們通過小波變換與傅里葉變換的詳細(xì)
發(fā)表于 01-13 11:02
?1.6w次閱讀
MATLAB學(xué)習(xí)筆記之傅里葉變換1
線性時不變系統(tǒng)** （輸入信號與輸出信號滿足線性關(guān)系，且系統(tǒng)參數(shù)不隨時間變換），無論你采用微分方程或者傳遞函數(shù)還是狀態(tài)空間，所以可以說正弦信號是系統(tǒng)的特征向量，當(dāng)然指數(shù)信號也是系統(tǒng)的特征向量，用于
傅里葉變換如何用于深度學(xué)習(xí)領(lǐng)域
機(jī)器學(xué)習(xí)和深度學(xué)習(xí)中的模型都是遵循數(shù)學(xué)函數(shù)的方式創(chuàng)建的。從數(shù)據(jù)分析到預(yù)測建模，一般情況下都會有數(shù)學(xué)原理的支撐，比如：歐幾里得距離用于檢測聚類中的聚類。 傅里葉變換是一種眾所周知的將函數(shù)從一個域轉(zhuǎn)換
深入淺出的學(xué)習(xí)傅里葉變換
學(xué)習(xí)傅里葉變換需要面對大量的數(shù)學(xué)公式，數(shù)學(xué)功底較差的同學(xué)聽到傅里葉變換就頭疼
Matlab實(shí)現(xiàn)傅里葉變換的步驟
傅里葉變換是將按時間或空間采樣的信號與按頻率采樣的相同信號進(jìn)行關(guān)聯(lián)的數(shù)學(xué)公式。
傅里葉變換基本性質(zhì) 傅里葉變換本質(zhì) 傅里葉變換的應(yīng)用
各個領(lǐng)域都有廣泛的應(yīng)用。本文將詳細(xì)介紹傅里葉變換的基本性質(zhì)、本質(zhì)和應(yīng)用。 一、傅里葉變換的基本性質(zhì) 1. 線性性：若f1(t)與f2(t)的傅里葉變換分別為F1(f)和F
傅里葉變換的定義 傅里葉變換的意義
傅里葉變換的定義 傅里葉變換的意義? 傅里葉變換，表示能將滿足一定條件的某個函數(shù)表示成三角函數(shù)（正弦和/或余弦函數(shù)）或者它們的積分的線性組合。 在不同的研究領(lǐng)域，傅里葉變換具有多種不同
傅里葉變換基本原理及在機(jī)器學(xué)習(xí)應(yīng)用
連續(xù)傅里葉變換（CFT）和離散傅里葉變換（DFT）是兩個常見的變體。CFT用于連續(xù)信號，而DFT應(yīng)用于離散信號，使其與數(shù)字?jǐn)?shù)據(jù)和機(jī)器學(xué)習(xí)任務(wù)更加相關(guān)。
發(fā)表于 03-20 11:15
?1349次閱讀
工商網(wǎng)監(jiān)
評論