【正文】 4 1/5 1/2 1/3 1/4 1/5 1/6 1/3 1/4 1/5 1/6 1/7 1/4 1/5 1/6 1/7 1/8 1/5 1/6 1/7 1/8 1/9 增加可讀性 format rat, n=input(‘n=’) for i=1:n for j=1:n, h(i, j)=1/(i+j1)。n=39。 end, x1 x1 =+308 （ 3） for語句 for k= 初值：增量：終值 語句組 A， end 將語句組 A反復執(zhí)行 N次，每次執(zhí)行時程序中的 k值不同。 while x~=inf, x1=x。 end, y1 y1 = 求 MATLAB相對精度， y不斷減小，直至 MATLAB 分不出 1+y與 1的差別為止。 y=1。 end, x1 x1 = +307 求 MATLAB中的最大實數(shù)。 while x~=inf, x1=x。 ? 修改為用戶無輸入時程序自動中止。, end ? n = [ ] a =empty ? 輸入數(shù) n，判斷奇偶性。, else a=39。, elseif rem(n,2)==0 a=39。), if isempty(n)==1 a=39。 n = 8 a =even ? n = [ ] a =odd ? n=input( 39。odd39。even39。n=39。 兩個對元素群運算的函數(shù)： 列中有一個元素為 0，即為 0 列中元素全為 1，才為 1 列中有一個元素為 1，即為 1 A =0 0 1 1 B =0 1 0 1 Aamp。B A|B ~A xor(A, B) G=magic(6) rem(G,3) 將邏輯運算用于元素群，得出同階的 01矩陣。 語句分割符，不顯示 = 賦值符 ‘ 轉置，引用 ! 操作系統(tǒng)命令 % 注釋符 == 關系相等符 關系大小符 ~= 關系不等符 amp。 A=magic(6) A =35 1 6 26 19 24 3 32 7 21 23 25 31 9 2 22 27 20 8 28 33 17 10 15 30 5 34 12 14 16 4 36 29 13 18 11 rem(A,3) ans = 2 1 0 2 1 0 0 2 1 0 2 1 1 0 2 1 0 2 2 1 0 2 1 0 0 2 1 0 2 1 1 0 2 1 0 2 p=(rem(A, 3)==0) p =0 0 1 0 0 1 1 0 0 1 0 0 0 1 0 0 1 0 0 0 1 0 0 1 1 0 0 1 0 0 0 1 0 0 1 0 lp=find (p)39。) ? disp([x, sin(x) cos(x) tan(x)]) ? 顯示 x sin(x) cos(x) tan(x) ? 0 0 0 ? ? ? ? ? ? ? 邏輯判斷及流程控制 1. 關系運算 a= 2+2==4 a=(2+2==4) a=(34) a=(43) a=(3=4) a=(4=3) a=(43) a=(34) 等于， a = 1 a = 1 小于， a = 1 a = 0 小于等于， a = 1 a = 0 大于， a =1 a =0 a=(4=3) a=(3=4) a=(3~=4) A=magic(6) rem(A,3) p=(rem(A, 3)==0) lp=find (p)39。 ? 基本函數(shù)庫（ elfun） ? x=[0: : pi/4]’ ? x = ? 0 ? ? ? ? ? ? ? ? disp(39。 區(qū)別左邊運算 輸入算式 d d^3 d.^3 輸出結果 1 4 7 8 5 2 3 6 0 627 636 510 804 957 516 486 612 441 1 64 343 512 125 8 27 216 0 輸入算式 3.^d 3^d 輸出結果 3 81 2187 6561 243 9 27 729 1 +005 * + + + + + + ? （ 3）元素群的函數(shù) ? 除矩陣運算的乘、右除、左除、冪指數(shù)（ / \ ^）、 sqrtm、 expm、logm函數(shù)外，基本函數(shù)庫中的常用函數(shù)都可用于元素群運算。 3 6 0] d^3 d.^3 3.^d 3^d 元素群的冪次運算是各個元素自行作冪次運算，對每個元素的這種運算和對標量運算一樣。（標量會自動擴展為同階矩陣參與運算） x=[1, 2, 3] y=[4, 5, 6] z=x.*y z=x.\y z=x.^y z=x.^2 z=2.^[x y] z =4 10 18 x*y不能成立 z = 元素群沒有左除右除之分 z =1 32 729 x^y 能成立嗎？ z =1 4 9 x^2能成立嗎？ z =2 4 8 16 32 64 2^[x y] 能成立嗎？ d=[1 4 7。運算符前加 .號，表示元素群運算。 t = 0 z = 10 7 4 1 2 5 k = 1 2 3 4 5 6 ? theta = 0 ? theta= linspace (0, 2*pi, 9) w=logspace (0, 1, 11) 用 linspace函數(shù) 格式： linspace（初值、終值、點數(shù)） logspace函數(shù)，自變量按等比級數(shù)賦值。 ? B8 =8 3 4 1 5 9 6 7 2 0 1 2 元素群運算 （ 1）數(shù)組及其賦值 數(shù)組是單行或單列的矩陣，一個 N階的數(shù)組可以表述為一個 N組向量。 fliplr矩陣左右翻轉 flipud矩陣上下翻轉 reshape階數(shù)重組（元素總數(shù)不變） B4=rot90(A) B5=diag(A) B6=tril(A) B7=triu(A) B8=A(: )39。 3 5 7 1。底數(shù)和指數(shù)不能同時為矩陣。 3 6 0] D^2 2.^D D^S 冪次運算：矩陣為底數(shù)，指數(shù)是標量，同矩陣乘法一樣，為保內階數(shù)相同，底數(shù)的矩陣必須是方陣。 3 4] D=[1 4 7。 7] ? X=A\B ? A = 6 3 4 ? 2 5 7 ? 8 4 3 ? B = 3 ? 4 ? 7 ? X = ? ? （ 3）矩陣的乘方和冪次函數(shù) MATLAB的運算符 *、 /、 \、 和 ^，指數(shù)函數(shù) expm、對數(shù)函數(shù) logm和開方函數(shù) sqrtm是對矩陣進行的，即把矩陣作為一個整體來運算。 8 4 3] ? B=[3。 ans = ? 0 ? ? 0 ? a/d ans = ? ? ? 解線性方程組 Ax=B ? 6x1+3x2+4x3=3 ? 2 x1+5 x2+7 x3=4 ? 8 x14 x23 x3=7 ? A=[6 3 4。*b ans = ? 6 16 20 ? 9 23 25 ? 12 30 30 ? a*b39。 8 5 2。4 5 6] ? b=[2 4 0。2] ? ax39。 7 8 9] ? x=[x1,x2,x3] ? b=[2。 ? a=[1 2 3。 未知矩陣在左 . D的逆陣右乘以 B，記作 /D 右除。 V存在條件： A的行列式不等于 0，det(A)≠0 V=A1 MATLAB內部函數(shù) inv，得出 A的逆陣 V。 I為 eye(n)。B(k,j)值為 A陣第 i行和 B陣第 j列對應元素乘積的和。 兩矩陣相乘的必要條件是 內階數(shù)相等 。 a = 1 2 3 ans = 1 2 3 ans = 1 2 3 4 5 6 4 5 6 4 5 6 7 8 9 7 8 9 7 8 9 矩陣 A n p階與 矩陣 B p m階的乘積 C是 n m階矩陣。 1 0 1 0 0 0 eye(3)*a 左、右乘結果不同，只有單位矩陣例外。 y =[2 1 0]。 X與 y內階數(shù)不同，將 y轉置 y’。 pi*x