Matlab-期末考試題庫(共12套卷)

第一套華東交通大學(xué)2015-2016學(xué)年第一學(xué)期Matlab期末考試一、填空題（4*6=24分）1、在MATLAB命令窗口中的“>>”標(biāo)志為MATLAB的命令行提示符，“│”標(biāo)志為輸入提示符。2、符號表達式sin(2*a+t)+m

中獨立的符號變量為_t_。3、在通常情況下，左除x=a\b是

的解，右除x=b/a是

的解，一般情況下，。4、為了使兩個plot的圖形在同一個坐標(biāo)顯示，可以使用___holdon__命令進行圖形保持；可以使用_gridon__命令為圖形添加網(wǎng)格。5、倘若要是對x進行賦值，從5到25，間隔是0.1,為x=5:0.1:25

；倘若要是對x進行賦值，從2到20，中間間隔100個點,為x=linspace(2,20,100)

6、A=[1,2;3,4];B=[1,0;0,1];

A*B=[1

2；3

4]_,

A.*B=_[1

0；0

4]_；3.a*x=b,x*a=b,a\b,b/a;二、選擇題（4*6=24分）7、如果x=1:

2

:

8,則x(1)和x(4)分別是(B)

A．1，8

B．1,

7

C．2,

8

D．2,

78、運行如下程序后，輸入９回車，命令窗口（command

windows）顯示的結(jié)果為(A)

c=input('請輸入一個字符','s');

if

c>='A'

&

c<='Z'

disp(setstr(abs(c)+abs('a')-abs('A')));

elseif

c>='a'&

c<='z'

disp(setstr(abs(c)-

abs('a')+abs('A')));

elseif

c>='0'&

c<='9'

disp(abs(c)-abs('0'));

else

disp(c);

end

A.

9

B.

8

C.

7

D.

69、MATLAB表達式2*2^3^2的結(jié)果是(A)

A．128

B．4096

C.

262144

D．25610、在循環(huán)結(jié)構(gòu)中跳出循環(huán)，執(zhí)行循環(huán)后面代碼的命令為(B)

(A)

return

(B)

break

(C)

continue

(D)

keyboard

11、在圖形指定位置加標(biāo)注命令是（C）

A.

title(x,y,?y=sin(x)?);

B.

xlabel(x,y,?y=sin(x)?);

C.

text(x,y,?y=sin(x)?);

D.

legend(x,y,?y=sin(x)?);12、下列哪個變量的定義是不合法的（A）。

(A)

abcd-3

(B)

xyz_3

(C)

abcdef

(D)

x3yz三、程序題（52分）13（10分）、在同一窗口繪制6個子圖，已知t=0.1:0.1:2*pi，每個子圖分別繪制（其中y1=sin(t),y2=cos(2t),y3=sin(t)+cos(t),y4=3t,y5=2ln(t),y6=et）；答案13.t=0.1:0.1:2*pi；

subplot(2,3,1)

y1=sin(t)；

plot(t,y1)

title('y1=sin(t)')

subplot(2,3,2)

y2=cos(2*t)；

plot(t,y2)

title('y2=cos(2t)')

subplot(2,3,3)

y3=sin(t)+cos(t)；

plot(t,y3)

title('y3=sin(t)+cos(t)')

subplot(2,3,4)

y4=3*t；

plot(t,y4)

title('y4=3t')

subplot(2,3,5)

y5=2*log(t)；

plot(t,y5)

title('y5=2ln(t)')

subplot(2,3,6)

y6=exp(t)；

plot(t,y6)

title('y6=e^t')14（10分）、某商場對顧客所購買的商品實行打折銷售，標(biāo)準(zhǔn)如下(商品價格用price來表示)：price<200沒有折扣200≤price<5003%折扣500≤price<10005%折扣1000≤price<25008%折扣2500≤price<500010%折扣5000≤price14%折扣輸入所售商品的價格，求其實際銷售價格。答案14.price=input('請輸入商品價格');switchfix(price/100)case{0,1}%價格小于200rate=0;case{2,3,4}%價格大于等于200但小于500rate=3/100;casenum2cell(5:9)%價格大于等于500但小于1000rate=5/100;casenum2cell(10:24)%價格大于等于1000但小于2500rate=8/100;casenum2cell(25:49)%價格大于等于2500但小于5000rate=10/100;otherwise%價格大于等于5000rate=14/100;endprice=price*(1-rate)%輸出商品實際銷售價格15（10分）、在同一圖形窗口中按照x=0.1的步長間隔分別繪制曲線y1=,0x；y1=sin(x)+cos(x),1≤x≤2。要求y1曲線為紅色虛線，數(shù)據(jù)點用圓圈標(biāo)識；y2曲線為藍色點劃線；給出圖例；標(biāo)記圖名為y1,y2。答案15.x1=0:0.1:1;

y1=x1.^2.*exp(-x1);x2=1:0.1:2;

y2=sin(x2)+cos(x2);plot(x1,y1,'r--o',x2,y2,'b-.')legend('y1','y2')title('y1,y2')

16（10分）、寫出執(zhí)行以下代碼后C，D，E的值

A=[1,2,3;4:6;7:9];

C=[A;[10,11,12]],

D=C(1:3,[2

3])

E=C(2,[1

2])

16.C

=

1

2

3

4

5

67

8

9

10

11

12

D

=

2

3

5

6

8

9

E

=

4

517（12分）、已知系統(tǒng)函數(shù)為試畫出零極點分布圖，并求系統(tǒng)的沖激響應(yīng)h(t)和頻率響應(yīng)H(jw)，并判斷系統(tǒng)是否穩(wěn)定。17.已知系統(tǒng)函數(shù)H（s），求系統(tǒng)的沖激響應(yīng)h(t)和頻率響應(yīng)H(jw)，可調(diào)用impulse函數(shù)和freqs函數(shù)。num=[1];den=[1221];sys=tf(num,den);poles=roots(den)figure(1);pzmap(sys);t=0:0.02:10;h=impulse(num,den,t);figure(2);plot(t,h)title('ImpulseRespone')[H,w]=freqs(num,den);figure(3);plot(w,abs(H))xlabel('\omega')title('MagnitudeRespone')第二套Matlab模擬試卷填空題（每空2分，總分20）。1.MATLAB常用操作界面包括命令窗口、工作空間窗口（瀏覽器）、命令歷史窗口、當(dāng)前目錄窗口、內(nèi)存數(shù)組編輯器、M文件編輯/調(diào)試器、幫助導(dǎo)航/瀏覽器、圖形窗口等。2.MATLAB中逗號主要用作要顯示計算結(jié)果的指令與其后指令的分隔；用作輸入量與輸入量之間的分隔符；用作數(shù)組元素分隔符號。3.工作空間瀏覽器主要用于內(nèi)存變量的查閱、保存和編輯。4.MATLAB實現(xiàn)將全下標(biāo)轉(zhuǎn)換為單下標(biāo)的指令為Sub2ind、據(jù)單下標(biāo)換算出全下標(biāo)的指令為Ind2sub。5.二維數(shù)組的標(biāo)識有“全下標(biāo)”標(biāo)識、“單下標(biāo)”標(biāo)識、“邏輯1”標(biāo)識。二、選擇題（每空2分，總共20分）1.下列哪個變量的定義是不合法的（A）(A)abcd-3(B)xyz_3(C)abcdef(D)x3yz2.下列哪條指令是求矩陣的行列式的值（C）(A)inv(B)diag(C)det(D)eig3.在循環(huán)結(jié)構(gòu)中跳出循環(huán)，執(zhí)行循環(huán)后面代碼的命令為（B）(A)return(B)break(C)continue(D)keyboard4.清空Matlab工作空間內(nèi)所有變量的指令是（C）(A)clc(B)cls(C)clear(D)clf5.用round函數(shù)四舍五入對數(shù)組[2.486.393.938.52]取整，結(jié)果為（C）(A)[2638](B)[2648](C)[2649](D)[3749]6.已知a=2:2:8,b=2:5，下面的運算表達式中，出錯的為（C）(A)a'*b(B)a.*b(C)a*b(D)a-b7.角度，計算其正弦函數(shù)的運算為（D）(A)SIN（deg2rad(x)）(B)SIN(x)(C)sin(x)(D)sin(deg2rad(x))8.下面的程序執(zhí)行后array的值為(A)fork=1:10ifk>6break;elsearray(k)=k;endend(A)array=[1,2,3,4,5,6](B)array=[1,2,3,4,5,6,7,8,9,10](C)array=6(D)array=10.9．i=2;a=2i;b=2*i;c=2*sqrt(-1);程序執(zhí)行后；a,b,c的值分別是多少？(A)a=4,b=4,c=2.0000i(C)(B)a=4,b=2.0000i,c=2.0000i(C)a=2.0000i,b=4,c=2.0000i(D)a=2.0000i,b=2.0000i,c=2.0000i10.求解方程x4-4x3+12x-9=0的所有解（A）(A)1.0000,3.0000,1.7321,-1.7321(B)1.0000,3.0000,1.7321i,-1.7321i(C)1.0000i,3.0000i,1.7321,-1.7321(D)-3.0000i,3.0000i,1.7321,-1.7321三、閱讀程序并回答問題（每題5分，共20分）1、寫出下列指令運行結(jié)果。A=zeros(2,4);A(:)=1:8;s=[235];A(s)Sa=[102030]'A(s)=Saans=235Sa=102030A=12030710468寫出下列指令運行結(jié)果。A=reshape(1:16,2,8)reshape(A,4,4)s=;A(s)=0答案：A=reshape(1:16,2,8)A=13579111315246810121416reshape(A,4,4)ans=15913261014371115481216s=;A(s)=0A=005700131524001012003、寫出下列指令運行結(jié)果。A=[1,2;3,4];B=[-1,-2;2;1];S=3;A.*BA*BS.*AS*B答案：ans=-1-464ans=305-2ans=36912ans=-3-6634、下面的函數(shù)主要完成什么功能？functionf=factor(n)ifn<=1f=1;elsef=factor(n-1)*n;End答：利用函數(shù)的遞歸調(diào)用，求n！。四、程序填空(每空3分，共15分)編寫MATLAB程序，在同一圖形窗口中繪制y1和y2兩條曲線，并對圖形坐標(biāo)軸進行控制。程序如下：x=-3:0.1:3;y1=2*x+5;y2=x.^2-3*x;plot(x,y1);%繪制曲線y1①;%保持曲線y1①holdonplot(x,y2);%繪制曲線y2m1=max([y1,y2]);m2=min([y1,y2]);②;%用axis調(diào)制坐標(biāo)軸，橫坐標(biāo)在[-33]之間，且完全顯示圖形③;%加入圖形注釋，橫坐標(biāo)為“時間”，縱坐標(biāo)為“幅度”編寫MATLAB程序，對于用戶輸入的向量或矩陣，計算其平均值和方差：④;%提示輸入數(shù)據(jù)并保存到x中，命令行顯示“請輸入數(shù)據(jù)”y=reshape(x,1,prod(size(x)));m=mean(y);%計算均值m⑤;%計算方差v①holdon②axis([-33m2m1])③xlabel(‘時間’);ylabel(‘幅度’);④x=input(‘請輸入數(shù)據(jù)\n’)或x=input(‘請輸入數(shù)據(jù)’)⑤v=var(y)或1/(length(y)-1)*(sum((y-m).^2))五、編程題（40分）分別在同一圖形窗的不同子圖繪制y=sin(t)sin(9t)和y=sin(t)sin(9t)及其包絡(luò)線。（5分）t=(0:pi/100:pi)';y1=sin(t)*[1,-1];y2=sin(t).*sin(9*t);t3=pi*(0:9)/9;y3=sin(t3).*sin(9*t3);subplot(1,2,1)plot(t,y1,'r:',t,y2,'b',t3,y3,'bo')subplot(1,2,2)plot(t,y2,'b')axis([0,pi,-1,1])編寫程序?qū)崿F(xiàn)f(n)=f(n-1)+f(n-2)（f(1)=1和f(2)=2）函數(shù)。（5分）functionf=fab(n)if(n==1)f=1;elseif(n==2)f=2;elsef=fab(n-1)+fab(n-2);End3、求解以下線性方程組，要求寫出程序代碼和運行結(jié)果（5分）2x1-3x2+x3+2x4=8x1+3x2+x4=6x1-x2+x3+8x4=17x1+x2-2x3+2x4=5>>symsabcx>>symsx1x2x3x4>>clear>>symsx1x2x3x4>>f=2*x1+3*x2+x3+2*x4-8;>>g=x1+3*x2+x4-6;>>h=x1-x2+x3+8*x4-1;>>i=7*x1+x2-2*x3+2*x4-5;>>[x1,x2,x3,x4]=solve(f,g,h,i)x1=13/17x2=320/187x3=212/187x4=19/1874、打印出所有的水仙花數(shù)。所謂“水仙花數(shù)”，是指一個三位數(shù)，其各位數(shù)字立方之和等于該數(shù)本身。（10分）fork=100:999 a=fix(k/100); b=rem(fix(k/10),10); c=rem(k,10); ifa.^3+b.^3+c.^3==k fprintf(“%u,\t\t”,k)

; endend第三套MATLAB試卷一、填空題（每題2分，共10題，共20分，錯一個得1分，錯2個得0分）MATLAB命令窗口中可用命令清除工作區(qū)中的變量；用命令清除命令窗口中的內(nèi)容。MATLAB中的運算包括、、。MATLAB中的M文件有、。MATLAB中的程序控制結(jié)構(gòu)包括、、。已知矩陣A[123;456;789],A(6)A（1,3）。MATLAB中用函數(shù)完成矩陣的求逆運算，用函數(shù)求矩陣的行列式。在MATLAB中，如果需要繪制出具有不同縱坐標(biāo)標(biāo)度的兩個圖形，可以使用函數(shù)。MATLAB提供了函數(shù)和函數(shù)來繪制三維曲面圖。用if語句判斷60≤x≤80，在MATLAB中if語言后面的判斷應(yīng)寫為。為了使兩個plot的圖形在同一坐標(biāo)顯示，可以使用命令進行圖形保持，可以使用命令為圖形添加網(wǎng)格。選擇題（每題2分，共10題，共20分）1.下列哪個變量的定義是不合法的(A)abcd-3(B)xyz_3(C)abcdef(D)x3yz2.在循環(huán)結(jié)構(gòu)中跳出循環(huán)，執(zhí)行循環(huán)后面代碼的命令為(A)return(B)break(C)continue(D)keyboard3、下列變量名中是合法的。(A)char_1;(B)x*y;(C)x\y;(D)end4.如果x=1:-1:-8,則x(1)和x(5)分別是

A．-1，-3B．-1,-7C．-1,-6D．1,-35.在MATLAB中下列數(shù)值的表示不正確的是

A．+25B．1.3e-5C．4-3*e^2D．3-2*pi6.用round函數(shù)四舍五入對數(shù)組[2.486.393.938.52]取整，結(jié)果為(A)[2638](B)[2648](C)[2649](D)[3749]

7.下列哪條指令是求矩陣的行列式的值(A)inv(B)diag(C)det(D)eig8.在圖形指定位置加標(biāo)注命令是A.title(x,y,’y=sin(x)’);B.xlabel(x,y,’y=sin(x)’);C.text(x,y,’y=sin(x)’);D.legend(x,y,’y=sin(x)’);9．i=2;a=2i;b=2*i;c=2*sqrt(-1);程序執(zhí)行后；a,b,c的值分別(A)a=4,b=4,c=2.0000i(B)a=4,b=2.0000i,c=2.0000i(C)a=2.0000i,b=4,c=2.0000i(D)a=2.0000i,b=2.0000i,c=2.0000i10、在循環(huán)結(jié)構(gòu)中跳出循環(huán)，但繼續(xù)下次循環(huán)的命令為(A)return;(B)break;(C)continue;(D)keyboard三、簡答題（每題7.5分，共4題，共30分）1.A=[12;34];B=[12;14];寫出下面幾條指令的運行結(jié)果：C=A*BD=A+BE=A.*B在MATLAB窗口中，用哪兩種方法啟動Simulink？3.判斷下列一段程序的錯誤，并進行改正。x=-2:0.2:2;y=x.^2;plot(x,y)str1=num2str(min(x));%數(shù)字轉(zhuǎn)化為字符串str2=num2str(max(x));%數(shù)字轉(zhuǎn)化為字符串out=[Valueofffrom'str1'to'str2'];xlabel(out);一個3位整數(shù)各位數(shù)字的立方和等于該數(shù)本身則稱該數(shù)為水仙花數(shù)。輸出100-999之間的全部水仙花數(shù)。請補全程序中所空缺的部分。Form=100:999m1=fix(m/100);m2=rem(fix(m/10),10);m3=rem(m,10);Ifm==()disp(m)endend編程題（每小題15分，共2題，共30分）1.編寫MATALAB程序，完成下列任務(wù)（將程序保存為test04.m文件）：（1）在區(qū)間[0,4*pi]上均勻地取20個點構(gòu)成向量；（2）分別計算函數(shù)y1=sin(t)與y2=2cos(2t)在向量t處的函數(shù)值；（3）在同一圖形窗口繪制曲線y1=sin(t)與y2=2cos(2t)，要求y1曲線為黑色點畫線，y2曲線為紅色虛線圓圈；并在圖中恰當(dāng)位置標(biāo)注兩條曲線的圖例；給圖形加上標(biāo)題“y1andy2”。2、Fibonacci數(shù)列的各元素為：1、1、2、3、5、8、……，滿足以下關(guān)系F1=1，F(xiàn)2=1，F(xiàn)n=Fn-1+Fn-2用while和for兩種方法實現(xiàn)。答案填空題1clearclc2算術(shù)運算關(guān)系運算邏輯運算3命令文件函數(shù)文件4順序結(jié)構(gòu)循環(huán)結(jié)構(gòu)選擇結(jié)構(gòu)5.83InvdetplotyymeshsurfIfx<=80&>=60holdongridon選擇題1-5ABADB6-10CCCCC簡答題C=310D=24E在MATLAB命令窗口中輸入“simulink”再回車，或者單擊工具欄中的SIMULINK圖標(biāo)3.out=['Valueofffrom'str1'to'str2];或out=['Valueofffromstr1tostr2’];4.m1*m1*m1+m2*m2*m2+m3*m3*m3編程題1.程序如下：（1）t=linspace(0,4*pi,20);（2）y1=sin(t);y2=2*cos(2*t);（3）plot(t,y1,'k-.');text(1.2,sin(1.2),'y1\leftarrow','FontSize',12);holdon;plot(t,y2,'r--o');text(6,1.5,'y2\leftarrow','FontSize',12);title('y1andy2')2.按M函數(shù)文件格式創(chuàng)建文件開頭functionf=chengxu（n)用while循環(huán)實現(xiàn)程序功能f(1)=1;f(2)=1;i=2;whilei<=nf(i+1)=f(i-1)+f(i);i=i+1;end用for循環(huán)實現(xiàn)f(1)=1;f(2)=1;fori=2:nf(i+1)=f(i-1)+f(i);end第四套MATLAB期末考試一、填空題（每題2分，共10）1：MATLAB的通用命令中若對函數(shù)功能和用法不熟悉，可用help命令；顯示變量、陣列或文本時，可用disp命令；清除工作空間中的變量時，可用clear命令；保存和恢復(fù)工作空間，或者工作空間中的特定變量時，可用save和load命令。2：在MATLAB圖形系統(tǒng)中，如果需要給已經(jīng)繪制好的圖形加上柵格線，可以使用gridon命令，如果需要在指定的坐標(biāo)位置標(biāo)注文本時使用的函數(shù)為text。3：MATLAB中用于控制不定次重復(fù)循環(huán)語句為while，若在循環(huán)執(zhí)行過程中需要終止該循環(huán)時采用的語句為break。4：用if判斷語句判斷80≥x>60,在MATLAB中if語言后的判斷應(yīng)寫為ifx<=80&x>60。5：為了使兩個plot的圖形在同一個坐標(biāo)顯示，可以使用holdon命令進行圖形保持；選擇題（每題2分，共20）1.下列哪個變量的定義是不合法的C(A)abcdef(B)xyz_3(C)abcd-3(D)x3yz2:已知a=2:2:8,b=2:5，下面的運算表達式中，出錯的為C(A)a'*b(B)a.*b(C)a*b(D)a-b3:下面的程序執(zhí)行后array的值為(D)fork=1:10ifk>6break;elsearray(k)=k;endend(A)array=10.(B)array=[1,2,3,4,5,6,7,8,9,10](C)array=6(D)array=[1,2,3,4,5,6]4：設(shè)A=[243;531;367]，則sum(A)，length(A)和size(A)的結(jié)果(A)A.[101311]3[33]B.[9916]3[33]C.[9916]9[33]D.[101311]9[33]5:下列關(guān)于腳本文件和函數(shù)文件的描述中不正確的是(B)A.去掉函數(shù)文件第一行的定義行可轉(zhuǎn)變成腳本文件;B.函數(shù)文件可以在命令窗口直接運行;C.腳本文件可以調(diào)用函數(shù)文件;D.函數(shù)文件中的第一行必須以function開始;6:．i=2;a=2i;b=2*i;c=2*sqrt(-1);程序執(zhí)行后；a,b,c的值分別是多少？(C)a=4,b=4,c=2.0000i(B)a=4,b=2.0000i,c=2.0000i(C)a=2.0000i,b=4,c=2.0000i(D)a=2.0000i,b=2.0000i,c=2.0000i7:下列哪條語句是正確的(

C

)

A=[1

2;3

4];B=[3;7];

A.*B;(B)

A=cell(2,3);

A(1,2)=[4;5]

(C)

A=[2

3

5

7;9

4

6

1;7

3

2

5];B=[1

7;0

5];A(2:end,2:2:end)

(D)

x=-5:5;y=-5:5;z=x.*x-y.*y;surf(x,y,z)

8:三、程序執(zhí)行題（共30分）1：寫出執(zhí)行以下代碼后，MATLAB命令窗口上顯示的x矩陣的值(6分)x=[0,1,0,2,0,3,0,4];fork=1:8ifx(k)==0x(k)=k;elsex(k)=2*k+1;endenddisp(x);答案為解下列線性方程組，請寫出MALAB語言程序。（8分）4x1+2x3+3x4=93x2+2x3+3x4=82x1+x4=10x3-2x2+3x4=7答案：A=[4023;0323;2001;0-213;];%3分b=[98107]';%2分x=A\b;或x=inv(A)*b;%3分3：下面的函數(shù)主要完成什么功能？（8分）functionf=factor(n)ifn<=1f=1;elsef=factor(n-1)*n;end答案：利用函數(shù)的遞歸調(diào)用，求n！。4：下面的程序完成功能是什么？（8分）t=0:pi/50:3*pi;y0=exp(-t/2);y=exp(-t/3).*sin(2*t);plot(t,y,'-r',t,y0,':b',t,-y0,':b')xlabel(‘\bf\itt’);ylabel(‘\bf\ity’);答案：在同一個圖中繪制兩個函數(shù)，這兩個函數(shù)分別是：y=e(-t/3)y0=e(-t/3)*sin(3t)四、編程題（每題10分，共20）1：求解以下線性方程組，要求寫出程序代碼和運行結(jié)果（10分）答案：>>symsabcx>>symsx1x2x3x4>>clear>>symsx1x2x3x4>>f=2*x1+3*x2+x3+2*x4-8;>>g=x1+3*x2+x4-6;>>h=x1-x2+x3+8*x4-1;>>i=7*x1+x2-2*x3+2*x4-5;>>[x1,x2,x3,x4]=solve(f,g,h,i)2：設(shè)x=sint,

y=sin(nt+a)，若a=1，令n

=1,2,3,4,在四個子圖中分別畫出其曲線。若n=2,取a=0，π/3，π/2，及π，在四個子圖中分別畫出其曲線。matlab源程序如下;t=-pi:0.1:pi;

%t的取值范圍a=1;n=1;%a=1，n=1

x=sin(t);%x的表達y=sin(n*t+a);

%y的表達式subplot(221),plot(x,y)%在子圖第一張，畫出圖像hold

on保持figure不關(guān)閉t=-pi:0.1:pi;

%t的取值范圍a=1;n=2;

%a=1，n=2

x=sin(t);%x的表達式y(tǒng)=sin(n*t+a);

%y的表達式

subplot(222),plot(x,y)

%在子圖的第二張畫出圖像

hold

on

%保持figure不關(guān)閉

t=-pi:0.1:pi;

%t的取值范圍

a=1;n=3;

%a=1，n=3

x=sin(t);

%x的表達式

y=sin(n*t+a);

%y的表達式

subplot(223),plot(x,y)；%在子圖的第三張，畫出圖像

hold

on

%保持figure不關(guān)閉

t=-pi:0.1:pi;

%t的取值范圍

a=1;n=4;

%a=1，n=4

x=sin(t);

%x的表達式y(tǒng)=sin(n*t+a);

%y的表達式subplot(224),plot(x,y)

%在子圖的第四張上五、程序填空和改錯(每空5分，共20)1：拋物線的方程為,輸入c后，繪制該曲線，尋找曲線的最小值和對應(yīng)的x值。小明編寫程序如下：clc;clear;closeall;c=input('c=');a=1;x=-2:0.1:2;xnum=length(x)orxnum=size(x,2);%(1)對變量xnum賦為數(shù)組x的個數(shù)y=zeros(1,xnum);y=a*x^2+c*x;→y=a*x.^2+c*x;%(2)按照函數(shù)形式對y賦值figure;plot(x,y) %(3)繪制曲線xlabel('x');ylabel('y');[ymin,Imin]=min(y);xmin=x(ymin);→xmin=x(Imin);%(4)找到最小函數(shù)值對應(yīng)的x值第五套matlab期末考試試卷考試時間：120分鐘考試分數(shù)：100分一、填空題（每空1分，20分）1、MATLAB常用操作界面包括、工作空間窗口、、、內(nèi)存數(shù)組編輯器、M文件編輯/調(diào)試器、幫助導(dǎo)航/瀏覽器、圖形窗口等。2、MATLAB中Inf或inf表示、NaN或nan表示、nargout表示。3、MATLAB中逗號主要用作；用作輸入量與輸入量之間的分隔符；。4、工作空間瀏覽器主要用于內(nèi)存變量的、和。5、MATLAB實現(xiàn)將全下標(biāo)轉(zhuǎn)換為單下標(biāo)的指令為、據(jù)單下標(biāo)換算出全下標(biāo)的指令為。6、二維數(shù)組的標(biāo)識有、、“邏輯1”標(biāo)識。7、在一個元胞數(shù)組A中尋訪第2行第3列元胞元素用；尋訪數(shù)組第2行第3列元胞中的內(nèi)容用。8、MATLAB中clf用于清除、clc用于清除、clear用于清除MATLAB。選擇題（每題2分，20分）1.在matlab中_______用于括住字符串。A.，B.；C.‘’D.“”2.下列___D____是合法變量。A.1_1B.a/bC.變量1D.Eps3.產(chǎn)生所有元素均為1的2行3列矩陣的命令是______。A.ones(2,3)B.ones(3,2)C.eye(2,3)D.eye(3,2)4.下列哪個變量的定義是不合法的A.abcd-3B.xyz_3C.abcdefD.x3yz5.已知數(shù)組a=,若要產(chǎn)生a的一個子矩陣，則下列命令中錯誤的是_____________。a([12],[34])B.a(1:2,3:end)C.a(linspace(1,2,2),3:4)D.a(:,:)6.下面的程序執(zhí)行后array的值為_____________。fork=1:10ifk>6break;elsearray(k)=k;endendA.array=[1,2,3,4,5,6]B).array=[1,2,3,4,5,6,7,8,9,10]C.array=6D.array=10.7.在循環(huán)結(jié)構(gòu)中跳出循環(huán)，但繼續(xù)下次循環(huán)的命令為_________。A.return;B.break;C.continue;D.keyboad8.圖形窗口（Figure）顯示網(wǎng)格的命令是_____________。A.axisonB.gridonC.boxonD.holdon9.已知a=2:2:8,b=2:5，下面的運算表達式中，出錯的為____。A.a'*bB.a.*bC.a*bD.a-b10.運行如下程序后,命令窗口（commandwindows）顯示的結(jié)果為__。A=[13,-56,78;25,63,-735;78,25,563;1,0,-1];y=max(max(A))A.y=564B.y=9C.y=-735D.y=563三、閱讀程序并回答問題（每題4分，共28分）1、寫出下列指令運行結(jié)果。A=[1,2,3;4:6;7:9];C=[A;[10,11,12]],D=C(1:3,[23])E=C(2,[12])2、寫出下列指令運行結(jié)果。A=reshape(1:16,2,8)reshape(A,4,4)s=;A(s)=03、寫出下列指令運行結(jié)果。A=[1,2;3,4];B=[-1,-2;2;1];S=3;A.*BA*BS.*AS*B4、下面的函數(shù)主要完成什么功能？functionf=factor(n)ifn<=1f=1;elsef=factor(n-1)*n;end5、寫出下列指令運行結(jié)果。ch=‘ABc123d4e56Fg9’;subch=ch(1:5)revch=ch(end:-1:1)k=find(ch>=‘a(chǎn)’&ch<=‘z’);ch(k)=ch(k)-(‘a(chǎn)’-‘A’);char(ch)6、寫出下列指令運行結(jié)果。A(1,1)={'thisiscell'};A{1,2}={[123;456]};A{2,1}=[1+2*i];A{2,2}=A{1,2}{1}+(A{1,2}{1}(1,1)+A{1,2}{1}(2,2));celldisp(A)7、下面的程序完成功能是什么？t=0:pi/50:4*pi;y0=exp(-t/3);y=exp(-t/3).*sin(3*t);plot(t,y,'-r',t,y0,':b',t,-y0,':b')xlabel(‘\bf\itt’);ylabel(‘\bf\ity’);gridon;四、編程題（32分）1、在同一圖上分別用紅色實線和綠色虛線繪制y1=sin(x)和y2=cos(x)在區(qū)間[0，4*pi]的曲線，并用星號*標(biāo)出兩條曲線的交點以及建立圖例。（6分）2、求解以下線性方程組，要求寫出程序代碼。（6分）2x1-3x2+x3+2x4=8x1+3x2+x4=6x1-x2+x3+8x4=17x1+x2-2x3+2x4=53、某商場對顧客所購買的商品實行打折銷售，標(biāo)準(zhǔn)如下(商品價格用price來表示)：price<200沒有折扣200≤price<5003%折扣500≤price<10005%折扣1000≤price<25008%折扣2500≤price<500010%折扣5000≤price14%折扣輸入所售商品的價格，求其實際銷售價格。（8分）4、編寫程序?qū)崿F(xiàn)f(n)=f(n-1)+f(n-2)（f(1)=1和f(2)=2）函數(shù)。（6分）5、編寫M函數(shù)文件,利用for循環(huán)或while循環(huán)完成計算函數(shù)y=sum(n)=QUOTE的任務(wù)，并利用該函數(shù)計算n=20時的和（將總程序保存為？？.m文件）。（6分）參考答案填空題1、命令窗口命令歷史窗口當(dāng)前目錄窗口。2、無窮大不是一個數(shù)函數(shù)輸出宗量數(shù)目。3、顯示計算結(jié)果的指令與其后指令的分隔用作數(shù)組元素分隔符號。4、查閱保存編輯。5、Sub2indInd2sub。6、“全下標(biāo)”標(biāo)識“單下標(biāo)”標(biāo)識。7、A(2,3)A{2,3}。8、清除圖形窗清除指令窗中顯示內(nèi)容清除MATLAB工作空間中保存的變量。二、選擇題（每題2分，共20分）1、C2、D3、A4、A5、D6、A7、C8、B9、C10、D三、閱讀程序并回答問題（每題4分，共28分）1、C=123456789101112D=235689E=452、A=13579111315246810121416ans=15913261014371115481216A=005700131524001012003、ans=-1-464ans=305-2ans=36912ans=-3-6634、利用函數(shù)的遞歸調(diào)用，求n！。5、subch=ABc12revch=9gF65e4d321cBAans=ABC123D4E56FG96、A{1,1}=thisiscellA{2,1}=1.0000+2.0000iA{1,2}{1}=123456A{2,2}=7891011127、在同一個圖中繪制兩個函數(shù)，這兩個函數(shù)分別是：y=e(-t/3)y0=e(-t/3)*sin(3t)其中y是用紅色的細實線，而y0是用蘭色的虛線繪制，t的區(qū)間是(0:4)，t的步長為/50，t軸和y軸分別用斜粗題指示，圖中有網(wǎng)格表示。四、編程題1、答案clfx=0:pi/200:2*pi;y1=sin(x);y2=cos(x);zz=x(find(abs(y1-y2)<0.005))z=min(zz)plot(x,y1,'r-',x,y2,'g-.')holdonplot(zz,sin(zz),'*')legend('sin','cos')2、答案t=(0:pi/100:pi)';y1=sin(t)*[1,-1];y2=sin(t).*sin(9*t);t3=pi*(0:9)/9;y3=sin(t3).*sin(9*t3);subplot(1,2,1)plot(t,y1,'r:',t,y2,'b',t3,y3,'bo')subplot(1,2,2)plot(t,y2,'b')axis([0,pi,-1,1])3、答案price=input('請輸入商品價格');switchfix(price/100)case{0,1}%價格小于200rate=0;case{2,3,4}%價格大于等于200但小于500rate=3/100;casenum2cell(5:9)%價格大于等于500但小于1000rate=5/100;casenum2cell(10:24)%價格大于等于1000但小于2500rate=8/100;casenum2cell(25:49)%價格大于等于2500但小于5000rate=10/100;otherwise%價格大于等于5000rate=14/100;endprice=price*(1-rate)%輸出商品實際銷售價格答案functionf=fab(n)if(n==1)f=1;elseif(n==2)f=2;elsef=fab(n-1)+fab(n-2);end5、functionsum=text4(n);sum=0;m=1;whilem<=n;sum=sum+m^mm=m+1;end當(dāng)n=20時用上述程序求得的和為：1.0688e+026第六套一、選擇題（2×10）1.變量名是以____開頭，后接字母、數(shù)字或下劃線的字符序列，最多63個字符。變量名區(qū)分大小寫。AA.字母B.數(shù)字C.下劃線D.大寫字母2.冒號表達式的一般格式為：BA.初始值:終止值:步長B.初始值:步長:終止值C.初始值:終止值:元素總數(shù)D.初始值:元素總數(shù):終止值3.linspace函數(shù)調(diào)用的格式為：CA.linspace(初始值:終止值:步長)B.linspace(初始值:步長:終止值)C.linspace(初始值:終止值:元素總數(shù))D.linespace(初始值:元素總數(shù):終止值)4.A\B等效于：AA.inv(A)*BB.inv(B)*AC.A*inv(B)D.B*inv(A)5.case子句后面的表達式不可以為：A.標(biāo)量B.字符串C.單元矩陣D.數(shù)值矩陣6.以下關(guān)于函數(shù)的說法錯誤的是：BA.在M文件中只有一個主函數(shù)B.主函數(shù)不一定在最上方C.主函數(shù)不能被其他文件的函數(shù)調(diào)用D.help和lookfor找不到子函數(shù)7.函數(shù)subplot（m,n,p），m,n,p中哪個表示區(qū)號，區(qū)號按行還是列優(yōu)先編號？DA.p列B.n列C.n行D.p行8.下列哪個函數(shù)用來繪制三維網(wǎng)格圖：BA.surf函數(shù)B.mesh函數(shù)C.plot3函數(shù)D.以上都不是9.gcf函數(shù)表示：AA.獲取當(dāng)前圖形窗口的句柄B.獲取當(dāng)前坐標(biāo)軸的句柄C.獲取最近被選中的圖形對象的句柄D.按照指定的屬性來獲取圖形對象的句柄10.N次多項式用長度為___的___向量表示。CA.N+1列B.N列C.N+1行D.N行二、填空題（2×10）1.注釋以_%開頭，后面是注釋的內(nèi)容.2.__clear__命令用于刪除MATLAB工作空間中的變量，_clc___命令用于清屏。3.矩陣元素按列（列/行）編號。4.點運算要求兩矩陣的維數(shù)相同。5.MATLAB實現(xiàn)將全下標(biāo)轉(zhuǎn)換為單下標(biāo)的指令為Sub2ind、據(jù)單下標(biāo)換算出全下標(biāo)的指令為Ind2sub。6.命令文件可以(不可以/可以)直接運行，函數(shù)文件不可以(不可以/可以)直接運行。7.函數(shù)文件定義的變量為局部變量（局部/全局）。三、程序閱讀（4×10）寫出程序運行結(jié)果。1.

>>A=[3579;0241;6448];

>>B=A([12],2:3)

答案：B=57242.>>a=[6

9

3;2

7

5];

>>

b=[2

4

1;4

6

8];

>>

c=a.*b

答案：c=12

36

3

8

42

40

3.>>A=[6448];>>B=[1110;1110;1110;0005];>>C=[A;B(:,1:3)’]答案：C=64481110111011104.x=[0,1,0,2,0,3,0,4];

for

k=1:8

if

x(k)==0

x(k)=k;

else

x(k)=2*k+1;

end

end

答案：x=1539513717四、編程（2×10）1.輸入三個數(shù)按從大到小的順序輸出。（使用if）答案：a=input('a=?');ifa(1)<a(2)temp=a(1);a(1)=a(2);a(2)=temp;endifa(1)<a(3)temp=a(1); a(1)=a(3);a(3)=temp;endifa(2)<a(3)temp=a(2);a(2)=a(3);a(3)=temp;enda2.請輸入星期幾的第一個字母來判斷一下是星期幾，如果第一個字母一樣，則繼續(xù)判斷第二個字母。f=input('輸入第一個字母（大寫字符）');switchfcase'S's=input('輸入第二個字母（小寫字符）');ifs=='a'disp('Saturday')elseifs=='u'disp('Sunday');elsedisp('dataerror');endbreak;case'F'disp('Friday');break;case'M'disp('Monday');break;case'T's=input('輸入第二個字母（小寫字符）’)；ifs=='u'disp('Tuesday');elseifs=='h'disp('Thursday');elsedisp('dataerror');endbreak;case'W'disp('Wednesday');break;otherwisedisp('dataerror');end第七套填空題（每題5分，共25分）rand(2,3)的含義為生成2×3隨機矩陣。MATLAB語言支持的3種最基本的控制結(jié)構(gòu)為順序結(jié)構(gòu)、選擇結(jié)構(gòu)和循環(huán)結(jié)構(gòu)。3）在循環(huán)結(jié)構(gòu)中跳出當(dāng)前循環(huán)，繼續(xù)下一次循環(huán)的命令為comtinue。4）MATLAB中逗號主要用作要顯示計算結(jié)果的指令與其后指令的分隔；用作輸入量與輸入量之間的分隔符；用作數(shù)組元素分隔符號。5）工作空間瀏覽器主要用于內(nèi)存變量的內(nèi)存變量的查閱、保存和編輯。二、單選題（每題5分，共25分）1.下列哪個變量的定義是不合法的A(A)abcd-3(B)xyz_3(C)abcdef(D)x3yz2.下列哪條指令是求矩陣的行列式的值C(A)inv(B)diag(C)det(D)eig3.用round函數(shù)四舍五入對數(shù)組[2.486.393.938.52]取整，結(jié)果為C(A)[2638](B)[2648](C)[2649](D)[3749]4.已知a=2:2:8,b=2:5，下面的運算表達式中，出錯的為C(A)a'*b(B)a.*b(C)a*b(D)a-b5.下面的程序執(zhí)行后array的值為(A)fork=1:10ifk>6break;elsearray(k)=k;endend(A)array=[1,2,3,4,5,6](B)array=[1,2,3,45,6,7,8,9,10](C)array=6(D)array=10簡答題（20分）1.求積分（5分）2.求微分方程的數(shù)值解（5分）dsolve('D2y*(x-1)+2*Dy+y=0','y(0)=1,Dy(0)=1','x')3.拋物線的方程為,輸入c后，繪制該曲線，尋找曲線的最小值和對應(yīng)的x值。（10分）小明編寫程序如下：clc;clear;closeall;c=input('c=');a=1;x=-2:0.1:2;; %(1)對變量xnum賦為數(shù)組x的個數(shù)y=zeros(1,xnum);y=a*x^2+c*x;→%(2)按照函數(shù)形式對y賦值figure;; %(3)繪制曲線xlabel('x');ylabel('y');[ymin,Imin]=min(y);xmin=x(ymin);→%(4)找到最小函數(shù)值對應(yīng)的x值答案：clc;clear;closeall;c=input('c=');a=1;x=-2:0.1:2;xnum=length(x)orxnum=size(x,2);%(1)對變量xnum賦為數(shù)組x的個數(shù)y=zeros(1,xnum);y=a*x^2+c*x;→y=a*x.^2+c*x;%(2)按照函數(shù)形式對y賦值figure;plot(x,y) %(3)繪制曲線xlabel('x');ylabel('y');[ymin,Imin]=min(y);xmin=x(ymin);→xmin=x(Imin);%(4)找到最小函數(shù)值對應(yīng)的x值編程題（30分）1.自己設(shè)計一個m函數(shù)文件，用來畫出橢圓x*x/16+y*y/9=1。（10分）答案：M文件的程序：%畫橢圓tt=0:0.02:2*pi;xx=4*cos(tt);yy=3*sin(tt);axisequalplot(xx,yy,'m*')title('橢圓x^2/16+y^2/9=1')2.在xy平面內(nèi)選擇區(qū)域[-8,8]×[-8,8]，利用mesh、meshc、meshz和surf繪制。要求用子圖在一個圖形窗口中繪制。答案：x=-8:0.1:8;y=x;[x,y]=meshgrid(x,y);m=sqrt(x.^2+y.^2);z=sin(m)./m;subplot(2,2,1);mesh(x,y,z);subplot(2,2,2);meshc(x,y,z);subplot(2,2,3);meshz(x,y,z);subplot(2,2,4);surf(x,y,z)3.利用menu函數(shù)輸入選擇參數(shù)len，輸入菜單的名稱為“selectannum”，根據(jù)len的取值分別給變量x賦不同的值。當(dāng)len=1時，x的值為的值為4行3列全0矩陣；當(dāng)len=2時，x的值為3行4列全為5的矩陣；當(dāng)len=3時，x的值為5行3列均值為0，方差為1的正態(tài)分布的隨機矩陣；當(dāng)len=4時，x的值為字符串“notanarray”；當(dāng)len為其它值時，則顯示“Error”。要求使用switch函數(shù)。答案：len=menu('selectanum','1','2','3','4','5');%3分switchlen%switch語句格式寫對2分case1y=zeros(4,3);%1分case2y=5.*ones(3,4);%1分case3y=5+randn(5,3);%1分case4y=’notanarray’;%1分otherwisedisp(‘Error’)%1分end第八套一、填空題（一題3分，共15分）1、標(biāo)點符號_______可以使命令行不顯示運算結(jié)果。答案：%2、x

為0

～4

，步長為0.1pi

的向量，使用命令_______創(chuàng)建。

答案：x=0:0.1*pi:4。3、輸入矩陣A=[123;7-58;9-1-7]，使用全下標(biāo)方式用A(2,2)

取出元素“-5

”，使用單下標(biāo)方式用_______取出元素“-5

”。答案：A（5）。在while

表達式,

語句體,

End

循環(huán)語句中，表達式的值__

__時表示循環(huán)條件為真，語句體將被執(zhí)行，否則跳出該循環(huán)語句。答案：非零。要從鍵盤讀入一個字符串并賦值給變量x，且給出提示“Hello”，應(yīng)使用命令_________.答案：x=input(‘Hello’,’s’)。選擇題（一題3分，共15分）1、

下列表示可以作為MATLAB的變量的是(D)

A．a(chǎn)bcd-2

B．xyz_2#

C．@h

D．X_1_a

2、如果x=1:

2

:

8,則x(1)和x(4)分別是(B)

A．1，8

B．1,

7

C．2,

8

D．2,

73、指出下列錯誤的指令（B）。

A.syms

a

b;

B.

syms

a,

b;

C.syms(‘a(chǎn)',

’b');

D.

syms(‘a(chǎn)','b','positive');4、下列哪條指令是求矩陣的行列式的值（C）。

A

.inv

B.

diag

C.

det

D.eig

5、清空

Matlab

工作空間內(nèi)所有變量的指令是（C）。

A.

clc

B.cls

C.clear

D.

clf

三、程序分析題（一題10分，共20分）

1、分析下面程序運行后s1,n1,s2,n2,m的值。

s1=0;s2=0;n1=0;n2=0；

x=[1,-4,-8,3,10,-9,7,-3,10,

8,-5,-2,2,0

];

m=length(x);

for

i=1:m

if

x(i)<0

s1=s1+x(i);

n1=n1+1;

else

s2=s2+x(i);

n2=n2+1;

end

end

s1,n1,s2,n2,m

請回答s1,n1,s2,n2,m的結(jié)果。答案：s1=-31,n1=6,s2=41,n2=8,m=14。2、求出mysum的值mysum=0；i=1；

while(i<=100)

mysum=mysum+i；

i=i+1；

end

mysum

答：

mysum

=

5050

程序改錯題（一題10分，共10分）

下面的語句用來判斷一個人的體溫是否處于危險狀態(tài)（華氏溫度），請找出錯誤并改正。

temp=input('Enter

temperature:')

if

temp<97.5

disp('Temperature

below

normal');

elseif

temp>97.5

改正：temp>97.5

&temp<=99.5

disp('Temperature

normal');

elseif

temp>99.5

改正：temp>99..5

&temp<=103.0

disp('Temperature

slightly

high');

else

temp>103.0

改正：去掉temp>103.0

disp('Temperature

dangerously

high');

end

程序編程題（共40分）1、用MATLAB編程求下列聯(lián)立方程的解

3x+4y-7z-12w=4

；5x-7y+4z+

2w=-3；

x

+8z-

5w=9

；-6x+5y-2z+10w=-8

。求系數(shù)矩陣的秩，求出方程組的解。（5分）答案：a=[3

4

-7

-12;

5

-7

4

2

;1

0

8

-5;

-6

5

-2

10];

c=[4;

-3;

9;-8];

b=rank(a)d=a\c

X=d(1)Y=d(2)Z=d(3)W=d(4)>>

b

=

4

x=-1.4841;y=

-0.6816;z=

0.5337;w=-1.24292、編寫函數(shù)文件求半徑為r的圓的面積和周長。（5分）答案：function[s,p]=fcircle(r)s=pi*r*r;p=2*pi*r;3、Fibonacci數(shù)列的各元素為：1、1、2、3、5、8、……，滿足以下關(guān)系F1=1，F(xiàn)2=1，F(xiàn)n=Fn-1+Fn-2，用M函數(shù)文件實現(xiàn)，數(shù)列的元素個數(shù)為輸入變量，且當(dāng)某個元素大于50時，退出循環(huán)結(jié)構(gòu)。（10分）答案:functionf=shiyan5(n)f(1)=1;f(2)=1;fori=2:nf(i+1)=f(i-1)+f(i);iff(i+1)>50break;endend4、從鍵盤輸入任意個正整數(shù)，以0結(jié)束，輸出那些正整數(shù)中的素數(shù)。（20分）答案：

clc;clear;

zzs(1)=input('請輸入正整數(shù)：

');k=1;

n=0;%素數(shù)個數(shù)

while

zzs(k)~=0

flag=0;%是否是素數(shù)，是則為1

for

yz=2:sqrt(zzs(k))%因子從2至此數(shù)平方根

if

mod(zzs(k),yz)==0

flag=1;break;%非素數(shù)跳出循環(huán)

end

end

if

flag==0&zzs(k)>1%忽略0和1的素數(shù)

n=n+1;sus(n)=zzs(k);

end

k=k+1;

zzs(k)=input('請輸入正整數(shù)：

');

end

disp(['你共輸入了'

num2str(k-1)

'個正整數(shù)。它們是：'])

disp(zzs(1:k-1))%不顯示最后一個數(shù)0

if

n==0

disp('這些數(shù)中沒有素數(shù)！')%無素數(shù)時顯示

else

disp('其中的素數(shù)是：')

disp(sus)

end第九套MATLAB試題一、填空題。（每空1分，共20分）1、MATLAB中(clc)用于清除指令窗中顯示的內(nèi)容，(clear)用于清除工作空間中保存的變量，(clf)用于清除圖形窗口。MATLAB中（help）命令只搜索出關(guān)鍵字完全匹配的結(jié)果，（lookfor）命令對搜索范圍內(nèi)的M文件進行關(guān)鍵字搜索。3、MATLAB變量名必須以（字母）開頭，MATLAB提供的標(biāo)準(zhǔn)函數(shù)名以及命令必須用（小寫字母）。4、t是0~4pi，步長為0.1pi的向量，使用冒號表達式建立向量為(t=0:0.1*pi:4*pi)，使用linspace函數(shù)建立向量為linspace(0,4*pi,10)。2-605、已知矩陣A=31-8,用全下標(biāo)方式(A(3,2))取得元素‘8’，用單586下標(biāo)方式(A(8))取得元素‘-8’.。6、在MATLAB矩陣除法中，A\B等效于(inv(A)*B)，B/A等效于（B*inv(A)）7、A=[12;34]，B=[10;01]，A*B=([12;34])，A.*B=([10;04])。8、為了使兩個plot的圖形在同一個坐標(biāo)顯示，可以用（holdon）命令進行圖形保持，一個圖形窗口上繪制多個圖形可以用（subplot()）函數(shù)進行分割