完善资料让更多小伙伴认识你，还能领取20积分哦， 立即完善>
3天内不再提示
电子发烧友论坛|
#include
#include "math.h" #include #include #define PI 3.14159265 #define f1 200 #define f2 1000 #define N 512 void gen_w_r2(float* w, int n); void bit_rev(float *x, int n); void divide(float* x, int n); void DSPF_sp_cfftr2_dit(float * x, float * w, short n); int i; float *w; float *x; float input[N],Rot[N]; void main() [ for(i=0;i<256;i++) [ input[2*i]=(sin(2*PI*f1*i/N)+sin(2*PI*f2*i/N)); input[2*i+1]=0; ] x=input; w=Rot; gen_w_r2(w, N); // Generate coefficient table bit_rev(w, N>>1); // Bit reverse coefficient table DSPF_sp_cfftr2_dit(x, w, N); // This is the radix 2 FFT benchmark from ti available input in normal order, output in coefficient table in bit-reversed order while ( 1 ); ] void gen_w_r2(float* w, int n) [ int i; float pi = 4.0*atan(1.0); float e = pi*2.0/n; for(i=0; i < ( n>>1 ); i++) [ w[2*i] = cos(i*e); w[2*i+1] = sin(i*e); ] ] //The follwoing C code is used to bit-reverse the output. void bit_rev(float *x, int n) [ int i, j, k; double rtemp, itemp; j = 0; for(i=1; i < (n-1); i++) [ k = n >> 1; while(k <= j) [ j -= k; k >>= 1; ] j += k; if(i < j) [ rtemp = x[j*2]; x[j*2] = x[i*2]; x[i*2] = rtemp; itemp = x[j*2+1]; x[j*2+1] = x[i*2+1]; x[i*2+1] = itemp; ] ] ] void divide(float* x, int n) [ int i; float inv = 1.0 / n; for(i=0; i < n; i++) [ x[2*i] = inv * x[2*i]; x[2*i+1] = inv * x[2*i+1]; ] ] 
0
|
|
相关推荐
4 个讨论
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
只有小组成员才能发言，加入小组>>
550 浏览 0 评论
1613 浏览 0 评论
2047 浏览 0 评论
为啥BQ7693003DBTR芯片在和BQ769X0盒子通讯时收不到信号？
1514 浏览 0 评论
DSP 28027F 开发板 XDS100v2调试探针诊断日志显示了 Error -150 (SC_ERR_FTDI_FAIL)如何解决
1338 浏览 0 评论
AT32F407在USART2 DMA发送数据时，接包接到了要发送的数据，程序还是处于等待传输完成的标识判断中，为什么？
1758浏览 29评论
2785浏览 23评论
请问下tpa3220实际测试引脚功能和官方资料不符，哪位大佬可以帮忙解答下
1725浏览 20评论
请教下关于TAS5825PEVM评估模块原理图中不太明白的地方，寻求答疑
1635浏览 14评论
两个TMP117传感器一个可以正常读取温度值，一个读取的值一直是0，为什么？
1646浏览 13评论
/9 
关注我们的微信
下载发烧友APP
电子发烧友观察
版权所有 © 湖南华秋数字科技有限公司
电子发烧友 （电路图） 湘公网安备 43011202000918 号 电信与信息服务业务经营许可证：合字B2-20210191
工商网监
湘ICP备2023018690号
1599
评论
淘帖
8764
19
