matlab線性代數(shù)試驗(yàn)

中國礦業(yè)大學(xué)線性代數(shù)實(shí)驗(yàn)指導(dǎo)書線性代數(shù)MATLAB實(shí)驗(yàn)指導(dǎo)書MATLAB是MatrixLaboratory的縮寫，是一個(gè)集數(shù)值計(jì)算、圖形處理、符號(hào)運(yùn)算、文字處理、數(shù)學(xué)建模、實(shí)時(shí)控制、動(dòng)態(tài)仿真和信號(hào)處理等功能為一體的數(shù)學(xué)應(yīng)用軟件，而且該系統(tǒng)的基本數(shù)據(jù)結(jié)構(gòu)是矩陣，又具有數(shù)量巨大的內(nèi)部函數(shù)和多個(gè)工具箱，使得該系統(tǒng)迅速普及到各個(gè)領(lǐng)域，尤其在大學(xué)校園里，許多學(xué)生借助它來學(xué)習(xí)大學(xué)數(shù)學(xué)和計(jì)算方法等課程，并用它做數(shù)值計(jì)算和圖形處理等工作。我們在這里介紹它的基本功能，并用它做與線性代數(shù)相關(guān)的數(shù)學(xué)實(shí)驗(yàn)。在正確完成安裝MATLAB軟件之后，直接雙擊系統(tǒng)桌面上的MATLAB圖標(biāo)，啟動(dòng)MATLAB，進(jìn)入MATLAB默認(rèn)的用戶主界面，界面有三個(gè)主要的窗口：命令窗口(CommendWindow),當(dāng)前目錄窗口(CurrentDirectory),工作間管理窗口(Workspace)。命令窗口是和Matlab編譯器連接的主要窗口，“>>”為運(yùn)算提示符，表示Matlab處于準(zhǔn)備狀態(tài)，當(dāng)在提示符后輸入一段正確的運(yùn)算式時(shí)，只需按Enter鍵，命令窗口中就會(huì)直接顯示運(yùn)算結(jié)果。實(shí)驗(yàn)1矩陣的運(yùn)算，行列式實(shí)驗(yàn)名稱：矩陣的運(yùn)算，行列式實(shí)驗(yàn)?zāi)康模簩W(xué)習(xí)在matlab中矩陣的輸入方法以及矩陣的相關(guān)運(yùn)算，行列式。實(shí)驗(yàn)原理：介紹相關(guān)的實(shí)驗(yàn)命令和原理(1)一般矩陣的輸入(2)特殊矩陣的生成(3)矩陣的代數(shù)運(yùn)算⑷矩陣的特征參數(shù)運(yùn)算(5)數(shù)字行列式和符號(hào)行列式的計(jì)算實(shí)驗(yàn)命令表1.1Matlab的基本運(yùn)算符號(hào)運(yùn)算符號(hào)=+一*\/A,.說明賦值加減乘左除右除幕運(yùn)算轉(zhuǎn)置群運(yùn)算1矩陣的輸入Matlab是以矩陣為基本變量單元的，因此矩陣的輸入非常方便。輸入時(shí)，矩陣的元素用方括號(hào)括起來，行內(nèi)元素用逗號(hào)分隔或空格分隔，各行之間用分號(hào)分隔或直接回車。f112]例1輸入矩陣A=—103，可以在命令窗口中輸入、4—56,>>A=[112;-103;4-56]A=112-1 0 34 -5 62特殊矩陣的生成某些特殊矩陣可以直接調(diào)用相應(yīng)的函數(shù)得到，例如：zeros(m,n)生成一個(gè)m行n列的零矩陣

中國礦業(yè)大學(xué)線性代數(shù)實(shí)驗(yàn)指導(dǎo)書ones(m,n) 生成一個(gè)m行n列元素都是1的矩陣eye(n) 生成一個(gè)n階的單位矩陣rand(m,n) 生成一個(gè)m行n列的隨機(jī)矩陣magic(n) 生成一個(gè)n階魔方矩陣?yán)?隨機(jī)生成一個(gè)3x2的矩陣。>>rand(3,2)ans=0.9501 0.48600.2311 0.89130.6068 0.76213矩陣的代數(shù)運(yùn)算、f1例3設(shè)A=、f1例3設(shè)A=0-3\2—1、 (-1A,A的轉(zhuǎn)置(或transpose(A))A+B加法k*A數(shù)k乘AA*B乘法inv(A)A的逆陣AAxA的x次方A\B左除A-1BA/B右除BA-1A.*B 矩陣元素符號(hào)前加“.”，其含義是矩陣元素的群運(yùn)算。01L22,求4，A+B，AB，A2，A-1B，A.*B51,程序設(shè)計(jì)結(jié)果如下：TOC\o"1-5"\h\z>>A=[12-1;012;-364]A=1 2 -1012-3 6 4>>B=[-101;022;351]B=-1 0 1022351>>A'ans=1 0 -3216-1 2 4>>A+Bans=0200340 11 5>>A*B中國礦業(yè)大學(xué)線性代數(shù)實(shí)驗(yàn)指導(dǎo)書TOC\o"1-5"\h\zans=-4 -1 46 12 415 32 13>>AA2ans=4 -2 -1-6 13 10-15 24 31>>inv(A)*Bans=-1.0000 0.1304 1.34780 0.3478 0.26090 0.8261 0.8696>>A.*Bans=-1 0 -1024-9 30 44矩陣的特征參數(shù)運(yùn)算在進(jìn)行科學(xué)運(yùn)算時(shí)，常常要用到矩陣的特征參數(shù)，如矩陣的行列式、秩、跡、條件數(shù)等，在Matlab可以用下述命令輕松地進(jìn)行這些運(yùn)算。det(A)A的行列式rank(A)A的秩trace(A)A的跡cond(A)A的條件數(shù)size(A)輸出A的行數(shù)和列數(shù)例4求向量組(0,-1,2,3)r,(1,4,0,-1)r,(3,1,4,2)t,(-2,2,-2,0)r的秩。程序運(yùn)行結(jié)果如下：>>A二[013-2;-1412;204-2;3-120]A=013-2-1412204-23-120>>rank(A)ans=3故可知向量組的秩為3.例5判斷向量組％=[1,1,2,3]t,a2=[1,-1,1,1]t,a3=[2,0,3,3]t,叱=[3,1,5,4]t是否線性相關(guān)？中國礦業(yè)大學(xué)線性代數(shù)實(shí)驗(yàn)指導(dǎo)書f1123、由a,a,a,a所組成的矩陣A二1-101，求出A的秩或者A的行列式12 3 421353V5134J即可判斷其線性相關(guān)性，因此，在matlab命令窗口下，鍵入：A=[1123;1-101;2135;3134];rank(A)ans=3即廠(A)=3<4，故a1,a2,a3,a4線性相關(guān).123例6計(jì)算D=456.789程序運(yùn)行結(jié)果如下：>>A=[123;456;789];>>D=det(A)ans=5含有符合變量的行列式定義符號(hào)變量symsx,y定義符號(hào)變量x,ydet(A(x,y)) 計(jì)算含有符號(hào)變量x,y的行列式e、I令1 1+x1 1例7計(jì)算1 1 1+j11 1 1 1+j程序運(yùn)行結(jié)果如下：symsxyA=[1+x111;11+x11;111+y1;1111+y];det(A)ans=2*x*y^2+2*x^2*y+x^2*y^2

中國礦業(yè)大學(xué)線性代數(shù)實(shí)驗(yàn)指導(dǎo)書實(shí)驗(yàn)2線性方程組的求解，特征向量與二次型實(shí)驗(yàn)名稱：線性方程組的求解，特征向量與二次型實(shí)驗(yàn)?zāi)康模簩W(xué)習(xí)在Matlab中行列式如何計(jì)算和線性方程組的求解問題；學(xué)習(xí)在Matlab中特征向量的計(jì)算和化二次型為標(biāo)準(zhǔn)形。實(shí)驗(yàn)原理：介紹相關(guān)的實(shí)驗(yàn)命令和原理(1)齊次線性方程組和非齊次線性方程組的求解(2)特征值與特征向量的計(jì)算；(3)化二次型為標(biāo)準(zhǔn)形的方法。1在Matlab中，求解線性方程組的方法有很多，本實(shí)驗(yàn)中介紹如下命令來直接求解。也可用程序設(shè)計(jì)的方法來求解線性方程組。rref(A) A的最簡行階梯形矩陣null(A,'r')求齊次方程組null(A,'r')求齊次方程組Ax=0的基礎(chǔ)解系A(chǔ)\y 求解Ax=y的一個(gè)特解(最小二乘問題)2方陣的特征值與特征向量在矩陣對角化中和微分方程組等問題中有著廣泛的應(yīng)用?？梢苑治鼍仃嚨膶腔瘑栴}，分析二次型的正定性等，在Matlab中與本實(shí)驗(yàn)相關(guān)的命令。P二Poly(A) 求A的特征多項(xiàng)式roots(P) 求多項(xiàng)式P的零點(diǎn)[V,U]=eig(A) A的特征值與特征向量1-1101、8-243化為最簡行階梯形矩陣。1121-1-15)程序運(yùn)行結(jié)果A二[71-1101;48-243;121-1-15];rref(A)ans=1.000000-2.20000.800001.00000-6.33331.500000 1.0000-31.7333 6.1000x+2x—x—2x=0例2求齊次方程組[2x1-j-++x:=0的基礎(chǔ)解系及全部解.3x+x-2x-x=0V1 2 3 4[1 2 -1-2]解：該方程組的矩陣表示形式為2-1-1 1X=0.、31-2-1,則在matlab命令窗口下，鍵入A二[12-1-2;2-1-11;31-2-1];

中國礦業(yè)大學(xué)線性代數(shù)實(shí)驗(yàn)指導(dǎo)書null(A,'r') %這里面的r表示矩陣的秩ans=0.6000 00.2000 1.00001.0000 00 1.0000即兩個(gè)基礎(chǔ)解系分別為n10.21.故原方程通解為y=k2+k2n2(kjk2為任意常數(shù)).x+3x-2x+4x+x=7., …、2x+6x+5x+2x=5TOC\o"1-5"\h\z例3求解線性方程組112 4 54x+11x+8x+5x=31 2 3 5x+3x+2x+x+x=-2解：程序結(jié)果如下： 12345>>B=[13-2417;260525;4118053;13211-2];rref(B)B=13-2417260525411805313211-2ans=1.000000-9.50004.000035.500001.000004.0000-1.0000-11.0000001.0000-0.75000-2.2500000000所以原方程組等價(jià)于方程組x-9.5x+4x=35.5<x+4x-x=-11x-0.75x=-2.25l3 4故方程組的通解為：中國礦業(yè)大學(xué)線性代數(shù)實(shí)驗(yàn)指導(dǎo)書]-11+-2.250例4求矩陣AJ-"OI的特征值與特征向量。程序運(yùn)行結(jié)果如下：>>A=[-110;-430;102]A=TOC\o"1-5"\h\z-1 1 0-4 3 01 0 2>>E=eig(A)E=211>>[V,D]=eig(A)V=0 0.4082 0.40820 0.8165 0.81651.0000 -0.4082 -0.4082D=200010001f5001例5將矩陣A二0