一個(gè)很好的運(yùn)算反漸開(kāi)線(xiàn)函數(shù)的循環(huán)公式

阿松

這是從一個(gè)德國(guó)人的程序里抄來(lái)的,，這個(gè)德國(guó)人的名字叫Emil Meyer

% c4 ]0 t h& L" B0 [( G/ N: G% R

inva=T

0 O9 \1 Q- h$ m- \

alpha_i+1=-1/(alpha_i+T)+alpha_i+1/tan(alpha_i)

/ N% J( [2 @$ Y: [' a# S: x- @

初始 alpha 值

 alpha_0=pi/2  (T>pi)

 alpha_0=T^0.25 (T<pi)

 

1 M; g$ V' P4 i# L# S$ o

經(jīng)過(guò)我的使用,，感覺(jué)很少的運(yùn)算次數(shù)就能達(dá)到很高的精度,。對(duì)減少計(jì)算機(jī)的負(fù)荷很有幫助,。

logxing

不錯(cuò),，試了一下，確實(shí)很快,，比二分法,，牛頓法都快
  R$ L6 `/ d0 ^  r; `, ~1 i' u原理思考中

齒輪工

能不能給個(gè)完整的求反漸開(kāi)線(xiàn)的函數(shù)源代碼?十分感謝,現(xiàn)在正在編個(gè)小程序,用到了這個(gè)函數(shù),結(jié)果不會(huì)求,能不能幫幫小弟啊?

guojinhuang

1# 阿松
, v! W6 J3 \" N. b7 N0 l! P/ F你好，能否把這個(gè)求解的源程序發(fā)給我,，急用,！

guojinhuang

2# logxing
你好，能否把求解的VB源程序發(fā)給我呢,！謝謝了,！

zhang5382

已經(jīng)給了“漁”，還要人家給“魚(yú)”,，估計(jì)很難有長(zhǎng)進(jìn)了,。

zhang5382

alpha_0一般可選為：(3*T)^(1/3)-T/3,，更快。

lalawo72

哇 ,，zhang5382 厲害啊  老書(shū)上我看到過(guò)