《数学软体70553673》-公开课件.pptVIP

複數 maple(‘abs’,a)在2009a版本會出現下列訊息： maple(abs,a) ??? Error using == maple at 52 Maple engine is not installed. 複數 如果不用MuPad 利用sym()函數直接將2+3i定義成為符號亦可： mathsoft * 例題 解聯立方程組 2a – b – c = 4 a – b + c = 0 -3a + 2b – 2c = -4 寫為 A*x=B , A=[2 -1 -1; 1 -1 1; -3 2 -1] x=[a b c]’ , B=[4 0 -4]’ 或 x*A=B , A=[2 1 -3; -1 -1 2; 1 1 -1] x=[a b c] , B=[4 0 -4] (註: x, B改為列向量) mathsoft * 矩陣的除法運算 mathsoft * 例題 解下面over-determined 方程組 令 A=[-1 0 2]’ , B=[-2 -1 1]’ 解A*s = B mathsoft * 例題 得到s=0.8, 表示對所有實數 s, s=0.8 使得||A*s-B||2為最小(最小平方差意義)。 mathsoft * 矩陣/向量的函數庫 前面提過的內建函數 abs, sqrt, cos, log, …. 對於矩陣是作用在每一個元素上 關於主對角、上三角或下三角的函數 diag(A) 取A矩陣主對角線上的元素,產生一個行向量。若A為向量,則產生一個主對角矩陣,其主對角線上的元素為A。 diag(A,k) 取A矩陣第k對角線上的元素,產生一個行向量。K0表示在主對角線上方, K0表示在下方。 triu(A) 取A矩陣主對角線以上(含主對角)的元素,其餘為零所產生的一個矩陣,其維度大小與A相同。 mathsoft * 矩陣/向量的函數庫 關於主對角(diagonal)、上三角或下三角(upper/lower triangular) 的函數 triu(A,k) 取A矩陣第k對角線以上(含對角線)的元素,其餘為零所產生的一個矩陣。K0表示在主對角線上方, K0表示在下方。 tril(A) 與triu(A)類似, 只是將主對角線以上, 改為以下。 tril(A,k) 與triu(A,k)類似, 只是將第k對角線以上, 改為以下。 mathsoft * 矩陣/向量的函數庫 關於計算矩陣/向量的函數 x向量與y 向量的內積 dot(x,y) norm(x) 向量 x的2-範數(2-norm)。 norm(x,p) 向量 x的p-範數(p-norm), ||x||p=(∑|xi|p)1/p 。經常使用的p有1 ,2, ∞。 norm(A) 矩陣 A的2-範數(2-norm) 。 norm(A,p) 矩陣 A的p-範數(p-norm), ||A||p=max||Ax||p/||x||p 。經常使用的p有1 ,2, ∞ 。 mathsoft * 矩陣/向量的函數庫 關於統計方面的函數, 原則以向量為主, 若是矩陣則逐一處理矩陣的行向量。 max(x) 向量x中的最大元素。 sum(x) 向量x中元素的總和。 prod(x) 向量x中元素的總乘積。 mean(x) 向量x中元素的平均數。 median(x) 向量x中元素的中位數。 mathsoft * 例題 mathsoft * 例題 mathsoft * 例題 tan(pi/2)→∞,使得tan(pi/4)=1或tan(pi)=0 相對小, 顯示為0. 改用format long e的格式比較能區別其差異. mathsoft * 例題 tan(?pi/2)→?∞,使得其他在(-pi/2 , pi/2)區間內 tan(x)的值相對小(正負不計), 顯示為0. mathsoft * 例題 將圖形x軸的範圍改為(-1.5, 1.5)就改善了顯示的樣貌。 Matlab Notebook 建立MS word的檔案，裡面可以執行Matlab的指令並顯示結果。 首先在Matlab Command Window執行Notebook setup Word 會自動產生setup.doc的檔案 鍵入文字 將要執行的指令標示(marked)為一個input cell, 按下Ctrl-Enter(或evaluate cell)即可執行。 mathsoft * Matlab Notebook mathsoft * Matlab Notebook mathsoft * 這些在一個中括弧內的綠色文字,代表input cell, 即是要執行的指令 Matlab Notebook mathsoft * Matlab Notebook mathsoft