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

Matlab基礎(chǔ)練習(xí)題

常量、變量、表達(dá)式

1、MAHAB中，下面哪些變量名是合法的？()

(A)_num(B)num(C)num-(D)-num

2、在MATLAB中，要給出一個(gè)復(fù)數(shù)z的模，應(yīng)當(dāng)運(yùn)用()函

數(shù)。

(A)mod(z)(B)abs(z)(C)double(z)(D)

angle(z)

3、下面屬于MATLAB的預(yù)定義特殊變量的是？()

|(八)eps(B)none(C)zero(D)exp批注［王士星1］：本題有問(wèn)題.出胞的本意應(yīng)當(dāng)是選擇

eps,但其實(shí)eps是一個(gè)函數(shù),并非變量。

4、推斷：在MATLAB的內(nèi)存工作區(qū)中，存放一個(gè)英文字符'a'須要

占用1個(gè)字節(jié)，存放一個(gè)中文字符‘啊'須要占用2個(gè)字節(jié)。(

錯(cuò)，都是2個(gè)字節(jié))

5、推斷：MATLAB中，i和j都是虛數(shù)單位Q,它們之間沒(méi)有什么

區(qū)分。(對(duì))

6、推斷：MATLAB||?,pi代表圓周率，它等于3.14。(錯(cuò)，后j

很多位小數(shù))

7、在MATLAB中，若想計(jì)算的.=型嗎型值，那么應(yīng)當(dāng)在MATLAB

的指令窗中輸入的MATLAB指令是

yl=2*sin(0.3*pi)/(1+sqrt⑸)0

8、在MATLAB中，a=1,b=i,貝I」a占_8一個(gè)字節(jié)，b占」6_個(gè)字

節(jié)，k占字節(jié)。批注［王士星2］:哪里來(lái)的c?

9、在MATLAB中，inf的含義是—無(wú)窮大nan的含義是—非數(shù)(結(jié)

果不定)—。

數(shù)組

]、在MATLAB中，X是一個(gè)一維數(shù)值數(shù)組，現(xiàn)在要把數(shù)組X中的全部

元素按原來(lái)次序的逆序排列輸出，應(yīng)當(dāng)運(yùn)用下面的()指

令。

(A)X[end:l](B)X[end：-l:l](C)X(end：-l:l)(D)

X(end:l)

2、在MATLAB中，A是一個(gè)字二維數(shù)組，要獲得A的行數(shù)和列數(shù)，應(yīng)

當(dāng)運(yùn)用的MATLAB的叮囑是()。

(A)class(A)(B)sizeof(A)(C)size(A)(D)isa(A)

3、在MATLAB中，用指令x=l:9生成數(shù)組x。現(xiàn)在要把x數(shù)組的其次

和第七個(gè)元素都賦值為0,應(yīng)當(dāng)在指令窗中輸入()

(A)x([27])=(00)(B)x([2,7])=[0,0]

(C)x[(2,7)]=[00](D)x[(27)]=(00)

4、在MATLAB中，依次執(zhí)行以下指令：clear;A=ones(3,4);

A(：)=[-6:5];這時(shí),若在指令窗中輸入指令b=A(:,2),,那么，

MATLAB輸出的結(jié)果應(yīng)當(dāng)是()

(A)b=-3-2-1(B)b=-2-101

(C)b=-5-13(D)b=-5-214

5、在MATLAB中，A=1:9,現(xiàn)在執(zhí)行如下指令LI=~(A〉5),則MATLAB

的執(zhí)行結(jié)果應(yīng)當(dāng)是LI=111110

000

6、在MATLAB中，要求在閉區(qū)間［0,5］上產(chǎn)生50個(gè)等距采樣的一維

數(shù)組b,請(qǐng)寫出具體的MATLAB指令1inspaco(0,5,50)°

7、在MATLAB中,A=[0:1/2:2]*pi,那么sin(A)=

-----0

8、在MATLAB中，A=[l,2,3;4,5,6;7,8,0],

B=［2,1,6;8,5,2;14,2,1］。寫出下面MATLAB語(yǔ)句執(zhí)行的結(jié)果:（為

節(jié)約篇幅，把矩陣寫成mal2str的形式）

⑴A==B00;010;000]

(2)A.*B218;322512;98160]

⑶A(:)'47258360]

(4)A(l,:)*B(:,3)13

9、在MATLAB中，寫出下面MATLAB語(yǔ)句執(zhí)行的結(jié)果:

(1)clear,A=ones(2,6)

A=

111111

111111

(2)A(:)=1:2:24

A=

59131721

3711151923

(3)A([l:3:7])

ans=

1713

(4)diag(diag(A))

ans=

10

07

(5)B=A(:,end：-l:1)

B二

211713951

2319151173

10.請(qǐng)編寫一段matlab程序，完成以下功能：

(1)生成一個(gè)100行，200列的二維隨機(jī)數(shù)組；

A=rand(100,200);

<2)找出數(shù)組A中全部大于0.49且小于0.51的元素的單下標(biāo)；

Idx=find(A(:)>0.49&A(:)<0.51)

(3)數(shù)組A中滿足(2)中的條件的元素有多少個(gè)？

length(Idx)

(4)求出數(shù)組A中滿足(2)中的條件的元素的和，并求出這些

元素的平均值；

sum(A(Idx)),m=mean(A(Idx))

(5)將(4)求出的平均值賦值給數(shù)組A中滿足(1)中的條件的

每個(gè)元素。

A(Idx)=m;

IL請(qǐng)編寫一段matlab程序，完成以下功能：

(1)找出100到200之間的全部質(zhì)數(shù)，將這些質(zhì)數(shù)存放在一個(gè)行

數(shù)組里；

X=100:200;

p=X(isprime(X))

<2)求出這些質(zhì)數(shù)之和；

sum(p)

(3)求出100到200之間的全部非質(zhì)數(shù)之和(包括100和200)。

sum(XCisprime(X)))

12.y=[o.7+-^#lsinx,編寫一段matlab程序，要求如下

L(l+x2)J

(1)在[0,2加區(qū)間，每隔0.01取一x數(shù)值，計(jì)算出相應(yīng)的)，的函

數(shù)值：

x=0:0.01:2*pi;

y=(0.7+2*cos(x)./(1+x.2)).*sin(x)

(2)依據(jù)MATLAB計(jì)算出的數(shù)據(jù).找出在[0,2m內(nèi)該函數(shù)的微小值

的坐標(biāo)。

[my,idx]=min(y)

x(idx)

數(shù)值計(jì)算

1、在MATLAB中，A是一個(gè)10X10數(shù)組，我們把該數(shù)組看成矩陣的話,

則此矩陣的行列式值二det(Q,此矩陣的逆矩陣(假如存

在的話)=inv(A)o(用MATLAB的函數(shù)表示)

2、求解下列方程組：

A=[l13-l;01-11;1122;1-11-1];

b=[2;l;4;0];

x=A\b

A=[34-7-12;5-742;l08-5;-65-210];

b=[4;-3;9;-8];

xyzw=A\b

3、求有理分式R=/的商多項(xiàng)式和余多項(xiàng)式

(X-+2X-2)(5A-3+2X-+1)

n=conv([3010],[3000.5]);

d=conv([l2-2],[5201]);

[q,r]=deconv(n,d)

4,一元多項(xiàng)式p=2/-3/+4x,寫出表示p的MATLAB語(yǔ)句p

40],求p=0的根的MATLAB語(yǔ)句是roots(p),求x=4.3

時(shí)P的數(shù)值的MATLAB語(yǔ)句是_r—(p,4.3)。

程序設(shè)計(jì)

】、M腳本文件和函數(shù)文件的有什么區(qū)分？

最主要的差別是，函數(shù)有自己的獨(dú)立工作區(qū)，和外界的信息交換

通過(guò)輸入輸出參數(shù)實(shí)現(xiàn)；而腳本沒(méi)有獨(dú)立工作區(qū)，其干脆運(yùn)用宿主程

序的工作區(qū)，并把結(jié)果留在宿主程序中。

2、駕馭下面程序設(shè)計(jì)中常用的函數(shù)及關(guān)鍵字:

inputpausewarningnarginnargout

forwhileifelseifelseswitchcase

otherwiseendtry

catchbreakcontinueend

3、分別用for和while循環(huán)語(yǔ)句計(jì)算K=*2,的程序。還請(qǐng)寫出一

1-0

種避開(kāi)循環(huán)的計(jì)算程序。(不要用symsum函數(shù))

K=0;

fori=0:63

K=K+i;

end

K=0;i=0;

whilei<=63

K=K+i;

i=i+l;

end

K=sum⑵'(0:63))

4、MATLAB供應(yīng)了兩種循環(huán)結(jié)構(gòu)，一種是for循環(huán)結(jié)構(gòu)，另一種是

[vhi],循環(huán)結(jié)構(gòu)。

5、編寫M函數(shù)求和s=l+2+3+

functions=mysum(n)

s=sum(l:n);

6、編寫M函數(shù)求積〃=lx2x3xxn

functionp=mysum(n)

p=prod(l:n);

7、編寫M函數(shù)列出從整數(shù)a到整數(shù)b之間能被3整除的奇數(shù)。

functionr=odds3(a,b)

r=a:b;

r=r(mod(r,2)==l&mod(r,3)==0)

8、編寫M腳本列出從100到200不能被3整除同時(shí)也不能被7整除、批注壓土星司：這個(gè)題和下個(gè)題的內(nèi)容「了

的數(shù)。

X=100:200;

X(mod(X,3)"=0&mod(X,7)~=0)

9、一張紙厚0.06mm且足夠大，試問(wèn)將紙對(duì)折多少次，其厚度將超

過(guò)10000m?

按題意，就是要求2%*0.06E-3>=10000,所以

n=cei1(log(10000/0.06e-3)/log(2))

io,編寫MATLAB腳本輸出“水仙花數(shù)”及水仙花數(shù)的個(gè)數(shù)。所謂“水

仙花數(shù)”是一個(gè)3位數(shù)，其各位數(shù)字的立方和等于該數(shù)本身。例

如153=「+53+3、

c=0;

forn=100:999

m=[fix(n/100)fix(mod(n,100)/10)mod(n,10)];

ifn==sum(m.”3),

fprintfC%i=%i"3+%/3+%i"3\n，,n,m)；

c=c+l;

end

end

fprintf('\n共%i個(gè)水仙花數(shù)\n',c)

輸出結(jié)果:

153=13+53+3%

370=33+73+03

371=33+73+3

407二43+03+3

共4個(gè)水仙花數(shù)

繪圖

1、MATLAB中,繪制三維曲面圖的函數(shù)是（)

(A)surf(B)plot(C)subplot(D)plot3

2、MATLAB中,要繪制三維空間曲線，應(yīng)當(dāng)運(yùn)用（)函數(shù)。

(A)polar(B)plot(C)subplot(D)plot3

3、在MATLAB中，A是一個(gè)1000行2列的二維數(shù)值數(shù)組，現(xiàn)在要把

A的第一列數(shù)據(jù)作為橫坐標(biāo)，把A的其次列數(shù)據(jù)作為縱坐標(biāo)，畫

出一條曲線，試寫出相應(yīng)的MATLAB語(yǔ)句

plol(A(:,1),A(:,2))

4、MATLAB繪圖指令中的—匈_指令允許用戶在同一個(gè)圖形窗

里布置幾個(gè)獨(dú)立的子圖。

5、編寫一段matlab程序，繪制出二元函數(shù)z=列B三維網(wǎng)線圖,

xy

要求如下:

(1)x?y的取值范圍為-9W%?9,-9<y<9:

(2)x,y每隔0.5取一個(gè)點(diǎn)；

(3)圖形的線型和顏色由MATLAB自動(dòng)設(shè)定

[x,y]=meshgrid((-9:0.5:9)+eps);%加eps避開(kāi)出現(xiàn)0/0

z=2*sin(x).*sin(y)./(x.*y);

mesh(x,y,z)

6、編寫一段matlab程序，繪制出函數(shù),，1=二皿3?2=§皿21)圖形的

X

MATLAB語(yǔ)句，要求如下：

(1)x的取值范圍為-3Vx?3;

(2)x每隔0.01取一個(gè)點(diǎn)；

(3)yl和y2的圖形要畫在同一幅圖里；

(4)圖形的線型和顏色由MATLAB自動(dòng)設(shè)定。

x=(-3:0.01:3)+eps;%加eps避開(kāi)出現(xiàn)除0警告

yl=x.*sin(l./x);

y2=sin(2*x);

plot(x,yl,x,y2)

符號(hào)計(jì)算

1、在MATLAB中，能正確的把x、y定義成符號(hào)變量的指令是(

)

(A)symxy(B)symx,y

(C)symsx,y(D)symsxy

2、要清除MATLAB工作空間中保存的變量，應(yīng)當(dāng)運(yùn)用leu指令。批注[王士星41：不是想說(shuō)清除符號(hào)運(yùn)算內(nèi)核中的變量

IE?

3、在MATLAB中,指令findsym(sym(*sin(w*t),),1)的執(zhí)行結(jié)果是

wo

4、在MATIAB中，寫出把x和y定義為符號(hào)變量的相應(yīng)語(yǔ)句—批注I王士星5]:上面不是剛考過(guò)嗎？

匕。

5、求解下列方程組的符號(hào)解：

A=[l13-l;01-11;1122;1-11-1];

b=[2;l;4;0];

x=sym(A)\b

A=[34-7-12;5-742;108-5;-65-210];

b=[4;-3;9;-8];

xyzw=sym(A)\b

6、求極限：

(l)求函數(shù)y=(l+x)，在x=0處的極限；

?symsx

?limit((1+x)"(l/x),0)

ans=

exp(l)

(2)求函數(shù)y=吧在x=o處的極限；

Ig5x

?symsx

?limit(sin(3*x)/tan(5*x),0)

ans=

3/5

(3)求函數(shù))，=與在趨向正無(wú)窮處的極限；

3

?symsxn

?limit(n*x/3x,inf)

ans=

0

(4)求函數(shù))，=嘩在趨向正無(wú)窮處的極限；

X'

?symsx

?limit(log(x)2/x3,inf)

ans=

0

7、求導(dǎo)數(shù)：

(l)求函數(shù)y=[-3x+3的50階導(dǎo)數(shù)；

廠

?symsx

?y=l/x"2-3*x+3;

?diff(y,50)

ans=

155111875328738228022424301646930321106325972023698

6112000000000000/x^52

<2)求函數(shù)y=asin侔'+6在f=匕處的3階導(dǎo)數(shù)；

symstabc

y=a*sin(b*exp(c"t)+t-a);

simple(subs(diff(y,t,3),t,b))

8、求不定積分：

?symsx

?int(l/sin(x)3)

ans=

-l/2/sin(x)-2*cos(x)+l/2*log(csc(x)-cot(x))

>>symsxa

?int(1/(a-2-x"2))

ans二

-l/2/a*log(a-x)+l/2/a*log(a+x)

>>symsx

>>int((sqrt(x-2-3)-sqr:(x~2+3))/sqrt(x"4-9))

ans=

(x^4-9)*(l/2)/(x*2-3)*(l/2)/(x*2+3)*(l/2)*asinh(l/3

*3F/2)*x)-l/(x-2+3廠(l/2)*(x-4-9)c(l/2)/(x-2-3廠(1/2

)*log(x+(x"2-3)\l/2))

9、求定積分及廣義積分

?symsxa

?int(sqrt(x*2+a),-2,2)

ans=

2*(4+a)~(l/2)+l/2*a*log(2+(4+a)~(l/2))-l/2*a*log(-2

+(4+a廠(1/2))

?symsx

?int(sin(x)-2*cos(x)2,-pi,pi)

a