matlab基礎(chǔ)練習(xí)題(帶答案)

1、For personal use only in study andresearch; not for commercial useMatlab基礎(chǔ)練習(xí)題常量、變量、表達式1、2、3、4、MATLAB中，下面哪些變量名是合法的？(A) _num (B) num_ ( C) num- 在MATLAB中，要給出一個復(fù)數(shù) z的模，B) abs(z)(C)的預(yù)定義特殊變量的是？none ( C) zero (D) exp 的內(nèi)存工作區(qū)中，存放一個英文字符 需要占用2個字節(jié)。(錯，都是(A) mod(z)下面屬于MATLAB (A) eps( B)判斷：在MATLAB 一個中文字符啊)(D) -num

2、應(yīng)該使用( double(z)函數(shù)。(D)angle(z)a2個字節(jié)需要占用1個字節(jié),)存放5、判斷：MATLAB中，i和j都是虛數(shù)單位匚1，它們之間沒有什么區(qū)別。( 對6、判斷：MATLAB中，pi代表圓周率，它等于 3.14。(錯，后面還有很多位小數(shù)7、在 MATLAB 中,若想計算的y1 = 2測(0嚴)值，那么應(yīng)該在 MATLAB的指令窗中1+V59、輸入的 MATLAB 在 MA TLAB 節(jié)。在 MA TLAB中,指令是仁2*sin(0.3*pi(1+sqrt(5) _。a = 1，b = i，貝U a占_8_個字節(jié)，b占_16_個字節(jié)，c占中,inf的含義是_無窮大_, nan的

3、含義是_非數(shù)(結(jié)果不定)數(shù)組1、2、X是一個一維數(shù)值數(shù)組, 的逆序排列輸出，應(yīng)該使用下面的(A) Xend:1( B) Xend:-1:1在MA TLAB中，A是一個字二維數(shù)組， 的命令是()。(A) class(A) ( B) sizeof(A) 在MATLAB 中，用指令 x=1:9 賦值為0,應(yīng)該在指令窗中輸入(A) x(2 7)=(0 0)(C) x(2,7)=0 0在 MATLAB中,現(xiàn)在要把數(shù)組X中的所有元素按原來次序)指令。(C) X (end:-1:1)( D) X(end:1)要獲取A的行數(shù)和列數(shù)，應(yīng)該使用的MATLAB(C) size(A)( D) isa(A)生成數(shù)組x。

4、現(xiàn)在要把x數(shù)組的第二和第七個元素都 )(B) x(2,7)=0,0(D) x(2 7)=(0 0)3、4、5、6、7、89、(A)b = -3-2 -1(B)b :=-2 -1 0 1(C)b = -5 -1 3(D)b :=-5 -2 1 4在 MA TLAB中，A =1:9，現(xiàn)在執(zhí)行如下指令L1 = (A5)，則MATLAB的執(zhí)行結(jié)果應(yīng)該是L1 =111110 00 0 。(2)A.*B(3)A(:)(4)A(1,:)*B(:,3)(1) A=B11 11(2)A(:) = 1:2:24A 3)A(1:3:7)17211923在 MATLAB 中，依次執(zhí)行以下指令

5、： clear;A=ones(3,4);A(:)=-6:5; 這時,若在指令窗中輸入指令 b=A(:,2)，那么，MATLAB輸出的結(jié)果應(yīng)該是()在MATLAB中，要求在閉區(qū)間0,5上產(chǎn)生50個等距采樣的一維數(shù)組b，請寫出具體的 MATLAB 指令 _linspace(0,5,50) _。在 MATLAB 中，A = 0:1/2:2*pi，那么 sin(A) = _ 0 1 0 -1 0。在 MA TLAB 中，A=1,2,3;4,5,6;7,8,0，B=2,1,6;8,5,2;14,2,1。寫出下面 MATLAB 語句執(zhí)行的結(jié)果：(為節(jié)省篇幅，把矩陣寫成mat2str的形式)0 0 0; 0

6、 1 0; 0 0 0 2 2 18;32 25 12;98 16 01 4 7 2 5 8 3 6 013在MATLAB中，寫出下面 MATLAB語句執(zhí)行的結(jié)果:(1) clear,A = on es(2,6)A =11111ans =1713(4)diag(diag(A)ans =1007(5) B = A(:,end:-1:1)B =21171392319151110、請編寫一段matlab程序，完成以下功能：(1) 生成一個100行，200列的二維隨機數(shù)組；A=ra nd(100,200);(2) 找出數(shù)組A中所有大于0.49且小于0.51的元素的單下標；ldx=find(A(:)0.

7、49 & A(:)0.51)(3) 數(shù)組A中滿足(2)中的條件的元素有多少個？len gth(ldx)(4) 求出數(shù)組A中滿足(2)中的條件的元素的和，并求出這些元素的平均值;sum(A(Idx), m=mea n( A(ldx)(5) 將(4)求出的平均值賦值給數(shù)組A中滿足(1)中的條件的每個元素。A(ldx)=m;11、請編寫一段matlab程序，完成以下功能：(1) 找出100到200之間的所有質(zhì)數(shù)，將這些質(zhì)數(shù)存放在一個行數(shù)組里;X=100:200;p=X(isprime(X)(2) 求出這些質(zhì)數(shù)之和；sum(p)(3) 求出100到200之間的所有非質(zhì)數(shù)之和(包括 100和200)。s

8、um(X(isprime(X)12、 y = 0.7 - 2cos2x sinx，編寫一段 matlab 程序，要求如下(1+x)(1) 在0,2二區(qū)間，每隔0.01取一 x數(shù)值，計算出相應(yīng)的 y的函數(shù)值;x=0:0.01:2*pi;y=(0.7+2*cos(x)./(1+x.A2).*s in(x)(2) 根據(jù)MATLAB計算出的數(shù)據(jù)，找出在0, 2二內(nèi)該函數(shù)的極小值的坐標。my,idx=mi n(y) x(idx)數(shù)值計算1、在MATLAB中，A是一個10X 10數(shù)組，我們把該數(shù)組看成矩陣的話，則此矩陣的行列式值=_det(A)_,此矩陣的逆矩陣(如果存在的話) =_inv(A) _。(用

9、MATLAB 的函數(shù)表示)2、求解下列方程組：，3x + 4y-7z-12w = 45x - 7y + 4z + 2w = -3x 8z -5w = 9-6x 5y - 2z 10w = -8玄 + x2 + 3x3 _ x4 = 2 嚴 _X3 +X4 =1x1 x2 2x3 2x4 二 4x1 - x2 x3 - x4 = 0A=1 1 3 -1;0 1 -1 1;1 1 2 2;1 -1 1 -1; b=2;1;4;0;x=AbA=3 4 -7 -12;5 -7 4 2;1 0 8 -5;-6 5 -2 10; b=4;-3;9;-8;xyzw=AbC 33 c L(3x +x 丫 x

10、+0.5 3、 求有理分式 R =32的商多項式和余多項式(x +2x -2 丫5x +2x +1 )n=co nv(3 0 1 0,3 0 0 0.5);d=co nv(1 2 -2,5 2 0 1);q,r=dec onv(n,d)424、 一元多項式 p=2x -3x +4x，寫出表示 p 的 MATLAB 語句 _p=2 0 -3 4 0_，求p=0的根的MATLAB語句是roots(p)_,求X = 43時p的數(shù)值的 MATLAB語句是 _ polyval(p,4.3)。程序設(shè)計1、M腳本文件和函數(shù)文件的有什么區(qū)別？最主要的差別是，函數(shù)有自己的獨立工作區(qū)，和外界的信息交換通過輸入輸出

11、參數(shù)實現(xiàn); 而腳本沒有獨立工作區(qū)，其直接使用宿主程序的工作區(qū)，并把結(jié)果留在宿主程序中。2、掌握下面程序設(shè)計中常用的函數(shù)及關(guān)鍵字：in putpausewarni ngnarginn argoutforwhileif elseifelse switch caseotherwiseendtrycatchbreakcon ti nueend633、分別用for和while循環(huán)語句計算 K = v 2i的程序。還請寫出一種避免循環(huán)的計算程i =9序。(不要用 symsum函數(shù))K=0;for i=0:63K=K+i;endK=0; i=0;while i= 10000，所以n = ceil(log(1

12、0000/0.06e-3)/log(2)10、編寫MATLAB腳本輸出“水仙花數(shù)”及水仙花數(shù)的個數(shù)。所謂“水仙花數(shù)”是一個3位數(shù)，其各位數(shù)字的立方和等于該數(shù)本身。例如15 13 53 33。c=0;for n=100:999m=fix(n/100) fix(mod(n,100)/10) mod(n,10);if n=sum(m.A3),fprin tf(%i = %iA3 + %jA3 + %iA3n, n,m)；c=c+1;endendfprintf(n 共 %i 個水仙花數(shù) n,c)輸出結(jié)果：153 = 1A3 + 5A3 + 3A3370 = 3A3 + 7A3 + 0A3371 = 3

13、A3 + 7A3 + 1A3407 = 4人3 + 0A3 + 7A3共4個水仙花數(shù)繪圖1、MATLAB中，繪制三維曲面圖的函數(shù)是()(A)surf( B) plot( C) subplot(D)plot32、MATLAB中，要繪制三維空間曲線，應(yīng)該使用()函數(shù)。(A) polar( B) plot( C)subplot(D)plot33、在MATLAB中，A是一個1000行2列的二維數(shù)值數(shù)組，現(xiàn)在要把A的第一列數(shù)據(jù)作為橫坐標，把 A的第二列數(shù)據(jù)作為縱坐標，畫出一條曲線，試寫出相應(yīng)的MATLAB語句plot(A(:,1),A(:,2) 。4、MATLAB繪圖指令中的_subplot_指令允許用

14、戶在同一個圖形窗里布置幾個獨立的子圖。5、 編寫一-段matlab程序，繪制出二元函數(shù)z二 一iy二維網(wǎng)線圖，要求如下xy(1)x，y的取值范圍為-9乞x乞9，-9 syms x limit(1+xF(1/x),0)ans =expsin 3x(2) 求函數(shù)y在X =0處的極限;tg5x syms x limit(si n(3*x)/ta n(5*x),0)ans =3/5(3)nx求函數(shù)y-在趨向正無窮處的極限;3x syms x n limit (n *x/3Ax,i nf)ans =0求函數(shù)ln2xy -x在趨向正無窮處的極限; syms x limit(log(x)A2/xA3, in

15、f)ans =07、求導(dǎo)數(shù)：1(1) 求函數(shù)y 2 _3x 3的50階導(dǎo)數(shù)；x syms x y=1/xA2-3*x+3; diff(y,50)ans =1551118753287382280224243016469303211063259720016986112000000000000/xA52cta(2) 求函數(shù)y二as in be t 在t = b處的3階導(dǎo)數(shù)；syms t a b cy=a*si n(b*exp(cAt)+tAa);simple(subs(diff(y,t,3),t,b)8、求不定積分：Ex2 dxa -xdx syms x in t(1/si n(xF3)ans =-

16、1/2/si n(x)A2*cos(x)+1/2*log(csc(x)-cot(x) syms x a in t(1/(aA2-xA2)ans =-1/2/a*log(a-x)+1/2/a*log(a+x) syms x in t(sqrt(xA2-3)-sqrt(xA2+3)/sqrt(xA4-9)ans =(xA4-9)A(1/2)/(xA2-3)A(1/2)/(xA2+3)A(1/2)*asi nh(1/3*3A(1/2)*x)-1/(xA2+3)A(1/2)*(xA4-9)A(1/2)/(xA2-3)A(1/2)*log(x+(xA2-3)A(1/2)9、求定積分及廣義積分f Jx2

17、+adx-2sin2xcos2 xdx-312 x2221 十 x y dydx1 syms x a in t(sqrt(xA2+a),-2,2)ans =2*(4+a)A(1/2)+1/2*a*log(2+(4+a)A(1/2)-1/2*a*log(-2+(4+a)A(1/2) syms x in t(si n(x)A2*cos(x)A2,-pi,pi)ans =1/4*pi syms x y in t(i nt(xA2+yA2,y,1,xA2),1,2)ans =1006/10510、求下面的積分，給出50位精度的數(shù)值:sin2 x sin2 y dydx syms x y J=i nt(

18、i nt(si n(x)A2+si n(y)A2,y,1,xA2),1,2); vpa(J,50)ans =11、級數(shù)求和:n +遲 n()2.1540459589705316265997501755762001048498664176916_n、3n 1 z -1n 42k 1 壬 2 i x -1 z.7 2k 1 x 1 syms z n symsum(z-1)A n/(n A2*2A n),n ,1,i nf)ans =(1/2*z-1/2)*hypergeom(1, 1, 1,2, 2,1/2*z-1/2) syms z n symsum(3* n+1)*(z-1)A n,n ,1,

19、i nf)ans =(4*z-4)*(-1/(z-2)+3/4/(z-2)A2*(z-1) syms z n symsum( n*(_1)A( n+1)*zA n,n ,1,i nf)ans =z/(z+1)A2 syms x positive syms k simple(symsum(2/(2*k+1)*(x-1)/(x+1)A(2*k+1),k,0,i nf)ans =log(-(1+(xA2-2*x+1)/(xA2+2*x+1)A(1/2)/(-1+(xA2-2*x+1)/(xA2+2*x+1)A(1/2)僅供個人用于學(xué)習(xí)、研究；不得用于商業(yè)用途For personal use only in study and research; not for commercial use.Nur f u r den pers?nlichen f u r Studien, Forschung, zu kommerziellen Zwecken ve