SCILAB第六章微积分应用.PDFVIP

SCILAB 第六章 微積分應用  撰文者 ：  1. 中央大學大氣系 楊善文  2.  ( 對本文有貢獻者 ，記得在此留名 ) Conférencier : Prénom Nom Date : Transparent 1 Fichier : /albafile/重要文件/網頁製作/SCILAB_Chap6.odp 微分運算  導函數代表著函數在某個時間與下個時間之間，變化量 的大小。當時間間隔極小時，即為瞬時變化率 。求解導 函數過程即稱為微分。定義 ：  Lim [ h-0, (f(t+h) – f(t)) / h ]  差分的定義 ：D = (f(b) - f(a)) / (b – a) ，差分在時間 間隔極小時 ，具有微分的性質。  Scilab 計算差分的指令是 diff();  diff 每對函數差分一次 ，函數矩陣內含值數量 -1 。 Conférencier : Prénom Nom Date : Transparent 2 Fichier : /albafile/重要文件/網頁製作/SCILAB_Chap6.odp diff 指令的使用  使用方法 ：diff( 函數 , 微分次數 ( 預設 1))  範例 ：  -- x=1:10;  -- y=x ^3;  -- diff(y)  7 19 37 6 1 9 1 127 169 217 271  ( 它使用的規則 ：y(x+1)-y(x) / 1)  -- diff(y,2)  12 18 24 30 36 42 48 54 == y=6x  -- diff(y,3)  6 6 6 6 6 6 6 == y=6 Conférencier : Prénom Nom Date : Transparent 3 Fichier : /albafile/重要文件/網頁製作/SCILAB_Chap6.odp diff 與原函數繪圖比較  範例 ：  clear  x=linspace(0,2*%pi,200);  y=sin(x);  y1=diff(sin(x));  y1(200)=y1(199);  // 此舉只是用來補其差分後會少一筆資料。  plot(x,y,color,[1 0 0]);  plot(x,y1,color,[0 1 0]);  // 紅色是原函數 ，綠色是差分後 。 Conférencier : Prénom Nom Date : Transparent 4 Fichier : /albafile/重要文件/網頁製作/SCILAB_Chap6.odp Conférencier : Prénom Nom Date : Transparent 5 Fichier : /albafile/重要文件/網頁製作/SCILAB_Chap6.odp 修正  各位應該可以發現到 ，得到的似乎不是我們想要的結果 。本來預期 diff(sin(x)) 應該是 cos(x) 才對 。卻發現它的振幅遠小於一 ，關鍵就 出在 diff() 的定義分母的部份 。不過各位讀者不用擔心 ，其實只要 寫一隻簡單的程式 ，就可以完