FFT算法在單片機中的使用&&LCD12864驅動
本次創新基金我是要做一個簡易的頻譜儀，核心就是要進行一個FFT運算。大家知道，如果采用DSP芯片效果那是相當好的。但由于項目資金以及時間不夠等情況，我采用的是ATMEL公司的AVR單片機，這款單片機的FLASH存儲和內存比51單片機犀利得多。
由于采用的是12864液晶，也就是一個橫128點豎64點的一個點陣，因而采用128點FFT運算已然夠了，因為即使得到再多的數據也無法在液晶上可視化顯示出來。本文是基于128點FFT運算。
程序如下：
#include
#include
#include
#define N 128
#define PI 3.141592653589
#define uchar unsigned char
#define uint unsigned int
typedef struct
{
int real;
int img;
}complex;
void initw(); //初始化旋轉因子
void bitReverse(); //比特反轉
void FFT();
complex x[N];
uchar vis[N];
void delayms(uint ms)
{
uint i,j;
for(i=0;i
{
for(j=0;j<3;j++);
}
}
void FFT()
{
int i,j,k,t,P,B,m;
complex up,down,product;
for (i=0;i<7;i++)
{
B=1<
for (j=0;j
{
t=1<<(6-i);
P=t*j;
for (k=j;k
{
complex product;
product.real=x[k+B].real*cos(2*PI*P/N)+x[k+B].img*sin(2*PI*P/N);
product.img=x[k+B].real*(-1)* sin(2*PI*P/N)+x[k+B].img*cos(2*PI*P/N);
x[k+B].real=x[k].real-product.real;
x[k+B].img=x[k].img-product.img;
x[k].real=x[k].real+product.real;
x[k].img=x[k].img+product.img;
}
}
}
}
void initw() //初始化旋轉因子
{
int i;
for (i=0;i
vis[i]=0;
}
void bitReverse() //比特反轉
{
int i,j=0;
int k=0;
int q=0;
complex tmp3;
for (i=0;i
{
int tmp=i,tmp2=0,j;
for(j=0;j<7;j++)
tmp2+=((tmp>>j)&1)*(1<<(6-j));
if(vis[i]==0)
{
tmp3=x[i];
x[i]=x[tmp2];
x[tmp2]=tmp3;
vis[i]=1;
vis[tmp2]=1;
}
}
}
void main()
{
uchar ii,y;
float tmp;
for (ii=0;ii<20;ii++)
{
x[ii].real=3;
x[ii].img=0;
}
for (ii=20;ii<128;ii++)
{
x[ii].real=0;
x[ii].img=0;
}
initw();
bitReverse();
FFT();
while(1);
}
上圖是8點FFT運算，按照上圖的流程所示，FFT運算主要有兩步，一步是比特反轉，就是右邊不是按照0、1、2、3……這樣順序進行計算的，而左邊是的，兩邊的關系就是進行一個比特反轉。可以看到右邊0對應二進制為000，左邊對應二進制為000，右邊1二進制001，左邊4對應二進制100，依次下去，可以清楚看到，對于8位FFT運算，對應二進制有三位，而左右兩邊的關系恰巧是按照中間位進行了個反轉。
FFT運算第二步就是乘以旋轉因子，注意的是這里是復數運算，虛部和實部都要加入運算。乘以旋轉因子后對進行加減運算得到新的值，依次下去得到最終解。
由于單片機內存的限制，因而對于傳統的FFT算法，我進行了些改進，原則就是盡量地少使用變量，一個變量可以重復的使用是最理想的了，大家可以在程序中看出。個人意見這是能節省變量最少的了，如果有好的方法，希望可以告訴我下，我的郵箱是albertvictordu@139.com，謝謝!
下面是12864液晶驅動程序的寫法：
LCD12864液晶，即像素為128*64的顯示液晶。它的每一行橫向一共有128個可顯示點，每一列縱向有64個，這些“點”其實也都是一個個發光二極管。它可以在一個16*16的點陣區域上顯示一個中文，也可以在一個8*16的點陣區域顯示一個非中文字符，一般稱為半寬字體。即一個中文字所占顯示面積是一個非中文字符的兩倍。
關于驅動函數的書寫，是液晶顯示的基礎，整個液晶驅動主要有四個函數組成：
1、寫命令函數;
2、寫數據函數;
3、讀狀態函數;
4、讀數據函數;
這四個函數并不是必須全部寫的，具體要看你實現的功能，如果只是單純的顯示漢字和字符，寫命令、寫數據、讀狀態這三個函數就夠了，如過你還需要進行一些繪圖的操作，那讀數據函數也必須書寫。
另外關于讀狀態函數，其實也就是用于判忙操作，原則上每次對控制器進行讀寫操作之前，都必須進行讀寫檢測，由于單片機的操作速度慢于液晶控制器的反應速度，因此可不進行讀寫檢測，或者只進行簡短的延時即可。因此，讀狀態函數也可以不寫，只用簡短的延時函數替換即可。
單片機用于控制LCD的管腳主要為RS、RW和E管腳，分別的功能是RS為0時，對應單片機訪問的是命令寄存器，為1時對應數據寄存器;RW為1時，對應單片機操作為讀操作，為0時對應單片機為寫操作;E是使能信號。
讀操作如下圖所示
寫操作如下圖所示
在12864液晶中，開發商將一些基本指令已經寫入到命令寄存器中，我們調用該指令就可以完成相應的功能。
LCD初始化
初始化操作如下：
1. 芯片上電;
2. 延時40ms以上;
3. 復位操作：RST出現一個上升沿(RST=1;RST=0;RST=1;);
4. 功能設定;
5. 延時100us以上;
6. 再次進行功能設定;
7. 延時37us;
8. 顯示開關控制;
9. 延時100us以上;
10. 清除顯示;
11. 延時10ms以上;
12. 進入點設置;
13. 初始化結束;
LCD液晶屏初始化過程如圖所示為：
打點函數
打點函數是創建GUI的基礎，打點函數的書寫分為以下幾個步驟：
1. 進入擴展模式
2. 寫入打點地址
3. 讀取該地址的數據
4. 修改該地址的數據
5. 將修改后的數據輸入LCD中
6. 進入普通模式
GDRAM地址分布情況，需要注意的是橫縱坐標的起始地址都是0x80，還有上下半屏的橫坐標是不一樣的，下半屏的橫坐標要加上0x08，而縱坐標跟對應的上半屏的縱坐標是一樣的。GDRAM地址分布圖，如圖所示。
下面的函數是12864與FFT算法的一個結合，里面設置了一個門函數，12864上顯示的結果則是一個sinc函數，證明結果是正確的。
#include
#include
#include
#define N 128
#define PI 3.141592653589
#define uchar unsigned char
#define uint unsigned int
#define RS (1<<4)
#define RW (1<<5)
#define EN (1<<6)
//
typedef struct
{
int real;
int img;
}complex;
void initw(); //初始化旋轉因子
void bitReverse(); //比特反轉
void FFT();
complex x[N];
uchar vis[N];
void delayms(uint ms)
{
uint i,j;
for(i=0;i
{
for(j=0;j<3;j++);
}
}
//此處定義字符串
//寫數據
void WriteDataLCM(unsigned char WDLCM) //寫數據函數
{
// ReadStatusLCM(); //檢測忙
delayms(1);
PORTA|=RS; //RS=1
delayms(1);
PORTA&=~RW; //RW=0
delayms(1);
PORTA|=EN; //EN=1
delayms(1);
PORTB=WDLCM; //輸出數據
delayms(1);
PORTA&=~EN; //EN=0
delayms(1);
}
//寫指令
void WriteCommandLCM(unsigned char WCLCM) //寫命令函數
{
// ReadStatusLCM(); //根據需要檢測忙
delayms(1);
PORTA&=~RS; //RS=0
delayms(1);
PORTA&=~RW; //RW=0
delayms(1);
PORTA|=EN; //EN=1
delayms(1);
PORTB=WCLCM; //輸出指令
delayms(1);
PORTA&=~EN; //EN=0
delayms(1);
}
//讀狀態:檢測忙
void ReadStatusLCM() //讀狀態函數
{
uchar temp;
uchar flag = 1;
while(flag==1)
{
PORTB=0xff;
delayms(1);
DDRB=0x00; //端口B改為輸入
delayms(1);
PORTA&=~RS; //RS=0
delayms(1);
PORTA|=RW; //RW=1
delayms(1);
PORTA|=EN; //EN=1
delayms(10);
temp = PINB; //讀端口B
delayms(10);
DDRB=0xff; //端口B改為
delayms(10);
PORTA&=~EN; //EN=0
delayms(1);
if(temp>>7==0)
flag = 0;
}
}
uchar read_data() //讀數據函數
{
uchar lcd_data;
PORTB=0xff;
DDRB=0x00;
PORTA|=RW;
PORTA|=RS;
delayms(10);
PORTA|=EN;
delayms(10);
lcd_data=PINB;
delayms(10);
PORTA&=~EN;
DDRB=0xff;
return(lcd_data) ;
}
void point(uchar x,uchar y) //打點函數，最重要的函數，GUI的基礎
{
uchar x_Dyte,x_byte;
uchar y_Dyte,y_byte;
uchar GDRAM_hbit,GDRAM_lbit;
WriteCommandLCM(0x36);
x_Dyte=x/16;
x_byte=x&0x0f;
y_Dyte=y/32;
y_byte=y&0x1f;
WriteCommandLCM(0x80+y_byte);
WriteCommandLCM(0x80+x_Dyte+8*y_Dyte);
read_data();
GDRAM_hbit=read_data();
GDRAM_lbit=read_data();
delayms(10);
WriteCommandLCM(0x80+y_byte);
WriteCommandLCM(0x80+x_Dyte+8*y_Dyte);
delayms(10);
if(x_byte<8)
{
WriteDataLCM(GDRAM_hbit|(0x01<<(7-x_byte)));
WriteDataLCM(GDRAM_lbit);
}
else
{
WriteDataLCM(GDRAM_hbit);
WriteDataLCM(GDRAM_lbit|(0x01<<(15-x_byte)));
}
WriteCommandLCM(0x30);
}
//LCM初始化
void LCMInit(void)
{
WriteCommandLCM(0x30); //三次顯示模式設置，不檢測忙信號
delayms(10);
WriteCommandLCM(0x30);
delayms(10);
WriteCommandLCM(0x30);
delayms(10);
WriteCommandLCM(0x30); //顯示模式設置,開始要求每次檢測忙信號
WriteCommandLCM(0x08); //關閉顯示
WriteCommandLCM(0x01); //顯示清屏
WriteCommandLCM(0x06); //顯示光標移動設置
WriteCommandLCM(0x0C); //顯示開及光標設置
}
void clear(uchar dat) //清屏函數
{
uchar i,j,k;
uchar addr=0x80;
for(i=0;i<2;i++)
{
for(j=0;j<32;j++)
{
for(k=0;k<8;k++)
{
WriteCommandLCM(0x36);
WriteCommandLCM(0x80+j);
WriteCommandLCM(addr+k);
WriteDataLCM(dat);
WriteDataLCM(dat);
}
}
addr=0x88;
}
WriteCommandLCM(0x36);
WriteCommandLCM(0x30);
}
void heng(uchar a)
{
uchar i;
for(i=0;i<127;i++)
point(i,a);
}
void su(uchar a)
{
uchar i;
for(i=0;i<63;i++)
point(a,i);
}
void FFT()
{
int i,j,k,t,P,B,m;
complex up,down,product;
for (i=0;i<7;i++)
{
B=1<
for (j=0;j
{
t=1<<(6-i);
P=t*j;
for (k=j;k
{
complex product;
product.real=x[k+B].real*cos(2*PI*P/N)+x[k+B].img*sin(2*PI*P/N);
product.img=x[k+B].real*(-1)* sin(2*PI*P/N)+x[k+B].img*cos(2*PI*P/N);
x[k+B].real=x[k].real-product.real;
x[k+B].img=x[k].img-product.img;
x[k].real=x[k].real+product.real;
x[k].img=x[k].img+product.img;
}
}
}
}
void initw() //初始化旋轉因子
{
int i;
for (i=0;i
vis[i]=0;
}
void bitReverse() //比特反轉
{
int i,j=0;
int k=0;
int q=0;
complex tmp3;
for (i=0;i
{
int tmp=i,tmp2=0,j;
for(j=0;j<7;j++)
tmp2+=((tmp>>j)&1)*(1<<(6-j));
if(vis[i]==0)
{
tmp3=x[i];
x[i]=x[tmp2];
x[tmp2]=tmp3;
vis[i]=1;
vis[tmp2]=1;
}
}
}
void xian(uchar x,uchar y)
{
uchar i;
for(i=63;i>=y;i--)
point(x,i);
}
//主函數
void main(void)
{
uchar ii,y;
float tmp;
//端口初始化
DDRA=0xff;
PORTA=0xff;
DDRB=0xff;
PORTB=0xff;
DDRD=0xff;
PORTD=0x00;
delayms(20);
delayms(20);
LCMInit(); //LCM初始化 //液晶初始化
delayms(100);
clear(0x00);
heng(0);
heng(63);
su(0);
su(127);
for (ii=0;ii<20;ii++)
{
x[ii].real=3;
x[ii].img=0;
}
for (ii=20;ii<128;ii++)
{
x[ii].real=0;
x[ii].img=0;
}
initw();
bitReverse();
FFT();
for(ii=64;ii<128;ii++)
{
tmp=sqrt((x[ii].real*x[ii].real)+(x[ii].img*x[ii].img));
y= 63-(int)tmp;
point(ii-64,y);
xian(ii-64,y);
}
for(ii=0;ii<64;ii++)
{
tmp=sqrt((x[ii].real*x[ii].real)+(x[ii].img*x[ii].img));
y= 63-(int)tmp;
point(ii+64,y);
xian(ii+64,y);
}
while(1);
}
得到的圖片：
MATLAB仿真圖形：
編輯：admin 最后修改時間：2018-05-19
