《數(shù)學實驗》期末總復習題

電子科技大學成都學院二零壹零至二零壹壹學年第一學期

數(shù)樊實驗課程考試題A(120分鐘)閉卷考試時間：

成績構(gòu)成：平時(50%)+期末(50%)=總成績(100%)

一二三合計

注意：請同學們將答案填寫在答題紙上，否則無效。

一、單項選擇題(共40分)

1、三階幻方又稱為九宮圖，提取三階幻方矩陣對角元并構(gòu)造對角陣用()

(A)diag(magic(3));(B)diag(magic);

(C)diag(diag(magic(3)));(D)diag(diag(magic))<>

2、符號計算與一般數(shù)值計算有很大區(qū)別，它將得到準確的符號表達式。在MATLAB命令窗口中鍵入命令

symsxyl=x0；y2=sqrt(x)；int(yl-y2,x,0J)屏幕顯示的結(jié)果是()

(A)yl=xA(l⑵(B)ans=1/3；(C)y2=xA2；(D)ans=-1/3

3、MATLAB命令symsx;f=sin(x);V=pi*int(f*f,x,0,pi)功能是()

(A)繪出函數(shù)f在[0,27]圖形；(B)計算函數(shù)f在[0,2%]的積分；

(C)計算旋轉(zhuǎn)曲面所圍的體積；(D)計算旋轉(zhuǎn)曲面的表面積。

4、用賦值語句給定x數(shù)據(jù)，計算47sin(3+2x)+e3In3對應的MATLAB表達式是()

(A)sqrt(7sin(3+2x)+exp(3)log(3))(B)sqrt(7*sin(3+2*x)+exp(3)*log(3))

(C)sqr(7*sin(3+2*x)+eA3*log(3))(D)sqr(7sin(3+2x)+eA31og(3))

5、MATLAB語句strcat(int2str(2010),年是1s,年)的功能是()

(A)將數(shù)據(jù)2010轉(zhuǎn)換為符號；(B)將數(shù)據(jù)2010與符號變量合并；

(C)將幾個符號變量合并為一個；(D)將符號變量轉(zhuǎn)換為數(shù)值變量；

6、十二屬相為“鼠?；⑼谬埳唏R羊猴雞狗豬"，命令k=mod(2008,12)+1的結(jié)果是()

(A)k指向第二動物牛；(B)k指向第三動物虎；

(C)k指向第四動物兔；(D)k指向第五動物龍。

7、下面有關MATLAB函數(shù)的說法，哪一個是錯誤的()

(A)函數(shù)文件的第一行必須由function開始，并有返回參數(shù)，函數(shù)名和輸入?yún)?shù)；

(B)MATLAB的函數(shù)可以有多個返回參數(shù)和多個輸入?yún)?shù)；

(C)在函數(shù)中可以用nargout檢測用戶調(diào)用函數(shù)時的輸入?yún)?shù)個數(shù)；

(D)如果函數(shù)文件內(nèi)有多個函數(shù)，則只有第一個函數(shù)可以供外部調(diào)用：

8、MATLAB命令x=rand(10,l)生成10個隨機數(shù)，將它們從大到小排序，使用()

(A)y=sort(x);z=y(10:1);(B)[y,II]=sort(x);z=y(II);

(C)y=sort(x);z=y(10:-1;1);(D)[y,II]=sort(x);z=x(II);

9、MATLAB繪三維曲面需要構(gòu)建網(wǎng)格數(shù)據(jù)，語句[x,y]=meshgrid(-2:2)返回數(shù)據(jù)中()

(A)x是行向量，y是列向量；(B)x是列向量，y是行向量；

(C)y是行元素相同的矩陣；(D)x是行向量相同的矩陣

10、設a,b,c表示三角形的三條邊，表達式a+b>c|a+c>b|b+c>a,()

(A)是三條邊構(gòu)成三角形的條件；(B)是三條邊不構(gòu)成三角形的條件；

(C)構(gòu)成三角形時邏輯值為真；(D)不構(gòu)成三角形時邏輯值為假

二、程序閱讀理解(24分)

1、用plot命令繪制平面曲線

v0=515;alpha=45*pi/l80;

T=2*v0*sin(alpha)/9.8;%第二行

t=T*(0:16)/16;

x=v0*t*cos(alpha);

y=vO*t*sin(alpha)-0.5*9.8*t.A2;

plot(x,y,x,y,，r*，)

Xmax=x(17)

8000

6000

4000

2000

0

00.511.522.53

(1)對下面有關程序的功能的說法確切的是()

(A)以515為初速度的拋射曲線的繪制；

(B)以515為初速度的拋射曲線的繪制以及計算射程；

(C)以515為初速度以45度為發(fā)射角的拋射曲線的繪制以及計算射程;

(D)以515為初速度以45度為發(fā)射角的拋射曲線的繪制。

(2)第二行語句的功能是()

(A)計算射程；(B)計算飛行時間；

(C)計算航點的橫坐標；(D)計算航點的縱坐標。

2、關于“牟合方蓋”的實驗程序如下

h=2*pi/100;t=0:h:2*pi;

r=0:0.05:1;x=r'*cos(t);y=r'*sin(t);

z=sqrt(l-x,A2);%第三行

meshz(x,y,z),axisoff

colormap([001])

view(-47,56),holdon

x1=cos(t);y1=sin(t);z1=abs(sin(t));

plot3(xl,yl,zl,'ro')；

(1)下面有關程序的功能的說法確切的是(

(A)繪圓柱面N+y2=i,"+z2=1的交線；

(B)繪圓柱面12+$=1,/+22=1所圍區(qū)域的邊界曲面；

(C)繪圓柱面N+y2=i,N+z2=l的交線及所圍區(qū)域的邊界曲面；

(D)繪圓柱面》2+$=i,#+z2=1的交線及所圍區(qū)域的邊界曲面的上半部分。

(2)關于第三行語句褚誤的解釋是()

(A)z是矩形域上曲頂柱面高度值；(B)z是與y同型的矩陣；

(C)z是圓域上曲頂柱面高度值；(D)z是與x同型的矩陣

3、非負函數(shù)y=f(x)在有限區(qū)間上的圖形為上半平面的一條曲線，曲線繞x軸旋轉(zhuǎn)時，產(chǎn)生以x為對稱軸

的旋轉(zhuǎn)曲面，其體積V=可："(x)F辦.

symsabx

f=exp(a*x)*sin(b*x);

fl=subs(f,a,-0.2);

f2=subs(fl,b,0.5);

V=pi*int(f2*f2,x,0,2*pi)

double(V)

t=(0:20)*pi/10;

theta=t;r=f2(t);

x=t'*ones(size(t));%第九行

y=r'*cos(theta);%第十行

z=r**sin(theta);%第^^一行

mesh(z,y,x)

colormap([000])

axisoff

view(-17,54)

(1)關于程序的功能確切的說法()

(A)計算曲線段f(x)=exp(ax)sin(bx),繞X軸旋轉(zhuǎn)的旋轉(zhuǎn)曲面體積

(B)計算曲線段f(x)=exp(-0.2x)sin(0.5x),0<x<2^繞X軸旋轉(zhuǎn)的旋轉(zhuǎn)曲面體積

(C)計算曲線段f(x)=exp(ax)sin(bx),()4％<2左繞X軸旋轉(zhuǎn)的旋轉(zhuǎn)曲面體積并繪圖

(D)計算曲線段f(x)=exp(-0.2x)sin(0.5x),0<x<27繞X軸旋轉(zhuǎn)的旋轉(zhuǎn)曲面體積

并繪圖

(2)由第9行至第十一行語句可得旋轉(zhuǎn)曲面的方程.)

(A)f(x)=Jy+z+r7(B)f(x)=Jy2+z2

(C)=f(\lx2+z2)(D)y=/(Vx2+z2+/2)

三、程序填空(36分)

1、中國農(nóng)歷60年一大輪回，按天干“甲乙丙丁戊己庚辛壬癸”和地支“子丑寅卯辰巳午未申酉成亥”循環(huán)排

列而成。已知20XX年是農(nóng)歷己丑年，通過簡單計算可以找出年份與天干/地支對應的規(guī)律。下面數(shù)學實驗

程序?qū)斎肽攴荩嬎悴⑤敵鲎址r(nóng)歷紀年。填空完善程序。(8分)

functioncalendar=year(year)

ifnargin==0,year=2009;end

Sl=?甲乙丙丁戊已庚辛壬癸，；

S2d子丑寅卯辰巳午未申酉戍亥，；

k1=mod(year-4,l0)+1;①;%定位天干序數(shù)

sl=Sl(kl);

k2=mod(year-4,12)+1②；%定位地支序數(shù)

s2=S2(k2);

calendar=strcat(int2str(year),'年是',s1,s2,'年')

2.對于二重積分dxdy,。是直線y=x-2和拋物線y2=x所圍區(qū)域。下面數(shù)學實驗程序的

D

功能是利用MATLAB符號計算方法將二重積分處理為二次積分進行計算(選擇先對X積分后對y積分的

秩序)，然后將計算結(jié)果的符號表達式轉(zhuǎn)換為雙精度實數(shù)。(8分)

symsxy;

f=x*y*y;

xl=y*y；x2=y+2①

Sl=int(f,x,xl,x2);

S2=int(Sl,y,-l,2);

S=②

3.某年A、B兩城鎮(zhèn)人各有1000人，以后每年A鎮(zhèn)人口10%遷往B鎮(zhèn);B鎮(zhèn)人口15%遷往A鎮(zhèn)。則有遷

移矩陣LL有兩個互異特征值;I],屬于兩個特征值的特征向量分別為：ai，。2,記初始人口分布

Xo=flOOO,lOOOf,則存在不全為零的數(shù)C1，62使得乂0=，1。1+。2。2(解此方程組求出口，C2)。下面

n

實驗程序主要功能是利用LX0=C\Z"a\+。2丸2”&2計算人口變化。(8分)

n=input(*inputn:-);

A=[0.9,0.15;0.1,0.85];

X0=[1000;1000];

[P,D]=eig(A);

C=double(S2)①

cl=C(I);c2=C(2);

alfal=P(:,l);

alfa2=P(:,2);

lamdal=D(l,l);

lamda2=D(2,2);

Xn=P\X0;②

4、長征三號甲運載火箭提供給探月衛(wèi)星的初始速度不足以將衛(wèi)星送往月球軌道。為提高到奔月速度，中國

航天工程師使用了衛(wèi)星變軌技術(shù)。數(shù)學實驗程序根據(jù)變軌中軌道周期和近地點距離數(shù)據(jù)，利用開普列第二

定律模擬計算計算衛(wèi)星飛行的最大速度。填空完善下面實驗程序。(12分)

R=6378;Time=[l6,15.63,23.3,50.5,225]*3600;

h=[200,600,600,600,600];H=[51000,51000,71000,128000,370000];

a=(h+H+2*R)/2;

c=c1*lamdalAn*alfal+c2*lamda2An*alfa2?

b=sqrt(a.*a-c.*c);

S=sqrt(a.*a-c.*c);

Vmax=2*S./(R+h)

Vmin=S=0.5*Vmax.*(R+h)③

電子科技大學成都學院二零壹零至二零壹壹學年第一學期

數(shù)學實驗課程考試題B(120分鐘)閉卷考試時間：

成績構(gòu)成：平時(50%)+期末(50%)=總成績(100%)

二三合計

注意：請同學們將答案填寫在答題紙上，否則無效。

一、單項選擇題(共40分)

1＞下面有關MATLAB變量名和函數(shù)名的說法中，錯誤的說法是()

(A)變量名的第一個字符必須是一個英文字母

(B)變量名可由英文字母、數(shù)字和下劃線混合組成

(C)變量名不得包含空格和標點，但可以有下連字符

(D)變量名和函數(shù)名對于英文的大小使用沒有區(qū)別

2、符號計算與一般數(shù)值計算有很大區(qū)別，它將得到準確的符號表達式。在MATLAB命令窗口中鍵入命令

symsxyl=xA2；y2=sqrt(x)；int(yl-y2,x,0,l)屏幕顯示的結(jié)果是()

(A)yl=xA(l⑵(B)ans=l/3；(C)y2=xA2；(D)ans=-l/3

3、MATLAB命令symsx;f=sin(x);V=pi*int(f*f,x,O,pi)功能是()

(A)繪出函數(shù)f在［0,27］圖形；(B)計算函數(shù)f在［0,27］的積分；

(C)計算旋轉(zhuǎn)曲面所圍的體積；(D)計算旋轉(zhuǎn)曲面的表面積。

4、用賦值語句給定x數(shù)據(jù)，計算gsin(3+2x)+e31n3對應的MATLAB表達式是()

(A)sqrt(7sin(3+2x)+exp(3)log(3))(B)sqrt(7*sin(3+2*x)+exp(3)*log(3))

(C)sqr(7*sin(3+2*x)+eA3*log(3))(D)sqr(7sin(3+2x)+eA310g(3))

5、MATLAB語句strcat(int2str(2010),'年是，,s,年')的功能是()

(A)將數(shù)據(jù)2010轉(zhuǎn)換為符號；(B)將數(shù)據(jù)2010與符號變量合并；

(C)將凡個符號變量合并為一個；(D)將符號變量轉(zhuǎn)換為數(shù)值變量；

6、十二屬相為“鼠?；⑼谬埳唏R羊猴雞狗豬"，命令k=mod(2008,12)+1的結(jié)果是()

(A)k指向第二動物牛；(B)k指向第三動物虎；

(C)k指向第四動物兔；(D)k指向第五動物龍。

7、MATLAB命令A=magic(3)創(chuàng)建3階幻方矩陣，求A的特征值絕對值最小用()

(A)min(abs(eig(A)));(B)min(eig(abs(A)));(C)min(eig(A));(D)min(abs(A));

8、MATLAB命令x=rand(10,1)生成10個隨機數(shù)，將它們從大到小排序，使用()

(A)y=sort(x);z=y(10:l);(B)[y,H]=sort(x);z=y(n);

(C)y=sort(x);z=y(10:-1;1);(D)[y,II]=sort(x);z=x(II);

9、MATLAB命令[x,y]=meshgrid(l:3);H=l./(x+y-l)產(chǎn)生的矩陣H是()

-11r-123-「123'■11/21/3

(A)222(B)123(C)233(D)1/21/31/4

3331233451/31/41/5

10、設a,b,c表示三角形的三條邊，表達式a+b>c|a+c>b|b+c>a,()

(A)是三條邊構(gòu)成三角形的條件；(B)是三條邊不構(gòu)成二角形的條件;

(C)構(gòu)成三角形時邏輯值為真；(D)不構(gòu)成三角形時邏輯值為假

二、程序閱讀理解(24分)

1、用plot命令繪制平面曲線

v0=515;alpha=45*pi/l80;

T=2*v0*sin(alpha)/9.8;%第二行

t=T*(0:16)/16;

x=v0*t*cos(alpha);

y=v0*t*sin(alpha)-0.5*9.8*t.A2;

plot(x,y,x,y/r*，)

Xmax=x(17)

8000

6000

4000

2000

0

00.511.522.53

(1)對下面有關程序的功能的說法確切的是()

(A)以515為初速度的拋射曲線的繪制；

(B)以515為初速度的拋射曲線的繪制以及計算射程；

(C)以515為初速度以45度為發(fā)射角的拋射曲線的繪制以及計算射程;

(D)以515為初速度以45度為發(fā)射角的拋射曲線的繪制。

(2)第二行語句的功能是()

(A)計算射程；(B)計算飛行時間；

(C)計算航點的橫坐標；(D)計算航點的縱坐標。

2、關于“牟合方蓋”的實驗程序如下

h=2*pi/100;t=0:h:2*pi;

1-0:0.05:1;x=r'*cos(t);y=r'*sin(t);

z=sqrt(l-x.A2);%第三行

meshz(x,y,z),axisoff

colormap([001])

view(-47,56),holdon

x1=cos(t);yl=sin(t);z1=abs(sin(t));

plot3(xl,yl,zl,,ro,)；

(1)下面有關程序的功能的說法確切的是()

(A)繪圓柱面N+y2=1,N+Z2=1的交線；

(B)繪圓柱面N+y2=I,x2+z2=1所圍區(qū)域的邊界曲面；

(C)繪圓柱面N+y2=I,x2+z2=1的交線及所圍區(qū)域的邊界曲面；

(D)繪圓柱面H+y2=i,H+z2=1的交線及所圍區(qū)域的邊界曲面的上半部分。

(2)關于第三行語句錯誤的解釋是()

(A)z是矩形域上曲頂柱面高度值；(B)z是與y同型的矩陣；

(C)z是圓域上曲頂柱面高度值；(D)z是與x同型的矩陣

3、非負函數(shù)y=f(x)在有限區(qū)間上的圖形為上半平面的一條曲線，曲線繞x軸旋轉(zhuǎn)時，產(chǎn)生以x為對稱軸

的旋轉(zhuǎn)曲面，其體積.=?r"(切2小

symsabx”

f=exp(a*x)*sin(b*x);

fl=subs(f,a,-0.2);

f2=subs(fl,b,0.5);

V=pi*int(f2*f2,x,0,2*pi)

double(V)

t=(0:20)*pi/10;

theta=t;r=f2(t);

x=f*ones(size(t));%第九行

y=rf*cos(theta);%第十行

z=rr*sin(theta);%第十一行

mesh(z,y,x)

colormapdO00])

axisoff

view(-17,54)

(3)關于程序的功能確切的說法()

(A)計算曲線段f(x)=exp(ax)sin(bx),0＜兀42萬繞X軸旋轉(zhuǎn)的旋轉(zhuǎn)曲面體積

(B)計算曲線段f(x)=exp(-0.2x)sin(0.5x),0MxM2萬繞X軸旋轉(zhuǎn)的旋轉(zhuǎn)曲面體積

(C)計算曲線段f(x)=exp(ax)sin(bx),0?左427繞X軸旋轉(zhuǎn)的旋轉(zhuǎn)曲面體積并繪圖

(D)計算曲線段f(x)=exp(-0.2x)sin(0.5x),04x42%繞X軸旋轉(zhuǎn)的旋轉(zhuǎn)曲面體積

并繪圖

(4)由第到行至第十一行語句可得旋轉(zhuǎn)曲面的方程.)

(A)/(x)=Jy2+z2+fr(B)f(x)=Jy2+z2

(C)y=/(V%2+z2)(D)y=/(A/X2+z2+t2)

三、程序填空(36分)

1、判斷閏年條件有二

①能被4整除，但不能被100整除:

②能被4整除，又能被400整除。

yeai^inputC'inputyear:=,);

nl=year/4;

n2=year/100;

n3=year/400;

ifnl===fix(nl)&n2~=fix(n2),，不是閏年，)

_____________________①

dispC是閏年)

elseifn1==fix(nI)&n3==fix(n3)

dispC是閏年)

else

disp②

end

2.對于二重積分JJAy2dxdy,。是直線y=x-2和拋物線/=》所圍區(qū)域。下面數(shù)學實驗程序的

D

功能是利用MATLAB符號計算方法將二重積分處理為二次積分進行計算

(選擇先對X積分后對y積分的秩序)，然后將計算結(jié)果的符號表達式轉(zhuǎn)換

為雙精度實數(shù)。(8分)

symsxy;

f=x*y*y;

x1=y*y;x2=y+2

Sl=int(f,x,xl,x2);

S2=int(Sl,y,-l,2);

S=double(S2)

__________________________②

3.某年A、B兩城鎮(zhèn)人各有1000人，以后每年A鎮(zhèn)人口10%遷往B鎮(zhèn);B鎮(zhèn)人口15%遷往A鎮(zhèn)。則有遷

移矩陣L,L有兩個互異特征值為，/U，屬于兩個特征值的特征向量分別為：ai，&2,記初始人口分布

r

X()=[1000,1000],則存在不全為零的數(shù)ci，C2使得X°=ciai+c2a2(解此方程組求出ci，c2)o下面

n

實驗程序主要功能是利用LX0=ci九"a1+C2/l2〃a2計算人口變化。(10分)

n=input(finputn:-);

A=[0.9,0J5;0.1,0.85];

X0=[1000;I000];

[P,D]=eig(A);

C=P\X0;①

cl=C(l);c2=C(2);

alfal=P(:,l);

alfa2=P(:,2);

lamdal=D(l,l);

lamda2=D(2,2);

Xn=cl*lamdalAn*alfal+c2*lamda2An*alfa2

_________________________②

4、為了進入地月轉(zhuǎn)移軌道，嫦娥一號衛(wèi)星進行了四次變軌調(diào)速度。第一次變軌從16小時初始軌道進入16

小時軌道，第二次衛(wèi)星進入24小時軌道，第三次衛(wèi)星進入48小時軌道，第四次衛(wèi)星進入116小時地月轉(zhuǎn)

移軌道。上面小時數(shù)并不是準確軌道周期，變軌目的是將速度從10.3(km/s)逐漸提高到約10.9(km/s)。下面

數(shù)學實驗程序是在區(qū)間[10.3,10.9]上插入線性等分點，即每個軌道的最大速度以等差數(shù)列出現(xiàn)，然后近似計

算出每個軌道的周期參數(shù)。(10分)

functionsatel1()

R=6378;

h=[200,600,600,600,600];

H=[51000,51000,71000,128000,370000];

a=(h+H+2*R)/2;

c=(H-h)/2;

b=(H-h)/2;①；％計算短半軸數(shù)據(jù)

E2=(c./a).A2;

L=2*pi*a.*(1-E2/4-3*E2.A2/64)

formatbank

Vmax=linspace(10.3,10.9,5)

S=pi*a*b/Times②;%根據(jù)最大速度計算每秒鐘掃過的面積

Times二a.*b.*pi./S;

myTimes=Times/3600

電子科技大學成都學院二零壹零至二零壹壹學年第二學期

數(shù)學實驗課程考試題A(120分鐘)閉卷考試時間：

成績構(gòu)成：平時(50%)+期末(50%)=總成績(100%)

二三合計

注意：請同學們將答案填寫在答題紙上，否則無效。

一、單項選擇題(共40分)

1、下面有關MATLAB變量名和函數(shù)名的說法中，正確的說法是()

(A)變量名的第一個字符必須是一個英文字母

(B)變量名不可由英文字母、數(shù)字和下劃線混合組成

(C)變量名可以包含空格、標點和下連字符

(D)變量名和函數(shù)名對于英文的大小使用沒有區(qū)別

2、MATLAB命令x=5:2:100將創(chuàng)建等差數(shù)列，該數(shù)列是()

(A)以5為初值的97個數(shù),；(B)以100為終值的97的個數(shù)；

(C)以99為終值的48個數(shù)；(D)以100為終值的48個數(shù)。

3、theta=linspace(0,2*pi,100);r=cos(5*theta);polar(theta,r；k，)功能是()

(A)繪四葉玫瑰線；(B)繪五葉玫瑰線；(C)繪心臟線；(D)繪十葉玫瑰線。

4、MATLAB繪三維曲面需要構(gòu)建網(wǎng)格數(shù)據(jù)，語句［x,y］=meshgrid(-2:2)返回數(shù)據(jù)中()

(A)x是行向量，y是列向量；(B)x是列向量，y是行向量；

(C)y是行元素相同的矩陣；(D)x是行向量相同的矩陣

5、十二屬相為“鼠?；⑼谬埳唏R羊猴雞狗豬"，命令k=mod(2011-4,12)+1的結(jié)果是()

(A)k指向第二動物牛；(B)k指向第三動物虎；

(C)k指向第四動物兔；(D)k指向第五動物龍。

6、MATLAB命令symsx;f=sin(x);int(f*f,x,0,pi)的功能是()

(A)計算函數(shù)f在［0,2*pi］圖形；(B)計算函數(shù)協(xié)2在［0,2*pi］的積分；

(C)計算旋轉(zhuǎn)曲面所圍的體積；(D)計算旋轉(zhuǎn)曲面的表面積。

7、己知A是一個nXn的可逆矩陣，下面正確的說法是()

(A)rank(A)求A的行列式；(B)inv(A)求A的秩；

(C)sum(A)求A的列和；(D)［P,D］=eig(A)求A的特征值。

8、MATLAB計算正態(tài)分布隨機變量分布函數(shù)的方法是()

(A)binocdf(x,n,p);(B)normpdf(x,mu,p);(C)binopdf(x,n,p);(D)normcdf(xzn,p)o

9、v0=515m/s;515*515*sin(2*alpha)"8的功能是計算關于拋射體問題的()

(A)發(fā)射角；(B)最大高度；(C)最大飛行時間；(D)最大射程。

10、繪制繞Z軸旋轉(zhuǎn)的曲面方程是()

(A)/(%)=\］y2+z2(B)f(y)=y/z2+x2

22

(C)f(z)=Jy2+x。(D)z=yly+x

二、程序閱讀理解（24分）

1.數(shù)學實驗程序如下：

h=439;H=2384;R=6400;

a=(h+H+2*R)/2;c=(H-h)/2;

el=c/a;b=sqrt(a*a-c*c);

symse2t

f=sqrt(1-e2*cos(t)A2);

ft=subs(f,e2,el*el);

S=int(ft,0,pi/2);

L=4*a*double(S)；

V=L/(114*60)

sl=pi*a*b/(l14*60);%第十行

Vmax=2*sl/(h+R)

Vmin=2*sl/(H+R)

（1）實驗程序的運行后，將顯示的數(shù)據(jù)是（)

（A）衛(wèi)星運行的近地速度、遠地速度；(B)衛(wèi)星運行的近地速度、遠地速度和平均速度;

（C）衛(wèi)星運行時向徑每秒掃過的面積：(D)衛(wèi)星運行的平均速度數(shù)據(jù)

（2）第十行語句的功能是（）

（A）計算衛(wèi)星運行的最小速度；(B)計算衛(wèi)星運行時向徑每秒掃過的面積；

（C）計算衛(wèi)星運行的最大速度；(D)計算衛(wèi)星運行軌道的平均速度

2.數(shù)學實驗程序如下：

xy=[-4-4;4-4;44;-44;-4-4];

A=[cos(pi/24)-sin(pi/24);sin(pi/24)cos(pi/24)];

x=xy(:,l);y=xy(:,2);

line(x,y),pause(l)

fork=l:30

xy=89*xy*A,;%第六行

x=xy（:,l）;y=xy（:,2）;

line（x,y）,pause（l）

end

（1）第六行語句中的矩陣A'的功能是（）

（A）正交矩陣；

（B）作用于矩陣xy時，其效果是將xy逆時針旋轉(zhuǎn)為pi?24；

（C）作用于矩陣xy時，其效果是將xy旋轉(zhuǎn)為pi>24,并將xy壓縮為0.89；

（D作用于矩陣xy時，其效果是將xy順時針旋轉(zhuǎn)為pi/24。

（2）對下面有關程序的功能的說法確切的是（）

（A）邊長為4以原點為中心的正方形旋轉(zhuǎn)pi/24重復30次并繪圖；

（B）邊長為8以原點為中心的正方形旋轉(zhuǎn)pi/24并將其邊長壓縮0.89,重復30次并繪圖；（

（C）邊長為4以原點為中心的正方形順時針旋轉(zhuǎn)pi/24并將其邊長壓縮0.89,重復30次并繪圖；

（D）邊長為8以原點為中心的正方形逆時針針旋轉(zhuǎn)pi/24并將其邊長壓縮0.89,重復30次并繪圖。

3、Viviani體是圓柱體（x-R/2）2+y2<R2/4被球面x2+y2+z2=R2所割立體。下面的數(shù)學實驗程序功自

是取R=2求體積上半部分，先利用符號處理重積分并轉(zhuǎn)換為數(shù)值數(shù)據(jù)，再用蒙特卡羅方法計算體做對比。

數(shù)學實驗程序如下:

symsxy;

f=sqrt(4-xA2-yA2);

y1=sqrt(2*x-xA2);y2=sqrt(2*x-xA2);

Sl=int(f,y,y1,y2);S2=int(S1,x,0,2)

V=double(S2)

P=rand(10000,3)；

X=2*P(:1);Y=2*P(:,2);Z=2*P(:,3);

II=find((X-l).A2+Y.A2<=1&Z<=sqrt(4-X.A2-Y.A2));

Vl=8*length(II)/10000

(1)符號計算所用的積分公式是()

⑻丫=同：然）4一/723

(C)V=\+^^-x2-y2dy

J-yj2x-x~

（2）蒙特卡羅方法選用的隨機點變化范圍的立方體區(qū)域是（）

(A)Q={(x,y,z)|xe(0,2),yG(0,2),zG(0,2)}；

(B)。={(x,y,z)|%￡(0,2),ye(-1,1),zG(0,2)}

(C)Q={(x,y,z)10<x<2,0<y<1,0<z<2}

(D)Q={(x,y,z)|—2<x<2,0<y<1,0<z<2}

三、程序填空（24分）

1、判斷閏年條件有二

①能被4整除，但不能被100整除;

②能被4整除，又能被400整除。

year=input（'inputyear:=*）;

nl=year/4;

n2=year/100;

n3=year/400;

if_______________________________

dispC是閏年）

elseif___________________________②

dispC是閏年》

else

dispC不是閏年）

end

2.對于二重積分“dxdy,。是由兩條拋物線y=&和拋物線y=V所圍區(qū)域。下面數(shù)學實驗程

D

序的功能是利用MATLAB符號計算方法將二重積分處理為二次積分進行計算（選擇先對x積分后對y積

分的秩序），然后將計算結(jié)果的符號表達式轉(zhuǎn)換為雙精度實數(shù)。

symsxy;

f=x*sqrt(y);

xl=yA2

x2=一①

Sl=int(f,x,xl,x2);

S2=____________②

S=double（S2）

3.汽車緊急剎車問題數(shù)據(jù)擬合實驗

V:2025303540455055606570

T:202841537293118149182221266

V表示剎車時汽車行駛速度(英里/小時),T表示剎車后汽車滑行距離(英尺)

分別做二次多項式和三次多項式數(shù)據(jù)擬合實驗，并繪出數(shù)據(jù)擬合曲線的圖，計算出殘差平方和，完成如下

實驗程序填空

v=[2025303540455055606570]*1.609;

T=[202841537293118149182221266]*.3048;

P2==polyfit(v,T,2)

T2=polyval(P2,v);

R2=sum((T-T2).A2)

figure(2),plot(v,T；*\v,T2)

P3=?

T3=polyval(P3,v);

R3=②

figure⑶,plot(v,T,'*',v,T3)

四、數(shù)學建模(12分)

每年A鎮(zhèn)的人口10%遷往B鎮(zhèn);B鎮(zhèn)的人口15%遷往A鎮(zhèn).假設某年A、B兩鎮(zhèn)人口各有110人和90人.假

設兩鎮(zhèn)人口流動正常的話，問兩年后兩鎮(zhèn)人口數(shù)量分布如何?尋找方案使每年人口正常流動而A鎮(zhèn)和B鎮(zhèn)

的人口數(shù)量不增不減。

電子科技大學成都學院二零壹零至二零壹壹學年第一學期

數(shù)學實驗課程考試題B(120分鐘)閉卷考試時間：

成績構(gòu)成：平時(50%)+期末(50%)=總成績(100%)

一二三合計

注意：請同學們將答案填寫在答題紙上，否則無效。

一、單項選擇題(共40分)

1、下面有關MATLAB變量名和函數(shù)名的說法中，錯誤的說法是()

(A)變量名的第一個字符必須是一個英文字母

(B)變量名可由英文字母、數(shù)字和下劃線混合組成

(C)變量名不得包含空格和標點，但可以有下連字符

(D)變量名和函數(shù)名對于英文的大小使用沒有區(qū)別

2、符號計算與一般數(shù)值計算有很大區(qū)別，它將得到準確的符號表達式。在MATLAB命令窗口中鍵入命令

symsxyl=xA2；y2=sqrt(x)；int(yl-y2,x,0,l)屏幕顯示的結(jié)果是()

(A)yl=xA(l/2)(B)ans=l/3；(C)y2=xA2；(D)ans=-l/3

3、MATLAB命令symsx;f=sin(x);V=pi*int(f*f,x,O,pi)功能是()

(A)繪出函數(shù)f在[0,2萬]圖形；(B)計算函數(shù)f在[0,2%]的積分；

(C)計算旋轉(zhuǎn)曲面所圍的體積；(D)計算旋轉(zhuǎn)曲面的表面積。

4、用賦值語句給定x數(shù)據(jù)，計算J；sin(3+2x)+e31n3對應的MATLAB表達式是()

(A)sqrt(7sin(3+2x)+exp(3)log(3))(B)sqrt(7*sin(3+2*x)+exp(3)*log(3))

(C)sqr(7*sin(3+2*x)+eA3*log(3))(D)sqr(7sin(3+2x)+eA31og(3))

5、MATLAB語句strcat(int2str(2010),'年是；s,,年，)的功能是()

(A)將數(shù)據(jù)2010轉(zhuǎn)換為符號；(B)將數(shù)據(jù)2010與符號變量合并；

(C)將幾個符號變量合并為一個；(D)將符號變量轉(zhuǎn)換為數(shù)值變量；

6、十二屬相為“鼠?；⑼谬埳唏R羊猴雞狗豬"，命令k=mod(2008,12)+1的結(jié)果是()

(A)k指向第二動物牛；(B)k指向第三動物虎；

(C)k指向第四動物兔；(D)k指向第五動物龍。

7、MATLAB命令A=magic(3)創(chuàng)建3階幻方矩陣，求A的特征值絕對值最小用()

(A)min(abs(eig(A)));(B)min(eig(abs(A)));(C)min(eig(A));(D)min(abs(A));

8、MATLAB命令x=rand(10,l)生成10個隨機數(shù)，將它們從大到小排序，使用()

(A)y=sort(x);z=y(10:l);(B)[y,II]=sort(x);z=y(II);

(C)y=sort(x);z=y(10:-1;1);(D)[y,II]=sort(x);z=x(H);

9、MATLAB命令[x,y]=meshgrid(l:3);H=l./(x+y-l)產(chǎn)生的矩陣H是()

-iir-123--123-■11/21/3

(A)222(B)123(C)233(D)1/21/31/4

3331233451/31/41/5

10、設a,b,c表示三角形的三條邊，表達式a+b>c|a+c>b|b+c>a,()

(A)是三條邊構(gòu)成三角形的條件；(B)是三條邊不構(gòu)成三角形的條件;

(C)構(gòu)成三角形時邏輯值為真；(D)不構(gòu)成三角形時邏輯值為假

二、程序閱讀理解(24分)

1、用plot命令繪制平面曲線

v0=515;alpha=45*pi/l80;

T=2*v0*sin(alpha)/9.8;%第二行

t=T*(0:16)/16;

x=v0*t*cos(alpha);

y=v0*sin(alpha)?0.5*9.8*t.八2;

plot(x,y,x,y,,r*>)

Xmax=x(17)

8000

6000

4000

2000

00.511.522.53

(1)對下面有關程序的功能的說法確切的是()

(A)以515為初速度的拋射曲線的繪制；

(B)以515為初速度的拋射曲線的繪制以及計算射程；

(C)以515為初速度以45度為發(fā)射角的拋射曲線的繪制以及計算射程;

(D)以515為初速度以45度為發(fā)射角的拋射曲線的繪制。

(2)第二行語句的功能是()

(A)計算射程；(B)計算飛行時間；

(C)計算航點的橫坐標；(D)計算航點的縱坐標。

2、關于“牟合方蓋”的實驗程序如下

h=2*pi/100;t=0:h:2*pi;

r=0:0.05:l;x=r'*cos(t);y=r'*sin(t);

z=sqrt(l-x.A2);%第三行

meshz(x,y,z),axisoff

colortnap([001])

view(-47,56),holdon

x1=cos(t);y1=sin(t);z1=abs(sin(t));

plot3(xl,yl,zl,'ro')；

(1)下面有關程序的功能的說法確切的是()

(A)繪圓柱面x2+y2=l,x2+z2=1的交線；

(B)繪圓柱面x2+y=i,x2+z2=1所圍區(qū)域的邊界曲面；

(C)繪圓柱面N+y2=i,%2+z2=1的交線及所圍區(qū)域的邊界曲面；

(D)繪圓柱面x2+y2=i,%2+z2=1的交線及所圍區(qū)域的邊界曲面的上半部分。

(2)關于第三行語句錯誤的解釋是()

(A)z是矩形域上曲頂柱面高度值；(B)z是與y同型的矩陣；

(C)z是圓域上曲頂柱面高度值；(D)z是與x同型的矩陣

3、非負函數(shù)y=f(x)在有限區(qū)間上的圖形為上半平面的一條曲線，曲線繞x軸旋轉(zhuǎn)時，產(chǎn)生以x為對稱軸

的旋轉(zhuǎn)曲面，其體積

symsabx

f=exp(a*x)*sin(b*x);

fl=subs(f,a,-0.2);

f2=subs(fl,b,0.5);

V=pi*int(f2*f2,x,0,2*pi)

double(V)

t=(0:20)*pi/10;

theta=t;r=f2(t);

x=tF*ones(size(t));%第九行

y=r**cos(theta);%第十行

z=r'*sin(lhela);%第^^一行

mesh(z,y,x)

colormap([000])

axisoff

view(-17,54)

(5)關于程序的功能確切的說法()

(A)計算曲線段f(x)=exp(ax)sin(bx),()《xw2乃繞X