找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 3849|回复: 10
打印 上一主题 下一主题
收起左侧

单片机+LCD12864液晶驱动的频谱显示

  [复制链接]
跳转到指定楼层
楼主
12864液晶驱动的频谱显示
电路原理图如下:


单片机源程序如下:
  1. #include<math.h>
  2. #include"stc12c5a60s2.h"
  3. #include"LCD12864.h"
  4. #define uint unsigned int
  5. #define uchar unsigned char
  6. struct compx                                                                           //定义复数结构体
  7. {
  8.         float real;
  9.         float imag;
  10. };
  11. xdata struct compx s[ 64 ];                                                 //FFT数据缓存放在XDATA空间
  12. struct compx EE(struct compx,struct compx);                 //复数乘法函数的声明
  13. void FFT(struct compx xin[],int N);                                  //FFT函数的声明
  14. struct compx EE(struct compx a1,struct compx b2)         //复数乘法
  15. {
  16.         struct compx b3;
  17.         b3.real=a1.real*b2.real-a1.imag*b2.imag;
  18.         b3.imag=a1.real*b2.imag+a1.imag*b2.real;
  19.         return(b3);
  20. }
  21. /*FFT函数*/
  22. void FFT(struct compx xin[],int N)                                   
  23. {
  24.         int f,m,nv2,nm1,i,k,j=1,l;
  25.         struct compx v,w,t;
  26.         nv2=N/2;
  27.         f=N;
  28.         for(m=1;(f=f/2)!=1;m++){;}
  29.         nm1=N-1;
  30.         for(i=0;i<nm1;i++)                                           //倒序操作
  31.         {
  32.                 if(i<j)
  33.                 {
  34.                         t=xin[j];
  35.                         xin[j]=xin[i];
  36.                         xin[i]=t;
  37.                 }
  38.                 k=nv2;                                                       //k为倒序中相应位置的权值
  39.                 while(k<j)
  40.                 {
  41.                         j=j-k;
  42.                         k=k/2;
  43.                 }
  44.                 j=j+k;
  45.         }
  46.         {
  47.                 int le,lei,ip;
  48.                   float pi;
  49.                   for(l=1;l<=m;l++)
  50.                    {
  51.                 le=pow(2,l);                                                 //乘方
  52.                     lei=le/2;
  53.                     pi=3.14159265;
  54.                     v.real=1.0;
  55.                 v.imag=0.0;  
  56.                    w.real=cos(pi/lei);                                           //旋转因子
  57.                     w.imag=-sin(pi/lei);
  58.                     
  59.                 for(j=1;j<=lei;j++)                                           //控制蝶形运算的级数
  60.                      {
  61.                         for(i=j-1;i<N;i=i+le)                                         //控制每级蝶形运算的次数
  62.                               {
  63.                                     ip=i+lei;
  64.                                        t=EE(xin[ ip ],v);
  65.                                        xin[ ip ].real=xin[ i ].real-t.real;   //蝶形计算
  66.                                        xin[ ip ].imag=xin[ i ].imag-t.imag;
  67.                                        xin[ i ].real=xin[ i ].real+t.real;
  68.                                        xin[ i ].imag=xin[ i ].imag+t.imag;
  69.                               }
  70.                               v=EE(v,w);   
  71.                       }     
  72.                    }
  73.           }
  74. }
  75. void showbar()
  76. {                                
  77.     xdata unsigned char dis_data[31][16];  //用于显示的数据
  78.     unsigned char i,j,high;
  79.         unsigned char dis_rdata[16];
  80.         for(i=0;i<16;i++)          //读取FFT转换数据
  81.         {
  82.                 float t0=0;
  83.                 t0=sqrt(pow((s[i  ].real+s[i+1].real),2)+pow((s[i  ].imag+s[i+1].imag),2))/2;
  84.                 dis_rdata[i]=(unsigned char)t0;
  85.         }
  86. /**************************************************************/
  87. /*****************写入下屏,i:0~16,表示16列**************/
  88. /***************************high:表示列的高度*********/
  89.       for (i=0;i<=32;i++) //在显示下屏显存必须清0,后面有操作是或运算
  90.           {
  91.         for (j=0;j<16;j++)
  92.                 {
  93.                          dis_data[i][j]=0;
  94.                 }         
  95.           }
  96.      for(i=0;i<16;i++)
  97.          {   
  98.                  high=dis_rdata[i];
  99.                  if(high>62)high=62;        //12864的高度只有64,判断高于64时,为62***/
  100.                  if(high>31)          //下屏和下屏的高度各处为32,大于32时,多于的在上屏显示
  101.                          {         
  102.                                  for(j=0;j<31;j++)//高度大于32时,下屏柱体全充满,多于的在上屏显示
  103.                         {
  104.                           dis_data[j][i]|=0xfe;  //  或的方法是数据结合12864的显示特性进行的算法
  105.                         }                                                  //0xfe使每条柱体之间有一条线的间隔
  106.                          }        
  107.                          else  //下屏和下屏的高度各处为32,小于32个,在下屏显示
  108.                          {
  109.                             for(j=31-high;j<31;j++)
  110.                         {
  111.                             dis_data[j][i]|=0xfe;  //  或的方法是数据结合12864的显示特性进行的算法
  112.                         }
  113.                          }                                                
  114.     }                  
  115.     for(i=0;i<=31;i++)//写入下屏图形数据时采用逐行扫描的方式
  116.     {  
  117.         write_com(0x80|i); //写入x和y的地址
  118.         write_com(0x88|0);                           
  119.         for(j=0;j<16;j++)write_data(dis_data[i][j]); //写入数据
  120.     }
  121. /**************************************************************/
  122. /*****************写入上屏,i:0~16,表示16列**************/
  123. /***************************high:表示列的高度*********/
  124.         for (i=0;i<=32;i++) //后面在上屏显示,显存必须清0,后面有操作或运算
  125.         {
  126.         for (j=0;j<16;j++)
  127.                 {
  128.                          dis_data[i][j]=0;
  129.                 }         
  130.         }
  131.          for(i=0;i<16;i++)
  132.          {   
  133.                  high=dis_rdata[i];
  134.                  if(high>62)high=62;
  135.               if(high>31)
  136.                  {         
  137.                         for(j=64-high;j<=31;j++)
  138.                 {
  139.                   dis_data[j][i]|=0xfe;  //  或的方法是数据结合12864的显示特性进行的算法
  140.                 }
  141.                  }
  142.                  else
  143.                  {
  144.                          for(j=64-high;j<=31;j++)
  145.                 {
  146.                   dis_data[j][i]|=0x00;  //  或的方法是数据结合12864的显示特性进行的算法
  147.                 }
  148.                  }                                                         
  149.     }                  
  150.     for(i=0;i<=31;i++)//写入图形数据时采用逐行扫描的方式
  151.     {  
  152.         write_com(0x80|i);
  153.         write_com(0x90|0);                           
  154.         for(j=0;j<16;j++)write_data(dis_data[i][j]);
  155.     }
  156.         
  157. }
  158. /*主函数*/
  159. void main()                                                                     
  160. {
  161.         int N=64,i;                                                             //变量初始化,64点FFT运算
  162.         float offset;
  163.     lcd_init();          //12864
  164.         P1ASF=0x01;                                                                      //P10口做AD 使用
  165.         P1M0 = 0x01;                                                //0000,0001用于A/D转换的P1.x口,先设为开漏
  166.         P1M1 = 0x01;                                                //0000,0001 P1.0先设为开漏。断开内部上拉电阻
  167.         ADC_CONTR=0xC8;                                                                 //40.96K采样率
  168.         while(!(ADC_CONTR&0x10));
  169.         offset=((float)ADC_RES*4+(float)(ADC_RESL%0x04)); //AD结果高8位左移2位,低2位不变,然后相加
  170.         while(1)
  171.         {                if(P3==(P3&0xFE))IAP_CONTR=0x60;

  172.                         for(i=0;i<N;i++)                                      //采集音频信号
  173.                           {
  174.                                 ADC_CONTR=0xC8;                                         //40.96K采样率
  175.                                 while(!(ADC_CONTR&0x10));
  176.                                    s[i].real=((float)ADC_RES*4+(float)(ADC_RESL%0x04)-offset)/4;//((((int)ADC_DATA-128)/2))*4;
  177.                                    s[i].imag=0;
  178.                           }  
  179.                           FFT(s,N);                                       //调用FFT函数进行变换
  180.                           showbar();                                                        //显示频谱                                                               
  181.         }
  182. }
复制代码

Keil5程序51hei下载:
12864.rar (168.18 KB, 下载次数: 140)

评分

参与人数 1黑币 +50 收起 理由
admin + 50 共享资料的黑币奖励!

查看全部评分

分享到:  QQ好友和群QQ好友和群 QQ空间QQ空间 腾讯微博腾讯微博 腾讯朋友腾讯朋友
收藏收藏4 分享淘帖 顶1 踩
回复

使用道具 举报

沙发
ID:13396 发表于 2020-2-26 19:45 | 只看该作者
谢谢分享。原理图是JPEG格式的就更好了。
回复

使用道具 举报

板凳
ID:688460 发表于 2020-5-2 11:30 | 只看该作者
好好学习!谢谢楼主!
回复

使用道具 举报

地板
ID:990486 发表于 2022-5-2 17:18 | 只看该作者
很不错!
回复

使用道具 举报

5#
ID:990486 发表于 2022-5-2 17:18 | 只看该作者
谢谢
好好学习!谢谢楼主!
回复

使用道具 举报

6#
ID:990486 发表于 2022-5-2 17:19 | 只看该作者
稍作修改,实验了一下,不错
回复

使用道具 举报

7#
ID:1043639 发表于 2022-12-11 21:43 | 只看该作者
全51hei就你注释的非常详细和清楚,点赞
回复

使用道具 举报

8#
ID:1066888 发表于 2023-3-16 21:52 | 只看该作者
proteus仿真能通过吗?有做过这个仿真的朋友吗?
回复

使用道具 举报

9#
ID:1066888 发表于 2023-4-10 10:46 | 只看该作者
谢谢楼主!下载收藏学习一下了
回复

使用道具 举报

10#
ID:963310 发表于 2023-5-13 10:05 | 只看该作者
有做过的吗?
回复

使用道具 举报

11#
ID:38792 发表于 2023-5-22 09:33 | 只看该作者
看看效果如何。
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

手机版|小黑屋|51黑电子论坛 |51黑电子论坛6群 QQ 管理员QQ:125739409;技术交流QQ群281945664

Powered by 单片机教程网

快速回复 返回顶部 返回列表