高級計算器切線迭代法的計算程序

摘要：目的：編制切線迭代法的計算程序，提供一種簡單的計算方法。方法：使用CASIO，X-3600P計算器，利用M、K寄存器輸入程序變量，直接按P.鍵運行程序。結果：第一次輸入變量進行運算得出結果后，程序內自動生成全部變量，下一輪計算變量零輸入，快速得到X值。結論：切線迭代法的計算程序，解決了煩瑣復雜的循環(huán)計算問題，滿足了平均發(fā)展速度計算工作的需要，具有操作簡單、運算方便、結果準確的優(yōu)點。

關鍵詞：計算器;切線迭代法;平均發(fā)展速度;程序

中圖分類號：TP311 文獻標識碼：A

文章編號：1009-3044（2020）09-0234-02

平均發(fā)展速度，是表示某種現(xiàn)象在一個較長時期中逐年平均發(fā)展變化程度的指標，是一種廣泛應用的統(tǒng)計方法。切線迭代法，就是先選定初始值Xo，將其代入迭代公式，不斷重復進行循環(huán)運算，直至最終兩次計算結果大體相近或完全相等，即為平均發(fā)展速度[1]。

CASIO_fX-3600 P及同類功能計算器，使用比較普遍。編程運算X值，降低了計算工作的強度，提高了工作效率。現(xiàn)介紹如下：

（限于篇幅公式從略）

1.3 退出與運算

MODE.

2 運算實例

2.1 根據(jù)公式手工計算出C值

2.2 使用P.區(qū)程序.進行迭代計算出近似值

循環(huán)運算結果顯示，Xs值與X6值完全相等，可以認為0.773487218（77.35%）即為淄博市10年法定報告?zhèn)魅静〉钠骄l(fā)展速度。

用平均發(fā)展速度減去1（100%），如果得到正值即為平均遞增速度。如果得到負值即為平均遞減速度。本例x= -22.65，說明淄博市法定報告?zhèn)魅静“l(fā)病率從1983年起，平均每年以22.65%的速度下降。

3 討論

1）程序運算前，首先將Xo=1存人K4寄存器備用，程序運算得出X1值后，再次存人K4寄存器，往復循環(huán)，省略輸入。

2）P1區(qū)寫入程序時，第2步至第8步完成了切線迭代公式分子的部分計算，第26步將C值再加入K1寄存器，切線迭代公式的分子計算完成。第9步至第15步完成了切線迭代公式的分母計算。第5步與第6步完成了X的指數(shù)運算。第17步將1累加入K3寄存器，第18步K5寄存器和第19步K6寄存器再減去1，做到了n、n-l的遞減。運算結束后，K5寄存器數(shù)據(jù)為一1，K6寄存器數(shù)據(jù)為0，第23步將n加入K5寄存器，第24步將n存人K6寄存器，為下一次運算輸入了變量。第20步讀取K3寄存器n值，用于判斷代人是否結束，然后按RUN鍵顯示X值[3]。

3）P1區(qū)寫入程序時，第21步暫停，第22步便減去K3寄存器數(shù)據(jù)，將K3寄存器清零，第28步和第31步的操作，同樣將K2寄存器和K1寄存器清零，這樣在程序運算時，不用再對計算器進行清零操作。

4）根據(jù)有關文獻資料評價，四種常用平均發(fā)展速度的計算方法，切線迭代法精度最高，準確性優(yōu)于幾何平均法、一次平均法、二次平均法。切線迭代法的計算程序，為基層衛(wèi)生工作人員計算平均發(fā)展速度，提供了一種簡單的計算方法。

參考文獻：

[1]董軍，胡松華，四種計算平均增長速度方法的比較[J].數(shù)理醫(yī)藥學雜志，1996（3）：216-218.

[2]鞏翠華，丁桂芬，劉佰燕，等.切線迭代法計算平均發(fā)展速度[J].中國公共衛(wèi)生，1994，10（5）：224-225.

[3]畢永.園型分布的計算程序[Jl.疾病監(jiān)測，1995，10（8）：246-248.

【通聯(lián)編輯：梁書】

作者簡介：畢永（1960-），男，衛(wèi)生檢驗師。