msc.51fudian.asm

51汇编浮点子程序 51汇编浮点子程序

 MOV @R0,A
 DEC R0
 POP ACC
 MOV @R0,A
 MOV A,R3 ；将返回地址压入堆栈
 PUSH ACC
 MOV A,R2
 PUSH ACC
 RET  ；返回主程序

    （１１）  标号：  ＦＣＭＰ    功能：浮点数代数值比较（不影响待比较操作数）

入口条件：待比较操作数分别在[R0]和[R1]中。
出口信息：若CY=1，则[R0] < [R1]，若CY=0且A=0则 [R0] = [R1]，否则[R0] > [R1]。
影响资源：A、B、PSW    堆栈需求：  ２字节

FCMP: MOV A,@R0 ；数符比较
 XRL A,@R1
 JNB ACC.7,CMP2
 MOV A,@R0 ；两数异号，以[R0]数符为准
 RLC A
 MOV A,#0FFH
 RET  
CMP2: MOV A,@R1 ；两数同号，准备比较阶码
 MOV C,ACC.6
 MOV ACC.7,C
 MOV B,A

MOV A,@R0
 MOV C,ACC.7
 MOV F0,C ；保存[R0]的数符
 MOV C,ACC.6
 MOV ACC.7,C
 CLR C ；比较阶码
 SUBB A,B
 JZ CMP6
 RLC A ；取阶码之差的符号
 JNB F0,CMP5
 CPL C ；[R0]为负时，结果取反
CMP5: MOV A,#0FFH ；两数不相等
 RET  
CMP6: INC R0 ；阶码相同时，准备比较尾数
 INC R0
 INC R1
 INC R1
 CLR C
 MOV A,@R0
 SUBB A,@R1
 MOV B,A ；保存部分差 
 DEC R0
 DEC R1
 MOV A,@R0
 SUBB A,@R1
 DEC R0
 DEC R1

ORL A,B ；生成是否相等信息
 JZ CMP7
 JNB F0,CMP7
 CPL C ；[R0]为负时，结果取反
CMP7: RET

    （１２）  标号：  ＦＡＢＳ    功能：浮点绝对值函数

入口条件：操作数在[R0]中。
出口信息：结果仍在[R0]中。
影响资源：A    堆栈需求：  ２字节

FABS: MOV A,@R0 ；读取操作数的阶码
 CLR ACC.7 ；清除数符
 MOV @R0,A ；回传阶码
 RET

    （１３）  标号：  ＦＳＧＮ    功能：浮点符号函数

入口条件：操作数在[R0]中。
出口信息：累加器 A=1 时为正数，A=0FFH时为负数，A=0 时为零。

影响资源：PSW、A    堆栈需求：  ２字节

FSGN: INC R0 ；读尾数
 MOV A,@R0
 INC R0
 ORL A,@R0
 DEC R0
 DEC R0
 JNZ SGN
 RET  ；尾数为零，结束
SGN: MOV A,@R0 ；读取操作数的阶码
 RLC A ；取数符
 MOV A,#1 ；按正数初始化
 JNC SGN1 ；是正数，结束
 MOV A,#0FFH ；是负数，改变标志
SGN1: RET

    （１４）  标号：  ＦＩＮＴ    功能：浮点取整函数

入口条件：操作数在[R0]中。
出口信息：结果仍在[R0]中。
影响资源：PSW、A、R2、R3、R4、位1FH    堆栈需求：  ６字节

FINT: LCALL MVR0 ；将[R0]传送到第一工作区中
 LCALL INT ；在工作寄存器中完成取整运算
 LJMP MOV0 ；将结果传回到[R0]中
INT: MOV A,R3
 ORL A,R4
 JNZ INTA
 CLR 1FH ；尾数为零，阶码也清零，结束取整
 MOV R2,#41H
 RET  
INTA: MOV A,R2
 JZ INTB ；阶码为零否？
 JB ACC.7,INTB；阶符为负否？
 CLR C
 SUBB A,#10H ；阶码小于１６否？
 JC INTD
 RET  ；阶码大于１６，已经是整数
INTB: CLR A   ；绝对值小于一，取整后正数为零，负数为负一

MOV R4,A
 MOV C,1FH
 RRC A
 MOV R3,A
 RL A
 MOV R2,A
 JNZ INTC
 MOV R2,#41H
INTC: RET
INTD: CLR F0 ；舍尾标志初始化
INTE: CLR C
 LCALL RR1 ；右规一次
 ORL C,F0 ；记忆舍尾情况
 MOV F0,C
 CJNE R2,#10H,INTE；阶码达到１６（尾数完全为整数）否？
 JNB F0,INTF ；舍去部分为零否？
 JNB 1FH,INTF；操作数为正数否？
 INC R4 ；对于带小数的负数，向下取整
 MOV A,R4
 JNZ INTF
 INC R3
INTF: LJMP RLN ；将结果规格化

    （１５）  标号：  ＦＲＣＰ    功能：浮点倒数函数

入口条件：操作数在[R0]中。
出口信息：OV=0时，结果仍在[R0]中，OV=1时，溢出。
影响资源：PSW、A、B、R2～R7、位1EH、1FH   堆栈需求： ５字节

FRCP: MOV A,@R0
 MOV C,ACC.7
 MOV 1FH,C ；保存数符
 MOV C,ACC.6 ；绝对值传送到第二工作区
 MOV ACC.7,C
 MOV R5,A
 INC R0
 MOV A,@R0
 MOV R6,A
 INC R0
 MOV A,@R0
 MOV R7,A
 DEC R0
 DEC R0
 ORL A,R6
 JNZ RCP
 SETB OV ；零不能求倒数，设立溢出标志
 RET  
RCP: MOV A,R6
 JB ACC.7,RCP2；操作数格式化否？

CLR C ；格式化之
 MOV A,R7
 RLC A
 MOV R7,A
 MOV A,R6
 RLC A
 MOV R6,A
 DEC R5
 SJMP RCP
RCP2: MOV R2,#1 ；将数值１．００传送到第一工作区
 MOV R3,#80H
 MOV R4,#0
 LCALL DIV3 ；调用工作区浮点除法，求得倒数
 LJMP MOV0 ；回传结果

    （１６）  标号：  ＦＳＱＵ    功能：浮点数平方

入口条件：操作数在[R0]中。
出口信息：OV=0时，平方值仍然在[R0]中，OV=1时溢出。
影响资源：PSW、A、B、R2～R7、位1EH、1FH      堆栈需求：  ９字节


FSQU: MOV A,R0 ；将操作数
 XCH A,R1 ；同时作为乘数
 PUSH ACC ；保存Ｒ１指针
 LCALL FMUL ；进行乘法运算
 POP ACC
 MOV R1,A ；恢复Ｒ１指针
 RET

    （１７）  标号：  ＦＳＱＲ    功能：浮点数开平方（快速逼近算法）

入口条件：操作数在[R0]中。
出口信息：OV=0时，平方根仍在[R0]中，OV=1时，负数开平方出错。
影响资源：PSW、A、B、R2～R7      堆栈需求：  ２字节

FSQR: MOV A,@R0
 JNB ACC.7,SQR
 SETB OV ；负数开平方，出错
 RET
SQR: INC R0

INC R0
 MOV A,@R0
 DEC R0
 ORL A,@R0
 DEC R0
 JNZ SQ
 MOV @R0,#41H；尾数为零，不必运算
 CLR OV
 RET
SQ: MOV A,@R0
 MOV C,ACC.6 ；将阶码扩展成８ｂｉｔ补码
 MOV ACC.7,C
 INC A ；加一
 CLR C
 RRC A ；除二
 MOV @R0,A ；得到平方根的阶码，回存之
 INC R0 ；指向被开方数尾数的高字节
 JC SQR0 ；原被开方数的阶码是奇数吗？
 MOV A,@R0 ；是奇数，尾数右规一次
 RRC A
 MOV @R0,A
 INC R0
 MOV A,@R0
 RRC A
 MOV @R0,A

DEC R0
SQR0: MOV A,@R0
 JZ SQR9 ；尾数为零，不必运算
 MOV R2,A ；将尾数传送到R2R3中
 INC R0
 MOV A,@R0
 MOV R3,A
 MOV A,R2 ；快速开方，参阅定点子程序说明
 ADD A,#57H
 JC SQR2
 ADD A,#45H
 JC SQR1
 ADD A,#24H
 MOV B,#0E3H
 MOV R4,#80H
 SJMP SQR3
SQR1: MOV B,#0B2H
 MOV R4,#0A0H
 SJMP SQR3
SQR2: MOV B,#8DH
 MOV R4,#0D0H
SQR3: MUL AB
 MOV A,B
 ADD A,R4
 MOV R4,A
 MOV B,A
 MUL AB
 XCH A,R3
 CLR C

SUBB A,R3
 MOV R3,A
 MOV A,B
 XCH A,R2
 SUBB A,R2
 MOV R2,A
SQR4: SETB C
 MOV A,R4
 RLC A
 MOV R6,A
 CLR A
 RLC A
 MOV R5,A
 MOV A,R3
 SUBB A,R6
 MOV B,A
 MOV A,R2
 SUBB A,R5
 JC SQR5
 INC R4
 MOV R2,A
 MOV R3,B
 SJMP SQR4
SQR5: MOV A,R4
 XCH A,R2
 RRC A
 MOV F0,C
 MOV A,R3
 MOV R5,A
 MOV R4,#8
SQR6: CLR C
 MOV A,R3
 RLC A
 MOV R3,A
 CLR C
 MOV A,R5
 SUBB A,R2
 JB F0,SQR7

JC SQR8
SQR7: MOV R5,A
 INC R3
SQR8: CLR C
 MOV A,R5
 RLC A
 MOV R5,A
 MOV F0,C
 DJNZ R4,SQR6
 MOV A,R3 ；将平方根的尾数回传到[R0]中
 MOV @R0,A
 DEC R0
 MOV A,R2
 MOV @R0,A
SQR9: DEC R0 ；数据指针回归原位
 CLR OV ；开方结果有效
 RET  
             
    （１８）  标号：  ＦＰＬＮ    功能：浮点数多项式计算

入口条件：自变量在[R0]中，多项式系数在调用指令之后，以４０Ｈ结束。
出口信息：OV=0时，结果仍在[R0]中，OV=1时，溢出。

影响资源：DPTR、PSW、A、B、R2～R7、位1EH、1FH    堆栈需求： ４字节

FPLN: POP DPH ；取出多项式系数存放地址
 POP DPL
 XCH A,R0    ；R0、R1交换角色，自变量在[R1]中
 XCH A,R1
 XCH A,R0
 CLR A ；清第一工作区
 MOV R2,A
 MOV R3,A
 MOV R4,A
 CLR 1FH
PLN1: CLR A ；读取一个系数，并装入第二工作区
 MOVC A,@A+DPTR
 MOV C,ACC.7
 MOV 1EH,C
 MOV C,ACC.6
 MOV ACC.7,C
 MOV R5,A
 INC DPTR
 CLR A
 MOVC A,@A+DPTR
 MOV R6,A
 INC DPTR

CLR A
 MOVC A,@A+DPTR
 MOV R7,A
 INC DPTR ；指向下一个系数
 MOV C,1EH ；比较两个数符
 RRC A
 XRL A,23H
 RLC A
 MOV F0,C ；保存比较结果
 LCALL AS1 ；进行代数加法运算
 CLR A ；读取下一个系数的第一个字节
 MOVC A,@A+DPTR
 CJNE A,#40H,PLN2；是结束标志吗？
 XCH A,R0 ；运算结束，恢复R0、R1原来的角色
 XCH A,R1
 XCH A,R0
 LCALL MOV0 ；将结果回传到[R0]中
 CLR A
 INC DPTR
 JMP @A+DPTR ；返回主程序
PLN2: MOV A,@R1 ；比较自变量和中间结果的符号

XRL A,23H
 RLC A
 MOV 1FH,C ；保存比较结果
 LCALL MUL0 ；进行乘法运算
 SJMP PLN1 ；继续下一项运算

    （１９）  标号：  ＦＬＯＧ    功能：以１０为底的浮点对数函数

入口条件：操作数在[R0]中。
出口信息：OV=0时，结果仍在[R0]中，OV=1时，负数或零求对数出错。
影响资源：DPTR、PSW、A、B、R2～R7、位1EH、1FH  堆栈需求：９字节

FLOG: LCALL FLN ；先以ｅ为底求对数
 JNB OV,LOG
 RET  ；如溢出则停止计算
LOG: MOV R5,#0FFH；系数０．４３４３０（１／Ｌｎ１０）

MOV R6,#0DEH
 MOV R7,#5CH
 LCALL MUL1 ；通过相乘来换底
 LJMP MOV0 ；传回结果

    （２０）  标号：  ＦＬＮ    功能：以ｅ为底的浮点对数函数

入口条件：操作数在[R0]中。
出口信息：OV=0时，结果仍在[R0]中，OV=1时，负数或零求对数出错。
影响资源：DPTR、PSW、A、B、R2～R7、位1EH、1FH  堆栈需求： ７字节

FLN: LCALL MVR0 ；将[R0]传送到第一工作区
 JB 1FH,LNOV；负数或零求对数，出错
 MOV A,R3
 ORL A,R4
 JNZ LN0
LNOV: SETB OV
 RET
LN0: CLR C

LCALL RL1 ；左规一次
 CLR A
 XCH A,R2 ；保存原阶码，清零工作区的阶码
 PUSH ACC
 LCALL RLN ；规格化
 LCALL MOV0 ；回传
 LCALL FPLN ；用多项式计算尾数的对数
 DB 7BH,0F4H,30H；０．０２９８０８
 DB 0FEH,85H,13H；－０．１２９９６
 DB 7FH,91H,51H；０．２８３８２
 DB 0FFH,0FAH,0BAH；－０．４８９７
 DB 0,0FFH,0CAH；０．９９９１８
 DB 70H,0C0H,0；１．１４４２×１０-5
 DB 40H ；结束
 POP ACC ；取出原阶码
 JNZ LN1
 RET  ；如为零，则结束

LN1: CLR 1EH ；清第二区数符
 MOV C,ACC.7
 MOV F0,C ；保存阶符
 JNC LN2
 CPL A ；当阶码为负时，求其绝对值
 INC A
LN2: MOV R2,A ；阶码的绝对值乘以０．６９３１５
 MOV B,#72H
 MUL AB
 XCH A,R2
 MOV R7,B
 MOV B,#0B1H
 MUL AB
 ADD A,R7
 MOV R7,A ；乘积的尾数在R6R7R2中
 CLR A
 ADDC A,B
 MOV R6,A
 MOV R5,#8 ；乘积的阶码初始化（整数部分为一字节）
LN3: JB ACC.7,LN4；乘积格式化
 MOV A,R2
 RLC A
 MOV R2,A
 MOV A,R7

RLC A
 MOV R7,A
 MOV A,R6
 RLC A
 MOV R6,A
 DEC R5
 SJMP LN3
LN4: MOV C,F0 ；取出阶符，作为乘积的数符
 MOV ACC.7,C
 LJMP ASN ；与尾数的对数合并，得原操作数的对数

    （２１）  标号：  ＦＥ１０    功能：以１０为底的浮点指数函数

入口条件：操作数在[R0]中。
出口信息：OV=0时，结果仍在[R0]中，OV=1时，溢出。
影响资源：DPTR、PSW、A、B、R2～R7、位1EH、1FH  堆栈需求：６字节

FE10: MOV R5,#2 ；加权系数为３．３２１９（Ｌｏｇ2１０）
 MOV R6,#0D4H

MOV R7,#9AH
 SJMP EXP ；先进行加权运算，后以２为底统一求幂

    （２２）  标号：  ＦＥＸＰ    功能：以ｅ为底的浮点指数函数

入口条件：操作数在[R0]中。
出口信息：OV=0时，结果仍在[R0]中，OV=1时，溢出。
影响资源：DPTR、PSW、A、B、R2～R7、位1EH、1FH  堆栈需求：６字节

FEXP: MOV R5,#1 ；加权系数为１．４４２７２（Ｌｎｇ2ｅ）
 MOV R6,#0B8H
 MOV R7,#0ABH
EXP: CLR 1EH ；加权系数为正数
 LCALL MVR0 ；将[R0]传送到第一工作区
 LCALL MUL1 ；进行加权运算

SJMP E20 ；以２为底统一求幂

    （２３）  标号：  ＦＥ２    功能：以２为底的浮点指数函数

入口条件：操作数在[R0]中。
出口信息：OV=0时，结果仍在[R0]中，OV=1时，溢出。
影响资源：DPTR、PSW、A、B、R2～R7、位1EH、1FH  堆栈需求：６字节

FE2: LCALL MVR0 ；将[R0]传送到第一工作区
E20: MOV A,R3
 ORL A,R4
 JZ EXP1 ；尾数为零
 MOV A,R2
 JB ACC.7,EXP2；阶符为负？
 SETB C
 SUBB A,#6 ；阶码大于６否？
 JC EXP2
 JB 1FH,EXP0；数符为负否？
 MOV @R0,#3FH；正指数过大，幂溢出

INC R0
 MOV @R0,#0FFH
 INC R0
 MOV @R0,#0FFH
 DEC R0
 DEC R0
 SETB OV
 RET
EXP0: MOV @R0,#41H；负指数过大，幂下溢，清零处理
 CLR A
 INC R0
 MOV @R0,A
 INC R0
 MOV @R0,A
 DEC R0
 DEC R0
 CLR OV
 RET  
EXP1: MOV @R0,#1 ；指数为零，幂为１．００
 INC R0
 MOV @R0,#80H
 INC R0
 MOV @R0,#0
 DEC R0
 DEC R0
 CLR OV
 RET  
EXP2: MOV A,R2 ；将指数复制到第二工作区
 MOV R5,A
 MOV A,R3