请问在编写DSP程序时，数学函数如：sprt这种开方的，需要定义什么头文件吗

如何使用STM32F4的DSP库

我们平常所使用的CPU为定点CPU，意思是进行整点数值运算的CPU。当遇到形如11+11的浮点数运算时，定点CPU就遇到大难题了。对于32位单片机，利用Q化处理能发挥他本身的性能，但是精度和速度仍然不会提高很多。

现在设计出了一个新的CPU，叫做FPU，这个芯片专门处理浮点数的运算，这样处理器就将整点数和浮点数分开来处理，整点数交由定点CPU处理而浮点数交由FPU处理。我们见到过TI的DSP，还有STM32F4系列的带有DSP功能的微控制器。前者笔者没有用过，不作评论，而后者如果需要用到FPU的浮点运算功能，必须要进行一些必要的设置。

首先，由于浮点运算在FPU中进行，所以首先应该使能FPU运行。在system_init()中，定义__FPU_PRESENT和__FPU_USED

/ FPU settings------------------------------------------------------------/

#if (__FPU_PRESENT == 1)&& (__FPU_USED == 1)

SCB->CPACR |= ((3UL<< 102)|(3UL << 112)); /set CP10 and CP11 Full Access /

#endif

这样就使能了FPU。

对于上述改变，当程序中出现这种简单的加减乘除运算FPU就起作用了。但是对于复杂的如三角运算、开方运算等，我们就需要加入mathh头文件。但是如果单纯的加入他，那么Keil会自动调用内部的mathh，该头文件是针对ARM处理器的，专门用于定点CPU和标准算法（IEEE-754）。对于使用了FPU的STM32F4是没有任何作用的。所以，需要将mathh换成ST的库，即arm_mathh。在该头文件中，涉及到另一个文件core_cmxh（x=0、3、4），当然了，如同STM32F1系列一样，在工程中加入core_cm4h即可。

到这里，算是全部设置完毕，之差最后一步，调用！但是别小看了这一步，因为如果调用的不正确，全面的设置就白费了。在使用三角函数如sin()、cos()时不要直接写如上形式，因为他们函数的名字来自于mathh，所以你调用的仍旧是Keil库中的标准mathh。要使用arm_mathh中的arm_sin_f32()函数（见Line5780，原函数见DSP_Lib\Source\FastMathFunctions），可以看到他利用的是三次样条插值法快速求值（见Line263 / Cubic interpolation process /）。

注意一下例外函数，sqrt()，在arm_mathh中为arm_sqrt_f32()。使用他的时候需要同时开启#if(__FPU_USED == 1) && defined ( __CC_ARM )才行，切记！还可以发现开方函数还有q15和q31之分，我想他们的区别就是精度的问题，但是他们没有应用FPU来计算，说白了就是利用0x5f3759df这个数进行快速开方

您好，1 首先是接口的预定义

----------------------------------------------

#define LCD_DATA (((volatile Uint16 )0x0070E0)) // GPIOA7-A0对应DB7-DB0

#define RS GpioDataRegsGPBDATbitGPIOB0

#define RW GpioDataRegsGPBDATbitGPIOB1 //别弄错0 1 2

#define EN GpioDataRegsGPBDATbitGPIOB2 // 实际接线要对应

void InitGpio(void)

{

EALLOW;

GpioMuxRegsGPAMUXbitPWM1_GPIOA0 = 0; // 设置为普通GPIO使用

GpioMuxRegsGPADIRbitGPIOA0 = 1; // 设置为输出

GpioMuxRegsGPAMUXbitPWM2_GPIOA1 = 0;

GpioMuxRegsGPADIRbitGPIOA1 = 1;

GpioMuxRegsGPAMUXbitPWM3_GPIOA2 = 0;

GpioMuxRegsGPADIRbitGPIOA2 = 1;

GpioMuxRegsGPAMUXbitPWM4_GPIOA3 = 0;

GpioMuxRegsGPADIRbitGPIOA3 = 1;

GpioMuxRegsGPAMUXbitPWM5_GPIOA4 = 0;

GpioMuxRegsGPADIRbitGPIOA4 = 1;

GpioMuxRegsGPAMUXbitPWM6_GPIOA5 = 0;

GpioMuxRegsGPADIRbitGPIOA5 = 1;

GpioMuxRegsGPAMUXbitT1PWM_GPIOA6 = 0;

GpioMuxRegsGPADIRbitGPIOA6 = 1;

GpioMuxRegsGPAMUXbitT2PWM_GPIOA7 = 0;

GpioMuxRegsGPADIRbitGPIOA7 = 1;

GpioMuxRegsGPBMUXbitPWM7_GPIOB0 = 0;

GpioMuxRegsGPBDIRbitGPIOB0 = 1;

GpioMuxRegsGPBMUXbitPWM8_GPIOB1 = 0;

GpioMuxRegsGPBDIRbitGPIOB1 = 1;

GpioMuxRegsGPBMUXbitPWM9_GPIOB2 = 0;

GpioMuxRegsGPBDIRbitGPIOB2 = 1;

EDIS;

}

----------------------------------------------

一般液晶的控制线是直接对I/O口的位进行 *** 作，数据线是按字进行 *** 作。在这容易出错的是：（1）数据线地址的对应。DSP的GPIO数据地址一般为16位一个地址（F28335有的是32个GPIO一组，给出了一个地址，实际上是有两个地址的，给出的那一个地址是低16位的）。需要注意的是，液晶数据线一般为8位，那么把八位数据送出的时候，实际给的是DSP的16位数据的低八位，所以接线上要接低八位的GPIO；如果接高八位的GPIO，软件上要用下面一行程序进行移位 dat = dat << 8; //左移8位，向高位移动。（2）在进行GPIO初始化和预定义的时候，一般都会复制，但是别忘记改一些0 1 2 3等数，接线上也要一一对应，仔细检查。

2 51程序移植到DSP的时序问题

----------------------------------------------

void Display_Data_All(uchar hz)

{

while(hz != '\0')

{

WriteData12864(hz);

hz++;

delay(20);//2就不够！！！！！！

}

}

----------------------------------------------

由于51单片机的晶振一般为110592MHz，而DSP等控制器的晶振为30MHz，实际执行起来最高有150MHz，而液晶为低速外设，所以移植后可能会不显示，显示乱码等情况。我在调试12864液晶的时候就出现过只显示乱码数字不显示汉字的情况，这不是字库损坏，而是因为写汉字的时间要比写数字的时间长，而程序中延时过短。上面程序中把delay(2)改为delay(20)就解决问题了。

实际上，真正造成影响的是，程序执行过快。它认为显示完一个字之后，又很快进入下一个字的 *** 作；实际上液晶要一定的时间才能写完（见液晶 *** 作时序图），所以写数据的程序中要加长延时。至于RS、RW、EN等控制引脚，延时与否影响不大。

3 240128液晶的调试

240128液晶有busy和int返回信号，实际上不需要接即可。程序中也可以不测忙。。程序中写控制指令两者中间也要加长延时，更不用说写数据之间的延时。

----------------------------------------------

void lcd_regwrite(Uint16 regname,Uint16 regdata) // 写控制指令

{

lcd_regwr(regname);

delay(10); // 加长延时

lcd_regwr(regdata);

}

void lcd_character(uchar cha,int count) // 显示中文或字符

{

int i;

for(i=0;i<count;i++)

{

delay(10); // 加长延时

lcd_datawrite(cha);

++cha;

}

}

将FLASH中的程序COPY到RAM中运行，通常的目的是加快程序的运行速度，通常有两种情况需要这样去 *** 作：

1、程序中对基要求比较高的函数，如中断；

2、程序需要对FLASH进行 *** 作，这时就要把程序先复制到RAM中运行然后才能对FLASH *** 作。

RamfuncsLoadStart、RamfuncsLoadEnd、RamfuncsRunStart这三个变量是在CMD文件中创建的，创建方式如下：

LOAD_START(RamfuncsLoadStart),

LOAD_END(RamfuncsLoadEnd),

RUN_START(RamfuncsRunStart),

分别表示了装载函数的首地址，装载函数的结束地址和装载函数的运行地址；

执行完MemCopy(&RamfuncsLoadStart, &RamfuncsLoadEnd, &RamfuncsRunStart);后，便将FLASH中相关的程序COPY到了RAM中，之后的程序运行时，只要调用FLASH中RamfuncsLoadStart地址开始的相关函数，系统都会自动地指向RAM中相应的函数入口地址运行。

要做好有效值运算，必须遵守相关的前提和约束。

1、dsp是对数字信号进行处理，因此，先要保证源头的数据的正确性和充分性。首先，电路带宽必须大于信号带宽，其次，采样频率必须满足采样定理要求，至少是信号带宽的2倍以上，建议采用10倍以上。

2、进行方均根运算时，关键在于“均”，也就是必须明白是多少数据的“均”。对于交流信号，严格讲，必须是信号周期的整倍数时间内的采样点的“均”。因此，运算前，必须知道信号的频率（周期）。信号频率除以采样频率再乘以任意自然数都可以作为运算的点数。

3、对于低频信号而言，若运算的信号周期数太大，数据更新时间太长，可设置较少的周期数。对于较高频率的信号，可以设置较多的周期数。

4、对于较高频率的信号，若运算的点数足够多，也就是信号的时间跨度远远大于信号周期，可以不考虑整数倍的问题。

5、有了上述基础，按照计算或设定的运算点数，提取最近采样的相应点数的数据，依次平方、求和、开方，即可得到有效值。

何编程语言是难学的,真正的工作中有时候学习一个新语法只有不到一周的时间,语法而言都一样,如果还停留在语法层学习的上面,编写程序的道路你就连门都没有入

但是后续的 数据结构和算法 稍微要看下书了,开始接触算法了,着和语法没有关系,之所以要先学好语法是为了能看懂用语法描叙的算法,学通了用任何语言来描叙都一样,到了这个阶段的就相当于抬起一条腿准备跨到起跑线一样,但还是没有入门

到编译原理 和 图灵机 再到 自己编写微型的 *** 作系统 就需要有个老师来引导了

我学软件开发就是从C开始一路由 ->国家数据库3级-> 程序员->高级程序员->系统分析师考上去然后通过近9年的工作,体会是如果你能够在市面上或者学校里精确买到或者学到的一门知识往往就表示你的层次还是不能让你得兴应手的在行业中创造你能想像的东西,真正的解决问题需要很多本书再加上80%以上的自己创造和理解才能做出来的时候才算真正懂得了你的工作

一个程序员需要这样的经历

楼上的估计最高学历不过是理工科的一个大学生,应该还没毕业,没有在计算机行业呆过,在做芯片处理的工作中比如DSP芯片设计需要你天天接触算法,但是这和C语言本身有关系吗傅立叶算法是数学家傅立叶设计出来的数学模型但是不适合做在计算机软件里面(运算量太大了比如离散的傅立叶变换等同于用序列Y(n×1列矢量)乘以n×n矩阵Fn，需要n×n次乘法。若n=1024,则是104，8576次乘法运算。什么概念呢？如果你选用的CPU单周期指令为25ns, 单周期也可以完成一次乘法运算，那么要计算1024点的傅立叶变换则需要262144ms,这还不包括加法或其它运算),我给出C算法如下:

void kkfft(double pr[], double pi[], int n, int k, double fr[], double fi[], int l, int il)

{

int it,m,is,i,j,nv,l0;

double p,q,s,vr,vi,poddr,poddi;

for (it=0; it<=n-1; it++)

{

m = it;

is = 0;

for(i=0; i<=k-1; i++)

{

j = m/2;

is = 2is+(m-2j);

m = j;

}

fr[it] = pr[is];

fi[it] = pi[is];

}

//----------------------------

pr[0] = 10;

pi[0] = 00;

p = 6283185306/(10n);

pr[1] = cos(p);

pi[1] = -sin(p);

if (l!=0)

pi[1]=-pi[1];

for (i=2; i<=n-1; i++)

{

p = pr[i-1]pr[1];

q = pi[i-1]pi[1];

s = (pr[i-1]+pi[i-1])(pr[1]+pi[1]);

pr[i] = p-q;

pi[i] = s-p-q;

}

for (it=0; it<=n-2; it=it+2)

{

vr = fr[it];

vi = fi[it];

fr[it] = vr+fr[it+1];

fi[it] = vi+fi[it+1];

fr[it+1] = vr-fr[it+1];

fi[it+1] = vi-fi[it+1];

}

m = n/2;

nv = 2;

for (l0=k-2; l0>=0; l0--)

{

m = m/2;

nv = 2nv;

for(it=0; it<=(m-1)nv; it=it+nv)

for (j=0; j<=(nv/2)-1; j++)

{

p = pr[mj]fr[it+j+nv/2];

q = pi[mj]fi[it+j+nv/2];

s = pr[mj]+pi[mj];

s = s(fr[it+j+nv/2]+fi[it+j+nv/2]);

poddr = p-q;

poddi = s-p-q;

fr[it+j+nv/2] = fr[it+j]-poddr;

fi[it+j+nv/2] = fi[it+j]-poddi;

fr[it+j] = fr[it+j]+poddr;

fi[it+j] = fi[it+j]+poddi;

}

}

if(l!=0)

for(i=0; i<=n-1; i++)

{

fr[i] = fr[i]/(10n);

fi[i] = fi[i]/(10n);

}

if(il!=0)

for(i=0; i<=n-1; i++)

{

pr[i] = sqrt(fr[i]fr[i]+fi[i]fi[i]);

if(fabs(fr[i])<0000001fabs(fi[i]))

{

if ((fi[i]fr[i])>0)

pi[i] = 900;

else

pi[i] = -900;

}

else

pi[i] = atan(fi[i]/fr[i])3600/6283185306;

}

return;

}

另外,团IDC网上有许多产品团购,便宜有口碑

由于正弦波是单一频率信号。可以采取下述措施： 1、可以在AD输入根据信号频率加一个带通滤波器和一个抗混叠的低通滤波器，或者共用一个低通滤波器。 2、如果滤波器对有用信息造成了不可忽视的衰减，可以根据滤波器的幅频响应特性及信号频率在dsp中做补偿。 3、输出加一个低通滤波器或积分器，消除DA量化时造成的“小台阶”，低通滤波器的截止频率应该高于信号频率，远远低于DA的转换频率。

以上就是关于请问在编写DSP程序时，数学函数如：sprt这种开方的，需要定义什么头文件吗全部的内容，包括:请问在编写DSP程序时，数学函数如：sprt这种开方的，需要定义什么头文件吗、DSP28335控制12864液晶，读取LCD数据，lcd_read_data程序应该怎么编写、dsp 如何使用自己写的.h文件等相关内容解答，如果想了解更多相关内容，可以关注我们，你们的支持是我们更新的动力！