请教 DSP FFT库到底怎么调用啊~
需要准备些什么？调用什么函数呢
只用DSPF_sp_bitrev_cplx 和DSPF_sp_cfftr2_dit 这个就可以了？
DSPF_sp_bitrev_cplx 中的index数组是怎么得到的啊？
谢谢
可是我用的DSPLIB 3.1.1中没有您说的那个函数啊
tw_genr2fft(wd,N); // 产生旋转因子
bit_rev1(wd, N>>1); // 把系数倒置，为下面函数输入准备
#include <math.h>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <ctype.h>
#include <DSPF_sp_cfftr2_dit.h>
//#include <utility.h>
# define PI 3.14159265358979323846
# define Numlength 32
void gen_w_r2(float* , int );
void bit_rev(float* , int );
void main()
{
int i;
float w[Numlength];
float Input[Numlength*2],Output[Numlength*2];
for(i=0;i<Numlength;i++) //输入Input
{
Input[2*i]=cos(PI*2*i/Numlength*3)*1024;
Input[2*i+1]=0;
}
gen_w_r2(w,Numlength); //产生旋转因子
// tw_genr2fft(w, Numlength); // 产生旋转因子
bit_rev(w,Numlength>>1); //把系数倒置，为下面函数输入准备
DSPF_sp_cfftr2_dit(Input,w,Numlength); ///计算该FFT，需要先调用上面两个函数
bit_rev(Input,Numlength); //输出是倒置的，所以重新倒置为正常
for(i=0;i<(Numlength*2);i++)
{
Output[i]=Input[i];
}
while(1);
}
void gen_w_r2(float* w, int n)
{
int i, j=1;
float pi = 4.0*atan(1.0);
float e = pi*2.0/n;
for(j=1; j < n; j <<= 1)
{
for(i=0; i < ( n>>1 ); i += j)
{
*w++ = cos(i*e);
*w++ = -sin(i*e);
}
}
}
void bit_rev(float* x, int n)
{
int i, j, k;
float 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;
}
}
}
东西有点多，刚开始上论坛，写的很杂，见谅，请给点指导，谢谢
penghui zhang 说：我用的是c6713，软件仿真，程序如下：#include <math.h>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <ctype.h>
#include <DSPF_sp_cfftr2_dit.h>
//#include <utility.h>
# define PI 3.14159265358979323846
# define Numlength 32void gen_w_r2(float* , int );
void bit_rev(float* , int );void main()
{
int i;
float w[Numlength];
float Input[Numlength*2],Output[Numlength*2];
for(i=0;i<Numlength;i++) //输入Input
{
Input[2*i]=cos(PI*2*i/Numlength*3)*1024;
Input[2*i+1]=0;
}gen_w_r2(w,Numlength); //产生旋转因子
// tw_genr2fft(w, Numlength); // 产生旋转因子
bit_rev(w,Numlength>>1); //把系数倒置，为下面函数输入准备
DSPF_sp_cfftr2_dit(Input,w,Numlength); ///计算该FFT，需要先调用上面两个函数
bit_rev(Input,Numlength); //输出是倒置的，所以重新倒置为正常
for(i=0;i<(Numlength*2);i++)
{
Output[i]=Input[i];
}
while(1);
}
void gen_w_r2(float* w, int n)
{
int i, j=1;
float pi = 4.0*atan(1.0);
float e = pi*2.0/n;
for(j=1; j < n; j <<= 1)
{
for(i=0; i < ( n>>1 ); i += j)
{
*w++ = cos(i*e);
*w++ = -sin(i*e);
}
}
}
void bit_rev(float* x, int n)
{
int i, j, k;
float 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;
}
}
}东西有点多，刚开始上论坛，写的很杂，见谅，请给点指导，谢谢
6713的例程里有没有要求变量要内存对齐呢？
类似这样的
#pragma DATA_ALIGN(CFFT_In,8);
使用#pragma预编译命令
我把我测试成功的DSP C6748的代码发上来你参考下哈
bin wang2 说：我用的是c6713，软件仿真，程序如下：#include <math.h>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <ctype.h>
#include <DSPF_sp_cfftr2_dit.h>
//#include <utility.h>
# define PI 3.14159265358979323846
# define Numlength 32void gen_w_r2(float* , int );
void bit_rev(float* , int );void main()
{
int i;
float w[Numlength];
float Input[Numlength*2],Output[Numlength*2];
for(i=0;i<Numlength;i++) //输入Input
{
Input[2*i]=cos(PI*2*i/Numlength*3)*1024;
Input[2*i+1]=0;
}gen_w_r2(w,Numlength); //产生旋转因子
// tw_genr2fft(w, Numlength); // 产生旋转因子
bit_rev(w,Numlength>>1); //把系数倒置，为下面函数输入准备
DSPF_sp_cfftr2_dit(Input,w,Numlength); ///计算该FFT，需要先调用上面两个函数
bit_rev(Input,Numlength); //输出是倒置的，所以重新倒置为正常
for(i=0;i<(Numlength*2);i++)
{
Output[i]=Input[i];
}
while(1);
}
void gen_w_r2(float* w, int n)
{
int i, j=1;
float pi = 4.0*atan(1.0);
float e = pi*2.0/n;
for(j=1; j < n; j <<= 1)
{
for(i=0; i < ( n>>1 ); i += j)
{
*w++ = cos(i*e);
*w++ = -sin(i*e);
}
}
}
void bit_rev(float* x, int n)
{
int i, j, k;
float 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;
}
}
}东西有点多，刚开始上论坛，写的很杂，见谅，请给点指导，谢谢
6713的例程里有没有要求变量要内存对齐呢？
类似这样的
#pragma DATA_ALIGN(CFFT_In,8);
使用#pragma预编译命令
我把我测试成功的DSP C6748的代码发上来你参考下哈
[/quote]
工程中引用的函数
原文件名 DSPF_sp_fftSPxSP.s
原文件名 DSPF_sp_ifftSPxSP.s
penghui zhang 说：关于变量对齐的问题我在C6713的例程里没有，您给的这两个.s文件是汇编程序，我不知道如何调用，我把函数DSPF_sp_cfftr2_dit()自己用C在主函数之前重新定义，运行可以通过，问题就出在用库文件的DSPF_sp_cfftr2_dit()运行出错，您能不能推荐几个资料我看看，以下是我现在在看的资料，另外TI给的英文资料TMS320C67x_DSP_Library_Programmer's_Reference_Guide.pdf。谢谢您给我之前的帮助
直接声明函数名就可以了~
你可以试试安装新版本的DSP库
资料的话直接看DSP库自带的说明文档就好了
不客气~