【正文】 矩陣而得、以 G 為子矩陣用復(fù)制函數(shù)(repｍａｔ)生成得大矩陣Ｈ。D］ E =１０00000８〉〉 F=(２：5，:)？?？ F=(2:5，：)Erroｒ: Exprｅｓｓion or statｅmenｔ is incorreｃｔ—pｏssiｂｌy unbalaｎcｅd(, ｛, or ［、＞〉 F=E（2：5,：)F =00００〉〉 G＝respace(E,3,４）？？? Undefiｎeｄ manｄ/fｕnｃtｉon “reｓpace’、〉G＝rｅspaｃe(F，３,4）？?? Ｕndｅfiｎeｄ manｄ／fuｎcｔion ”reｓｐaｃe“、〉 G=ｒeｓhapｅ(Ｆ,3,4）G =00１６00〉〉 Ｈ=repmat（G,2,2)Ｈ =00０５０８０0７0０0３１０１０0１000０2)矩陣得運(yùn)算 ? 矩陣得加減、數(shù)乘與乘法已知矩陣: A=[1 ２－１], A =１ 〉 B=［1 ０2] Ｂ =0〉 A+B ans =０４１ 〉〉 2＊A ａns =４2 〉 2*A－３*B ans =8 〉〉 A＊Ｂ aｎｓ =４４—2 ? 矩陣得逆矩陣 〉〉 foｒｍat ｒat。3 1 4；2 2 １]，b=[１ １ 2],d=b’ a =１３２b ＝１d =１１ c１=b*ｉnv(ａ）,c2＝b/a ｃ1 =6／７3/7－4/7c２ ＝6/7３/７－4／7 c3=inv（a)*d, c4=ab c3 ＝２/74/7？？? Erｒor usinｇ ==〉 ｍｌｄivｉｄｅ Matrｉx dimensions ｍust agree、〉 ｃ３＝inv(ａ)*ｄ，ｃ4=ad ｃ3 =2/7－4/7c4 =2/74/７練習(xí): 按下列要求求出各種得矩陣運(yùn)算得值求矩陣得秩、特征值與特征向量、矩陣得乘冪與開方。4 5 6。Ａ2＝ｒeshaｐｅ([10:18],3，3)Ａ2 ＝1０１３１61４1７１215〉〉 Ｔ1(:，:，1）＝ｏneｓ(3)；T1（：,：,2)=ｚeros(3)T1（:,：,１）=１１１Ｔ１(:，:,2)=０000000０0 〉 T２=ones（3，3,2）Ｔ2（：,:，1)=１１１ T2（:,:，２)=１１１ T３＝ｃat(3,A1，Ａ2)，T4=repmaｔ（A1,［1,1,２]）Ｔ3(：,:,１)=Ｔ3(:,：,２)=1０1３1６1７１5T4(：,:，1）=２６８T4（:,:,2)＝７８2）多維數(shù)組得創(chuàng)建數(shù)組運(yùn)算用小圓點(diǎn)加在運(yùn)算符得前面表示,以區(qū)分矩陣得運(yùn)算?！怠?A＝［１：6]。 Ｃ１＝A＋B，C2=Ａ—B Ｃ１ =２Ｃ2 =0３ C３=A、＊B,C4＝B、／A，C5=A、B C3 =６ C4 =１、0０0０0、5000０、3３3３0、25000、200０0、１667 C５ ＝ 0000０、５0000、3３325０00、200００、１6６7關(guān)系運(yùn)算或邏輯運(yùn)算得結(jié)果都就是邏輯值.＞〉 I=A〉3，C６＝A（I）I =000１１ C６ =４６ 〉〉 A1=Ａ3,I2=A1＆A A1 =—2—10２I2 =１１0１〉〉 I3＝~I I3 =１0０0 1)字符串得創(chuàng)建＞〉 S1=”Iｌike MAＴＬＡB’ S1 = Ilike MATLＡＢ ＞＞ Ｓ2=“I＇’m a ｓｔuenｔ、” S２ = Ｉ“ｍ a stuｅnt、〉 S３＝［S2,”ａnd’,Ｓ１] Ｓ３ ＝ I“m a sｔｕent、anｄＩliｋe MATＬAB 2）求字符串長(zhǎng)度 〉〉 length(S１)aｎｓ = siｚe（Ｓ１)aｎs =13）字符串與一維數(shù)值數(shù)組得相互轉(zhuǎn)換 CS1=abs（S1)CS1 =1０1３26５８46６ 〉 ＣS２=doublｅ(Ｓ１)CS2 =１０8１05１01３27７6５７6 chａｒ(CS2）anｓ = Iliｋe MＡＴLＡB setsｔr（CS2）ans = Ｉlikｅ MAＴLAB 練習(xí)：用ｃｈaｒ()與向量生成得方法創(chuàng)建如下字符串ＡaBbCｃ、XxＹｙZz、〉 S1=65：90?！怠?C=[S1。C=C(：）’。char(S3)ａｎs = ＡaＢbCcDdEeＦfGgHhIiJjKkLlMｍNnOoＰpQqRrＳsTtＵｕVvWｗXxＹyＺz 實(shí)驗(yàn)三ＭAL ＴAＢ Ｂ 數(shù)值運(yùn)算 一、實(shí)驗(yàn)?zāi)康谜莆?MATLAB 得數(shù)值運(yùn)算及其運(yùn)算中所用到得函數(shù)，掌握結(jié)構(gòu)數(shù)組與細(xì)胞數(shù)組得操作。2）多項(xiàng)式插值與擬合。4）結(jié)構(gòu)數(shù)組與細(xì)胞數(shù)組。在 MＡTLAB 中，: 在 MATLＡB 中表示為 s=[1 3５ ０ 9] 2）多項(xiàng)式得加減法相當(dāng)于向量得加減法,但必須注意階次要相同。如多項(xiàng)式與多項(xiàng)式相加?！担?s1=[2 ３ 1１] ＞ s2=［１ ３ —5 4 7］ 〉 s３=cｏnv（ｓ1，ｓ2） s４=ｄecoｎv（s3,s１)答?！?ｓ1＝[2 4 ２］ roots(s1）答。s１ =２3 aｎs ＝９0 x ＝１６９1０ y =Colｕmns １ thrｏugh ８19３３5２５798８６1285Cｏlｕｍns 9 throuｇｈ 1０ 練習(xí):求得“商”及余數(shù).〉 s1=[１ ０ 1]；s２=［1 3］。附表 1－1 X 1 2 3 4 5 ６ 7 8 9 10 Ｙ １6 ３2 71分別用擬合(二階至三階)與插值（線性與三次樣條)得方法來估算Ｘ=5 時(shí) Y 得值?！?x=1：10 〉 y＝［１6 32 70 1４2 ２6０ 4３６ 682 1０10 1342 1960］ 〉 p１=ｐｏｌyｆｉt（x,ｙ,1)％一階擬合 y1=pｏｌｙval(p1，5）%計(jì)算多項(xiàng)式 P1 在ｘ＝９、5 得值 答。＞〉 x=1:2:９ 〉〉 ｄifｆ(x)答?！?x=linspace(0，2＊pi，1００)。〉〉 x=ones(1,1０）＞ cumsum（x)＞ x＝linspaｃe(0,pｉ，100）； ＞ y=sin(x)； 〉〉 trａpｚ(x,y)〉 ｐ＝cumsum（y）； p（100)*ｐi／(100—1)答。地圖得數(shù)據(jù)見附表 1—２(單位ｍm).附表 1—２ X 7 ５ 13 1７、5 34 ４0、５ 4５ 48 56 61 65 75 8０、５ ９1 Y1 ４4 ４5 47 534 ４1 ４5 ４6 Y２ 4４ 59 7１８ 1 １8 續(xù)表 X 9６ 101 104 １05 115 118 125 1 ３5 1４2 146 150 157 158 Ｙ1 ４３ ３7 3３ 28 32 6５ 55 ５4 5２ ５0 ６6 66 ６8 Ｙ２ 121 124 １ 2１ １ 2１ 12１ 116 1 ２２ 83 81 82 86 85 ６8 提示：由高等數(shù)學(xué)得知識(shí),一條曲線得定積分就是它與 x 軸所圍成得面積,那么兩條曲線所圍成得面積可由兩條曲線得定積分相減得到。〉〉 stuｄent、ｎａｍｅ=’Jack”?！怠?ｓtuｄent（2)、ｎamｅ=’Lucｙ“； stｕｄent 或者用 sｔruct 函數(shù)創(chuàng)建。答?！?student(１)、subjｅｃt=[]%添加 subject 域并賦予空值 〉〉 ｓtuｄｅnt(1）、score=［］ 〉(studeｎg)〉fｉeldnａｍes(sｔｕdeng)〉fiｅldnaｍｅs(sｔudent）〉〉 getfielｄ(stｕdｅnt,{2｝，＇name＇)〉〉 stuｄenｔ=ｒｍfｉeld(stuｄｅnt,”subｊect’)%刪除 subｊeｃt 域 〉〉 student＝setｆiｅlｄ（student,{1},“scoｒｅ’，90）。結(jié)構(gòu)數(shù)組得大小為 22?！?A＝{’Ｈｏw are ｙou！”,oｎｅs(3)。%直接創(chuàng)建〉〉 B（1,１)={’Hellｏ woｒld”}；%由各個(gè)細(xì)胞元素創(chuàng)建 B(1，2）＝{mａgｉｃ(3)｝。答或者用ｃell 函數(shù)先創(chuàng)建空得細(xì)胞數(shù)組，然后再給各個(gè)元素賦值 c=ceｌl（１，２)； ｃ（1，1)＝{’Hｅｌlｏ world’}?！?ans1＝A(１,1)＞ ans2＝A(１,１)〉 whos ans1 ａnｓ2 〉〉 cｅｌlｄisｐ（Ａ)〉 a１=Ａ{2,1}（1,2)〉 [a２ a3］=ｄeａｌ(A{1:2}）答；ａns1 ＝’How arｅ ｙou！’ ans2 =’Ｈｏw are yｏu!’NａｍeSizeBｙteｓClasｓAttributｅｓａｎｓ11x１８4celｌａｎs２１x18４ｃellＡ{1，1｝ = Hｏw are ｙoｕ！A｛2，1} =A｛１，2} =１Ａ{２，2}{1｝ = ceｌl a1 =a2 = Ｈow arｅ you!a3 =１２實(shí)驗(yàn)四MA LT ＡB B 符號(hào)運(yùn)算一、實(shí)驗(yàn)?zāi)康谜莆辗?hào)變量與符號(hào)表達(dá)式得創(chuàng)建,掌握MＡLＴAB得symbｏl工具箱得一些基本運(yùn)用。2）符號(hào)微積分運(yùn)算．3））、實(shí)驗(yàn)步驟符號(hào)運(yùn)算得引入在數(shù)值運(yùn)算中如果求，則可以不斷讓得讓ｘ趨近0，一球得表達(dá)式趨近什么數(shù),但終究不能令 x=０，因?yàn)樵跀?shù)值運(yùn)算中 0 不就是能作除數(shù)得。輸入如下命令: 〉 ｆ=sym(’sｉn(pi*x)/x＇） lｉｍｉｔ(f,’x“,0)答；f ＝ sin(ｐi*ｘ)/x ans = pi 2 符號(hào)常量、符號(hào)變量、符號(hào)表達(dá)式得創(chuàng)建1）使用 sym（）創(chuàng)建 輸入以下命令,觀察 Workspace 中 A、B、ｆ就是什么內(nèi)性得數(shù)據(jù),占用多少字節(jié)得內(nèi)存空間。A ＝ 1 B ＝ x f ＝ 2*x＾2＋3＊x１ ｆ1 = 1＋２ f2 = 3 f2 = ３ f4 = 2*ｘ+3 x =f4 = 5 通過瞧 MＡTLAB 得幫助可知,sym()得參數(shù)可以使字符串或就是數(shù)值類型,無論就是哪種類型都會(huì)生成符號(hào)類型數(shù)據(jù)。x = x y = y z ＝ z f１ ＝ ｘ^2+2*ｘ＋１ ｆ2 ＝ ｅｘp（ｙ)+eｘｐ(z）^2 f3 = x＾2＋2＊x＋1+exp(y)+ｅｘp（z)^２ 3 符號(hào)矩陣創(chuàng)建〉 syms a1 a2 ａ3 a４ 〉〉 A=[ａ１ a2。A = [ a１，a2］ [ a３, a4］ ａns ＝ a1 ans = ａ2 4 符號(hào)算術(shù)運(yùn)算1）符號(hào)向量相乘、相除 符號(hào)量相成與數(shù)值量相乘一樣,分成矩陣乘與數(shù)組乘。B＝sｙｍ([3 4 5］)。ｂ 3]。〉＞ Ｃ1＝A*B，C２＝Ａ、*B 〉 C3＝AB,C4=A、/B 答；c1 =３5 ｃ2 = 5/7 Ｃ1 = ［ 15, 20, 2５] Ｃ２ ＝ [ 3／5, 4/5，1] C1 ＝ [ 10*a＋2*a＊b,5*b+ａ^2］ ［2＊a*ｂ+６＊b,ｂ＾２+3＊a] C2 = [10*a，a＊b］ [ 2*b^2，3*a] C３ ＝ [２*a＊(ｂ—3)/（—15+a＊b)，（a^2３*b)／(—１5+a*b)] [2*ｂ*(a5）/(１５+a＊ｂ), —（5*aｂ^2）/(15+a＊b)] C4 = ［ ５/2/a，ａ／b] [1／2, 3／a] 2）符號(hào)數(shù)值任意精度控制與運(yùn)算 任意精度得 VPA 運(yùn)算可以使用命令 digitｓ(設(shè)定默認(rèn)得精度）與ｖpa(對(duì)指定對(duì)象以新得精度進(jìn)行計(jì)算)來實(shí)現(xiàn)。3）符號(hào)類型與數(shù)值類型得轉(zhuǎn)換 使用命令sｙm可以把數(shù)值型對(duì)象轉(zhuǎn)換成有理數(shù)性符號(hào)對(duì)象,命令vpa可以講數(shù)值型對(duì)象轉(zhuǎn)換為任意精度得 VPA 型符號(hào)對(duì)象．使用 dｏublｅ，nｕmeric 函數(shù)可以將有理數(shù)型與 VＰA 型符號(hào)對(duì)象轉(zhuǎn)換成數(shù)值對(duì)象、〉＞ clear 〉 a1=sym(＇2*ｓqrt(5)+pi’)