AVR汇编百例 - 浮点程序库

.ORG $A00

EXCH: MOV R5,R8 ;两浮点数交换子程序 MOV R8,R12 MOV R12,R5

EXCH1: MOV R5,R9 ;尾数交换 MOV R9,R13 MOV R13,R5

MOV R5,R10 ;双字节交换 MOV R10,R14 MOV R14,R5 MOV R5,R11 MOV R11,R15 MOV R15,R5 RET

DP: ANDI R16,$7F ;处理积/商数符，计算积/商阶码子程序 SBRC R9,7 SUBI R16,$80 SBRC R13,7

SUBI R16,$80 ;积/商符号放在r16,7

ADD R12,R8 ;移码相加（除数阶码已求补） LDI R17,$80 BRCC DP1

ADD R12,R17 ;移码求和有进位，将和再加上$80,再有进位为溢出 RET

DP1: SUB R12,R17 ;移码求和无进位，将和减去$80,有借位 RET ;或差为0也为溢出 NEG3: COM R15 ;3字节数据求补 COM R14 ;先求反后加1 COM R13

INC3: LDI R17,255

SUB R15,R17 ;以减去-1代替加1 SBC R14,R17 SBC R13,R17 RET

;范例63 ;浮点数比较大小子程序 X1为被减数 X2为减数 FPCP: SBRC R9,7 ;X1为正，跳行 RJMP CP1

SBRC R13,7 ;X2为正，跳行 RJMP CP2 ;X1，X2异号

FPCP1: CP R11,R15 ;X1，X2皆为正，以尾数低位字节，中位字节，高位字节和 CPC R10,R14 ;阶码的顺序（按无符号数）进行比较

CPC R9,R13 ;不等，阶码大者浮点数值也大；只有阶码和尾数对应相等， CPC R8,R12 ;两浮点数才相等

RET ;比较结果：Z=1时X1=X2，否则C=0时X1>X2,C=1时X1RJMP CP3 ;两负数比较,转

CP2: CP R13,R9 ;正数与负数比较,只比较尾数高位字节即可 RET

CP3: CP R15,R11 ;X1,X2皆为负，以尾数低位字节，中位字节，高位字节和 CPC R14,R10 ;阶码的顺序（按无符号数）进行比较

CPC R13,R9 ;但要将X1、X2交换位置后按正数比较过程进行 CPC R12,R8

CP4: RET ;比较结果：Z=1时X1=X2，否则C=0时X1>X2,C=1时X1;范例64

FPSU: LDI R17,$80 ;浮点减法子程序

SUB R13,R17 ;减数数符求反后作为加数 FPAD: TST R8 ;浮点加法子程序

BREQ DON1 ;被加数为0，加数为和 TST R12

BRNE FPLAD ;加数为0，取被加数为和 SAV0: MOV R12,R8 ;传送被加数取代加数 MOV R13,R9 MOV R14,R10

MOV R15,R11 DON1: RET

FPLAD: ANDI R16,$3f ;清除被加数，加数数符 SBRC R9,7

ORI R16,$80 ;被加数数符取到(R16,7) SBRC R13,7

ORI R16,$40 ;加数数符取到(R16,6) LDI R17,$80 OR R9,R17

OR R13,R17 ;恢复尾数最高位 MOV R17,R12

SUB R17,R8 ;计算阶差

BREQ GOON ;两阶相等，转 BRCC NX3

NEG R17 ;不够减求补 CPI R17,24

BRCC EXAD ;|阶差|>24,取被加数为和 NX2A: LSR R13 ROR R14 ROR R15

DEC R17

BRNE NX2A ;加数阶小，右移加数对阶 MOV R12,R8 ;取被加数阶为和之阶 BRCC GOON

RCALL INC3 ;舍入移出位 RJMP GOON NX3: CPI R17,24

BRCC COM1 ;阶差>24,取加数为和 LOOP: LSR R9 ROR R10 ROR R11 DEC R17

BRNE LOOP ;加数阶大，右移被加数对阶 BRCC GOON

RCALL INC3A ;舍入移出位 GOON: SBRC R16,6 SUBI R16,$80

SBRS R16,7 ;判别两数是否同号 RJMP SAMS ;同号转

SUB R15,R11 ;异号，执行减法，加数为被减数 SBC R14,R10 SBC R13,R9

BRCC NOM ;够减转

SUBI R16,$40 ;否则被减数数符求反为和之数符 RCALL NEG3 ;并将差求补 NOM: MOV R17,R13 OR R17,R14 OR R17,R15

BREQ DON0 ;差为0转 NMLOP: SBRC R13,7 RJMP COM1 LSL R15 ROL R14

ROL R13 DEC R12

BRNE NMLOP ;规格化

OV1: SEV ;阶码变为0，下溢（可取为0，不算溢出） RET

SAMS: ADD R15,R11 ADC R14,R10

ADC R13,R9 ;两数同号，执行加法 BRCC COM1 ROR R13 ROR R14

ROR R15

INC R12 ;有进位时右规1次($7F+1=$80溢出) BREQ OV1 ;阶码增1后变为0为上溢 BRNC COM1 RCALL INC3 COM1: CLV

SBRC R16,6 RET

COMA: LDI R17,$7F

AND R13,R17 ;正数数符为0 DON: RET

EXAD: RCALL SAV0 ;取被加数为和 SBRS R16,7

RJMP COMA ;配置数符 RET

DON0: CLR R12 ;浮点数为0 RET

;范例65 ;浮点乘法子程序 FPMU: TST R8

BREQ M0 ;被乘数为0，积为0 TST R12

BRNE M1 ;乘数为0，积也为0 M0: RJMP G0

M1: RCALL DP ;处理积符号，计算积之阶码 BRCS OV2

BREQ OV2 ;判断溢出 LDI R17,$80 OR R9,R17

OR R13,R17 ;恢复尾数最高位 MOV R5,R13 MOV R6,R14

MOV R7,R15 ;乘数转入R5，R6，R7 LDI R17,25 ;设右移部分积次数 CLR R13 CLR R14

CLR R15 ;r13r14r15清除,存放积 CLC

LOOP1: BRCC M2 ADD R15,R11 ADC R14,R10

ADC R13,R9 ;乘数右移移出位为1,被乘数加入部分积1次 M2: ROR R13

ROR R14 ROR R15 ROR R5 ROR R6

ROR R7 ;部分积连同乘数右移1位 DEC R17

BRNE LOOP1 ;尾数相乘计算完成？ SBRC R13,7

RJMP M3 ;乘积最高位为1 转 ROL R5 ROL R15 ROL R14

ROL R13 ;乘积最高位为0，高4位字节左移1位 SBRS R5,7 RJMP M5

RCALL INC3 ;末位字节舍入 BRNE M5

SEC ;舍入后R13变为0 ROR R13 ;将其改为$80(即0.5) RJMP COM2

M5: DEC R12 ;舍入后R13不为0 BRNE COM2 ;阶码减1

OV2: SEV ;变为0为溢出 RET

M3: SBRC R5,7

RCALL INC3 ;乘积低3位字节舍入 COM2: LDI R17,$7F SBRS R16,7

AND R13,R17 ;正数将符号位请除 DON2: CLV RET

;范例66

FPDI: TST R12 ;浮点除法子程序 BREQ OV3 ;除数为0,溢出 TST R8 BRNE D1

RJMP G0 ;被除数为0,商为0

D1: NEG R12 ;除数阶码求补，以加补码代替减原码 RCALL DP ;处理商符号，计算商之阶码 BRCS OV3

BREQ OV3 ;判断溢出 LDI R17,$80

OR R9,R17

OR R13,R17 ;恢复尾数最高位

FPD3: LDI R17,25 ;左移相减试商25次，最后1次舍入 SUB R11,R15 SBC R10,R14 SBC R9,R13

BRCS D2 ;第一次尾数相减试商 INC R12 ;够减，商阶增1 SEC

BRNE D3 ;商阶增1后不为0，转计商；否则为溢出 OV3: SEV RET

D2: ADD R11,R15 ADC R10,R14

ADC R9,R13 ;不够减则恢复被除数 LOOP2: LSL R11 ROL R10

ROL R9 ;被除数算术左移

BRCS D4 ;进位位为1，够减，本位商1 SUB R11,R15 SBC R10,R14

SBC R9,R13 ;否则相减试商 BRCS D2A SEC

RJMP D3 ;够减，本位商1

D2A: ADD R11,R15 ;不够减，恢复被除数 ADC R10,R14 ADC R9,R13

CLC ;本位商0 RJMP D3

D4: SUB R11,R15 SBC R10,R14

SBC R9,R13 ;被除数减去除数 D3: DEC R17

BRNE D5 ;除法未完成，循环(1-1=0,不溢出) MOV R13,R5 MOV R14,R6

MOV R15,R7 ;取回商 BRCC COM3

RCALL INC3 ;第25位商舍入（$800000-$FFFFFF不溢出，故INC3不会溢出！） COM3: LDI R17,$7F SBRS R16,7

AND R13,R17 ;配置商数符 DON3: RET

D5: ROL R7 ;在R5，R6，R7中记商（不必预先清除） ROL R6

ROL R5 ;商数左移1位并记商

RJMP LOOP2

;范例67

FPSQ: ANDI R16,$7F ;模拟手算开平方子程序 SBRC R13,7

ORI R16，$80 ;负数 建虚根标志 FPS0: TST R12

BREQ DON4 ;0的平方根为0 LDI R17,$80

OR R13,R17 ;恢复尾数最高位 LSR R12 ;阶码算术右移1位 BRCC FSQ2

INC R12 ;移出位舍入

RCALL INC3 ;先将位数增1（提前舍入） BRCS FSQ1 ;C=1，不够减 SEC

ROR R13 ;若尾数变为0将其改为0.5($80-->r13) RJMP FSQ2 FSQ1: LSR R13 ROR R14

ROR R15 ;否则将为数算术右移

FSQ2: LDI R17,25 ;开出25位根,末位舍入 MOV R8,R17 LDI R17,$40

ADD R12,R17 ;根恢复为移码 CLR R5 CLR R6

CLR R7 ;根扩展区清除 CLR R9 CLR R10

CLR R11 ;根存储区清除 FSQ3: SUB R13,R17 SBC R7,R11 SBC R6,R10

SBC R5,R9 ;试根 BRCS FSQ3A SEC

RJMP FSQ4 ;够减，本位根1 FSQ3A: ADD R13,R17 ADC R7,R11

ADC R6,R10

ADC R5,R9 ;否则恢复开平方数之尾数 CLC ;本位商0 FSQ4: DEC R8

BRNE FSQ5 ;开出第25位根? FQDON: MOV R13,R9 MOV R14,R10

MOV R15,R11 ;回送根尾数 BRCC COM4

RCALL INC3 ;第25位根舍入 COM4: LDI R17,$7F

AND R13,R17 ;根尾数为正数 DON4: RET

FSQ5: ROL R11 ROL R10

ROL R9 ;根尾数带进位左移,记根

LSL R15 ROL R14 ROL R13 ROL R7 ROL R6 ROL R5 LSL R15 ROL R14 ROL R13 ROL R7 ROL R6

ROL R5 ;开平方数之尾数连同扩展区左移2位 BRCC FSQ3 ;未产生进位，循环

RJMP FQDON ;否则进位为第25位根（不须试,并结束子程序）！

;范例68 ;牛顿迭代开平方子程序 FSQR: TST R12

BREQ SQRT ;0的平方根为0 ANDI R16,$7E SBRC R13,7

ORI R16,$80 ;虚根标志 SBRC R12,0

INC R16 ;阶码为奇数 LDI R17,$7F

AND R13,R17 ;尾数变为正数 LSR R12

LDI R17,$40

ADC R12,R17 ;得到根之移码 PUSH R12 ;暂存 LDI R17,$80 MOV R12,R17 SBRC R16,0

INC R12 ;得到X1的阶码(0.5≤X1<2) RCALL LD1 ;存 X1 LSR R13 ROR R14 ROR R15 LDI R17,$40

SBRS R16,0 ;阶码为奇数时算术右移尾数即得到X1之尾数；否则将其最 ;高位字节加上$40 OR R13,R17 ;得到首次根r0=(1+x1)/2 LDI R17,3

MOV R0,R17 ;迭代3次 FSQLP: RCALL LD2 RCALL GET1 RCALL FPDI RCALL GET2 RCALL FPAD

DEC R12 ;计算r(i+1)=(x1/ri+ri)/2 DEC R0

BRNE FSQLP ;r3的尾数为根之尾数 POP R12 ;取回根之阶码

SQRT: RET ;r16,7=1 为虚数根

;范例69 ;基本运算程序的演示程序 DMST1: .EQU SPL=$3D .EQU SPH=$3E

LDI R16,2 ;high(ramend) OUT SPH,R16

LDI R16,$5F ;low(ramend) OUT SPL,R16

LDS R11,$60 ;r11,7:数符 r11,6 ：阶符 r11,5--0:阶(最大为38) LDS R12,$61 ;r12-r15:尾数 LDS R13,$62 LDS R14,$63

LDS R15,$64 ;尾数共8位BCD码 RCALL DTOB ;转为二进制浮点数 RCALL LD2 ;暂存

LDS R11,$65 ;r11,7:数符 r11,阶符 r11,5--0:阶(最大为38)

LDS R12,$66 ;r12-r15:尾数 LDS R13,$67 LDS R14,$68

LDS R15,$69

RCALL DTOB ;转为二进制浮点数 RCALL GET2 ;取第一操作数

RCALL FPAD ;调基本运算子程序之一(FPSU/FPMU/FPDI) RCALL BTOD ;转回十进制浮点数 DMRET: RJMP DMRET

;范例70 ;辅助子程序

KP2: MOV R8,R12 ;复制第二操作数 MOV R9,R13 MOV R10,R14 MOV R11,R15 RET

LD1: STS $70,R12 ;存浮点数 STS $71,R13 STS $72,R14 SYS $73,R15 RET

LD2: STS $74,R12 ;存浮点数 STS $75,R13 STS $76,R14 STS $77,R15 RET

LD3: STS $78,R12 ;存浮点数 STS $79,R13 STS $7A,R14 STS $7B,R15 RET

GET1: LDS R8,$70 ;取浮点数 LDS R9,$71 LDS R10,$72 LDS R11,$73 RET

GET2: LDS R8,$74 ;取浮点数 LDS R9,$75 LDS R10,$76 LDS R11,$77 RET

GET3: LDS R8,$78 ;取浮点数 LDS R9,$79

LDS R10,$7A LDS R11,$7B RET

INVPI: LDI R17,$86 ;取浮点数180/л MOV R8,R17 LDI R17,$65 MOV R9,R17 LDI R17,$2E MOV R10,R17 LDI R17,$E1 MOV R11,R17 RET

G90: LDI R17,$87 ;取浮点数90 MOV R8,R17 LDI R17,$34 MOV R9,R17 CLR R10 CLR R11 RET

DTOR: RCALL PI18 ;角度化为弧度 RJMP FPMU

RTOD: RCALL INVPI ;弧度化为角度 RJMP FPMU

GHPI: LDI R17,$81 ;取浮点数л/2 MOV R8,R17 LDI R17,$49 MOV R9,R17 LDI R17,$0f MOV R10,R17 LDI R17,$DB MOV R11,R17 RET

G01: LDI R17,$7D ;取浮点数0.1 MOV R8,R17 LDI R17,$4C MOV R9,R17 LDI R17,$CC MOV R10,R17 LDI R17,$CD MOV R11,R17 RET

G1: LDI R17,$81 ;取浮点数1 MOV R8,R17 CLR R9

CLR R10 CLR R11 RET

PI18: LDI R17,$7B ;取浮点数л/180 MOV R8,R17 LDI R17,$0E MOV R9,R17 LDI R17,$FA MOV R10,R17 LDI R17,$35 MOV R11,R17 RET

GINT: LDI R17,R12 ;浮点数取整 CPI R17,$81 BRCC GINT1

RCALL G0 ;阶码<$81,结果为0 RJMP KP2

GINT1: ANDI R16,$DD SBRC R13,7

ORI R16,2 ;记数符 CPI R17,$98

BRCC GOVER ;阶码>$97,溢出

RCALL BRK ;分解出整数部分（在R9 R10 R11） SBRS R16,1

RET ;正数返回

NEG3A: COM R11 ;负数求(r9 r10 r11)之补 COM R10 COM R9

INC3A: LDI R17,255 SUBI R11,R17 SBCI R10,R17

SBCI R9,R17 ;求反后加1 RET

GOVER: ORI R16,$20 ;设整数部分超过23位标志 RET

BRK: ANDI R16,$DF ;将正浮点数分解为整数/小数两部分 LDI R17,$80

OR R13,R17 ;恢复尾数最高位 CLR R9 CLR R10 CLR R11

MOV R17,R12 SUBI R17,$80 BREQ BRKRT

BRCS LOOPT

CPI R17,$19 ;整数部分超过24位 BRCC GOVER ;为溢出 LOOP4: LSL R15 ROL R14 ROL R13 ROL R11 ROL R10 ROL R9 DEC R17

BRNE LOOPT ;左移位数为阶码-$80,整数部分进入r9-r11中 BRKRT: RET

LOOPT: LSR R13 ;只有小数部分右移尾数（$80-阶码）位 ROR R14 ROR R15 INC R17

BRNE LOOPT RET

NRML: ANDI R16,$BF ;1字节正整数（在R13中）规格化为浮点数 CLR R14 CLR R15

LDI R12,$88 ;设阶码 RJMP NMLOP

G10: LDI R17,$84 ;取浮点数10 MOV R8,R17 LDI R17,$20 MOV R9,R17 CLR R10 CLR R11 RET

GLN2: LDI R17,$80 ;取浮点数ln2(=0.6931471806) MOV R8,R17 LDI R17,$31 MOV R9,R17 LDI R17,$72 MOV R10,R17 LDI R17,$18 MOV R11,R17 RET

GLN10: LDI R17,$82 ;取浮点数ln10(=2.302585093) MOV R8,R17 LDI R17,$13 MOV R9,R17 LDI R17,$5D

MOV R10,R17 LDI R17,$8E MOV R11,R17 RET

INVX: TST R12 ;计算1/X, X=0时溢出 BRNE INV OV4: SEV RET

INV: RCALL G1 ;取1 RJMP FPDI

;范例71 ;用荷纳法计算多项式值子程序

FPLN1: ORI R16,$10 ;设计算奇函数(lnx,sinx,arcsinx 等)标志 RCALL LD3 ;存X RCALL KP2

RCALL FPMU ;计算X2

RJMP FLN0 ;

FPLN2: ANDI R16,$EF ;设计算偶函数(EXP,COSX等)标志 FLN0: RCALL LD1 ;存T,T=X 或T=X2 POP R30

POP R31 ;系数表数据地址进入Z LSL R30

ROL R31 ;由按字取数变为按字节取数 LPM ;r0<--(z)取阶码 MOV R8,R0

ADIW R30,1 ;指针增1

LPM ;取尾数高位字节 MOV R9,R0

ADIW R30,1 ;z+1

LPM ;取尾数中位字节 MOV R10,R0

ADIW R30,1 ;z+1

LPM ;取尾数低位字节

MOV R11,R0 ;取浮点数到r8 r9 r10&r11 ADIW R30,1 ;z+1

PLN: RCALL M1 ;计算(....((An*T+A(n-1))*T+A(n-2))*T+....+Ai)*T LPM

MOV R8,R0 ADIW R30,1 LPM

MOV R9,R0 ADIW R30,1 LPM

MOV R10,R0 ADIW R30,1 LPM

MOV R11,R0 ;取A(i-1) ADIW R30,1

RCALL FPLAD ;计算(....((An*T+A(n-1))*T+A(n-2))*T+....+Ai)*T+A(i-1) LPM

RCALL GET1 DEC R0

BRNE PLN ;1为停止符号；否则继续计算 PEND: SBRS R16,4 RJMP REND

RCALL GET3 ;奇函数 取出自变量

RCALL M1 ;自变量乘以计算结果才是函数值 REND: LSR R31

ROR R30 ;Z指针折半后

ADIW R30,1 ；增1为后继指令地址 IJMP ;转到该地址去执行

;范例72

LNX: TST R12 ;对数函数子程序 BREQ OV5 SBRS R13,7 RJMP LN1

OV5: SEV ;求负数或0的对数为错误 RET

LN1: ANDI R16,$7E ;R16,7:（T-1）/(T+1)或(2T-1)/(2T+1)之符号 R16，0:p之符号 ; m

MOV R0,R12 ;设X=2 *T, 则LnX=m*Ln2+LnT,存入p=m LDI R17,$F3 CP R15,R17 LDI R17,$04 CPC R14,R17 LDI R17,$35

CPC R13,R17 ; _

BRCC LN5 ;T>√2/2时跳转

DEC R0 ;取p=m-1 LnX=(m-1)*Ln2+LN(2T) MOV R17,R15 OR R17,R14 OR R17,R13 MOV R12,R17

BREQ LN5A ;2T-1=0 只须计算(m-1)Ln2 RCALL KP2 ;R12 NOUSED!

LSL R9 ROL R10

ROL R11 ;(2T-1) LSR R13 ROR R14 ROR R15 LDI R17,$80

OR R13,R17 ;2T+1 LDI R17,$7E

MOV R12,R17 ;取1/(2T+1)的阶码 RJMP LNTLP

LN5: ORI R16,$80 ;(T-1)为负，数符位改为1 RCALL KP2 RCALL NEG3A LDI R17,$80

ADD R9,R17 ;计算(T-1) LSR R13 ROR R14 ROR R15 LDI R17,$C0 OR R13,R17 LDI R17,$7F

MOV R12,R17 ;取1/(T+1)的阶码 LNTLP: LSL R11 ROL R10

ROL R9 ;(2T-1)或(T-1)规格化

DEC R12 ;调整(2T-1)/(2T+1))或(T-1)/(T+1)的阶码 SBRS R9,7 RJMP LNTLP

RCALL FPD3 ;计算(2T-1)/(2T+1)或(T-1)/(T+1) 位r16,7为商之数符 PUSH R0

RCALL FPLN1 ;计算LnT或Ln(2T)

.DB $7E,$12,$49,$25 ;0.14285714 ;er.total<0.000000029! .DB $7E,$4C,$CC,$CD ;0.2

.DB $7F,$2A,$AA,$AB ;0.33333333 .DB $81,$00,$00,$00 ;1 .DB $01,$00 ;结束符 INC R12 POP R0

LN5A: LDI R17,$80 ADD R0,R17

BREQ LN53 ;p=$80结束 BRCS LN51 NEG R0

INC R16 ;p为负数

LN51: RCALL LD1 ;存LnT或Ln(2T) MOV R13,R0

RCALL NRML ;|P|规格化 RCALL GLN2 ;取ln2

RCALL FPMU ;计算|p|*ln2

RCALL GET1 ;取LnT或Ln(2T) SBRS R16,0 RJMP LN52

RCALL FPSU ;p<0 计算lnT-|p|*ln2或Ln(2T)-|p|*ln2 RET

LN52: RCALL FPAD ;p>0 计算lnT+|p|*ln2或Ln(2T)+|p|*ln2 LN53: RET

;范例73 ;对数衍生函数子程序 LGX: RCALL LNX ;计算lnx RCALL GLN10 ;取ln10 RCALL EXCH

RJMP FPDI ;转计算lgx=lnx/ln10 LGAX: RCALL LD2 ;存a RCALL EXCH

RCALL LNX ;计算lnx RCALL GET2 ;取a RCALL LD2 ;存lnx RCALL EXCH

RCALL LNX ;计算lna

RCALL GET2 ;转计算logax=lnx/lna

RJMP FPDI

;范例74

EXP: MOV R17,R12 ;指数函数子程序 CPI R17,$68 ;X之阶<$68 E1: BRCC E2 RCALL G0

ROR R12 ;(R12)=$80 INC R12 ;取exp=1 RET

E2: ANDI R16,$3F ;r16,6:数符 SBRC R13,7

ORI R16,$40 ;负数 LDI R17,$7F

AND R13,R17 ;取正（取|X|）

LDI R17,$33 CP R15,R17 LDI R17,$0F CPC R14,R17 LDI R17,$30 CPC R13,R17 LDI R17,$87

CPC R12,R17 ;|X|与88.02969 比较 BRCS E3 ;|X|<88.02969 转 SBRS R16,6 RJMP OV6

G0: CLR R12 ;若x<-88.02969 CLR R13

CLR R14 ;Exp=0 CLR R15 RET

OV6: SEV ;x>88.02969,Exp溢出 RET

E3: CLR R0 ;X整数部分予清除 LDI R17,$81 MOV R8,R17 LDI R17,$38 MOV R9,R17 LDI R17,$AA MOV R10,R17 LDI R17,$3B

MOV R11,R17 ;取log2e(=1/ln2) RCALL FPMU ;计算X/ln2 LDI R17,$80 SBRC R16,6 OR R13,R17 MOV R17,R12 CPI R17,$81

BRCS E6 ;X/ln2整数部分为0 转

RCALL BRK ;否则分解该数为整数I（在R11），小数F两部分 LDI R17,$80 MOV R12,R17

RCALL NOM ;小数部分规格化为浮点数 SBRC R16,6

NEG R11 ;整数部分求补 MOV R0,R11 ; E6: PUSH R0

RCALL FPLN2 ;计算EXP(F*ln2)

.DB $69,$5A,$92,$9F ;0.10178086 E-6 ;er.total<0.000000024

.DB $6D,$31,$60,$11 ;0.13215487 E-5 .DB $70,$7F,$E5,$FE ;0.15252734 E-4 .DB $74,$21,$84,$89 ;0.15403530 E-3 .DB $77,$2E,$C3,$FF ;0.13333558 E-2 .DB $7A,$1D,$95,$5B ;0.96181291 E-2 .DB $7C,$63,$58,$47 ;0.55504109 E-1 .DB $7E,$75,$FD,$F0 ;0.24022651 .DB $80,$31,$72,$18 ;0.69314718 .DB $81,$00,$00,$00 ;1 .DB $01,$00 ;结束符 POP R0

ADD R12,R0 ;整数部分I 加入阶码中 RET

;范例75 ;指数衍生函数子程序 DXP: RCALL GLN10 ;取ln10

RJMP EXP0 ;转计算EXP(X*ln10) AXP: RCALL LD2 ;存X RCALL EXCH

RCALL LNX ;计算lna RCALL GET2 ;取出x EXP0: RCALL FPMU

RJMP EXP ;转计算EXP(X*lna)

;范例76 ;双曲函数和反双曲函数子程序 SHX: RCALL SUB11 ;计算双曲正弦 RCALL FPSU BRNE NX48 RET

CHX: RCALL SUB11 ;计算双曲余弦 RCALL FPAD NX48: DEC R12 RET

SUB11: RCALL EXP RCALL LD2 RCALL INVX RJMP GET2

ASHX: RCALL SUB2 ;计算反双曲正弦 RCALL FPAD ASH: RCALL FPSQ RCALL GET2 RCALL FPAD

RJMP LNX

ACHX: RCALL SUB2 ;计算反双曲余弦 RCALL EXCH RCALL FPSU RJMP ASH

SUB2: RCALL LD2 ;存X RCALL KP2

RCALL FPMU ;得到X2 RJMP G1 ;取浮点数1

; 范例77 ;正弦函数子程序

SINX: RCALL RTOD ;弧度化为角度 SINX1: CLR R16 ;X1为角度 SBRC R13,7

INC R16 ;存数符 LDI R17,$7F ;X1-->|X1| AND R13,R17 NX30: RCALL G90 INC R8

INC R8 ;取360°

RCALL FPCP1 ;|X1|与360°比较 BREQ GE0 ;相等,转出

BRCC NX31 ;|X1|<360° 转出 RCALL EXCH

RCALL FPSU ;否则|X1|-360°-->|X1| RJMP NX30 ;循环 NX31: DEC R8

RCALL FPCP1 ;|X1|与180°比较 BREQ GE0 ;相等,转出

BRCC NX32 ;|X1|<180°，转 RCALL EXCH

RCALL FPSU ;否则|X1|-180°-->|X1| INC R16 ;将数符求反 NX32: RCALL G90

RCALL FPCP1 ;|X1|与90°比较 BRCC NX36 INC R8

RCALL FPSU ;|X1|>90°,取180°-|x1|-->|x1| RJMP NX36

GE0: RJMP G0 ;|X1|=0 则sinX=0 NX36: RCALL DTOR ;变回弧度X MOV R17,R12

CPI R17,$79 ;阶码<$79,sinX=X

BRCS PP2

RCALL FPLN1 ;计算sin|X|

.DB $60,$30,$92,$32 ; 0.16059044 E-9 er.total<0.0000000071 .DB $67,$D7,$32,$2A ;-0.25052108 E-7 .DB $6E,$38,$EF,$1C ; 0.27557319 E-5 .DB $74,$D0,$0D,$01 ;-0.19841270 E-3 .DB $7A,$08,$88,$88 ; 0.83333333 E-2 .DB $7E,$AA,$AA,$AA ;-0.16666667 .DB $81,$00,$00,$00 ;1 .DB $01,$00 ;结束符 PP2: LDI R17,$80 SBRC R16,0

PP3: OR R13,R17 ;配置数符 DON6: RET

;范例78 ;衍生三角函数子程序

CTNX: RCALL RTOD ;弧度化为角度 CTNX1: RCALL TANX1 ;计算tgX RJMP INVX ;取倒数为ctgX

TANX: RCALL RTOD ;弧度化为角度 TANX1: RCALL LD2 ;存X RCALL SINX1 ;计算sinX RCALL GET2 ;取X RCALL LD2 ;存sinX RCALL EXCH

RCALL COSX1 ;计算cosX BRNE NX39 OV7: SEV

RET ;cosX=0,溢出

NX39: RCALL GET2 ;取sinX RJMP FPDI ;tgX=sinX/cosX

COSX: RCALL RTOD ;弧度化为角度 COSX1: RCALL G90 ;取浮点数90° RCALL FPSU

RJMP SINX1 ;cosX=sin(90-X)

;范例79 ;反正弦函数子程序 ASINX: MOV R17,R12 CPI R17,$78

BRCS DON6 ;X阶码<$78,acrsinX=X

ANDI R16,8 ;清除数符和|X|>0.5标志，保留计算acosx标志（R16，3） SBRC R13,7

INC R16 ;记数符 LDI R17,$7F

AND R13,R17 ;取绝对值 X-->|X| RCALL G1 RCALL FPCP1 BREQ AA BRCC AA1

OV8: SEV ;|X>1,溢出 RET

AA: RCALL GHPI RCALL EXCH

RJMP PP2 ;|X|=1,arcsinX=±л/2 AA1: MOV R17,R12 CPI R17,$80

BRNE AS1 ;|X|<0.5,y=|x| MOV R17,R13 OR R17,R14 OR R17,R15

BREQ AS1 ;X=0.5,y=|x| ORI R16,$20 ;X>0.5，建标 RCALL NEG3 LDI R17,$80 ADD R13,R17 LDI R17,$7F

MOV R12,R17 ;((1-|x|)/2)方根之阶最大为$7F NRMLP: LSL R15 ROL R14 ROL R13 DEC R12 SBRS R13,7

RJMP NRMLP ; __________ RCALL FPS0 ;√(1-|X|)/2-->y AS1: RCALL FPLN1 ;计算arcsiny

.DB $7A,$3D,$43,$C4 ;0.11551801 E-1 er. total<0.0000000245 .DB $7A,$64,$CC,$CD ;0.13964844 E-1 .DB $7B,$0E,$27,$62 ;0.17352764 E-1 .DB $7B,$37,$45,$D1 ;0.22372159 E-1 .DB $7B,$78,$E3,$8E ;0.30381944 E-1 .DB $7C,$36,$DB,$6E ;0.44642857 E-1 .DB $7D,$19,$99,$9A ;0.075

.DB $7E,$2A,$AA,$AA ;0.16666667 .DB $81,$00,$00,$00 ;1 .DB $01,$00 ;结束符 SBRS R16,5

RJMP PP2 ;|x|≤0.5 转配置数符，有acsin|x|=acsiny INC R12 ;否则取2arcsiny(=arccosx) SBRC R16,3 ;测试计算ARCCOSX的标志 RJMP ACSRT ;有计算ARCCOSX标志，转清除该标志（其余计算在ACOSX子程序中完成）

RCALL GHPI ;否则取л/2

AS2: RCALL FPSU ;|X|>0.5时,arcsin|X|=л/2 -2arcsiny PP20: RJMP PP2 ;转去配置数符

;范例80 ;函数值为弧度的反三角函数子程序 ACOSX: ORI R16,8 ;设计算arccosx标志 RCALL ASINX ;调反正弦函数子程序 RCALL GHPI ;取л/2

SBRC R16,3 ;计算ARCCOS|X|标志未被清除？ RJMP AS3 ;是，转计算arccosx=л/2-arcsinx SBRS R16,0 ;x>0且x>0.5 RJMP ACSRT ;有arccosx=2arcsiny!

INC R8 ;否则取л；即当x<0且|X|>0.5时,有arccosX=л-2arcsiny AS3: RCALL FPSU

ACSRT: ANDI R16,$F7 ;清除计算arccosx标志 RET

ATANX: MOV R17,R12 ;反正切函数子程序 CPI R17,$98 BRCS AT1

RCALL GHPI ;X阶码大于$98,取л/2 RCALL EXCH ROL R9 BRCC AT2 LDI R17,$80

OR R13,R17 ;arctgx=л/2 AT2: RET

AT1: MOV R17,R12

CPI R17,$74 ;X阶码小于$74,arctgX=X BRCS AT2 RCALL KP2

RCALL LD1 ;存X RCALL FPMU RCALL G1

RCALL FPAD ; _______ RCALL FPSQ ;计算 √(1+X2) RCALL GET1 RCALL FPDI

RJMP ASINX ;转计算arctgx=arcsin(X/√<(1+X2)

ACTNX: RCALL ATANX ;反余切函数子程序 RCALL GHPI

RJMP FPSU ;arcctgX=л/2-arctgx

;范例81 ;函数值为角度的反三角函数子程序

ASNX: RCALL ASINX ;反正弦函数子程序,结果以角度表示 RJMP RTOD

ACSX: RCALL ACOSX ;反余弦函数子程序,结果以角度表示 RJMP RTOD

ATNX: RCALL ATANX ;反正切函数子程序,结果以角度表示 RJMP RTOD

ACNX: RCALL ACTNX ;反余切函数子程序,结果以角度表示 RJMP RTOD

;范例82 ;函数计算子程序演示程序 DMST2: LDI R16,2 OUT SPH,R16

LDI R16,$5F ;堆栈指针初始化 OUT SPL,R16

LDS R11,$65 ;取操作数（自变量X） LDS R12,$66 ;r11,7:数符 r11,6：阶符 LDS R13,$67 ;r11,5--0:阶(最大为38) LDS R14,$68

LDS R15,$69 ;r12-r15:十进制尾数（8位BCD码） RCALL DTOB ;翻为二进制浮点数 RCALL LNX ;调函数子程序之一

RCALL BTOD ;将函数值转为十进制浮点数 DMHER: RJMP DMHER

;范例83 ;阶乘子程序

NP: RCALL G1 ;取浮点数1

MOV R17,R12 ;二进制整数N在R12中 CPI R17,2 ;N<2,N！=1 BRCS GG CPI R17,34 BRCS NX59

OV9: SEV ;N>33,溢出 RET

GG: RJMP SAV0 ;取N!=1 NX59: MOV R0,R12 ;存N DEC R0 ;N-1

LDI R17,1

PUSH R17 ;取T=1 并存入 LDI R17,$81 STS $70,R17 CLR R17 STS $71,R17 STS $72,R17

STS $73,R17 ;存储浮点数1 L43: POP R13 ;取T INC R13 ;T+1-->T PUSH R13 ;存T

RCALL NRML ;T规格化

RCALL GET1 ;取阶段阶乘结果 RCALL FPMU ;得到当前T! RCALL LD1 DEC R0

BRNE L43 ;T=N时得到N! POP R0 RET

;范例84 ;长整数(r9,r10,r11,r12)规格化为浮点数 LINOM: BST R9,7 ;数符存于T BRTC LI10

CLR R16 ;负数求补 SUB R16,R12 MOV R12,R16 CLR R16 SBC R16,R11 MOV R11,R16 CLR R16 SBC R16,R10 MOV R10,R16 CLR R16 SUB R16,R9 MOV R9,R16

LI10: LDI R16,$A0 ;取阶32（长整数共32位） LP10: SBRC R9,7

RJMP NX63 ;最高位为1，已规格化 LSL R12 ROL R11 ROL R10

ROL R9 ;否则左规1位 DEC R16 ;阶码减1

CPI R16,$80 BRNE LP10

RJMP G0 ;左规达32次，浮点数为0

;范例85 ;定点十进制数翻为二进制浮点数 DTOB1: RCALL LD1 ;存入十进制小数 RCALL CONV2 ;定点整数十翻二 RCALL GET1 ;取出十进制小数 RCALL LD1 RCALL SAV0

RCALL CONV4 ;定点小数十翻二

RCALL GET1 ;取出二进制定点整数 LDI R16,$98 ;予设阶码 LP11: SBRC R9,7

RJMP NX63 ;最高位为1,已规格化 LSL R15 ROL R14 ROL R13 ROL R12 ROL R11 ROL R10

ROL R9 ;整数和小数部分左移一位 DEC R16 ;阶码减1 CPI R16,$60 BRNE LP11

RET ;得到浮点数0 NX63: MOV R13,R9 MOV R14,R10

MOV R15,R11 ;尾数取到r13-r15 SBRS R12,7 RJMP PP6

RCALL INC3 ;尾数截去部分舍入 BRNE PP6

INC R16 ;尾数变为0将阶码增1 SEC

ROR R13 ;$80-->r13，即将尾数变为0.5 PP6: MOV R12,R16 ;取回阶码 BLD R13,7 ;装入数符(T-->R13,7) RET

;范例86 ;浮点数十翻二

DTOB: ANDI R16,$FC ;r11,7：数符 r11,6：阶符 r11，5--0：阶(最大为38)

SBRC R11,6 ;R12---R15；8BCD码尾数 INC R16 ;阶符存于R16，0 SBRC R11,7

ORI R16,2 ;数符存于R16，1 MOV R17,R11

ANDI R17,$3F ;取阶 MOV R0,R17 ;存于R0 MOV R8,R12 OR R8,R13 OR R8,R14 OR R8,R15

BREQ PP8 ;十进制浮点数尾数为0，取二进制浮点数0 PUSH R16

RCALL CONV4 ;十进制浮点数尾数翻为二进制定点小数 MOV R16,R15 MOV R15,R14 MOV R14,R13

MOV R13,R12 ;二进制定点小数转入r13r14r15r16 LDI R17,$80 ;予设阶码 MOV R12,R17 LP14: SBRC R13,7 RJMP NX67 LSL R16 ROL R15 ROL R14 ROL R13 DEC R12

RJMP LP14 ;二进制定点小数规格化为浮点数 NX67: SBRS R16,7 RJMP NX66

RCALL INC3 ;调整 BRNE NX66 INC R12

SEC ;调整后结果为0将其改为0.5 ROR R13 ;即$80-->r13 NX66: LDI R17,$7F POP R16 SBRS R16,1

AND R13,R17 ;配置数符 SBRS R16,0

RJMP DBL4 ;正阶转 DBL1: LDI R17,$10 SUB R0,R17 BRCS DBL2

RCALL INVDP ;

RCALL FPMU ;阶码减10, X*10ˉ1o -->X RJMP DBL1

DBL2: ADD R0,R17 ;不够减则恢复阶 BREQ PP8

DBL3: RCALL G01 ;取0.1 RCALL FPMU

DEC R0 ;X*0.1-->X,阶减1 BRNE DBl3 RET

DBL4: LDI R17,$10

SUB R0,R17 ;阶减10 BRCS DBL5

RCALL DDP ;X*101o -->X RCALL FPMU RJMP DBL4

DBL5: ADD R0,R17 ;不够减则恢复阶 BREQ PP8

DBL6: RCALL G10 RCALL FPMU

DEC R0 ;X*10-->X,阶减1 BRNE BDL6

PP8: RET

INVDP: LDI R17,$5F ;取浮点数10ˉ1o MOV R8,R17 LDI R17,$5B MOV R9,R17 LDI R17,$E6 MOV R10,R17 LDI R17,$FF MOV R11,R17

RET ;

DDP: LDI R17,$A2 ;取浮点数101o MOV R8,R17 LDI R17,$15 MOV R9,R17 LDI R17,$02 MOV R10,R17 LDI R17,$F9 MOV R11,R17 RET

;范例87 ;浮点数二翻十

BTOD: TST R12

BREQ PP4 ;转取十进制浮点数0

ANDI R16,$FC ;予清十进制浮点数数符及阶符(r16,1&0) CLR R0 ;予清十进制浮点数之阶 SBRC R13,7

ORI R16,2 ;取数符 LDI R17,$7F

AND R13,R17 ;取绝对值 BTA: RCALL DDP

RCALL FPCP1 ;|X|与101o 比较 BREQ BTB

BRCC BTC ;|X|<101o 转 BTB: RCALL INVDP

RCALL FPMU ;|X|*10ˉ1o-->|X| LDI R17,$10

ADD R0,R17 ;十进制浮点数阶加10 RJMP BTA

BTC: RCALL INVDP ;

RCALL FPCP1 ;|X|与10ˉ1o 比较 BREQ BTC1 ;

BRCS BT0 ;|X|>10ˉ1o 转 BTE: RCALL DDP ;

RCALL FPMU ;|X|*101o -->|X| LDI R17,$10

ADD R0,R17 ;十进制浮点数阶加10 ORI R16,1 ;置负阶 RJMP BTC

BTC1: LDI R17,9 ;|X|=10ˉ1o 特别处理 ADD R0,R17 ; -9 ORI R16,1 ;取0.1*10 SJMP BT4 BT0: RCALL G1

RCALL FPCP1 ;|X|与1比较 BREQ BT1

BRCC BT2 ;|X|<1转 BT1: RCALL G01

RCALL FPMU ;|X|*0.1-->|X| INC R0 ;十进制浮点数阶加1 RJMP BT0

BT2: RCALL G01

RCALL FPCP1 ;|X|与0.1比较 BREQ BT4

BRCS BDS ;|X|≤0.1转出 BT3: RCALL G10

RCALL FPMU ;|X|*10--->|X| INC R0 ;十进制浮点数阶加1 ORI R16,1 ;置负阶 RJMP BT2

PP4: RJMP KP2 ;十进制浮点数取为0 BT4: LDI R17,$10 MOV R9,R17 CLR R10 CLR R11

CLR R12 ;十进制浮点数尾数取为0.10000000 BT6: MOV R8,R0 ;取十进制浮点数阶 SBRS R8,3 RJMP BT7 SBRC R8,1

SUBI R8,$FA ;对产生非法BCD调整（加6） BT7: LDI R17,$40 SBRC R16,0

OR R8,R17 ;配置阶符(r8,6) LSL R17 SBRC R16,1

OR R8,R17 ;配置阶浮(r8,7) RET

BDS: RCALL BT6 ;BT6将十进制浮点数阶，阶符和数符配置到R8 LDI R17,$80

OR R13,R17 ;恢复尾数最高位 LDI R17,$98

SUB R17,R12 ;右移次数为($98-阶码)

RJMP CONV31 ;调CONV31子程序完成尾数二翻十，结果在(r9r10r11r12)

;范例88 ;二进制浮点数快速翻为定点十进制数，整数在r9,r10,r11中，小数在r13,r14,r15中

FBTOD: RCALL BRK ;二进制浮点数分解为整数和小数两部分 SBRC R16,5

RET ;整数部分多于24位,溢出 MOV R0,R13 MOV R5,R14

MOV R8,R15 ;小数部分转入R0R5R8

RCALL CONV1 ;定点整数二翻十，结果在R12,R13,R14,R15 RCALL LD1 ;十进制整数-->RAM MOV R15,R8 MOV R14,R5

MOV R13,R0 ;取回二进制小数

RCALL CONV3 ;定点小数二翻十，结果在r9,r10,r11,r12

RCALL EXCH1 ;十进制定点小数转入r13,r14,r15,r12

RCALL GET1 ;取出十进制定点整数r8,r9,r10,r11)/小数在r13,r14,r15,r12 CLR R16 ;清除无用的标志! RET

;范例89

.ORG $E80 ;最小二乘法拟和直线子程序

.EQU NUMB=10 ;取10点,即十对浮点数，按增地址存放Y1,X1,Y2,X2,..Yn,Xn .EQU TABLA=$9000 ;数据表，第一个浮点数为Y1 STRT: LDI R28,$70

CLR R29 ;POINT TO $0070

LP51: ST Y+,R29 ;累加和或暂存区清除（LD1,LD2,LD3,LD4和LD5子程序工作区）

CPI R28,$84 BRNE LP51

LDI R16,NUMB ;取拟合点数 MOV R0,R16

LDI R29,$90

CLR R28 ;参加拟合数据首地址$9000 IN R16,MCUCR,7

SBR R16,$C0 ;片外RAM，选一个读写等待周期 OUT MCUCR,R16

LOOP3: RCALL GETA ;取浮点数Yi 占4字节 即Yi0,Yi1,Yi2,Yi3 RCALL INVX ;计算1/Yi RCALL LD6 ;暂存

RCALL GET1 ;取累加和 n

RCALL FPAD ;1/Yi加入累加和(∑1/Yi是 ∑1/Yi 简写形式，下同) RCALL LD1 ; i=1 RCALL GET6 ;取1/Yi PUSH R28

PUSH R29 ;保护堆栈指针

RCALL GETA ;取浮点数Xi（Xi0,Xi1,Xi2,Xi3）占4字节 POP R29

POP R28 ;恢复堆栈指针，仍指向Xi RCALL FPMU ;计算Xi/Yi RCALL LD7 ;暂存 RCALL GET2

RCALL FPAD ;Xi/Yi加入累加和∑(Xi/Yi) RCALL LD2

RCALL GET7 ;取出Xi/Yi RCALL SAV0 ;

RCALL FPMU ;计算(Xi/Yi)2 RCALL GET3

RCALL FPAD ;(Xi/Yi)2加入累加和∑(Xi/Yi)2 RCALL LD3

RCALL GET6 ;取1/Yi RCALL SAV0 ;

RCALL FPMU ;计算1/Yi2 RCALL LD6 ;暂存 RCALL GET4

RCALL FPAD ;1/Yi2 加入累加和∑1/Yi2 RCALL LD4

RCALL GET6 ;取出1/Yi2 RCALL GETA ;再取Xi RCALL FPMU ;计算Xi/Yi2 RCALL GET5

RCALL FPAD ;Xi/Yi2加入累加和∑Xi/Yi2 RCALL LD5

DEC R0 ;点数减1

BRNE LOOP3 ;未到总点数n，循环 RCALL GET4 RCALL SAV0

RCALL GET3

RCALL FPMU ;计算(∑1/Yi2)*(∑(Xi/Yi)2) 并存入 RCALL LD6

RCALL GET5 ;取出∑Xi/Yi2 RCALL SAV0

RCALL FPMU ;计算(∑Xi/Yi2)2 RCALL GET6

RCALL FPSU ;计算c=(∑1/Yi2)*(∑(Xi/Yi)2-(∑Xi/Yi2)2 RCALL LD6 ;存入 RCALL GET2 RCALL SAV0

RCALL GET4

RCALL FPMU ;计算(∑(Xi/Yi)*(∑1/Yi2)并存入 RCALL LD7 RCALL GET1 RCALL SAV0

RCALL GET5

RCALL FPMU ;计算(∑1/Yi)*(∑(Xi/Yi2) 并存入 RCALL GET7

RCALL FPSU ;计算d=(∑Xi/Yi)*(∑1/Yi2)-(∑1/Yi)*(∑Xi/Yi2)) RCALL GET6 ;取c RCALL EXCH

RCALL FPDI ;计算b=d/c并存入 RCALL LD7

RCALL GET5 ;取 ∑Xi/Yi2

RCALL FPMU ;计算(∑Xi/Yi2)*b RCALL GET1

RCALL FPSU ;计算(∑1/Yi)-(∑Xi/Yi2)*b RCALL GET4 ;取 ∑1/Yi2 RCALL EXCH

RCALL FPDI ;计算a=(∑1/Yi-(∑Xi/Yi2)*b)/∑1/Yi2 RCALL LD6 ;结果a在$84-$87中,b在$88-$8b中 RER

GETA: LD R12,Y+ LD R13,Y+ LD R14,Y+

LD R15,Y+ ;从外部SRAM中取浮点数到R12-R15 RET

LD4: STS $7C,R12 ;存浮点数 STS $7D,R13 STS $7E,R14 STS $7F,R15 RET

LD5: STS $80,R12 ;计算∑Xi/Yi2的存储单元 STS $81,R13 STS $82,R14 STS $83,R15 RET

LD6: STS $84,R12 ;暂存1/Yi，1/Yi2等浮点数 STS $85,R13 STS $86,R14 STS $87,R15 RET

LD7: STS $88,R12 ;暂存Xi/Yi等浮点数 STS $89,R13 STS $8A,R14 STS $8B,R15 RET

GET4: LDS R8,$7C ;取浮点数 LDS R9,$7D LDS R10,$7E LDS R11,$7F

RET ;

GET5: LDS R8,$80 ;取∑Xi/Yi2或中间结果 LDS R9,$81 LDS R10,$82 LDS R11,$83

RET

GET6: LDS R8,$84 ;取浮点数1/Yi，1/Yi2等 LDS R9,$85 LDS R10,$86 LDS R11,$87 RET

GET7: LDS R8,$88 ;取浮点数Xi/Yi等 LDS R9,$89 LDS R10,$8A LDS R11,$8B RET

;范例90

GETAD: LDI R17,0Bxxx01110;PC0&PC4输入/PC1-PC3输出&PC3(CAL) OUT DDRC,R17 ; CBI PORTC，1

GAD1: SBI PORTC,4

SBIB PINC,4 ;查DRDY

RJMP GAD1 ;低为数据准备好 GAD2: SBI PORTC,4

SBIC PINC,4 ;PINC:$13/PORTB:$15 RJMP GAD2 ;DRDY低有效 CBI PORTC,2 ;置片选有效 LDI R16,16 ;16位数据 GETL0: CLC ;予清除C SBI PORTC,0

SBIC PINC,0 ;接收一位数据 SEC

ROL R14 ;数据高位在前

ROL R13 ;在R13R14里带进位左移 SBI PORTC,1

CBI PORTC,1 ;发出时钟，下降沿读出数据 DEC R16

BRNE GETL0

SBI PORTC,2 ;置片选无效 MOV R4,R14 ; MOV R3,R13 ;保存

GADCOM: CLR R15 ;3字节小数r13r14r15(0)规格化为浮点数 LDI R16,$80

MOV R12,R16 ;阶码为$80 GAD: SBRC R13,7 RJMP GETL2 LSL R14

ROL R13 ;尾数左移，阶码递减 DEC R12 BRNE GAD

RET ;如果(r12)=0 得到0浮点数 GETL2: LDI R16,$7F AND R13,R16 ;正数

LDI R16,$82 ;取浮点数2.5(基准源为2.5v) MOV R8,R16 LDI R16,$20 MOV R9,R16 CLR R10

CLR R11 RCALL FPMU RET

;相乘 ;(r12)不为0