1、2407FFT 128 点 FFT的sin 和 cos值存储表TMS320LF2407上实现快速傅里叶变换(FFT)快速傅立叶变换(FFT)的原理FFT的程序代码(1)主程序#include f2407_c.h#include math.h#define N 32 / FFT变换的点数 extern void fft(void); extern void resave(void); interrupt void phantom(void); void sysinit(void);extern int input2*N; / 输入数据的存储数组 int indatiN=0; / -/ 128 点
2、 FFT所需的数据/ 采样函数: x=1/4+1/4cos(3*2*pi*f*t)+1/4cos(6*2*pi*f*t)+1/4cos(9*2*pi*f*t);/ f=50Hz/ -/* int indatrN=16394,15871,14425,12398, 10276 ,8584 , 7767 , 8088 ,9557 , 11913 , 14660 ,17155 , 18724 , 18802 , 17044 , 13411 , 8197 , 1995 , -4389 ,-10071, -14231 , -16255 ,-15844, -13057 ,-8309 , -2296, 412
3、5 , 10079, 14819, 17843 , 18969, 18342 , 16394 ,13739 , 11055 , 8950 , 7848, 7921 , 9070 , 10961, 13110 , 14992 , 16159 , 16334 ,5479 , 13792 , 11675 , 9640 , 8197 , 7741, 8457, 10264 , 12815 , 15554 , 17812 ,18939 , 18429 , 16034 , 11825 , 6203 , -156, -6405, -11662 , -15165 , -16394 , -15165 ,-116
4、62 , -6405 , -156 , 6203 , 11825 , 16034, 18429, 18939 , 17812 , 15554 , 12815 ,10264 , 8457 , 7741 , 8197, 9640 , 11675, 13792, 15479 , 16334 , 16159, 14992 ,13110 , 10961 , 9070 ,7921 , 7848 , 8950 , 11055 , 13739 , 16393 , 18342 , 18969 ,17843 , 14819 , 10079 ,4125 , -2296 , -8309 , -13057 , -158
5、44 , -16255 , -14231 , -10071 ,-4389 , 1995 , 8197 , 13411 , 17044 , 18802 , 18724 , 17155 , 14660 , 11913 , 9557 ,8088 , 7767, 8584 , 10276 , 12398 , 14425 , 15871 ,;*/ -/ 32点FFT所需的数据 / 采样函数: x=1/4+1/4cos(3*2*pi*f*t)+1/4cos(6*2*pi*f*t)+1/4cos(9*2*pi*f*t);/ f=50Hz ;pi=;/ -/* int indatrN=16384, 10270
6、, 9551 , 18713 , 8192 , -14222 , -8304 , 14810 , 16384 ,7843 , 13102 ,15469, 8192 , 12807 , 18418 , -0156 ,-16384 , -0156 , 18418 , 12807 , 8192 , 15469,13102 , 7843 , 16383 , 14810 , -8304 , -14222 , 8192 , 18713 , 9551 , 10270 ,;*/int indatrN=0x07ff, 0x07ff, 0x07ff , 0x07ff , 0x07ff , 0x07ff , 0x0
7、7ff , 0x07ff , 0x07ff ,0x07ff , 0x07ff ,0x07ff, 0x07ff , 0x07ff , 0x07ff , 0x07ff ,0x0F801 , 0x0F801 , 0x0F801 , 0x0F801 , 0x0F801 , 0x0F801,0x0F801 , 0x0F801, 0x0F801 , 0x0F801 , 0x0F801 , 0x0F801 , 0x0F801 , 0x0F801 , 0x0F801 , 0x0F801 ,; / -/ 64点 FFT所需的数据/ 采样函数: x=1/4+1/4cos(3*2*pi*f*t)+1/4cos(6*
8、2*pi*f*t)+1/4cos(9*2*pi*f*t);/ f=50Hz/ -/* int indatrN=16384 , 14416, 10270, 7762, 9551 , 14651 , 18713, 17034, 8192 , -4387 , -14222,-15834 , -8304 , 4123 , 14810 , 18957 , 16384 , 11049 , 7843 , 9064 , 13102 , 16149 ,15469 , 11668 , 8192 , 8452 , 12807 , 17801 , 18418 , 11818, -0156 , -11655 , -16
9、384 , -11655 , -0156 , 11818 , 18418 , 17801, 12807 , 8452 , 8192, 11668 , 15469 , 16149 ,13102 , 9064 , 7843 , 11049 , 16383 , 18957 , 14810 , 4123 , -8304 , -15834 , -14222 ,-4387 , 8192 , 17034 , 18713 , 14651 , 9551 , 7762 , 10270 , 14416 , ;*/ -/ 128 点 FFT的sin 和 cos值存储表/ -/* const int sintabN=0
10、x07fff,0x0,0x07fd9,0x0f9b9,0x07f62,0x0f375,0x07e9d,0x0ed38, 0x07d8a,0x0e708,0x07c2a,0x0e0e7,0x07a7d,0x0dad8,0x07885,0x0d4e1, 0x07642,0x0cf05,0x073b6,0x0c946,0x070e3,0x0c3aa,0x06dca,0x0be32, 0x06A6C,0x0B8E4,0x066CE,0x0B3C1,0x062F1,0x0AECD,0x05ED6,0x0AA0C, 0x05A81,0x0A57F,0x055F4,0x0A12A,0x05133,0x09D
11、0F,0x04C3F,0x09932, 0x0471C,0x09594,0x041CD,0x09237,0x03C56,0x08F1F,0x036B9,0x08C4B, 0x030FB,0x089C0,0x02B1E,0x0877D,0x02527,0x08584,0x01F19,0x083D7, 0x018F8,0x08277,0x012C7,0x08164,0x00C8B,0x0809F,0x00647,0x08029, 0x00000,0x08001,0x0F9B9,0x08029,0x0F375,0x0809F,0x0ED39,0x08164, 0x0E708,0x08277,0x0E
12、0E7,0x083D7,0x0DAD9,0x08584,0x0D4E2,0x0877D, 0x0CF05,0x089C0,0x0C947,0x08C4B,0x0C3AA,0x08F1F,0x0BE33,0x09237, 0x0B8E4,0x09594,0x0B3C1,0x09932,0x0AECD,0x09D0F,0x0AA0C,0x0A12A, 0x0A57F,0x0A57F,0x0A12A,0x0AA0C,0x09D0F,0x0AECD,0x09932,0x0B3C1, 0x09594,0x0B8E4,0x09237,0x0BE33,0x08F1F,0x0C3AA,0x08C4B,0x0C
13、947, 0x089C0,0x0CF05,0x0877D,0x0D4E2,0x08584,0x0DAD9,0x083D7,0x0E0E7, 0x08277,0x0E708,0x08164,0x0ED39,0x0809F,0x0F375,0x08029,0x0F9B9;*/ -/ 64点 FFT的sin 和 cos值存储表/ - /* const int sintabN= 0x7FFF,0x0000,0x7F61 ,0xF375, 0x7D89 , 0xE708 , 0x7A7C , 0xDAD9,0x7640, 0xCF05 ,0x70E1 ,0xC3AA ,0x6A6C , 0xB8E4 ,
14、 0x62F1 , 0xAECD ,0x5A81, 0xA57F, 0x5133, 0x9D0F, 0x471C , 0x9594, 0x3C56, 0x8F1F,0x30FB, 0x89C0 ,0x2527, 0x8584 ,0x18F8 , 0x8277 , 0x0C8B , 0x809F,0x0000 ,0x8001, 0xF375, 0x809F, 0xE708 , 0x8277 , 0xDAD9 , 0x8584,0xCF05 ,0x89C0 ,0xC3AA, 0x8F1F, 0xB8E4 , 0x9594 , 0xAECD , 0x9D0F,0xA57F, 0xA57F, 0x9D
15、0F, 0xAECD, 0x9594 , 0xB8E4 , 0x8F1F , 0xC3AA,0x89C0, 0xCF05, 0x8584 ,0xDAD9 ,0x8277 , 0xE708 , 0x809F , 0xF375, ;*/ - / 32 点 FFT的sin 和 cos值存储表/ -const int sintabN= 0x7FFF,0x0000,0x7D89,0xE708,0x7640,0xCF05,0x6A6C,0xB8E4, 0x5A81,0xA57F,0x471C,0x9594,0x30FB,0x89C0,0x18F8,0x8277, 0x0000,0x8001,0xE708,
16、0x8277,0xCF05,0x89C0,0xB8E4,0x9594, 0xA57F,0xA57F,0x9594,0xB8E4,0x89C0,0xCF05,0x8277,0xE708,;extern int table128;extern int nom; / 当nom=1时, FFT 需要归一化处理main() int i; double x=0,y; nom=1; / 需要归一化处理 sysinit(); for(i=0;iRi ,ARP:AR2, / AR2:ID, AR3:Y, AR5=M, AR6:input LACC *,1 SACL *+ / ID=ID*2,ARP:AR2,AR2:IW,/ AR3:Y, AR5=M, AR6:input LACC *,15 SACH * / IW=IW/2 LACC *-,15,AR3 SACH *+,AR2 / C2=IW/2/ 堆栈分布情况:ADDRESS/AR6/AR7/AR0/N/M/ID/IW/C2/Y/ ARP:AR2, AR2
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1