Matlab數(shù)學(xué)實(shí)驗(yàn)一2015答案版

1、Matlab數(shù)學(xué)實(shí)驗(yàn)一matlab初體驗(yàn)、實(shí)驗(yàn)?zāi)康募耙饬x1熟悉MATLA歆件的用戶環(huán)境；2 了解MATLA歆件的一般目的命令；3掌握MATLA敢組操作與運(yùn)算函數(shù)；通過(guò)該實(shí)驗(yàn)的學(xué)習(xí)，使學(xué)生能熟悉 matlab的基礎(chǔ)應(yīng)用，初步應(yīng)用MATLA瞅件解決一些簡(jiǎn)單問(wèn)題。、實(shí)驗(yàn)內(nèi)容1 .認(rèn)識(shí)matlab的界面與基本操作2 .了解matlab的數(shù)據(jù)輸出方式(format)3 、 MATLAB軟件的數(shù)組(矩陣)操作及運(yùn)算練習(xí)；三、實(shí)驗(yàn)任務(wù)根據(jù)實(shí)驗(yàn)內(nèi)容與步驟，完成以下具體實(shí)驗(yàn)，要求寫出實(shí)驗(yàn)報(bào)告(實(shí)驗(yàn)?zāi)康囊粏?wèn)題一原理一算法與編程一 計(jì)算結(jié)果或圖形一心得體會(huì))完成如下題目，并按照實(shí)驗(yàn)報(bào)告格式與要求填寫實(shí)驗(yàn)報(bào)告1.在c

2、ommand window中分別輸入如下值，瞧它們的值等于多少，并用matlab的help中查詢這些缺省預(yù)定義 變量的含義，用中文寫出它們的意義。i j eps inf nan pi realmax realmin2 .分別輸入一個(gè)分?jǐn)?shù)、整數(shù)、小數(shù)等 ，(如:a=1/9),觀察顯示結(jié)果，并使用format函數(shù)控制數(shù)據(jù)的顯示格式，如: 分另1J輸入 format short、 format long、 format short e、 format long g、 format bank、format hex 等，然后再 在命令窗口中輸入 a,顯示a的值的不同形式，并理解這些格式的含義。3 .測(cè)試

3、函數(shù)clear、clc的含義及所帶參數(shù)的含義 (利用matlab的help功能)。4、寫出在命令窗口中的計(jì)算步驟與運(yùn)行結(jié)果。計(jì)算9一也一也；81 (log(pi)+log(pi)/10g(10)-exp(1、2)A2/81ans =0、 0348(2) x=2;y=4; z=xA2+exp(x+y)-y*10g(x)-3 z =401、 65621 3(3)輸入變重a 5.3,b,在工作仝間中使用who,whos,并用save命令將變重存入D:exe01、2 5mat”文件。測(cè)試 clear命令，然后用load命令將保存的D:exe01、mat”文件載入 a=5、35、 3000 b=1 3;

4、 2 5 b =1 32 5 whoYour variables are:a b whosBytes Class8 double array32 double arrayName Sizea1x1b2x2Grand total is 5 elements using 40 bytes save D:exe01 clear 清除內(nèi)存中在全部變量 load D:exe015、對(duì)矩陣，求其行列式（det）、逆矩陣（inv）、矩陣的特征值與特征向量（eig）、矩陣的秩（rank）、矩陣的行最簡(jiǎn)形 （rref） 、以該矩陣為系數(shù)矩陣的線性方程組Ax=0 的通解 （null）;42213已知 A 305

5、, B 2 015321(1) 計(jì)算矩陣A 的行列式的值det(A) A=4,-2,2;-3,0,5;1,5,3; det(A)43 ,在 MATLAB 命令窗口中建立A、 B 矩陣并對(duì)其進(jìn)行以下操作17-70-#0130115ans =-15811 2T(2) 分別計(jì)算下列各式:2A B,A* B, A.* B,AB 1,A 1B,A2,AT A=4,-2,2;-3,0,5;1,5,3;B=1,3,4;-2,0,-3;2,-1,1; 2*A-Bans = A*B ans =1210247-14-7-30-8 A、 *Bans =4-6860-152-53 A*inv(B)ans =4-31-

6、205253-0、 0000-0、 00002、-2、 7143-8、 0000-8、2、 42863、 00002、000014292857 inv(A)*Bans =0、 48730、 41141、0、 3671-0、 43040、-0、 10760、 24680、000000000000 A*Aans =2424-7319-81336 Aans =在MATLAB中分別利用矩陣的初等變換及函數(shù)rank、函數(shù)inv求下列矩陣的秩1632(1) A 354 0 , 求 rank(A)=?111 2 4 A=1,-6,3,2;3,-5,4,0;-1,-11,2,4; rank(A) ans =3

7、350112001(2) B, 求 B1。1020 1202 B=3,5,0,1;1,2,0,0;1,0,2,0;1,2,0,2 inv(B)ans =2、 0000-1 、 0000-1 、 0000 0-4、 00002、 50002、 0000-0、 5000-0、 00000、 00000、 5000 0-1 、 00000、 50000、 50000、 5000在 MATLAB 中判斷下列向量組就是否線性相關(guān),并找出向量組1 (113 2)T ,2 ( 11 13)T, 3 (5 289)t, 4 ( 1317) T中的一個(gè)最大線性無(wú)關(guān)組。 a1=1 1 3 2a2=-1 1 -1

8、3a3=5 -2 8 9a4=-1 3 1 7A= a1, a2 ,a3 ,a4 ;R jb=rref(A)a1 =1132-11-13 a3 =5-289 a4 =-13171 、 0000001、090901、 000001、7879001 、 0000-0、06060000jb =123 A(:,jb)ans =1-1511-23-18239在 MATLAB 中判斷下列方程組解的情況,若有多個(gè)解,寫出通解。x1 x2 4x3 2x4 0(1) x1 x2 x3 2x4 03x1 x2 7x3 2x4 0x1 3x2 12x3 6x40 A=1,-1,4,2;1,-1,-1,2;3,1,7

9、,-2;1,-3,-12,6; rank(A)ans =3 rref(A)ans =1000010-200100000 A=1,-1,4,2;1,-1,-1,2;3,1,7,-2;1,-3,-12,6; format ratn=4;RA=rank(A)RA =3 if(RA=n)fprintf(% 方程只有零解)elseb=null(A,r)endb =0201 syms k X=k*b2*k2x13x2x34x12x24x353x18x22x3134x1x29x36(2) A=2 3 1;1 -2 4;3 8 -2;4 -1 9;b=4 -5 13 -6;B=A b; n=3; RA=ran

10、k(A)RA =2 RB=rank(B)RB =2rref(B)ans =102-101-1200000000 format ratif RA=RB&RA=n % 判斷有唯一解X=Abelseif RA=RB&RA a1=inv(A)a1 =-3/21/21/201/20-21/21 P,R=eig(A)P =-985/1393-528/2177379/125700379/419-985/1393-2112/2177379/1257-100020002A 的三個(gè)特征值就是:r1=-1,r2=2,r3=2 。三個(gè)特征值分別對(duì)應(yīng)的特征向量就是P1=1 0 1;p2=1 0 4;p3=1 3 122

11、2化方陣A 254 為對(duì)角陣。 A=2 2 -2;2 5 -4;-2 -4 5;P,D=eig(A)-0、29810、 89440、 3333-0、5963-0、 44720、 6667-0、74540-0、 66671、00000001、 000000010、 0000 B=inv(P)*A*PB =1 、0000-0、00000、00000、00001 、00000、0000-0、 0000010、 0000程序說(shuō)明:,D 與 A 相似。所求得的特征值矩陣D 即為矩陣A 對(duì)角化后的對(duì)角矩陣求一個(gè)正交變換,將二次型f2225x1 5x2 3x3 2x1 x2 6x1 x36X2X3化為標(biāo)準(zhǔn)型。 A=5 -1 3;-1 5 -3;3 -3 3; syms y1 y2 y3y=y1;y2;y3;P,D=eig(A)P =881/2158985/1393-780/1351-881/2158985/1393780/1351-881/10790-780/1351 x=P*y x =(6A(1/2)*y1)/