MATLAB學(xué)習(xí)筆記之傅里葉變換2
1.3.3 源代碼分析
（1）創(chuàng)建74HC595驅(qū)動(dòng)
<1>創(chuàng)建74hc595.h并輸入以下代碼。
#ifndef _74HC595_H_
#define _74HC595_H_
#include "sys.h"
sbit HC595_CLK_X = P0^3 ; //X軸595時(shí)鐘
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時(shí)鐘
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:寫(xiě)入的數(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:寫(xiě)入的數(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采樣計(jì)算(128位)
#endif
<2>創(chuàng)建fft.c并輸入以下代碼。
/*********************************************************************************************************
FFT 算 法 驅(qū) 動(dòng) 代 碼
*********************************************************************************************************/
#include "fft.h"
#include "intrins.h"
#include "math.h"
#include "74hc595.h"
//采樣存儲(chǔ)序列表
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 ) ;
//對(duì)稱(chēng)幅值計(jì)算
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 ; //計(jì)算實(shí)部平方
ulImage = FftImage[ i+1 ] ; //獲取虛部
ulImage *= ulImage ; //計(jì)算虛部平方
result[ i ] = sqrt( ulReal+ulImage )*4 ; //計(jì)算模值
}
}
/********************************************************
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 :定時(shí)器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)容及配圖由入駐作者撰寫(xiě)或者入駐合作網(wǎng)站授權(quán)轉(zhuǎn)載。文章觀點(diǎn)僅代表作者本人，不代表電子發(fā)燒友網(wǎng)立場(chǎng)。文章及其配圖僅供工程師學(xué)習(xí)之用，如有內(nèi)容侵權(quán)或者其他違規(guī)問(wèn)題，請(qǐng)聯(lián)系本站處理。
舉報(bào)投訴
-
matlab
+關(guān)注
關(guān)注
182文章
2963瀏覽量
230195 -
線(xiàn)性
+關(guān)注
關(guān)注
0文章
196瀏覽量
25128 -
傅里葉變換
+關(guān)注
關(guān)注
6文章
438瀏覽量
42566
發(fā)布評(píng)論請(qǐng)先 登錄
相關(guān)推薦
圖像頻率域分析之傅里葉變換
文章目錄傅里葉變換基礎(chǔ)傅里葉級(jí)數(shù)傅里葉積分傅里葉變換一維連續(xù)傅里葉變換一維離散傅里葉變換二維離散傅里葉變換正
發(fā)表于 05-22 07:41
學(xué)習(xí)傅里葉變換意義和方法
學(xué)習(xí)傅里葉變換需要面對(duì)大量的數(shù)學(xué)公式，數(shù)學(xué)功底較差的同學(xué)聽(tīng)到傅里葉變換就頭疼。事實(shí)上，許多數(shù)學(xué)功底好的數(shù)字信號(hào)處理專(zhuān)業(yè)的同學(xué)也不一定理解傅里葉變換的真實(shí)含義，不能做到學(xué)以致用！事實(shí)上，
發(fā)表于 06-28 07:31
DSP變換運(yùn)算-傅里葉變換
第24章 DSP變換運(yùn)算-傅里葉變換本章節(jié)開(kāi)始進(jìn)入此教程最重要的知識(shí)點(diǎn)之一傅里葉變換。關(guān)于傅里葉變換，本章主要是把傅里葉相關(guān)的基礎(chǔ)知識(shí)進(jìn)行必要的介紹，沒(méi)有這些基礎(chǔ)知識(shí)的話(huà)，后面
發(fā)表于 08-03 06:14
深入淺出的學(xué)習(xí)傅里葉變換
見(jiàn)過(guò)的介紹傅里葉變換的很好的文章，通俗易懂，轉(zhuǎn)發(fā)的，學(xué)習(xí)。
發(fā)表于 04-29 14:12
?10次下載
小波變換比傅里葉變換好在哪里_小波變換與傅里葉變換詳解
小波變換與傅里葉變換有什么區(qū)別嗎？小波變換與傅里葉變換哪個(gè)好？我們通過(guò)小波變換與傅里葉變換的詳細(xì)
發(fā)表于 01-13 11:02
?1.6w次閱讀
MATLAB學(xué)習(xí)筆記之傅里葉變換1
線(xiàn)性時(shí)不變系統(tǒng)** （輸入信號(hào)與輸出信號(hào)滿(mǎn)足線(xiàn)性關(guān)系，且系統(tǒng)參數(shù)不隨時(shí)間變換），無(wú)論你采用微分方程或者傳遞函數(shù)還是狀態(tài)空間，所以可以說(shuō)正弦信號(hào)是系統(tǒng)的特征向量，當(dāng)然指數(shù)信號(hào)也是系統(tǒng)的特征向量，用于
傅里葉變換如何用于深度學(xué)習(xí)領(lǐng)域
機(jī)器學(xué)習(xí)和深度學(xué)習(xí)中的模型都是遵循數(shù)學(xué)函數(shù)的方式創(chuàng)建的。從數(shù)據(jù)分析到預(yù)測(cè)建模，一般情況下都會(huì)有數(shù)學(xué)原理的支撐，比如：歐幾里得距離用于檢測(cè)聚類(lèi)中的聚類(lèi)。 傅里葉變換是一種眾所周知的將函數(shù)從一個(gè)域轉(zhuǎn)換
深入淺出的學(xué)習(xí)傅里葉變換
學(xué)習(xí)傅里葉變換需要面對(duì)大量的數(shù)學(xué)公式，數(shù)學(xué)功底較差的同學(xué)聽(tīng)到傅里葉變換就頭疼
Matlab實(shí)現(xiàn)傅里葉變換的步驟
傅里葉變換是將按時(shí)間或空間采樣的信號(hào)與按頻率采樣的相同信號(hào)進(jìn)行關(guān)聯(lián)的數(shù)學(xué)公式。
傅里葉變換基本性質(zhì) 傅里葉變換本質(zhì) 傅里葉變換的應(yīng)用
各個(gè)領(lǐng)域都有廣泛的應(yīng)用。本文將詳細(xì)介紹傅里葉變換的基本性質(zhì)、本質(zhì)和應(yīng)用。 一、傅里葉變換的基本性質(zhì) 1. 線(xiàn)性性：若f1(t)與f2(t)的傅里葉變換分別為F1(f)和F
傅里葉變換的定義 傅里葉變換的意義
傅里葉變換的定義 傅里葉變換的意義? 傅里葉變換，表示能將滿(mǎn)足一定條件的某個(gè)函數(shù)表示成三角函數(shù)（正弦和/或余弦函數(shù)）或者它們的積分的線(xiàn)性組合。 在不同的研究領(lǐng)域，傅里葉變換具有多種不同
傅里葉變換基本原理及在機(jī)器學(xué)習(xí)應(yīng)用
連續(xù)傅里葉變換（CFT）和離散傅里葉變換（DFT）是兩個(gè)常見(jiàn)的變體。CFT用于連續(xù)信號(hào)，而DFT應(yīng)用于離散信號(hào)，使其與數(shù)字?jǐn)?shù)據(jù)和機(jī)器學(xué)習(xí)任務(wù)更加相關(guān)。
發(fā)表于 03-20 11:15
?864次閱讀
工商網(wǎng)監(jiān)
評(píng)論