【文章內(nèi)容簡介】 函數(shù)都是對矩陣中的元素分別進行 . C=[1,2。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)于其他許多語言的一個特色 。 MATLAB的矩陣和數(shù)組運算 矩陣和數(shù)組 是 MATLAB中的基本單位 矩陣運算 數(shù)組運算 +加 減 *乘 ^冪 \左除 /右除 ’轉(zhuǎn)置 .+加 .減 .*乘 .^冪 .\左除 ./右除 .’轉(zhuǎn)置 以矩陣運算法則進行運算 以相應(yīng)位置元素進行運算 一、數(shù)組及其賦值 ? 數(shù)組通常是指單行或單列的矩陣，也稱為向量。 ? 其賦值方法： –（ 1） x=[初值：增量：終值 ] –（ 2） 線性分割函數(shù) x=linspace（起點 ,終點 ,點數(shù)） –（ 3） 對數(shù)分割函數(shù) logspace(起點對數(shù) ,終點對數(shù) ,點數(shù) ) 二、元素群的四則和冪次運算 ? 元素群運算：就是把兩矩陣按逐個元素進行運算。為了與矩陣作為整體的運算符號相區(qū)別，要在運算符 “ *、／、＼、 ^”前加一點符號 “ ． ” 。 元素群運算的英文原文是 ‘ 數(shù)組運算 ’ (Array Algorithm),實際上它適用于任何階的矩陣，故取名 ‘ 元素群運算 ’ 更為確切。 三、元素群的函數(shù)運算 ? 所有的 MATLAB函數(shù)都適用于作元素群運算，只有專門說明的幾個除外。 就是 *、 /、 \、 ^ 運算符和 sqrtm、 expm、logm三個函數(shù)。 ? 實例，做一個三角函數(shù)表： x=[0::pi/4]39。 [x,sin(x),cos(x),tan(x)] x=[0::pi/4]39。 disp(39。 x sin(x) cos(x) tan(x)39。) [x,sin(x),cos(x),tan(x)] x sin(x) cos(x) tan(x) ans = 0 0 0 disp([x,sin(x),cos(x),tan(x)]) 0 0 0 disp([x sin(x) cos(x) tan(x)]) 0 0 0 關(guān)系與邏輯運算 ? 關(guān)系運算符 ? 邏輯運算符 = = == ~= amp。 與 | 或 ~ 非 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。 B=0。 Aamp。B ans = 0 A|B ans = 1 ~A ans = 0 xor(A,B) ans = 1 A=[1 1 0。1 1 1。0 0 0。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。x10 %判斷條件 2是否滿足 y=2*x1, %滿足條件 2的結(jié)果 else %其他情況 ,即條件 3 y=3*x11, %滿足條件 3的結(jié)果 end %程序結(jié)束 請輸入一個值 x=9 x = 9 y = 9 保存文件 function y=yx(x)。 if x1 y=x, elseif x1amp。x10 y=2*x1, else y=3*x11, end 主窗口中輸入 : yx(5) y = 9 ans = 9 yx(8)。 y = 15 ? switch 語句 switch a %讀入一個語句 case A1 %情形 1 B1 %命令 1 case A2 %情形 2 B2 %命令 2 case … … otherwise Bn %最后一個命令 判斷輸入數(shù) n的奇、偶、空的程序： clear% n=input(39。n=39。)。 switch mod(n,2) case 1 A=‘奇 39。 case 0 A=‘偶 39。 otherwise A=‘空 39。 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。 循環(huán)結(jié)構(gòu) ? for 語句 for 變量 =表達式 命令 1 命令 2 … end ? while 語句 for 用于循環(huán)次數(shù)已知 while 表達式 命令 end while 用于循環(huán)次數(shù)未知 break 語句用來跳出循環(huán) ,繼續(xù)執(zhí)行循環(huán)語句的下一語句。 continue語句跳過循環(huán)體中所有剩下的語句，繼續(xù)下一次循環(huán)。 計算階乘的程序 (): clear prd=1。 %放置乘積的變量 ,且設(shè)置初值為 1 n=10。 %設(shè)置循環(huán)的次數(shù) for k=1:n %循環(huán)語句 prd=prd*k。 %計算階乘 end prd %輸出結(jié)果 主窗口輸出結(jié)果 : prd = 3628800 也可用 factorial(n)來計算： factorial(10) ans = 3628800 求和的程序 (一 ):求 1到 20 的階乘和 clear% sum=0。 for i=1:20 prd=1。 for k=1:i prd=prd*k。 end sum=sum+prd。 end Sum 主窗口輸出結(jié)果 : sum = +018 求和的程序 (二 ):求 1到 100的和 clear% sum=0。 i=1。 while i=100 sum=sum+i。 i=i+1。 end sum 主窗口輸出結(jié)果 : sum = 5050 例：小猴吃桃問題。設(shè)第 k天的桃子數(shù)為 ，則桃子數(shù)的變化規(guī)律為： 121 1 ?? ?kk ppkp給定 k, ，可算出第一天的桃子數(shù)。 kpp(10)=1。 for k=10:1:2。 p(k1)=2*(p(k)+1)。 end p(1) ans =1534 例：小猴吃桃問題。設(shè)第 天的桃子數(shù)為 ，則桃子數(shù)的變化規(guī)律為：給定 ，可算出第一天的桃子數(shù)。例 平面簡諧波： x=::。%設(shè)定 x的取值范圍 for i=0:2%用循環(huán)語句令 i分別取 0,1,2 t=*i。 y=*cos(pi*5*x200*pi*t)。%此為時刻 t的波函數(shù) if i==0。%用選擇語句分別用不同的顏色線型畫不同時刻的波形圖 plot(x,y,39。k39。)%用黑色實線畫 t=0時刻的波形圖 hold on%保存圖形的命令，否則后一幅圖會覆蓋前一幅 grid on%繪制網(wǎng)格 elseif i==1。 plot(x,y,39。r39。)%用紅色虛線畫 t= elseif i==2。 plot(x,y,39。b.39。)%用藍色虛點線畫 t= end%結(jié)束 if語句 end%結(jié)束 for語句 ? 直角坐標中的兩維