文号
t41511
阅读
17030
回复
16
收藏
10
求问:如何用单片机实现求对数和反对数
2012/03/07
IP:四川
有个问题解决不了,请编程算法高手指点:

如何用单片机求一个值的对数值(以10为底),以及已知对数值,如何求出原值?

情况是这样的,库函数里面有求对数和反对数的函数,但是每求一次,至少要耗费8千个机器周期。而我们要求每秒钟做10万到100万次对数和反对数运算,且这些运算只能占用机器资源的5%。用这个库函数的话,把CPU扩大10倍都办不到。有没有简单一些的算法,尤其是求解析解的算法?
14年10个月前IP:未同步
374419
1楼
输入数据范围?输出范围?精度?CPU主频?
引用
评论
14年10个月前IP:未同步
374420
2楼

attachment icon log.pdf 186.31KB PDF 2665次下载 预览
+50
科创币
虎哥
2012/03/07
嗯。
引用
评论
14年10个月前IP:未同步
374442
3楼
新手出点馊主意:
模拟电路求对数反对数是很简单的事。
实在不行的话,就D/A,模拟求对数,A/D吧。
+10
科创币
乖乖乖
2012/03/08
我也这么想
引用
评论
14年10个月前IP:未同步
374457
4楼
每秒钟做80000000000个机器周期-------
很简单吗?[s:178]
引用
评论
虎哥作者
14年10个月前IP:未同步
374462
5楼
引用第1楼warmonkey于2012-03-07 20:26发表的  :
输入数据范围?输出范围?精度?CPU主频?


误差以分贝计,不超过0.01dB。范围100dB。
引用
评论
14年10个月前IP:未同步
374522
6楼
求次方到好~  直接乘法指令乘乘乘乘乘乘  (指整数~)    求对数嘛。。。[s:275]  
这个LOG要是小数更难处理  单片鸡的CPU必须要有强大的浮点运算单元专用的浮点运算指令   能缩短整数CPU模拟小数运算的程序执行时间     换个猛点的单片鸡吧~
对单片鸡要求太高~  要不虎哥开发一款一个时钟周器1000条指令的CPU  [s:274]
引用
评论
14年10个月前IP:未同步
374534
7楼
引用第1楼warmonkey于2012-03-07 20:26发表的  :
输入数据范围?输出范围?精度?CPU主频?


行家一出手,就知有没有。
引用
评论
14年10个月前IP:未同步
374535
8楼
引用第5楼虎哥于2012-03-07 22:19发表的  :


误差以分贝计,不超过0.01dB。范围100dB。


你分段用线性近似就行了。
引用
评论
14年10个月前IP:未同步
374536
9楼
引用第5楼虎哥于2012-03-07 22:19发表的  :


误差以分贝计,不超过0.01dB。范围100dB。


通常还可以跟前后的计算总和优化,不要死跟着公式走。
引用
评论
14年10个月前IP:未同步
374537
10楼
引用第8楼婺源寻芳于2012-03-08 08:23发表的  :


你分段用线性近似就行了。



如果是连续检测到的资料,还可以用差分法,用高次多项式近似。
引用
评论
14年10个月前IP:未同步
374572
11楼
分段近似的方法很好。还可以扩展下。如,计算机里只保留[0.1,1]之间的数的分段近似对数值表。
对于0.1到1之间的数,用查表+线性插值运算就可以算出来。
对于其他区间的值,通过乘以10的n次幂挪到[0.1,1]之间再求。求得的对数值再加n就可以了。
引用
评论
14年10个月前IP:未同步
374575
12楼
所有运算,可以变成二进制,尽量不要用十进制。
引用
评论
14年10个月前IP:未同步
374578
13楼
对于计算机来说,由于是二进制运算,因此算以二为底的对数更好。相应的区间是[1,2]。
其他整数,只要左移或右移就可以挪到这个区间了。
其他浮点数,只要看表示幂次的位区间就可以了。
至于以2为底的对数如何转变成以10为底的对数,只要按如下公式就可以了:
以2为底的x 等于 以10为底的x 除以 以10为底的2。
引用
评论
14年10个月前IP:未同步
374584
14楼
能否先在PC上算出结果,给MCU一个表,~~~~~查表 套表的数据结果出来用,,,,就不要在MCU运算了............................我承认我在胡扯~~~~不要扣我分[s:213]
引用
评论
14年10个月前IP:未同步
374585
15楼
而且,事实上[1,2]区间比起[0.1,1]区间,线性度更好。
引用
评论
笔记
{{note.content}}
{{n.user.username}}
{{fromNow(n.toc)}} {{n.status === noteStatus.disabled ? "已屏蔽" : ""}} {{n.status === noteStatus.unknown ? "正在审核" : ""}} {{n.status === noteStatus.deleted ? '已删除' : ''}}
  • 编辑
  • 删除
  • {{n.status === 'disabled' ? "解除屏蔽" : "屏蔽" }}
我也是有底线的