IAR开发stm8系列,C语言实现16位乘法器和32位除法器函数

stm8是8位单片机,在ADC采样采用12bit采样值进行定点整型运算的时候,为了保证精度需要通过16位乘法器 进行扩大,通过32位除法器缩小运算。但是用c语言直接用"*"和"/"计算是无法实现。c语言的math函数库也没有提供这样的计算函数,只好自己编写了一段。开始想用c语言嵌套汇编的方式来做,但是发现调试通过的函数,放到大程序中执行总是计算错误,只好放弃,完全用c语言来实现,过程有点复杂,这里共享出来,给有需要的人。

1、16位乘法器的实现:采样2进制乘法的基础计算方法,用乘数的0-15位分别与被乘数相乘,结果放到x0(低16位),x1(高16位)中。y0,y1用来放置每一位乘数与被乘数的乘积,然后把结果和x0,x1累加。

2、因为c语言不能直接访问cpu_ccr_c来判断加法溢出,需要通过z1=x0+y0预加一次,看结果是否比两个加数小,如果z1比两个加数都小说明是溢出了,那么高位y1就加1。if((z1<x0)&(z1<y0)==1){x1=x1+1;}//判断x0+y0是否会溢出

3、低16位向高16位移位时,需要取出低16位数的高8位进行判断(8位单片机无法直接判断16位数的高8位数据结果),方法是对这个数右移8位,再赋值给一个8位的变量如y_h=y0>>8;然后用if((y_h&0x80)==0x80){y1=y1+1;}判断最高位是否为1,并进行相应的处理。

4、32位除法器也是类似的细节处理,基本原理采用的也是二进制计算的基础原理,从高位往低位除,把多次除法的商拼接起来,上一次除法的余数作为下一次除法的被除数一部分和低16位剩余的数进行移位拼接再做除法,循环往复。

前几日发布的算法,编译时总有问题,我也不清楚错误出在哪里,有高手请指点一下。

现在实现的算法也是基于过去做51单片机时的汇编语言基础思想做的,看起来啰嗦和繁琐,如果有更好更简洁的方法,还请指教。

例程可以参考如下:

cpp 复制代码
  //r1*r2=>mul_rs[0]是低16位,mul_rs[1]是高16位
int Mul_16Bit(uint16_t r1,uint16_t r2,uint16_t *mul_r ){
  uint16_t x0,x1,y0,y1,z0,z1;
  uint8_t i,y_h;
  x0=0;
  x1=0;
  y0=r1;
  y1=0;
  z0=r2;
  for(i=0;i<=15;i++){
    if(z0&0x1==0x1){
      z1=x0+y0;
      if((z1<x0)&(z1<y0)==1){x1=x1+1;}//判断x0+y0是否会溢出
      x0=x0+y0;
      x1=x1+y1;
    }
    y1=y1<<1;
    y_h=y0>>8;
    y0=y0<<1;
    if((y_h&0x80)==0x80){y1=y1+1;}
    if (z0==0){break;}
    z0=z0>>1;
  } 
  mul_r[0]=x0;
  mul_r[1]=x1;  
  return 0;
}

  //r2r1/r3=>mul_r[0]是低16位,mul_rs[1]是高16位
int Div_32Bit(uint16_t r2,uint16_t r1,uint16_t r3,uint16_t *mul_r ){
  uint16_t x0,x1,y1;
  uint16_t z;
  uint8_t i,sh,rh;
  i=0;
  x0=0;
  x1=r2/r3;//高位商数
  y1=r2%r3;//高位余数
  if (y1==0){
    x0=r1/r3;//低位商数
  }
  else{
    while(i<=15){
      sh=y1>>8;//取出y1的高8位
      while((sh&0x80)!=0x80){
        x0=x0<<1;
        y1=y1<<1;
        rh=r1>>8;//取出r1的高8位
        if((rh&0x80)==0x80){y1=y1+1;}
        r1=r1<<1;
        sh=y1>>8;//取出y1的高8位
        i++;
        if(i>15){break;}
      }
      z=y1/r3;
      x0=x0+z;
      y1=y1%r3;
    }
  }
  mul_r[0]=x0;
  mul_r[1]=x1;  
  return 0;
}
相关推荐
侃侃_天下9 小时前
最终的信号类
开发语言·c++·算法
echoarts10 小时前
Rayon Rust中的数据并行库入门教程
开发语言·其他·算法·rust
Aomnitrix10 小时前
知识管理新范式——cpolar+Wiki.js打造企业级分布式知识库
开发语言·javascript·分布式
每天回答3个问题10 小时前
UE5C++编译遇到MSB3073
开发语言·c++·ue5
伍哥的传说11 小时前
Vite Plugin PWA – 零配置构建现代渐进式Web应用
开发语言·前端·javascript·web app·pwa·service worker·workbox
晶振厂家-晶发电子11 小时前
晶振在5G时代的角色:高精度时钟的核心支撑
单片机·嵌入式硬件·5g·晶振·电子元器件·晶振知识
F1372980155711 小时前
WD5030A 芯片,12V降5V,输出电流12A,电路设计
stm32·单片机·嵌入式硬件·汽车·51单片机
小莞尔11 小时前
【51单片机】【protues仿真】基于51单片机的篮球计时计分器系统
c语言·stm32·单片机·嵌入式硬件·51单片机
三佛科技-1873661339711 小时前
分享机械键盘MCU解决方案
单片机·嵌入式硬件·计算机外设
李永奉11 小时前
51单片机-使用IIC通信协议实现EEPROM模块教程
单片机·嵌入式硬件·51单片机