一個很好的運算反漸開線函數的循環公式

阿松

這是從一個德國人的程序里抄來的，這個德國人的名字叫Emil Meyer

* a+ ?7 U4 |' ]3 s' S, |) E

inva=T

* @- I. C9 Q: V' T$ _

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

6 J s& y; J* |9 X

初始 alpha 值

 alpha_0=pi/2  (T>pi)

% i/ J/ l" [. k" H. ]& \. X

 alpha_0=T^0.25 (T<pi)

3 j" Z* r# e1 K8 g g

 

( x" Q" n8 ]9 X/ g0 G, e6 e5 M

經過我的使用，感覺很少的運算次數就能達到很高的精度。對減少計算機的負荷很有幫助。

logxing

不錯，試了一下，確實很快，比二分法，牛頓法都快
! B( `8 i, C; I, H$ U$ i& ~* y2 N- l, e原理思考中

齒輪工

能不能給個完整的求反漸開線的函數源代碼?十分感謝,現在正在編個小程序,用到了這個函數,結果不會求,能不能幫幫小弟啊?

guojinhuang

1# 阿松
1 K; p( z: l& b1 }你好，能否把這個求解的源程序發給我，急用！

guojinhuang

2# logxing
# ~- ?0 G$ Y' \& a1 c你好，能否把求解的VB源程序發給我呢！謝謝了！

zhang5382

已經給了“漁”，還要人家給“魚”，估計很難有長進了。

zhang5382

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

lalawo72

) h: M/ \  W3 f& q' _; p, O0 T
4 D+ T6 b/ s3 K, S; [" C  b, N4 }哇 ，zhang5382 厲害啊  老書上我看到過