g****t 发帖数: 31659 | 1 【 以下文字转载自 Mathematics 讨论区 】
发信人: guvest (我爱你老婆Anna), 信区: Mathematics
标 题: 请教个exp(x)的估计问题
发信站: BBS 未名空间站 (Wed Mar 31 21:44:29 2010, 美东)
对给定c/d,0
求最佳整数估计,即整数0
使得exp(c/d)-a/b 最小. |
s*****o 发帖数: 22187 | 2 Are (a,b,c,d) only positive integers or all integers that <2^16?
【在 g****t 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 【 以下文字转载自 Mathematics 讨论区 】 : 发信人: guvest (我爱你老婆Anna), 信区: Mathematics : 标 题: 请教个exp(x)的估计问题 : 发信站: BBS 未名空间站 (Wed Mar 31 21:44:29 2010, 美东) : 对给定c/d,0: 求最佳整数估计,即整数0: 使得exp(c/d)-a/b 最小.
|
g****t 发帖数: 31659 | 3 都正整数.
【在 s*****o 的大作中提到】![](/moin_static193/solenoid/img/up.png) : Are (a,b,c,d) only positive integers or all integers that <2^16?
|
s*****o 发帖数: 22187 | 4 可以用穷举法编个程序吗?我知道很笨,但应该可行吧。瞎写个算法,跟定有错,你看
看能不能改改再用吧。
y=exp(c/d);
m=int(y);
n=int(y)+1;
min=2^16;
for i = 1 : 2^16-1
{ for j = i*m : i*n
{ delta = y-i/j;
if delta < min
{ min = delta;
a=i; b=j;
}
}
}
reture a,b;
【在 g****t 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 都正整数.
|
ET 发帖数: 10701 | 5 赞爱动手。。
【在 s*****o 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 可以用穷举法编个程序吗?我知道很笨,但应该可行吧。瞎写个算法,跟定有错,你看 : 看能不能改改再用吧。 : y=exp(c/d); : m=int(y); : n=int(y)+1; : min=2^16; : for i = 1 : 2^16-1 : { for j = i*m : i*n : { delta = y-i/j; : if delta < min
|
p***o 发帖数: 1252 | 6 有理逼近用连分数比较靠谱。
【在 s*****o 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 可以用穷举法编个程序吗?我知道很笨,但应该可行吧。瞎写个算法,跟定有错,你看 : 看能不能改改再用吧。 : y=exp(c/d); : m=int(y); : n=int(y)+1; : min=2^16; : for i = 1 : 2^16-1 : { for j = i*m : i*n : { delta = y-i/j; : if delta < min
|
s*****o 发帖数: 22187 | 7 偶笨,抛砖引玉:D
请问啥是连分数?
【在 p***o 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 有理逼近用连分数比较靠谱。
|
p***o 发帖数: 1252 | 8
你那个肯定能找着,就是慢了点。连分数wiki一下,不难。
【在 s*****o 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 偶笨,抛砖引玉:D : 请问啥是连分数?
|
s*****o 发帖数: 22187 | 9 恩,谢谢。很长知识,不过好像不太容易界定a,b<2^16。
【在 p***o 的大作中提到】![](/moin_static193/solenoid/img/up.png) : : 你那个肯定能找着,就是慢了点。连分数wiki一下,不难。
|
g****t 发帖数: 31659 | 10 需要的是整数之比.
因为要做firmware,real time应用.
device的CPU只能处理单字节整数.
偶笨,抛砖引玉:D
请问啥是连分数?
【在 s*****o 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 偶笨,抛砖引玉:D : 请问啥是连分数?
|
g****t 发帖数: 31659 | 11 太慢了,没法用 :-(
【在 s*****o 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 可以用穷举法编个程序吗?我知道很笨,但应该可行吧。瞎写个算法,跟定有错,你看 : 看能不能改改再用吧。 : y=exp(c/d); : m=int(y); : n=int(y)+1; : min=2^16; : for i = 1 : 2^16-1 : { for j = i*m : i*n : { delta = y-i/j; : if delta < min
|
j********l 发帖数: 551 | 12 I think what you are looking for is a rational representation of an
irrational number, or rational approximation. One example is how to
approximate pi to a rational number like 22/7. Google it and there are quite
a few papers on this topic.There is a recursive algorithm for it but I don'
t quite remember the name.
【在 g****t 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 太慢了,没法用 :-(
|