MATLAB的程式档：以微分的计算为例.pdfVIP

MATLAB 的程式檔: 以微分的計算為例 last modified March 30, 2009 MATLAB作為一個數學計算的工具, 不只是提供許多計算的指令而已, 它還是一個電 腦語言 (接近 C 語言), 不但充分運用了一般電腦語言的邏輯能力, 還結合數學相關的 指令成為一個具高度計算能力的程式語言, 可以解決複雜程度較高的數學問題。 本章 以微分的計算初窺 MATLAB 解決數學問題的典型與能力, 也是從學習 MATLAB 的指令跳脫到 MATLAB 程式的第一步。 本章將學到關於程式設計 MATLAB的程式結構、 程式的邏輯概念與程式檔案的管理。 本章關於 MATLAB 的指令與語法 操作元 (operators): .∧ ./ 指令:polyval, format long, semilogx, gtext, diff, polyder, log, fprintf, input 語法: 註解的使用時機與方式。 1 1 背景介紹 微分根據其定義: df (x) f (x + h) − f (x) f (x) = = lim ( ) (1) dx h→0 h 代表函數在某個位置 (x) 的變化情形 (速率), 是陡峭的還是平緩的, 在不同的應用領 域上有其特殊的意義。 這種瞬間的變化在函數上以極限 (lim) 的觀念來表達; 當函數 從 x 移動到 x + h時, 函數值的變化從 f (x) 到 f (x + h), 這個變化速率表示為 f (x + h) − f (x) (2) (x + h) − x h 愈小代表愈短距離的速率, 當 h 小到幾乎等於 0 時, 稱為函數在 x 的瞬間變化率 (一般也稱為 「斜率」), 代表函數在這個點上的動態表現。 圖 1表示出函數在兩個點上 的變化情形。 圖 1: 微分示意圖 配合函數的型態, 已推展出許多計算極限的微分技巧, 這些技巧為人所用所熟悉, 但是 電腦可沒這麼聰明, 雖有人工智慧的輔助, 還是不及人腦來得靈光。 而且技巧總是有 2 限, 碰到一些 「不友善」 的函數時, 往往不夠用。 這時電腦可派上用場了, 除了可以迅 速並大量計算微分問題外, 函數的 「長相」 通常它是認不得的。 撇開像 Mathematica 或 Maple 及 MATLAB Simulink(骨子裡也是 Maple) 這類可以根據數學符號就能 計算微分的軟體, 以電腦程式來計算數值微分, 通常會從上述的定義來下手。 以下列舉 幾個範例來練習, 看看電腦如何計算函數 f (x) 在任一點的導數 (或稱切線斜率)。 1 2 練習 2 範例1: 假設f (x) = x + 3x + 1 , 求f (3)。 • 定義上的極限lim 在數值計算上如何表示? • h值該給多少才算趨近於0? 試試看當 h = 0.1, 0.01, 0.001, 0.0001時, 導數值 分別為多少? 與實際值相差多少? 哪一個選擇比較好? 導數值由下列公式計算: f (3 + h) − f (3) f (3) = h 上述的公式需要計算兩個多項函數值, 輸入時稍嫌麻煩, 此時可以使用 MATLAB 提 供計算多項式