CCS3.3 5402 浮点运算 诡异问题求解

gongxd

请教高手们一个问题
我写了一段程序
用ccs3.3 dsp5402
仿真的时候是没有错误的
脱离仿真器运行就错误，程序也是在运行的，就是 有浮点运算的地方有问题

比如
x=11.0;
peakError=0.0001;
sampleRate=882000;
peakx=(float)x+peakError;
peakfrequency=(float)(sampleRate)/(float)(128.0) ;

peakfrequency的结果却是负数


还有 用double的 地方 也出现错误
仿真没问题 运行就出错
这个问题已经困扰我好几天了 现在头大了 请大家帮帮我  谢谢

zhangsusu

看到你到处发帖 5402本身就是一款定点的dsp 浮点算法都是通过C模拟的 效率很低
你需要查54x的C programming 手册 看看里面浮点算法的说明
不过 我觉得一个切实可行的方法就是描述清楚你的问题
包括 输入值范围 输出值范围 计算表达式（算法） 大家一起看能否定点（Q化）完成
定点算法本身就是dsp的特色 要是都去操作外设 dsp和51 也没啥大的区别了
论坛里讨论定点算法的比较少 可以借此活跃人气 呵呵

gongxd

回复【1楼】zhangsusu
-----------------------------------------------------------------------

谢谢 zhangsusu，呵呵

我现在疑惑的是 为什么 仿真 浮点运算是正确的 下载运行的就不对呢
当然通过定点可以解决这个问题
但是根本原因还是没找到啊
期待高人指点

程序如下






//计算阻尼比
void cacldamp(Int16 *sounddataorg)
{
//最大值点数组
Uint16 peakpointindex[peakpointcount];
Uint16 peakpoint[peakpointcount];
Uint16 peakpointindextmp;
//最大峰值点索引
Uint16 Maxpeakpointindex;
Uint16 i;
double damptmp;

//找出波形所有的峰值
peakpointindextmp=0;
for (i=SonicPM.SoundDecay.WindowStart;i<sounddataorgmax_1;i++)
{

        if (sounddataorg>=sounddataorg[i+1] && sounddataorg>sounddataorg[i-1])
        if (sounddataorg>=sounddataorg[i+2] && sounddataorg>sounddataorg[i-2])
        {
        //取过极大值点
                peakpointindex[peakpointindextmp]=i;
                peakpoint[peakpointindextmp]=sounddataorg;
                peakpointindextmp++;
        }
        //如果采集够指定的波形数就跳出循环
        if (peakpointindextmp>peakpointcount)
        {
                break;
        }

}//for

//如果极大值点个数大于0为正常，继续处理
if (peakpointindextmp>0)
{
        //for (i=0;i<peakpointindextmp;i++)
        //{
                if( peakpointindextmp>SonicPM.SoundDecay.WindowLength)
                {
                Maxpeakpointindex=getpeakpositonUint16(peakpoint,peakpointindextmp-SonicPM.SoundDecay.WindowLength);
                }
                else
                {
                Maxpeakpointindex=getpeakpositonUint16(peakpoint,peakpointindextmp);
                }


        //}

        //获得点数长度方便显示
        SonicPM.SoundDecay.WindowStartPoint=peakpointindex[Maxpeakpointindex];

        SonicPM.SoundDecay.WindowLengthPoint=peakpointindex[Maxpeakpointindex+SonicPM.SoundDecay.WindowLength];
        //计算峰值频率
        damptmp=100000*log(
        (        (double)peakpoint[Maxpeakpointindex]/
        (double)peakpoint[Maxpeakpointindex+SonicPM.SoundDecay.WindowLength]   )
        );
//阻尼比计算公式 1/n  * ln (An/A1)

/*
        damptmp=damptmp/
        ((double)peakpointindex[Maxpeakpointindex+SonicPM.SoundDecay.WindowLength]
        -
        (double)peakpointindex[Maxpeakpointindex]
        );
*/

        damptmp=damptmp/(double)(SonicPM.SoundDecay.WindowLength);
        damptmp=damptmp/2.0;
         damptmp=damptmp/3.1415926;

        SonicPM.DetDamp=(Int16)damptmp;
}
else
{
//计算峰值频率
SonicPM.DetDamp=-1;
}


其中 sounddataorg 是采集的声音数据 16位int
格式的

函数的功能是 找出波形的最大的幅值 然后间隔n个波形的最大值

计算波形的衰减系数

算式是  衰减=1/n * Ln(A1/An)

A1 An 都不超过 32768  也进行了强制转换 为double 了

仿真计算结果正确  下载后 计算结果为 0x7FFFFH

haeha

能用串口输出吗，看看从哪部开始计算结果不对了。
另外，用sizeof看下double的长度，我记得ccs3.3里,double和float精度是一样的

gongxd

能输出数据 可以运行程序后 用仿真器调试  看到执行 浮点乘除法时错误
加减没问题

但是用仿真器下载进程序去就没有问题

zhangsusu

回复【3楼】haeha  
-----------------------------------------------------------------------
用仿真器+watch 就行了 这么复杂的算法 串口搞不定

zhangsusu

回复【2楼】gongxd  
-----------------------------------------------------------------------

你仿真使用的simulator吗？
下载是用得load 还是 烧到flash上了？
对于
sampleRate=882000;
peakfrequency=(float)(sampleRate)/(float)(128.0) ;

何不写成
// 对于54x Int32 means long
Int32 sampleRate = 882000;
Int32 peakfrequency;

peakfrequency = sampleRate >> 7; // means sampleRate/128

gongxd

回复【6楼】zhangsusu
-----------------------------------------------------------------------

看了一天汇编，似乎找到问题了
我的程序 是通过单片机的 HPI 接口 boot的

上段程序中
damptmp=damptmp/2.0;
damptmp=damptmp/3.1415926;

用到了 2.0，3.1415926 的常数，编译到.const 了
汇编调试的时候，看到 这两个数据并没有初始化，所以计算错误


程序中还用到了 很多其他的常数，但都是整形的，似乎已经初始化了

我把

damptmp=damptmp/2.0;
改为
uint16 ddd=2;
damptmp=damptmp/(double)ddd;

就可以了


请问是这个样子么，那.const段要编写代码手动初始化么


如下的方法 可以么



设置.const段的装载地址和运行地址。其表达式如下。
.const : LOAD = FLASH PAGE 0 , RUN = EXT PAGE 1
{
__const_run = .;      /* 获取运行地址 */
*(.c_mark)        /* 标志装载地址 */
*(.const)          /* 分配.const段 */
__const_length = .- __const_run;  /* 计算段长 */
}
     采用了装载地址与运行地址分离的方式后，为了使程序正常运行，在初始化时，需要将.const段的内容从装载地址拷贝到运行地址内。这段程序可以在编译时由编译器自动生成。这还需要对软件所使用的库文件作简单的修改。该库名称即是rts.lib(表示不同类型的DSP，有2xx，25，50等等)。修改该库的方法是将源文件从库中提取出来，作修改，编译后再替代原有的文件。具体的操作如下：先执行DOS命令：
dspar -x rts.src boot.asm
这句的功能是从rts.src文件中提取出boot.asm文件。这个rts.src即是rts.lib的源文件。在boot.asm文件中能找到CONST_COPY这个标志量，为了实现所需要的功能它应被赋值为1。对boot.asm文件的编辑完成之后，就可以将其编译生成目标文件，执行语句：
dspa -v boot.asm
其中对于不同的DSP需要使用不同的参数，对于240xA来说，应该使用2xx来代替“”。语句执行完后会生成boot.obj文件。再执行语句：
dspar -r rts.lib boot.obj
这时它就替换了库里的同名文件。在编译时编译器就会自动增加拷贝.const段到数据空间的语句。这种方法不必修改程序，代价是牺牲了一定数据存储空间，时间开销主要出现在初始化中。这是最经济实用的方法之一。
使用代码工具
各种各样的代码生成工具可以大大增加编程的灵活性。除了前面提到的dspar与dspa两个指令外，C文法分析程序dspac.exe对.C文件进行文法分析，生成.IF中间文件，dspopt.exe对.IF文件进行优化，生成.OPT文件。此外，还有其它多种功能不同的代码工具。对于应用环境苛刻的程序，它们有利于分析如何对其进行优化。  





正在试验，dsp不太熟悉啊，呵呵



请问是这个样子么，那.const段要编写代码手动初始化么

gongxd

回复【6楼】zhangsusu
-----------------------------------------------------------------------

问题 解决 是这个地方的问题 谢谢 zhangsusu 的帮助

gongxd

又出现 新问题了

我的程序中用到了 FFT 函数 其中有正弦函数表
     .sintab

在段分配中
         .sintab :>       DARAM1   PAGE 1

放到 ram中 无法运行
也应该是没初始化

请问有什么好的方法 初始化 sintab么

gongxd

点击此处下载 ourdev_620059C3DH9W.rar(文件大小:232K) (原文件名:rtslib5000.rar)

HPIBoot cmd文件

birenpeter

回复【10楼】gongxd
-----------------------------------------------------------------------

太好了啊，我也遇到这个问题