【文章內(nèi)容簡(jiǎn)介】 函數(shù)都是對(duì)矩陣中的元素分別進(jìn)行 . 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 = 元素群運(yùn)算 ? 數(shù)組及其賦值 ? 元素群的四則運(yùn)算 ? 元素群的冪次運(yùn)算 ? 元素群的函數(shù) 元素群運(yùn)算能大大簡(jiǎn)化編程 ， 提高運(yùn)算的效率 ， 是 MATLAB優(yōu)于其他許多語(yǔ)言的一個(gè)特色 。 MATLAB的矩陣和數(shù)組運(yùn)算 矩陣和數(shù)組 是 MATLAB中的基本單位 矩陣運(yùn)算 數(shù)組運(yùn)算 +加 減 *乘 ^冪 \左除 /右除 ’轉(zhuǎn)置 .+加 .減 .*乘 .^冪 .\左除 ./右除 .’轉(zhuǎn)置 以矩陣運(yùn)算法則進(jìn)行運(yùn)算 以相應(yīng)位置元素進(jìn)行運(yùn)算 一、數(shù)組及其賦值 ? 數(shù)組通常是指單行或單列的矩陣，也稱為向量。 ? 其賦值方法： –（ 1） x=[初值：增量：終值 ] –（ 2） 線性分割函數(shù) x=linspace（起點(diǎn) ,終點(diǎn) ,點(diǎn)數(shù)） –（ 3） 對(duì)數(shù)分割函數(shù) logspace(起點(diǎn)對(duì)數(shù) ,終點(diǎn)對(duì)數(shù) ,點(diǎn)數(shù) ) 二、元素群的四則和冪次運(yùn)算 ? 元素群運(yùn)算：就是把兩矩陣按逐個(gè)元素進(jìn)行運(yùn)算。為了與矩陣作為整體的運(yùn)算符號(hào)相區(qū)別，要在運(yùn)算符 “ *、／、＼、 ^”前加一點(diǎn)符號(hào) “ ． ” 。 元素群運(yùn)算的英文原文是 ‘ 數(shù)組運(yùn)算 ’ (Array Algorithm),實(shí)際上它適用于任何階的矩陣，故取名 ‘ 元素群運(yùn)算 ’ 更為確切。 三、元素群的函數(shù)運(yùn)算 ? 所有的 MATLAB函數(shù)都適用于作元素群運(yùn)算，只有專門說(shuō)明的幾個(gè)除外。 就是 *、 /、 \、 ^ 運(yùn)算符和 sqrtm、 expm、logm三個(gè)函數(shù)。 ? 實(shí)例，做一個(gè)三角函數(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)系與邏輯運(yùn)算 ? 關(guān)系運(yùn)算符 ? 邏輯運(yùn)算符 = = == ~= amp。 與 | 或 ~ 非 xor異或 ? 函數(shù) all (x) 檢查列向量中元素是否全為非零 any(x) 檢查列向量中有無(wú)非零元素 邏輯判斷及流程控制 小于 大于 小于等于 大于等于 等于 不等于 關(guān)系運(yù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 邏輯運(yùn)算例子： 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 流程控制語(yǔ)句 if A1 % 表達(dá)式 1 B1 %命令 1 elseif A2 % 表達(dá)式 2 B2 %命令 2 else B3 %命令 3 end ? If 語(yǔ)句 表達(dá)式A語(yǔ)句組假真ifend表達(dá)式A語(yǔ)句組else假真B語(yǔ)句組ifend1表達(dá)式A語(yǔ)句組elseif假真B語(yǔ)句組ifendelse假C語(yǔ)句組2表達(dá)式????????????)10( 113 )101( 12)1( xxxxxxy例：計(jì)算下列函數(shù)的值 clear %清除變量 x=input(‘請(qǐng)輸入一個(gè)值 x=’) %任意輸入一個(gè)值 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é)束 請(qǐng)輸入一個(gè)值 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 語(yǔ)句 switch a %讀入一個(gè)語(yǔ)句 case A1 %情形 1 B1 %命令 1 case A2 %情形 2 B2 %命令 2 case … … otherwise Bn %最后一個(gè)命令 判斷輸入數(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 語(yǔ)句組 1 catch 語(yǔ)句組 2 end ? try語(yǔ)句 try語(yǔ)句先試探性執(zhí)行語(yǔ)句組 1，如果語(yǔ)句組 1在執(zhí)行過(guò)程中出現(xiàn)錯(cuò)誤，則將錯(cuò)誤信息賦給保留的 lasterr變量，并轉(zhuǎn)去執(zhí)行語(yǔ)句組 2。 循環(huán)結(jié)構(gòu) ? for 語(yǔ)句 for 變量 =表達(dá)式 命令 1 命令 2 … end ? while 語(yǔ)句 for 用于循環(huán)次數(shù)已知 while 表達(dá)式 命令 end while 用于循環(huán)次數(shù)未知 break 語(yǔ)句用來(lái)跳出循環(huán) ,繼續(xù)執(zhí)行循環(huán)語(yǔ)句的下一語(yǔ)句。 continue語(yǔ)句跳過(guò)循環(huán)體中所有剩下的語(yǔ)句，繼續(xù)下一次循環(huán)。 計(jì)算階乘的程序 (): clear prd=1。 %放置乘積的變量 ,且設(shè)置初值為 1 n=10。 %設(shè)置循環(huán)的次數(shù) for k=1:n %循環(huán)語(yǔ)句 prd=prd*k。 %計(jì)算階乘 end prd %輸出結(jié)果 主窗口輸出結(jié)果 : prd = 3628800 也可用 factorial(n)來(lái)計(jì)算： 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 例：小猴吃桃問(wèn)題。設(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 例：小猴吃桃問(wèn)題。設(shè)第 天的桃子數(shù)為 ，則桃子數(shù)的變化規(guī)律為：給定 ，可算出第一天的桃子數(shù)。例 平面簡(jiǎn)諧波： x=::。%設(shè)定 x的取值范圍 for i=0:2%用循環(huán)語(yǔ)句令 i分別取 0,1,2 t=*i。 y=*cos(pi*5*x200*pi*t)。%此為時(shí)刻 t的波函數(shù) if i==0。%用選擇語(yǔ)句分別用不同的顏色線型畫(huà)不同時(shí)刻的波形圖 plot(x,y,39。k39。)%用黑色實(shí)線畫(huà) t=0時(shí)刻的波形圖 hold on%保存圖形的命令，否則后一幅圖會(huì)覆蓋前一幅 grid on%繪制網(wǎng)格 elseif i==1。 plot(x,y,39。r39。)%用紅色虛線畫(huà) t= elseif i==2。 plot(x,y,39。b.39。)%用藍(lán)色虛點(diǎn)線畫(huà) t= end%結(jié)束 if語(yǔ)句 end%結(jié)束 for語(yǔ)句 ? 直角坐標(biāo)中的兩維