欧美极品高清xxxxhd,国产日产欧美最新,无码AV国产东京热AV无码,国产精品人与动性XXX,国产传媒亚洲综合一区二区,四库影院永久国产精品,毛片免费免费高清视频,福利所导航夜趣136

標題: 51單片機ADC采樣傅里葉變換頻譜輸出 [打印本頁]
作者: 51hei小林    時間: 2016-9-25 10:56
標題: 51單片機ADC采樣傅里葉變換頻譜輸出
  1. //完成顯示數據的轉換工作，包括定點下落速度控制以及定點更新定格
  2. // 對fft數據進行處理，得到各個頻率短的電壓增幅值
  3. //++++++++++++++++++++++++++++++++++++++//
  4. //復數乘法
  5. void ee(struct compx b1,uchar data b2)
  6. {
  7. temp.real=b1.real*iw[2*b2]-b1.imag*iw[2*b2+1];
  8. temp.imag=b1.real*iw[2*b2+1]+b1.imag*iw[2*b2];
  9. }
  10. //++++++++++++++++++++++++++++++++++++++//  
  11. //乘方函數
  12. uint mypow(uchar data nbottom,uchar data ntop)
  13. {
  14.     uint data result=1;
  15.     uchar data t;   
  16.     for(t=0;t<ntop;t++)result*=nbottom;
  17.     return result;
  18. }
  19. //++++++++++++++++++++++++++++++++++++++//
  20. //快速傅里葉變換
  21. void fft(struct compx *xin,uchar data N)
  22. {
  23. uchar data  fftnum,i,j,k,l,m,n,disbuff,dispos,dissec;
  24. data struct compx t;
  25. fftnum=N;                         //傅里葉變換點數
  26. for(m=1;(fftnum=fftnum/2)!=1;m++);//求得M的值
  27. for(k=0;k<=N-1;k++)               //碼位倒置
  28. {
  29. n=k;
  30. j=0;
  31. for(i=m;i>0;i--)             //倒置
  32. {
  33. j=j+((n%2)<<(i-1));
  34. n=n/2;
  35. }
  36. if(k<j){t=xin[1+j];xin[1+j]=xin[1+k];xin[1+k]=t;}//交換數據
  37. }  
  38. for(l=1;l<=m;l++)                //FFT運算
  39. {
  40. disbuff=mypow(2,l);          //求碟間距離
  41. dispos=disbuff/2;            //求碟間兩點之間距離
  42. for(j=1;j<=dispos;j++)
  43. for(i=j;i<N;i=i+disbuff) //遍歷M級所有碟形
  44. {
  45. dissec=i+dispos;     //求第二點位置
  46. ee(xin[dissec],(uint)(j-1)*(uint)N/disbuff);//復數乘法
  47. t=temp;
  48. xin[dissec].real=xin[i].real-t.real;
  49. xin[dissec].imag=xin[i].imag-t.imag;
  50. xin[i].real=xin[i].real+t.real;
  51. xin[i].imag=xin[i].imag+t.imag;
  52. }
  53. }
  54. }














  69. void processfft(void)
  70. {
  71.     uchar data pt=0,tmp;

  73.     for(pt=1;pt<65;pt++)
  74.     {
  75.         dd[pt].imag=0;                //清零需部
  76.     }
  77.     fft(dd,64);                       //對當前數據進行傅里葉變換

  79.     for(pt=1;pt<65;pt++)
  80.     {                                      
  81.         dd[pt].real=sqrt(dd[pt].real*dd[pt].real+dd[pt].imag*dd[pt].imag);//取均方根
  82.     }

  84. if(Menu==1)
  85. {
  86. for(pt=2;pt<34;pt+=2)
  87.    {      
  88. for(tmp=(dd[pt].real/32)+1,LEDBuf[pt]=0xFF;tmp>=1;tmp--)//tmp>1;不保留最低位那一行長亮
  89. {
  90. LEDBuf[pt]<<=1;
  91. }
  92. LEDBuf[pt]=~(LEDBuf[pt]);

  94. }
  95. }
  96. else if(Menu==2)
  97. {

  99. //ÏÂÂä¸D
  100. for(pt=2;pt<34;pt+=2)
  101.    {      
  102. tmp=(dd[pt].real/32)+1;
  103. if(refreshflag[pt]<tmp)            //刷新數據取較大高度值 存儲顯示
  104. {
  105. for(LEDBuf[pt]=0xFF;tmp>1;tmp--)//tmp>1;不保留  最低位那一行長亮；{}一次也不執行
  106. {
  107. LEDBuf[pt]<<=1;
  108. }
  109. refreshflag[pt]=(dd[pt].real/32)+1;
  110. }
  111. else
  112. {
  113. if(refreshflag[pt]>1)refreshflag[pt]--;                //定點下落速度控制  改變值可以改變下落速度
  114. for(LEDBuf[pt]=0xFF,tmp=refreshflag[pt];tmp>1;tmp--)   //tmp>1;不保留  最低位 那一行長亮
  115. {
  116. LEDBuf[pt]<<=1;
  117. }
  118. }
  119. }
  120. }
  121. else if(Menu==0)
  122. {

  124. //ÏÂÂä¸DÆμÆ×
  125. for(pt=2;pt<34;pt+=2)
  126.    {      
  127. tmp=(dd[pt].real/32)+1;
  128. if(refreshflag[pt]<tmp)          //刷新數據取較大高度值 存儲顯示
  129. {
  130. for(LEDBuf[pt]=0xFF;tmp>=1;tmp--)//tmp>1;不保留  最低位那一行長亮；{}一次也不執行
  131. {
  132. LEDBuf[pt]<<=1;
  133. }
  134. refreshflag[pt]=(dd[pt].real/32)+1;
  135. }
  136. else
  137. {
  138. if(refreshflag[pt]>1)refreshflag[pt]--;                        
  139. for(LEDBuf[pt]=0xFF,tmp=refreshflag[pt];tmp>=1;tmp--)
  140. {
  141. LEDBuf[pt]<<=1;
  142. }
  143. }
  144. LEDBuf[pt]=~(LEDBuf[pt]);
  145. }
  146. }
  147. }  
復制代碼
作者: yxf0915    時間: 2017-10-12 09:18
請問這里 “Menu” = 0，1，2 是什么意思？
作者: 51黑電子名字    時間: 2018-3-12 15:50
yxf0915 發表于 2017-10-12 09:18
**** 作者被禁止或刪除 內容自動屏蔽 ****

我也搞不懂這個，同問
作者: 磊雷雷    時間: 2018-7-29 12:13
可用嗎
作者: zouzichun    時間: 2019-9-2 15:14
這代碼亂寫的，沒啥用



歡迎光臨 (http://www.raoushi.com/bbs/) Powered by Discuz! X3.1