完善资料让更多小伙伴认识你，还能领取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 个讨论
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
只有小组成员才能发言，加入小组>>
351 浏览 1 评论
546 浏览 2 评论
NA555DR VCC最低电压需要在5V供电，为什么用3.3V供电搭了个单稳态触发器也使用正常?
790 浏览 3 评论
MSP430F249TPMR出现高温存储后失效了的情况，怎么解决？
660 浏览 1 评论
对于多级放大电路板，在PCB布局中，电源摆放的位置应该注意什么？
1140 浏览 1 评论
AT32F407在USART2 DMA发送数据时，接包接到了要发送的数据，程序还是处于等待传输完成的标识判断中，为什么？
85浏览 29评论
310浏览 23评论
请问下tpa3220实际测试引脚功能和官方资料不符，哪位大佬可以帮忙解答下
262浏览 20评论
请教下关于TAS5825PEVM评估模块原理图中不太明白的地方，寻求答疑
210浏览 14评论
两个TMP117传感器一个可以正常读取温度值，一个读取的值一直是0，为什么？
69浏览 13评论
/7 
关注我们的微信
下载发烧友APP
电子发烧友观察
版权所有 © 湖南华秋数字科技有限公司
电子发烧友 （电路图） 湘公网安备 43011202000918 号 电信与信息服务业务经营许可证：合字B2-20210191
工商网监
湘ICP备2023018690号
1388
评论
淘帖
8250
17
