【正文】 定位出最大．最小值的位置 [mx,indx] = min(count) mx = 7 9 7 indx = 2 23 24 試試看，你能看懂下面的命令是干什么的嗎？ [n,p] = size(count) e = ones(n,1) x = count – e*mu 點這 看看答案！ 下面這句命令則找出了整個矩陣的最小值 : min(count(:)) ans = 7 協(xié)方差及相關(guān)系數(shù) 下面 ,我們來看看第一列的方差 : cov(count(:,1)) ans = cov()函數(shù)作用于 矩陣 ,則會計算其協(xié)方差矩陣 . corrcoef()用于計算相關(guān)系數(shù) ,如 : corrcoef(count) ans = 數(shù)據(jù)的預(yù)處理 未知數(shù)據(jù) NaN(Not a Number不是一個數(shù) )被定義為未經(jīng)定義的算式的結(jié)果 ,如 0/ ,NaN常用來表示未知數(shù)據(jù)或未能獲得的數(shù)據(jù) .所有與 NaN 有關(guān)的運算其結(jié)果都是 NaN. a = magic(3)。,0) xlabel(39。,39。 set(0,39。 4:察看查找目錄下是否有文件 :。 save([file int2str(i)],’j’)。 長命令行 如想另起一行輸入命令 ,在末尾加上 ...即可 ,如 : s = 1 – 1/2 + 1/3 – 1/4 + 1/5 – 1/6 + 1/7 ... – 1/8 + 1/9 – 1/10 + 1/11 – 1/12。 實時工作室？允許你產(chǎn)生來自你 的圖表塊的 C 代碼，使之能用于各種實時系統(tǒng)。 5. MATLAB 應(yīng)用程序接口（ API）。 datafun 數(shù)據(jù)分析和傅立葉變換。 elmat 初步矩陣，和矩陣操作。 graphics 制圖法。它也為包括低階的指令在內(nèi) ,允許你建造完整的圖形用戶界面（ GUIs）， MATLAB應(yīng)用。 2. MATLAB 工作環(huán)境。 strfun 字符串。 MATLAB 系統(tǒng) MATLAB 系統(tǒng)由 5 主要的部分構(gòu)成： 1. MATLAB 語言。 這是高階的矩陣 /數(shù)組語言 ,帶控制流動陳述，函數(shù)，數(shù)據(jù)結(jié)構(gòu)，輸入 /輸出，而且面向目標的編程特點。 iofun 輸入 /輸出。 這是你作為 MATLAB用戶或程序編制員的一套工具和設(shè)施。制圖法功能在 MATLAB工具箱中被組織成 5文件夾： graph2d 2的維數(shù)上的圖表。 uitools 圖形用戶界面工具。 elfun 初步的數(shù)學函數(shù)。 polyfun 插入物，并且多項式。 這是允許你寫 C、 Fortran 語言與 MATLAB交互。 關(guān)于工具箱 工具箱是為了解答特別種類的問題擴展 MATLAB環(huán)境的 MATLAB函數(shù)的綜合的（ M文件）收集 MatLab 工作環(huán)境 命令窗口 若輸入 A = [1 2 3。 MatLab 工作間 你可用 who 或 whos來察看當前工作間中有哪些變量 ,如 : whos Name Size Bytes Class A 4x4 128 double array D 3x5 120 double array M 10x1 40 cell array S 1x3 628 struct array h 1x11 22 char array n 1x1 8 double array s 1x5 10 char array v 1x14 28 char array Grand total is 93 elements using 984 bytes 若要從工作間中刪除所有的變量 ,用 clear 保存、重載工作間 你可以將工作間保存為一個二進 制的 M 文件，以后還可以恢復回來： save june10 也可只保存工作間中的部分變量值 save june10 x y z 重載時只需輸入 load june10 文件名保存在字符串中 這樣可以像調(diào)用函數(shù)一樣調(diào)用工作間 save(’myfile’,’VAR1’,’VAR2’) A = ’myfile’。 end 查找路徑 當你輸入“ yourpig時發(fā)生了什么呢 ? 1:察看是否是變量 。 對于查找路徑中的文件 ,what顯示當前目錄下的文件 ,加上路徑后可顯示輸入的路徑下所有的 MatLab文件 .如 : what matlab/elfun 以下二命令分別顯示、編輯 m 文件 type rank edit rank 圖像窗口 下面的命令產(chǎn)生一個與命令窗口隔離的圖形窗口 , figure plot函數(shù)則會在新的窗口中繪制圖形 ,如 t = 0:pi/100:2*pi。defaultaxeslinestyleorder’,’||.’) set(0,39。Location 239。Time39。 a(2,2) = NaN a = 8 1 6 3 NaN 7 4 9 2 sum(a) ans = 15 NaN 15 在做統(tǒng)計時，常需要將 NaN 轉(zhuǎn)化為可計算的數(shù)字或去掉 ,以下是幾種方法 : 注 :判斷一個值是否為 NaN,只能用 isnan(),而不可用 x==NaN。 消掉 NaN X(any(isnan(X)’ ),:) = []。 nout = sum(outliers) nout = 1 0 0 count(any(outliers39。 plot(t,y,’o’), grid on 多項式回歸 由圖可以看出應(yīng)該可以用多項式來表達： y=a0+a1*t+a2*t^2 系數(shù) a0,a1,a2可以由最小平方擬合來確定，這一步可由反除號＂＼＂來完成 解下面的三元方程組可得： X = [ones(size(t)) t t.^2] X = 0 0 a = X\y a = – a 即為待求的系數(shù)，畫圖比較可得 T = (0::)’。,t,y,39。 Y = [ones(size(T)) exp(– T) (– T)]*a。o39。,cdate,pop,39。 plot(cdate,pop1,39。) res1 = pop – pop1。 pop2 = polyval(p,sdate)。+39。 plot(cdate,pop4,39。) res4 = pop – pop4。 logpred1 = 10.^polyval(logp1,sdate)。+39。+39。–39。 grid on r = pop – 10.^(polyval(logp2,sdate))。 [pop2,del2] = polyval(p2,sdate,S2)。g–39。r:39。 y = filter(b,a,x)。,t,y,39。,39。 plot(year,wolfer) title(’Sunspot Data’) 現(xiàn)在來看看其 FFT: Y = fft(wolfer)。 nyquist = 1/2。 period(index) ans = 變換后的幅度和相位 abs()和 angle()是用來計算幅度和相位的 先創(chuàng)建一信號 ,再進行分析 ,unwarp()把相位大于π的變換為 2π的補角 : t = 0:1/99:1。 p = unwrap(angle(y))。Abs. Magnitude39。Frequency [Hertz]39。 X = 8 3 4 1 5 9 6 7 2 如ｘ與ｙ均是列向量，則 x*y 無解 ,但下二表達式卻可以 : x39。 若函數(shù)形式是 :y(t)=c1+c2*exp(t)。–39。 讓我們臨時把 A 變一變 : A = pascal(6) A = 1 1 1 1 1 1 1 2 3 4 5 6 1 3 6 10 15 21 1 4 10 20 35 56 1 5 15 35 70 126 1 6 21 56 126 252 A 是二項式系數(shù) ,每一項是其左方與上方系數(shù)之和 ,求其 Cholesky 分解系數(shù)有 : R = chol(A) R = 1 1 1 1 1 1 0 1 2 3 4 5 0 0 1 3 6 10 0 0 0 1 4 10 0 0 0 0 1 5 0 0 0 0 0 1 Ｒ認識二項式系數(shù)． 這樣對于線性方程便可化簡： A*x = b R39。Q = I 正交陣的好處在于 ,她保持了原陣的長度 ,角度 ,并且在計算的過程中不會擴大誤差 . RQ 分解如下 : A = Q R 或 A P = Q R 其中 ,Q 是正交陣 ,R 是上三角陣 . 矩陣的冪與指數(shù) 若 A 是方陣 ,p 是正數(shù) ,則 X = A^2 X = 3 6 10 6 14 25 10 25 46 若Ａ是方陣，且是非奇異的，則 X=A^(P)將 inv(A) P 次方 ,如 : Y = B^(–3) Y = — — — — 分數(shù)詞冪將由Ａ的特征值決定．