后使用快捷导航没有帐号？
查看: 8570|回复: 17
请教关于fft的问题。
在线时间1 小时
TA的帖子TA的资源
一粒金砂（初级）, 积分 0, 距离下一级还需 5 积分
一粒金砂（初级）, 积分 0, 距离下一级还需 5 积分
本帖最后由 dontium 于
13:39 编辑 各位老师好，我想用2812做一个fft的程序，但不知如何下手，我想向各位请教一下，应该如何利用TI公司自带的fft库函数来实现，我现在把需要计算的小数值存放在外部ram中，不知道应该怎么处理，而且我知道2812是一个定点芯片，对于小数又应当如何处理呢，谢谢各位老师了。
在线时间0 小时
TA的帖子TA的资源
一粒金砂（初级）, 积分 2, 距离下一级还需 3 积分
一粒金砂（初级）, 积分 2, 距离下一级还需 3 积分
本帖最后由 dontium 于
13:39 编辑 对于小数可以用定标的方法
在线时间1 小时
TA的帖子TA的资源
一粒金砂（初级）, 积分 2, 距离下一级还需 3 积分
一粒金砂（初级）, 积分 2, 距离下一级还需 3 积分
那我通过的AD转换得到的16位数，如果不转换成小数，能不能直接进行FFT变换呢，如果可以的话，应该怎么实现呢，如果不行，那定标识怎么实现的，我不懂啊。
在线时间1 小时
TA的帖子TA的资源
一粒金砂（初级）, 积分 2, 距离下一级还需 3 积分
一粒金砂（初级）, 积分 2, 距离下一级还需 3 积分
那我再进行fft变换时，必须得进行定标处理吗，我如果直接转换成小数存在外部ram里，不知道这样做可以吗？
在线时间0 小时
TA的帖子TA的资源
一粒金砂（初级）, 积分 2, 距离下一级还需 3 积分
一粒金砂（初级）, 积分 2, 距离下一级还需 3 积分
看你的算法决定是否使用小数，数据放在哪都无所谓
在线时间0 小时
TA的帖子TA的资源
一粒金砂（初级）, 积分 2, 距离下一级还需 3 积分
一粒金砂（初级）, 积分 2, 距离下一级还需 3 积分
我是这么想的，我想把AD采到的数据进行放大，然后进行FFT变换，再将转换后的结果同时缩小相应的倍数，然后显示在液晶上，这样就不用进行格式的转化，不知道这样行不？
在线时间0 小时
TA的帖子TA的资源
一粒金砂（初级）, 积分 2, 距离下一级还需 3 积分
一粒金砂（初级）, 积分 2, 距离下一级还需 3 积分
很多书上都说用C语言来编写fft算法，速度很慢，是不是我不能使用fft库函数，而要使用汇编编写，然后用C语言调用才可以呢，我不会用汇编来编写2812的程序，那我应该怎么办呢？
在线时间1 小时
TA的帖子TA的资源
一粒金砂（初级）, 积分 2, 距离下一级还需 3 积分
一粒金砂（初级）, 积分 2, 距离下一级还需 3 积分
各位大侠，我还有一个问题要请教，当我用fft库函数作变换时，我想在ccs的图形显示区观察所得到的频谱图，然而利用fft库函数得到的函数值很多，我不知道应该观察那个函数，以及怎样将该函数数组以图形的形式显示出来，希望各位大侠指点一下，还有上面我所问的问题，都指点一下好吗，谢谢大家了。
在线时间0 小时
TA的帖子TA的资源
一粒金砂（初级）, 积分 2, 距离下一级还需 3 积分
一粒金砂（初级）, 积分 2, 距离下一级还需 3 积分
参考刘和平的2407dsp C语言开发应用 上面有FFT的算法介绍，不过是用汇编的
在线时间1 小时
TA的帖子TA的资源
一粒金砂（初级）, 积分 2, 距离下一级还需 3 积分
一粒金砂（初级）, 积分 2, 距离下一级还需 3 积分
可我现在想用TI公司的fft库函数直接完成FFT变换，就不需要再编写函数了，因为时间比较急，所以我只能出此下策了，可是我不会用这个库函数，希望大家指点一下好吗？
在线时间0 小时
TA的帖子TA的资源
一粒金砂（初级）, 积分 2, 距离下一级还需 3 积分
一粒金砂（初级）, 积分 2, 距离下一级还需 3 积分
最近我都要被fft折磨得快要崩溃了，光看资料实在是看不懂，不知道说的是什么意思，还有就是fft库函数里的那些函数是怎么使用的，我一个也看不懂，因此，我十分恳请EEWORLD的老师们能不吝赐教，不胜感激，谢谢各位了。
在线时间0 小时
TA的帖子TA的资源
一粒金砂（初级）, 积分 2, 距离下一级还需 3 积分
一粒金砂（初级）, 积分 2, 距离下一级还需 3 积分
网上找的资料，据说调通了。（我没试过）
由于2812时定点运算，sin值得大小不能太小！否则计算结果为零！2812的fft的错误在于fft.init1(),fft.init2()这两个函数没有这函数库里面定义，而是该用fft.init(&fft),而且由于是定点运算，输入的sin值不能太小，要不然计算结果为零！
#include \&DSP281x_Device.h\& // DSP281x Headerfile Include File
#include \&DSP281x_Examples.h\& // DSP281x Examples Include File
#include \&fft.h\&
#define BUF_SIZE 128
#pragma DATA_SECTION(ipcb, \&FFTipcb\"
#pragma DATA_SECTION(mag, \&FFTmag\"
// Prototype statements for functions found within this file.
// Global variables used in this example:
Uint16 LoopC
Uint16 ConversionC
Uint16 fftflag=0;//FFT转换启动标志:0不启动；1:启动
long Voltage1[BUF_SIZE];
long Voltage2[BUF_SIZE];
long ipcb[BUF_SIZE];
RFFT32 fft=RFFT32_128P_DEFAULTS;
//CFFT32 fft=CFFT32_128P_DEFAULTS;
long mag[BUF_SIZE/2+1];
const long win[BUF_SIZE/2]=HANNING128;
double j=0.12345;
// Step 1. Initialize System Control:
// PLL, WatchDog, enable Peripheral Clocks
// This example function is found in the DSP281x_SysCtrl.c file.
InitSysCtrl();
// For this example, set HSPCLK to SYSCLKOUT / 6 (25Mhz assuming 150Mhz SYSCLKOUT)
SysCtrlRegs.HISPCP.all = 0x3; // HSPCLK = SYSCLKOUT/6
// Step 2. Initialize GPI
// This example function is found in the DSP281x_Gpio.c file and
// illustrates how to set the GPIO to it\'s default state.
// InitGpio(); // Skipped for this example
// Step 3. Clear all interrupts and initialize PIE vector table:
// Disable CPU interrupts
// Initialize the PIE control registers to their default state.
// The default state is all PIE interrupts disabled and flags
// are cleared.
// This function is found in the DSP281x_PieCtrl.c file.
InitPieCtrl();
// Disable CPU interrupts and clear all CPU interrupt flags:
IER = 0x0000;
IFR = 0x0000;
// Initialize the PIE vector table with pointers to the shell Interrupt
// Service Routines (ISR).
// This will populate the entire table, even if the interrupt
// is not used in this example. This is useful for debug purposes.
// The shell ISR routines are found in DSP281x_DefaultIsr.c.
// This function is found in DSP281x_PieVect.c.
InitPieVectTable();
// Interrupts that are used in this example are re-mapped to
// ISR functions found within this file.
// EALLOW; // This is needed to write to EALLOW protected register
// PieVectTable.ADCINT = &adc_
// EDIS; // This is needed to disable write to EALLOW protected registers
// Step 4. Initialize all the Device Peripherals:
// This function is found in DSP281x_InitPeripherals.c
// InitPeripherals(); // Not required for this example
InitAdc(); // For this example, init the ADC
// Step 5. User specific code, enable interrupts:
PieCtrlRegs.PIEIER1.bit.INTx6 = 1;
IER |= M_INT1; // Enable CPU Interrupt 1
EINT; // Enable Global interrupt INTM
ERTM; // Enable Global realtime interrupt DBGM
LoopCount = 0;
ConversionCount = 0;
for(i=0;i {
// Voltage1=(long)1000*sin((2*3.14159*i)/25);
ipcb=(long)*sin((2*3.14159*i)/12);
// RFFT32_brev(Voltage1,ipcb,BUF_SIZE);
//initialize FFT module
fft.ipcbptr=
fft.magptr=//存放幅度
// fft.magptr=
// fft.winptr=(long *)//加窗
fft.init(&fft);
// fft.win(&fft);
RFFT32_brev(ipcb,ipcb,BUF_SIZE);
// RFFT32_brev(ipcb,ipcb,16);
fft.calc(&fft);
fft.split(&fft);
fft.mag(&fft);
for(i=0;i {
mag=sqrt(mag);
//###########################################################################
// FILE: F2812_EzDSP_RAM_lnk.cmd
// TITLE: Linker Command File For F2812 eZdsp examples that run out of RAM
// This linker file assumes the user is booting up in Jump to H0 mode
//###########################################################################
// Ver | dd mmm yyyy | Who | Description of changes
// =====|=============|======|===============================================
// 1.00| 11 Sep 2003 | L.H. | Changes since previous version (v.58 Alpha)
// | | | Added BEGIN section to the start of H0
// | | | Removed .bss, .const and .sysmem
// | | | These are for a small memory model. All examples
// | | | use the large model.
// | | | Added .esysmem section
// | | | Changed ramfuncs section to load and run from RAM
// | | | (previously this was type DSECT)
// | | | Moved peripheral register files to DSP28_Headers_BIOS.cmd
// | | | and DSP28_Headers_nonBIOS.cmd
// | | | Added CSM_RSVD memory section in FLASHA - this region
// | | | should be programmed with all 0x0000 when using the CSM
// -----|-------------|------|-----------------------------------------------
//###########################################################################
/* ======================================================
// For Code Composer Studio V2.2 and later
// ---------------------------------------
// In addition to this memory linker command file,
// add the header linker command file directly to the project.
// The header linker command file is required to link the
// peripheral structures to the proper locations within
// the memory map.
// The header linker files are found in \\DSP281x_Headers\\cmd
// For BIOS applications add: DSP281x_Headers_nonBIOS.cmd
// For nonBIOS applications add: DSP281x_Headers_nonBIOS.cmd
========================================================= */
/* ======================================================
// For Code Composer Studio prior to V2.2
// --------------------------------------
// 1) Use one of the following -l statements to include the
// header linker command file in the project. The header linker
// file is required to link the peripheral structures to the proper
// locations within the memory map */
/* Uncomment this line to include file only for non-BIOS applications */
/* -l DSP281x_Headers_nonBIOS.cmd */
/* Uncomment this line to include file only for BIOS applications */
/* -l DSP281x_Headers_BIOS.cmd */
/* 2) In your project add the path to \\DSP281x_headers\\cmd to the
library search path under project-&build options, linker tab,
library search path (-i).
/*========================================================= */
/* For this example, H0 is split between PAGE 0 and PAGE 1 */
/* BEGIN is used for the \&boot to HO\& bootloader mode */
/* RESET is loaded with the reset vector only if */
/* the boot is from XINTF Zone 7. Otherwise reset vector */
/* is fetched from boot ROM. See .reset section below */
RAMM0 : origin = 0x000000, length = 0x000400
BEGIN : origin = 0x3F8000, length = 0x000002
PRAMH0 : origin = 0x3F8002, length = 0x000FFE
RESET : origin = 0x3FFFC0, length = 0x000002
/* For this example, H0 is split between PAGE 0 and PAGE 1 */
RAMM1 : origin = 0x000400, length = 0x000400
DRAMH0 : origin = 0x3f9000, length = 0x001000
/* Setup for \&boot to H0\& mode:
The codestart section (found in DSP28_CodeStartBranch.asm)
re-directs execution to the start of user code.
Place this section at the start of H0 */
codestart : & BEGIN, PAGE = 0
ramfuncs : & PRAMH0 PAGE = 0
.text : & PRAMH0, PAGE = 0
.cinit : & PRAMH0, PAGE = 0
.pinit : & PRAMH0, PAGE = 0
.switch : & RAMM0, PAGE = 0
.reset : & RESET, PAGE = 0, TYPE = DSECT /* not used, */
FFTtf & PRAMH0, PAGE = 0
FFTipcb ALIGN(256): {}& DRAMH0, PAGE = 1
FFTmag & DRAMH0, PAGE = 1
.stack : & RAMM1, PAGE = 1
.ebss : & DRAMH0, PAGE = 1
.econst : & DRAMH0, PAGE = 1
.esysmem : & DRAMH0, PAGE = 1
在线时间1 小时
TA的帖子TA的资源
一粒金砂（初级）, 积分 2, 距离下一级还需 3 积分
一粒金砂（初级）, 积分 2, 距离下一级还需 3 积分
上午测试过，没问题。
主程序中：
#define N 128&&
#pragma DATA_SECTION(ipcb, \&FFTipcb\"
#pragma DATA_SECTION(mag, \&FFTmag\"
RFFT32 fft=RFFT32_128P_DEFAULTS;
long ipcb[N+2];
long mag[N/2+1];
for(i=0;i&N;i++)
&&ipcb=INDATA*65536;//赋值，128点的采样值&&
//y=0.25cos(314*t)+0.25cos(314*3*t)/3+0.25cos(314*9*t)/9& && &&&
fft.ipcbptr=
fft.magptr=
fft.init(&fft);
RFFT32_brev(ipcb,ipcb,N);
fft.calc(&fft);
fft.split(&fft);
fft.mag(&fft);
for(i=0;i&N;i++)
mag=sqrt(mag);
在CMD文件中添加\&FFTipcb\&、 \&FFTmag\&空间，方法见TI文档
还要包含fft.h、fft.lib
最好看一下TI提供的FFT的文档，说明的很仔细。
在线时间1 小时
TA的帖子TA的资源
一粒金砂（初级）, 积分 2, 距离下一级还需 3 积分
一粒金砂（初级）, 积分 2, 距离下一级还需 3 积分
对，只能是Q30或者是Q31
在线时间0 小时
TA的帖子TA的资源
一粒金砂（初级）, 积分 2, 距离下一级还需 3 积分
一粒金砂（初级）, 积分 2, 距离下一级还需 3 积分
那你给的这段程序怎样才能看出你用的是那种Q格式呢？我不大理解。
在线时间0 小时
TA的帖子TA的资源
一粒金砂（初级）, 积分 2, 距离下一级还需 3 积分
一粒金砂（初级）, 积分 2, 距离下一级还需 3 积分
我现在是外扩得AD，是16位宽度的，而Q30是32位的，怎样转换才可以呢？
在线时间1 小时
TA的帖子TA的资源
一粒金砂（初级）, 积分 2, 距离下一级还需 3 积分
一粒金砂（初级）, 积分 2, 距离下一级还需 3 积分
我一直用的54X系列的，我调用的是DSPLIB，优化好了的程序例如FFT，添加#include &displib.h&
#include &tms320.h&就可以调用优化后的程序了
在线时间0 小时
TA的帖子TA的资源
一粒金砂（初级）, 积分 2, 距离下一级还需 3 积分
一粒金砂（初级）, 积分 2, 距离下一级还需 3 积分
我现在还是不大懂，为何我输入的是1024点，输出的是512点的频谱呢，帮忙解释一下阿。
Powered by
逛了这许久，何不进去瞧瞧？Copyright &copy
600万学术达人首选