【正文】
本算法 22 實(shí)驗(yàn)一 MATLAB中矩陣與多項(xiàng)式的基本運(yùn)算實(shí)驗(yàn)任務(wù)1.了解MATLAB命令窗口和程序文件的調(diào)用。2.熟悉如下MATLAB的基本運(yùn)算:① 矩陣的產(chǎn)生、數(shù)據(jù)的輸入、相關(guān)元素的顯示;② 矩陣的加法、乘法、左除、右除;③ 特殊矩陣:?jiǎn)挝痪仃?、?”矩陣、“0”矩陣、對(duì)角陣、隨機(jī)矩陣的產(chǎn)生和運(yùn)算;④ 多項(xiàng)式的運(yùn)算:多項(xiàng)式求根、多項(xiàng)式之間的乘除。4 5 6。...B=[1 3 5。4 6 8]。另外A\B和inv(A)*B 應(yīng)該也是沒(méi)有定義的,但還是有結(jié)果,因此,我們有必要檢查運(yùn)算結(jié)果,鑒于軟件在精度方面不可避免的誤差。p=[3,0,4,5,7,5]。 whoYour variables are:A B Q X ans f p r x whos Name Size Bytes Class Attributes A 3x3 72 double B 3x3 72 double Q 1x2 16 double X 5x1 80 double plex ans 3x3 72 double f 1x1 170 sym p 1x6 48 double r 1x4 32 double x 1x1 126 sym A=fix(10*rand(1,3)),B=fix(20*rand(3)),C=fix(10*rand(3,2)),D=[3]diag(A),diag(B),diag(C),diag(D)A = 8 0 3B = 5 18 2 16 3 2 8 5 17C = 5 8 5 6 1 3D = 3ans = 8 0 0 0 0 0 0 0 3ans = 5 3 17ans = 5 6 ans = 3 disp(39。)。分號(hào)”影響,size輸出矩陣的行數(shù)和列數(shù),length輸出行列數(shù)較大值實(shí)驗(yàn)二 MATLAB繪圖命令實(shí)驗(yàn)任務(wù)熟悉MATLAB基本繪圖命令,掌握如下繪圖方法:1.坐標(biāo)系的選擇、圖形的繪制;2.圖形注解(題目、標(biāo)號(hào)、說(shuō)明、分格線)的加入;3.圖形線型、符號(hào)、顏色的選取。y1=x.*x。subplot(2,2,1)。39。r39。y=x^{2}39。y=x^{3}39?!rid off。loglog(x,y1,39。,x,y2,39。),title(39。),grid on。semilogx(x,y1,39。,x,y2,39。),xlabel(39。),ylabel(39。)。semilogy(x,y1,39。,x,y2,39。),text(50,8000,39。),text(40,100000,39。)【結(jié)果截圖】【程序】t=0::2*pi。subplot(2,2,1)。bar(t,r,39。),axis([0,1,]),title(39。)subplot(2,1,2)。c39。stairs39。subplot(1,2,1),mesh(x,y,z)subplot(1,2,2),contour(x,y,z)【結(jié)果截圖】實(shí)驗(yàn)三 MATLAB程序設(shè)計(jì)實(shí)驗(yàn)任務(wù)用一個(gè)MATLAB語(yǔ)言編寫(xiě)一個(gè)程序:輸入一個(gè)自然數(shù),判斷它是否是素?cái)?shù),如果是,輸出“It is one prime”,如果不是,輸出“It is not one prime.”。最好能具有如下功能:①設(shè)計(jì)較好的人機(jī)對(duì)話(huà)界面,程序中含有提示性的輸入輸出語(yǔ)句。如果操作者希望停止這種判斷,則可以退出程序。例:輸入 6, 因?yàn)?不是素?cái)?shù)。 【程序】:a=input(39。)。1 is not a prime or a posite number.39。%d is a prime.\n39。m=length(f)。%d is not a prime,%d=%d39。*%d39。 end fprintf(39。)。請(qǐng)輸入一個(gè)整數(shù)(end in 0):39。end【運(yùn)行結(jié)果】:請(qǐng)輸入一個(gè)整數(shù)(end in 0):11 is not a prime or a posite number.請(qǐng)輸入一個(gè)整數(shù)(end in 0):22 is a prime.請(qǐng)輸入一個(gè)整數(shù)(end in 0):3535 is not a prime,35=5*7請(qǐng)輸入一個(gè)整數(shù)(end in 0):0實(shí)驗(yàn)四 MATLAB的符號(hào)計(jì)算與SIMULINK的使用實(shí)驗(yàn)任務(wù)1.掌握MATLAB符號(hào)計(jì)算的特點(diǎn)和常用基本命令;2.掌握SIMULINK的使用。c d]。2,9]。x^3+5*x^29*x+1739。 den=[1 , 1 , 4] 。tp=spline(y , t , max(y)) %計(jì)算峰值時(shí)間ymaxmax(y) %計(jì)算峰值