【正文】 ,2 t=*i。設(shè)第 天的桃子數(shù)為 ，則桃子數(shù)的變化規(guī)律為：給定 ，可算出第一天的桃子數(shù)。 p(k1)=2*(p(k)+1)。 kpp(10)=1。 end sum 主窗口輸出結(jié)果 : sum = 5050 例：小猴吃桃問題。 while i=100 sum=sum+i。 end Sum 主窗口輸出結(jié)果 : sum = +018 求和的程序 (二 ):求 1到 100的和 clear% sum=0。 for k=1:i prd=prd*k。 %計算階乘 end prd %輸出結(jié)果 主窗口輸出結(jié)果 : prd = 3628800 也可用 factorial(n)來計算： factorial(10) ans = 3628800 求和的程序 (一 ):求 1到 20 的階乘和 clear% sum=0。 %放置乘積的變量 ,且設(shè)置初值為 1 n=10。 continue語句跳過循環(huán)體中所有剩下的語句，繼續(xù)下一次循環(huán)。 end n=13 A = 奇 n=46 A = 偶 n= A = 空 try 語句組 1 catch 語句組 2 end ? try語句 try語句先試探性執(zhí)行語句組 1，如果語句組 1在執(zhí)行過程中出現(xiàn)錯誤，則將錯誤信息賦給保留的 lasterr變量，并轉(zhuǎn)去執(zhí)行語句組 2。 case 0 A=‘偶 39。)。 y = 15 ? switch 語句 switch a %讀入一個語句 case A1 %情形 1 B1 %命令 1 case A2 %情形 2 B2 %命令 2 case … … otherwise Bn %最后一個命令 判斷輸入數(shù) n的奇、偶、空的程序： clear% n=input(39。 if x1 y=x, elseif x1amp。0 1 0] A = 1 1 0 1 1 1 0 0 0 0 1 0 all(A) ans = 0 0 0 any(A) ans = 1 1 1 流程控制語句 if A1 % 表達式 1 B1 %命令 1 elseif A2 % 表達式 2 B2 %命令 2 else B3 %命令 3 end ? If 語句 表達式A語句組假真ifend表達式A語句組else假真B語句組ifend1表達式A語句組elseif假真B語句組ifendelse假C語句組2表達式????????????)10( 113 )101( 12)1( xxxxxxy例：計算下列函數(shù)的值 clear %清除變量 x=input(‘請輸入一個值 x=’) %任意輸入一個值 if x1 %判斷條件 1是否滿足 y=x, %滿足條件 1的結(jié)果 elseif x=1amp。1 1 1。 Aamp。 與 | 或 ~ 非 xor異或 ? 函數(shù) all (x) 檢查列向量中元素是否全為非零 any(x) 檢查列向量中有無非零元素 邏輯判斷及流程控制 小于 大于 小于等于 大于等于 等于 不等于 關(guān)系運算例子： A=2+3==4 A = 0 A=2+3==5 A = 1 A=2+3=5 A = 1 A=2+3~=9 A = 1 邏輯運算例子： A=1。 x sin(x) cos(x) tan(x)39。 [x,sin(x),cos(x),tan(x)] x=[0::pi/4]39。 就是 *、 /、 \、 ^ 運算符和 sqrtm、 expm、logm三個函數(shù)。 元素群運算的英文原文是 ‘ 數(shù)組運算 ’ (Array Algorithm),實際上它適用于任何階的矩陣，故取名 ‘ 元素群運算 ’ 更為確切。 ? 其賦值方法： –（ 1） x=[初值：增量：終值 ] –（ 2） 線性分割函數(shù) x=linspace（起點 ,終點 ,點數(shù)） –（ 3） 對數(shù)分割函數(shù) logspace(起點對數(shù) ,終點對數(shù) ,點數(shù) ) 二、元素群的四則和冪次運算 ? 元素群運算：就是把兩矩陣按逐個元素進行運算。3 4] C = 1 2 3 4 C^2 ans = 7 10 15 22 expm(C) ans = \ logm(C) ns = + + sqrtm(C) ans = + + C.^2 ans = 1 4 9 16 exp(C) ans = log(C) ans = 0 sqrt(C) ans = 元素群運算 ? 數(shù)組及其賦值 ? 元素群的四則運算 ? 元素群的冪次運算 ? 元素群的函數(shù) 元素群運算能大大簡化編程 ， 提高運算的效率 ， 是 MATLAB優(yōu)于其他許多語言的一個特色 。 X = A\B X = 七、矩陣整體的冪次運算 ? MATLAB的乘冪函數(shù) ” ^”、指數(shù)函數(shù) expm、對數(shù)函數(shù) logm、和開方函數(shù) sqrtm是對矩陣進行的，另外，還有矩陣乘法和除法也是把矩陣作為一個整體來運算。4。 8,4,3]。 x=A\b x = 六、用矩陣除法解線性方程組 解方程組 6 x1 + 3 x2 + 4 x3 = 3 2 x1 + 5 x2 + 7 x3 = 4 8 x1 4 x2 3 x3 = 7 的程序為： A = [6,3,4。0。7 8 9]。x3]=A\b 1231 2 3 2* 3 5 4 07 8 9 2xA X x bx??? ? ? ???? ? ? ?? ? ? ???? ? ? ???? ? ? ?? ? ? ??? A=[1 2 3。4 5 6] B = 1 2 3 4 5 6 D/B ans = 五、線性方程組表為矩陣相乘 x1 + 2 x2 + 3 x3 = 2 3 x1 5 x2 + 4 x3 = 0 7 x1 + 8 x2 + 9 x3 = 2 可以表為 故有 X=[x1。4 5 6。 則其解為 X=B1*inv(D1)=B1/D1 符號 ‘ /’稱為右除。3 4。4 5 6。5 6 7 8]%2行 4列 B = 1 2 3 4 5 6 7 8 C=A*B %3行 4列 C = 11 14 17 20 23 30 37 44 35 46 57 6 CC=B*A %矩陣的內(nèi)階數(shù)不同不能相乘 ??? Error using == mtimes Inner matrix dimensions must agree. 三、矩陣 ‘ 除法 ’ 線性方程組 D*X=B，如果 D非奇異，即它的逆矩陣 inv(D)存在； 則其解用 MATLAB表為 ： X=inv(D)*B=D\B 符號 ‘ \’稱為左除，即分母放在左邊。3 4。3 4] A = 1 2 3 4 B=[5 6 。 ? 不符合交換律： A*B≠B*A，左乘不等于右乘； ? 若一個變量是標(biāo)量，則按標(biāo)量乘法。設(shè) 矩陣 A 為n p階，矩陣 B為 p m階，則 C是 n m階。3 4] A = 1 2 3 4 B=[5 6 。 size(A)=size(B) 如果其中一個變量是標(biāo)量，則會自動擴展成各元素均取此標(biāo)量值的同階次矩陣。 ? 元素大小差別很大時，不宜用矩陣顯示。 ? 一個矩陣中所有元素用同一顯示格式。 y=sin(x).*cos(x)+sin(3*x).*cos(x).^2。5 7]+[2 4。 A=3+ A = + B=5+ B = + C=[1+2i,3+4i。 ? conj(x)表示共軛，即把 x的虛部反號。另外，也可以將空矩陣賦給一個變量，如： B=[ ] 三、復(fù)數(shù) ? 虛數(shù)符號 MATLAB啟動時定為 i,j，可以不用乘號，連寫在數(shù)字后面。 B=find(A) B = [ ] 這里 [ ]是空矩陣的符號， B=find(A)表示列出矩陣 A中值大于 。例如 A=[ 。 ? 空矩陣 在 MATLAB里，把行數(shù)、列數(shù)為零的矩陣定義為空矩陣。 zeros是 MATLAB內(nèi)部函數(shù) ， 使用格式如下： zeros(m)：產(chǎn)生 m*m階零矩陣； zeros(m,n)：產(chǎn)生 m*n階零矩陣，當(dāng) m=n時等同 于 zeros(m)； zeros(size(A))：產(chǎn)生與矩陣 A同樣大小的零矩陣。 A=diag(V) A= 5 0 0 0 7 0 0 0 2 》 eye(2) ans= 1 0 0 1 》 zeros(2) ans= 0 0 0 0 》 ones(2) ans= 1 1 1 1 pan 友矩陣函數(shù) magic 魔方矩陣 hankel Hankel矩陣 rosser 對稱特征值測試矩陣 hilb Hilbert矩陣 pascal Pascal矩陣 invhilb 反 Hilbert矩陣 vander 范德蒙矩陣 … … （ 7）其他特殊矩陣 ? 零矩陣 :所有元素值為零的矩陣稱為零矩陣 。 ones(m) 對角矩陣： 對角元素向量 V=[a1,a2,…,an] A=diag(V) 隨機矩陣： rand(m,n)產(chǎn)生一個 m n的均勻分別的隨機 矩陣 （ 6）常用的特殊矩陣 如果已知 A為方陣，則 V=diag(A)可以提取 A的對角元素構(gòu)成向量 V。 eye(m) 零 矩 陣： zeros(m,n)。11:15]