debug正常，51运算错误

机械绿藤

debug正常，51运算错误 [复制链接]

已经用其他程序试验过了，不是xdata的问题，xdata里的数据能正常使用。我在用51做一个32点的fft。
部分代码如下：
#include
#include
#define  uchar  unsigned  char
#define  uint unsigned int
#define  N 32 //fft点数
float xdata re[32]={0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31};//实部
float xdata im[32]=0;//虚部
uint x[32];
float code RSin[16] = {1.0000,0.9808,0.9239,0.8315,0.7071,0.5556,0.3827,0.1951,0,-0.1951,-0.3827,-0.5556,-0.7071,-0.8315,-0.9239,-0.9808};
float code RCos[16] = {0.0000,0.1951,0.3827,0.5556,0.7071,0.8315,0.9239,0.9808,1.0000,0.9808,0.9239,0.8315,0.7071,0.5556,0.3827,0.1951};
float tempRe;//实部缓存
float tempIm;//虚部缓存
float sinVal = 0;//正弦值
float cosVal = 0;//余弦值
uint s=0;//
uint group_intf = 0;//
uint num =0;//
uint i = 0;//
uint indexA  =0;//
uint indexB = 0;//
uint index_Rotate = 0;//
uint group = 0;//
uint stage = 0;//
uint g =0;//
char TempData[3];
上面是程序用到的主要变量的定义
***************************************************************************************************
下面是主体
main()
{
//倒位序算法
tempRe=re[1];re[1]=re[16];re[16]=tempRe;
tempRe=re[2];re[2]=re[8];re[8]=tempRe;
tempRe=re[3];re[3]=re[24];re[24]=tempRe;
tempRe=re[5];re[5]=re[20];re[20]=tempRe;
tempRe=re[6];re[6]=re[12];re[12]=tempRe;
tempRe=re[7];re[7]=re[28];re[28]=tempRe;
tempRe=re[9];re[9]=re[18];re[18]=tempRe;
tempRe=re[11];re[11]=re[26];re[26]=tempRe;
tempRe=re[13];re[13]=re[22];re[22]=tempRe;
tempRe=re[15];re[15]=re[30];re[30]=tempRe;
tempRe=re[19];re[19]=re[25];re[25]=tempRe;
tempRe=re[23];re[23]=re[29];re[29]=tempRe;

//*****************fft算法*************************//
    stage = 5;
for(s = 0; s < stage; s++)  //级
{
  group = N/2/(int)(0x01<   group_intf = N / group;  //组间下标间隔

  for(g = 0; g < group; g++)
  {
   num = (int)(0x01 << s);  //组内蝶形数
   for (i = 0; i < num; i++)
   {
    indexA = g * group_intf + i;
    indexB = indexA + (int)(0x01<     index_Rotate = i * (int)(0x01<<(stage-s-1));  //旋转因子下标
    sinVal = RSin[index_Rotate];
    cosVal = RCos[index_Rotate];
   
    tempRe = (re[indexB] * sinVal - im[indexB] * cosVal);
    tempIm = (re[indexB] * cosVal + im[indexB] * sinVal);

    re[indexB] = re[indexA] - tempRe;
    im[indexB] = im[indexA] - tempIm;
    re[indexA] = re[indexA] + tempRe;
    im[indexA] = im[indexA] + tempIm;

    }
  }
}

//求实部与虚部的平方和，再开根号

   for(i=0;i<32;i++)
   {
     tempRe=sqrt(re*re+im*im);x=tempRe;
   }
//将结果显示在1602上，用其他程序，1602显示正常
LCD_Init();
LCD_Clear();//清屏
while (1)
{  
for(i=0;i<16;i=i+4)
  {
   TempData[0]=x[i/4]/100;//分解显示信息
   TempData[1]=(x[i/4]%100)/10;
   TempData[2]=(x[i/4]%100)%10;
   LCD_Write_Char(i,0,TempData[0]+48);
   LCD_Write_Char(i+1,0,TempData[1]+48);
   LCD_Write_Char(i+2,0,TempData[2]+48);
   LCD_Write_Char(i+3,0,'#');
  }
  for(i=16;i<32;i=i+4)
  {
   TempData[0]=x[i/4]/100;//分解显示信息
   TempData[1]=(x[i/4]%100)/10;
   TempData[2]=(x[i/4]%100)%10;
   LCD_Write_Char(i-16,1,TempData[0]+48);
   LCD_Write_Char(i-16+1,1,TempData[1]+48);
   LCD_Write_Char(i-16+2,1,TempData[2]+48);
   LCD_Write_Char(i-16+3,1,'#');
  }
  while(1);
}
}
************************************
最后显示结果“y35#y35#y35#y35#”“y35#y35#y35#y35#”。
怀疑原因在 for(i=0;i<32;i++){tempRe=sqrt(re*re+im*im);x=tempRe;}部分。
若直接用x=re（正常结果，re=469,-15,-15...,-16）,则469能够正常显示，后面的显示“y35”。

机械绿藤

有没有人能告诉我。51做数据运算应该注意什么？

zqjqq88

尽量使用加减法移位等方式代替乘除法，尤其是除法！

机械绿藤

内容那我的程序错在哪呢？

辛昕

先直接在显示到1602之前，用一组明确的数据填充到tempdata里先，确定1602显示对不对？即使不是1602底层有问题，也可以把问题范围缩小到显示那一部分的代码

辛昕

两个疑问，进入显示死循环之前，求虚实部平方和，为什么好端端要循环32次？

辛昕

2，已经进入死循环，为什么在最后一句还要死循环？

辛昕

还有，为什么1602显示部分，都是显示一个两位数，为什么下半部分加了48，后面就不加了？

机械绿藤

1602显示已测试，如果跳过算法部分，显示正常
之所以循环32次是因为那是一个数组，32个元素。
死循环确实多余了。
上下部分都加了啊。那是为了将数据转换成asci码

[ 本帖最后由 机械绿藤 于 2013-4-7 19:27 编辑 ]

ltbytyn

这样的程序能执行通过，我表示怀疑。

//求实部与虚部的平方和，再开根号

   for(i=0;i<32;i++)
   {
     tempRe=sqrt(re*re+im*im);x=tempRe;   
   }

看看这一句，x=tempRe;   x在什么地方定义了？还有循环32次求模值后，x最终求的是31次模值

苏莎莎

感觉到大家的热心解答了，找到自己的问题了吗？要不要帮你再问问其他网友？

机械绿藤

原程序里，有中括号，不知道复制过来怎么都没有了。循环次数确实错了，但应该不影响运算吧？
我是用手机上论坛，不能针对每位解答我的朋友回复，请大家见谅。

ltbytyn

fft计算部分没有问题,你的问题在显示部分。

看看你的循环，一个模值占4位（3位数字+“#”），16个模值占多少位？1602一行能显示几位数字？太狗血了
for(i=0;i<16;i=i+4)
  {
   TempData[0]=x[i/4]/100;//分解显示信息
   TempData[1]=(x[i/4]%100)/10;
   TempData[2]=(x[i/4]%100)%10;
   LCD_Write_Char(i,0,TempData[0]+48);
   LCD_Write_Char(i+1,0,TempData[1]+48);
   LCD_Write_Char(i+2,0,TempData[2]+48);
   LCD_Write_Char(i+3,0,'#');
  }

机械绿藤

我只让他显示前八个结果，以检验结果是否正确。

机械绿藤

循环次数没错

naga568

ltbytyn

for(i=0;i<16;i=i+4)
  {
   TempData[0]=x[i/4]/100;//分解显示信息
   TempData[1]=(x[i/4]%100)/10;
   TempData[2]=(x[i/4]%100)%10;
   LCD_Write_Char(i,0,TempData[0]+48);
   LCD_Write_Char(i+1,0,TempData[1]+48);
   LCD_Write_Char(i+2,0,TempData[2]+48);
   LCD_Write_Char(i+3,0,'#');
  }

i+3  应该是x坐标吧。1602一行显示16个字符，所以i+3<16是必须的。你看看自己程序

机械绿藤

i＝0，4，8，12，不会超出显示范围，显示部分没有错误，已经测试过了。我怀疑是浮点数与整形转换的问题，但就是查不出来。