上周上机课题目数值分析.pptVIP

上週上機課題目：數值分析 請利用下列公式算出 f ’(8.1) 、f ’(8.3) 的近似值 ? f ’(x0)=1/2h [-3f(x0)+4f(x0+h)-f(x0+2h)] 這裡的 f(x) 為 x lnx Note：注意！乘的運算為 * a=-3*8.1*log(8.1) ; %-3f(x0) b=4*8.2*log(8.2) ; %4f(x0+h) c=-8.3*log(8.3); %-f(x0+2h) 1/(2*0.1)*(a+b+c) %(1/2h)(a+b+c) ans = 3.09191394208330 Note： ? ? 猜猜！1/(2*0.1)與1/2*0.1有差別嗎？ 5 0.05 事實上，可以簡化成 a=8.1;b=0.1; z=1/(2*b)*(-3*a*log(a)+4*(a+b)*log(a+b)-… (a+2*b)*log(a+2*b)) z= 3.09191394208330 這樣有個好處： 可以隨意的修改a,b的值來調整z的值 數學軟體簡介 PART IIMatlab 介紹 今天進度：Matlab基本環境操作 重複命令與邏輯命令 M檔案 基本二維空間繪圖 for 迴圈 最簡單的重複命令是for迴圈 其基本形式為：  for 變數 = a:b； 運算式； end 其中變數的值會被依次設定為a~b，來執行介於for和end之間的運算式。因此，若無意外情況，運算式執行的次數會等於b-a+1。  for 迴圈 例如： format long %算至小數點下14位數 b = zeros(1,10); %b是一個1x10的零矩陣 for i = 1:10; b(i) = 1/(10^i); end b 0.1 0.01 0.001 0.0001 0.00001 0.000001 0.0000001 0 0.000000001 0.0000000001 for 迴圈 例如： for i=1:10; % i從 1到 10 a=8.1; h=1/(10^i); %從0.1到0.0000000001 z(i)=1/(2*h)*(-3*a*log(a)+4*(a+h)*… log(a+h)-(a+2*h)*log(a+2*h)); end z z = 3.09191394208334 3.09186456879100 3.09186406675011 3.09186406168394 3.09186406113326 3.09186405367257 3.09186402702721 3.09186507507775 3.09186809488438 3.09183789681811 for迴圈也可以是多層的，下例產生一個6 x 6的Hilbert矩陣h，其中Ai,j為 ：1/(i+j-1) format rat for i = 1:6; for j = 1:6; h(i,j) = 1/(i+j-1); end end 另一個重複命令是while迴圈 基本形式為：  while 條件式； 運算式； end 也就是說，只要條件是成立，運算式就會一再被執行。右例是產生一數列 x： [ 1 1/2 1/3 1/4 1/5 1/6 ] 邏輯命令 最簡單的邏輯命令是if, ..., end 其基本形式為： if 條件式； 運算式； end if min(x) 0.02 ; disp(Very Big! ); % disp　印出 else disp(So tiny! ); end Very Good! M檔案 Matlab可以將一連串命令貯存在一個檔案中，以後再將這群指令以你訂的檔案名稱為新指令，集體同時執行。這種可執行的檔案稱為M-files，因為其檔案名稱之字尾均以 『.m』附檔名。 例如； myfirst.m 即表示其檔案名稱為『 myfirst』，檔案型式為 M-file，你只要將目錄指向存放該檔案之檔案夾，就可在指令窗下打入 myf