數(shù)學(xué)實(shí)驗(yàn)11091

1、下載文檔收藏?cái)?shù)學(xué)實(shí)驗(yàn)試題(2009)A 期末復(fù)習(xí)題 期末復(fù)習(xí)題<<隱藏 窗體頂端窗體底端電子科技大學(xué)二零零八到二零零九學(xué)年第二學(xué)期期末考試數(shù)學(xué)實(shí)驗(yàn)課程考試題 A 卷 (120 分鐘) 考試形式：閉卷 考試日期：2009 年 7 月 8 日 課程成績構(gòu)成：平時(shí) 10 分，期中 0 分，實(shí)驗(yàn) 30 分，期末 60 分(本試卷滿分 100 分) 所有答案一律寫在答題紙上，寫在試卷上無效。 一、單項(xiàng)選擇題(20 分) 1、三階幻方又稱為九宮圖，提取三階幻方矩陣對角元并構(gòu)造對角陣用( C ) (A) diag(magic(3); (B) diag(magic); (C) diag(diag(

2、magic(3); (D) diag(diag(magic)。 2、MATLAB 命令 P=pascal(3)將創(chuàng)建三階帕斯卡矩陣，max(P)的計(jì)算結(jié)果是( D ) (A) 1 2 3 (B) 1 2 1 (C) 3 6 10 (D) 1 3 6 3、命令 J=1;1;1*1,2,3;A=j+j-1 將創(chuàng)建矩陣( A ) ?1 2 3? ? ? (A) 2 3 4 ; ? ? ?3 4 5? ? ? ?2 3 4? ? ? (B) 3 4 5 ? ? ?4 5 6? ? ? ?1 2 3? ? ? (C) 1 2 3 ? ? ?1 2 3? ? ? ?1 1 1 ? ? ? (D) 2 2

3、2 ? ? ?3 3 3? ? ? 4、data=rand(1000,2);x=data(:,1);y=data(:,2);II=find(y<sqrt(x)&y>x.2);的功能是( B ) (A) 統(tǒng)計(jì) 2000 個(gè)隨機(jī)點(diǎn)中落入特殊區(qū)域的點(diǎn)的索引值； (B) 統(tǒng)計(jì) 1000 個(gè)隨機(jī)點(diǎn)落入特殊區(qū)域的點(diǎn)的索引值； (C) 模擬 2000 個(gè)隨機(jī)點(diǎn)落入特殊區(qū)域的過程； (D) 模擬 1000 個(gè)隨機(jī)點(diǎn)落入特殊區(qū)域的過程。 5、MATLAB 計(jì)算二項(xiàng)分布隨機(jī)變量分布律的方法是( C ) (A) binocdf(x,n,p); (B) normpdf(x,mu,s); (C)bi

4、nopdf(x,n,p); (D) binornd(x,n,p)。 6、MATLAB 命令 syms e2;f=sqrt(1-e2*cos(t)2);S=int(f,t,0,pi/2)功能是(D) (A) 計(jì)算 f(x)在0,pi/2上的積分； (B) 計(jì)算 f(t)不定積分符號結(jié)果； (C) 計(jì)算 f(x)積分的數(shù)值結(jié)果； (D) 計(jì)算 f(t)定積分的符號結(jié)果。 7、y=dsolve(Dy=1/(1+x2)-2*y2,y(0)=0,x);ezplot(y)的功能是( A ) (A) 求微分方程特解并繪圖； (B) 解代數(shù)方程 (C) 求定積分； (D)求微分方程通解。 8、X=10000

5、;0.5*asin(9.8*X/(5152)的功能是計(jì)算關(guān)于拋射體問題的(A) (A) 十公里發(fā)射角； (B) 十公里飛行時(shí)間； (C)最大飛行時(shí)間； (D)最大射程。 9、theta=linspace(0,2*pi,100) ;r=cos(4*theta) ;polar(theta,r,k)功能是(D) (A) 繪四葉玫瑰線； (B)繪三葉玫瑰線； (C)繪心臟線； (D) 繪八葉玫瑰線。 10、北京和紐約的經(jīng)度分別是：東經(jīng) 118 和西經(jīng) 76，根據(jù)經(jīng)度差計(jì)算時(shí)差用(D) (A) fai1=118;fai2=-76;Dfai=(fai1+fai2)/24; (B) fai1=118;fai

6、2=-76;Dfai=(fai1+fai2)/15; (C) fai1=118;fai2=-76;Dfai=(fai1-fai2)/24; (D) fai1=118;fai2=-76;Dfai=(fai1-fai2)/15。 二、程序閱讀題 (40 分) 1、直方圖功能是將數(shù)據(jù)分為 n 個(gè)類，統(tǒng)計(jì)各個(gè)類的數(shù)據(jù)量并繪圖。借用現(xiàn)有的直方圖命令 hist，編寫新直方圖程序如下。 function m=myhist(data,n) if nargin=1,n=7;end Xmin=min(data);Xmax=max(data);h=(Xmax-Xmin)/n; m=hist(data,n)/leng

7、th(data)/h; t=linspace(Xmin,Xmax,n+1); II=1:4:4*n-3;JJ=1:n; x(II)=t(JJ);y(II)=zeros(1,n); x(II+1)=t(JJ);y(II+1)=m; x(II+2)=t(JJ+1);y(II+2)=m; x(II+3)=t(JJ+1);y(II+3)=zeros(1,n); plot(x,y,'k') (1) 變量 data 存放了 1000 個(gè)數(shù)據(jù)，在命令窗口調(diào)用 myhist(data)的結(jié)果是( C ) (A) 只繪數(shù)據(jù)的直方圖而不顯示被分類后各類的數(shù)據(jù)量； (B) 只顯示被分類后各類的數(shù)據(jù)量

8、而不繪數(shù)據(jù)的直方圖； (C) 既繪數(shù)據(jù)直方圖也顯示被分類后各類的數(shù)據(jù)量； (D) 根據(jù)默認(rèn)值在數(shù)據(jù)范圍內(nèi)插入七等分點(diǎn)繪直方圖。 (2) 關(guān)于新直方圖繪圖程序下面說法不正確的是( D ) (A)h 是 n 等分直方圖中小區(qū)間長度； (B) 修改程序最后一行可繪紅色直方圖； (C) 直方圖中所有小矩形面積之和為 1； (D) 直方圖中所有小矩形的高度和為 1。 2、3n+1 問題反映一個(gè)數(shù)學(xué)猜想：對任一自然數(shù) n，按如下法則進(jìn)行運(yùn)算：若 n 為偶數(shù)，則 將 n 除 2，若 n 為奇數(shù)，則將 n 乘 3 加 1。重復(fù)這種操作，結(jié)果終會(huì)為 1。實(shí)驗(yàn)程序如下。 function k,N=threeN(n

9、) if nargin=0,n=5;end k=1;N=n; while n=1 r=rem(n,2); if r=0 n=n/2; else n=3*n+1; end N=N,n;k=k+1; end (1)在 MATLAB 命令窗口中直接調(diào)用 threeN 運(yùn)行結(jié)果為( A ) %5 16 8 4 2 1 (A)只顯示 k 的最后數(shù)值為 6； (B) 只顯示 k 的最后數(shù)值 5； (C) 同時(shí)顯示 k 和 N 的數(shù)據(jù)； (D) 僅顯示 N 的所有數(shù)據(jù)。 (2)實(shí)驗(yàn)程序運(yùn)行過程中( B ) (A) 輸入變量 n 不發(fā)生改變； (B)N 是記錄數(shù)據(jù)變化的一維數(shù)組； (C) N 記錄每次數(shù)據(jù)變化

10、的單個(gè)數(shù)據(jù)； (D)n 是記錄數(shù)據(jù)變化的一維數(shù)組。 3、將半徑為 r 的球體(密度 < 1 )置入水中，球體將浮出水面一定高度。程序如下： function h,Rou=highNu(r) if nargin=0,r=10;end Rou=0.3:0.1:1; N=length(Rou); for k=1:N rouk=Rou(k); P=1,-3*r,0,4*r3*rouk; x=roots(P); II=find(x<2*r&x>0); h(k)=2*r-x(II); end (1)在 MATLAB 命令窗口省略輸入調(diào)用函數(shù) highNu 將顯示( A ) (A)

11、 球體浮出水面的高度數(shù)據(jù)； (B) 球體的 8 個(gè)不同的密度數(shù)據(jù)； (C) 球體沉入水下的深度數(shù)據(jù)； (D) 深度數(shù)據(jù)和密度數(shù)據(jù)。 (2) 程序中變量 x 存入如下方程的根( A ) (A) x ? 3rx + 4r 3 2 3 = 0 ； (B) x 3 ? 3r = 4r 3 ； (C) x 3 ? 3r = 0 ； (D) 1 ? 3rx 3 = 0 -9/(-9+4 exp(-t) 2 1.8 1.6 1.4 1.2 1 0.8 4、一階常微分方程確定一個(gè)平面向量場，初值條件確定了向量場中一條曲線。程序如下： x,y=meshgrid(0:.25:6,0:.05:2); k=y.*(1

12、-y); d=sqrt(1+k.2); px=1./d;py=k./d; quiver(x,y,px,py),hold on u=dsolve('Du=u*(1-u)','u(0)=.2'); v=dsolve('Dv=v*(1-v)','v(0)=1.8'); ezplot(u,0,6) ezplot(v,0,6) 0.6 0.4 0.2 0 0 1 2 3 4 5 6 t 圖 1 向量場圖 (1) 程序中所繪向量場對應(yīng)的一階常微分方程是( B 2 (A) y = 1 + y ； (B) y = y (1 ? y ) ； ) (

13、D) y = 1 / y (1 ? y ) 。 2 (C) y = 1 / 1 + y ； (2) 關(guān)于實(shí)驗(yàn)程序下面說法錯(cuò)誤的是( A ) (A) 程序中第一個(gè)初值條件所對應(yīng)的解曲線在圖 1 中上方； (B) 程序中第二個(gè)初值條件所對應(yīng)的解曲線在圖 1 中上方； (C) 程序繪圖原理是根據(jù)每一點(diǎn)處曲線切線的單位向量繪圖； (D) 當(dāng)初值數(shù)據(jù)大于 1 時(shí)解曲線單調(diào)減少，當(dāng)初值數(shù)據(jù)小于 1 時(shí)解曲線單調(diào)增加。 5、維維安尼體由柱面切割球體所得。下面程序的功能是演示柱面切割球體的過程。 function viviani(dt) if nargin=0,dt=10;end N=fix(360/dt);

14、 X,Y,Z=sphere(N); mesh(X,Y,Z),hold on x,y,z=cylinder(1,1,N); y=.5*y;x=.5*(1-x);z(1,:)=-ones(1,N+1); for p=10:N+1 II=1:p;u=x(:,II); v=y(:,II);w=z(:,II); mesh(u,v,w),pause(.5) end (1) 根據(jù)程序中語句，所繪圖形中( D ) (A) 圓柱的半徑為 1； (B) 圓柱的高度為 1； (C) 圓柱以 Z 軸對稱； (D) 圓柱的高度為 2。 (2) 關(guān)于實(shí)驗(yàn)程序以下錯(cuò)誤的說法是( D ) (A) 程序中輸入變量 dt 大則球

15、面網(wǎng)格線?。?(B) 程序正常運(yùn)行時(shí)球面圖形保持不變； (C) 程序繪圖時(shí)每半秒種圖形變動(dòng)一次； (D) 每循環(huán)一次只加繪柱面一條母線。 三、程序填空(40 分) 1、中國農(nóng)歷 60 年一大輪回，按天干“甲乙丙丁戊已庚辛壬癸”和地支“子丑寅卯辰巳午未申 酉戍亥”循環(huán)排列而成。已知 2009 年是農(nóng)歷已丑年，通過簡單計(jì)算可以找出年份與天干/地 支對應(yīng)的規(guī)律。 下面數(shù)學(xué)實(shí)驗(yàn)程序?qū)斎肽攴荩?計(jì)算并輸出字符串農(nóng)歷紀(jì)年。 填空完善程序。 function calendar=year(year) if nargin=0, year=2009;end S1= 甲乙丙丁戊已庚辛壬癸; S2=子丑寅卯辰巳午未

16、申酉戍亥; k1= mod(S1-4,10)+1 ; %定位天干序數(shù)%2010 庚寅 s1=S1(k1); k2= mod(S2-4,12) +1 ; %定位地支序數(shù) s2=S2(k2); calendar=strcat(int2str(year),年是,s1,s2,年) 2、紅、綠兩隊(duì)從相距 100 公里的地點(diǎn)同時(shí)出發(fā)相向行軍。紅隊(duì)速度為 10(公里/小時(shí))，綠隊(duì) 速度為 8(公里/小時(shí))。開始時(shí)，通訊員騎摩托從紅隊(duì)出發(fā)為行進(jìn)中的兩隊(duì)傳遞消息。摩托車 的速度為 60(公里/小時(shí))往返于兩隊(duì)之間。每遇一隊(duì)，立即回駛向另一隊(duì)。當(dāng)兩隊(duì)距離小于 0.2 公里時(shí)，摩托車停止，下面數(shù)學(xué)實(shí)驗(yàn)程序模擬計(jì)算摩

17、托車跑了多少趟。請?zhí)羁胀晟瞥绦颉?function k=moto(A,B) if nargin=0,A=0;B=100;end va=10;vb=8;vc=60; f=1;k=0; while (B-A)>0.2 if f=1 tk=(B-A)/(vb+vc); else tk= (B-A)/(vc+va); ; %計(jì)算 A 與 C 相遇時(shí)間 end A= A+va*tk ; %計(jì)算 A 點(diǎn)位置 B= B-vb*tk ; %計(jì)算 B 點(diǎn)位置 f=-f; k=k+1; end 3、為了進(jìn)入地月轉(zhuǎn)移軌道，嫦娥一號衛(wèi)星進(jìn)行了四次變軌調(diào)速度。第一次變軌從 16 小時(shí)初 始軌道進(jìn)入 16 小時(shí)軌道

18、，第二次衛(wèi)星進(jìn)入 24 小時(shí)軌道，第三次衛(wèi)星進(jìn)入 48 小時(shí)軌道，第 四次衛(wèi)星進(jìn)入 116 小時(shí)地月轉(zhuǎn)移軌道。 上面小時(shí)數(shù)并不是準(zhǔn)確軌道周期， 變軌目的是將速度 從 10.3(km/s)逐漸提高到約 10.9(km/s)。下面數(shù)學(xué)實(shí)驗(yàn)程序是在區(qū)間10.3,10.9上插入線性等 分點(diǎn)，即每個(gè)軌道的最大速度以等差數(shù)列出現(xiàn)，然后近似計(jì)算出每個(gè)軌道的周期參數(shù)。填空 完善程序。 function satel1() 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= s

19、qrt(a.2-c.2) E2=(c./a).2; L=2*pi*a.*(1-E2/4-3*E2.2/64) format bank Vmax=linspace(10.3,10.9,5) S= 0.5*Vmax.*(R+h) 鐘掃過的面積 Times=a.*b.*pi./S; myTimes=Times/3600 ; %計(jì)算短半軸數(shù)據(jù) ; %根據(jù)最大速度計(jì)算每秒 4、冰淇淋錐的下部為圓錐面，上部為半球面。計(jì)算體積的蒙特卡羅方法是在包含冰淇淋的 六面體內(nèi)產(chǎn)生 N 個(gè)均勻分布的隨機(jī)點(diǎn)， 并統(tǒng)計(jì)落入錐體內(nèi)的隨機(jī)點(diǎn)的數(shù)目 m。 根據(jù)比值 m/N 和六面體體積數(shù)據(jù)計(jì)算出錐體體積數(shù)據(jù)，這種隨機(jī)統(tǒng)計(jì)方法會(huì)產(chǎn)

20、生誤差，根據(jù)大數(shù)定律，誤 差變量服從正態(tài)分布。 下面數(shù)學(xué)實(shí)驗(yàn)程序使用上面二題中第 1 小題繪出誤差直方圖與正太分 布密度函數(shù)比較，填空完善程序。 function mu,sagma=monterror(L) if nargin=0,L=1000;end for k=1:L P=rand(2000,3); x=2*P(:,1)-1; y=2*P(:,2)-1; z= 2*P(:,3) ; %計(jì)算隨機(jī)點(diǎn)Z坐標(biāo)數(shù)據(jù) R2=x.2+y.2; R= sqrt(R2) ; %計(jì)算隨機(jī)點(diǎn)到坐標(biāo)原點(diǎn)距離 II=find(z>=R&z<=1+sqrt(1-R2); m=length(II);

21、 q(k)=8*m/2000; end X=q-pi; mu=mean(X); sagma=sqrt(sum(X-mu).2)/(L-1); myhist(X,7);hold on x=linspace(-3*sagma,3*sagma,50); y= normpdf(x,mu,samga) ; %計(jì)算正態(tài)分布密度函數(shù) 值 plot(x,y,'r') 5 4.5 4 3.5 3 2.5 2 1.5 1 0.5 0 -0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 附參考答案： 一、單項(xiàng)選擇題(每小題 2 分共 20 分) CDAB CDAA DD 二、程序閱讀

22、題(每小題 4 分共 40 分) 1、CD 2、AB 3、AA 4、BB 5、DD 三、程序填空(每小題 4 分共 40 分) 1、mod(year-4,10)+1; mod(year-4,12)+1; 2、(B-A)/(vc+va); A+va*tk; B-vb*tk; 3、sqrt(a.*a-c.*c); (R+h).*Vmax/2; 4、2*P(:,3); sqrt(R2); normpdf(x,mu,sagma) 或 exp(-(x-mu).2/(2*sagma2)/(sqrt(2*pi)*sagma) ? 1 e 2 ( x? )2 2 2 下載文檔收藏?cái)?shù)學(xué)實(shí)驗(yàn)試題A(2010) 數(shù)

23、學(xué)實(shí)驗(yàn)考題 數(shù)學(xué)實(shí)驗(yàn)考題<<隱藏 窗體頂端窗體底端學(xué)院 姓名 學(xué)號 任課老師 選課號 密封線以內(nèi)答題無效 電子科技大學(xué)成都學(xué)院二零零 至二零一零 電子科技大學(xué)成都學(xué)院二零零 九 至二零一零 學(xué)年第 二 學(xué)期期 末 考試 成都學(xué)院數(shù)學(xué)實(shí)驗(yàn) 課程考試題 A 卷（120 分鐘） 考試形式： 閉 考試日期 數(shù)學(xué)實(shí)驗(yàn)?zāi)?月 日 課程成績構(gòu)成：平時(shí) 10 分， 期中 0 分， 實(shí)驗(yàn) 30 分，期末 60 分（本試卷 100 分） 一 二 三 四 合計(jì) 一、單項(xiàng)選擇題（共 30 分） 1符號計(jì)算與一般數(shù)值計(jì)算有很大區(qū)別，它將得到準(zhǔn)確的符號表達(dá)式。在 MATLAB 命令窗口中鍵 入命令 syms

24、x，y1=x2；y2=sqrt(x);；int(y1-y2,x,0,1)，屏幕顯示的結(jié)果是（ ） （A）y1 =x(1/2) （B）ans= 1/3； （C）y2 =x2； （D）ans= -1/3 2在 MATLAB 命令窗口中鍵入命令 B=8,1,6;3,5,7;4,9,2；B*B(:,1)。結(jié)果是（ ） （A）ans= （B）ans= （C）ans= （D）ans= 91 67 67 67 67 91 67 67 67 67 91 67 3在 MATLAB 命令窗口中，鍵入命令 syms x,F=1/(2+cos(2x);ezplot(diff(F)，結(jié)果是（ ） （A）繪出函數(shù) F 在

25、0， 的圖形； （B）繪出函數(shù) F 在2，2 的圖形； （C）繪函數(shù) F 的導(dǎo)函數(shù)在0， 的圖形； （D）繪函數(shù) F 的導(dǎo)函數(shù)在2，2 的圖形 4在 MATLAB 命令窗口中，鍵入命令 binopdf(3,6,1/2),結(jié)果是 （A）ans= (B) ans= (C) ans= (D) ans= 15/64 5/16 3/32 1 5用賦值語句給定 x 數(shù)據(jù)，計(jì)算 7 sin(3 + 2 x ) + e ln 3 對應(yīng)的 MATLAB 表達(dá)式是（ 3 ） （A）sqrt(7sin(3+2x)+exp(3)log(3) （B）sqrt(7*sin(3+2*x)+exp(3)*log(3) （C）

26、sqr(7*sin(3+2*x)+e3*log(3) （D）sqr(7sin(3+2x)+ e3log(3) 6MATLAB 語句 strcat(int2str(2010),'年是', s,'年')的功能是（ ） （A）將數(shù)據(jù) 2010 轉(zhuǎn)換為符號； （B）將數(shù)據(jù) 2010 與符號變量合并； （C）將幾個(gè)符號變量合并為一個(gè)； （D）將符號變量轉(zhuǎn)換為數(shù)值變量； 7語句 L=sqrt(pi); x=fix(1000*L)/1000 的功能是（ ） （A）將無理數(shù) 取三位近似； （B）將無理數(shù) （C）將無理數(shù) 取三位近似數(shù)； 8MATLAB 繪三維曲面需要構(gòu)建網(wǎng)格數(shù)據(jù)

27、，語句x,y=meshgrid(-2:2)返回?cái)?shù)據(jù)中（ （A）x 是行向量，y 是列向量； （B）x 是列向量，y 是行向量； （C）x 是行元素相同的矩陣； （D）x 是行向量相同的矩陣 9下面有關(guān) MATLAB 函數(shù)的說法，哪一個(gè)是錯(cuò)誤的（ ） （A）函數(shù)文件的第一行必須由 function 開始，并有返回參數(shù)，函數(shù)名和輸入?yún)?shù)； （B）MATLAB 的函數(shù)可以有多個(gè)返回參數(shù)和多個(gè)輸入?yún)?shù)； （C）在函數(shù)中可以用 nargout 檢測用戶調(diào)用函數(shù)時(shí)的輸入?yún)?shù)個(gè)數(shù)； （D）如果函數(shù)文件內(nèi)有多個(gè)函數(shù)，則只有第一個(gè)函數(shù)可以供外部調(diào)用； 10設(shè) a,b,c 表示三角形的三條邊，表達(dá)式 a+b>

28、;c|a+c>b|b+c>a， （ ） 取四位近似數(shù) （D）將有理數(shù) 取四位近似數(shù) ） 第 1 頁 共 5頁 學(xué)院 姓名 學(xué)號 任課老師 選課號 密封線以內(nèi)答題無效 （A）是三條邊構(gòu)成三角形的條件； （C）構(gòu)成三角形時(shí)邏輯值為真； （B）是三條邊不構(gòu)成三角形的條件； （D）不構(gòu)成三角形時(shí)邏輯值為假 二、程序閱讀理解（24 分） 1用 plot 命令繪制平面曲線 v0=515;alpha=45*pi/180; T=2*v0*sin(alpha)/9.8; %第二行 t= T*(0:16) /16; x=v0 *t *cos(alpha); y=v0 *t *sin(alpha) -0

29、.5*9.8*t.2; plot(x,y,x,y,r*) Xmax=x(17) 8000 6000 4000 2000 0 0 0.5 1 1.5 2 2.5 3 （1）對下面有關(guān)程序的功能的說法確切的是（ ） （A）以 515 為初速度的拋射曲線的繪制； （B）以 515 為初速度的拋射曲線的繪制以及計(jì)算射程； （C）以 515 為初速度以 45 度為發(fā)射角的拋射曲線的繪制以及計(jì)算射程； （D）以 515 為初速度以 45 度為發(fā)射角的拋射曲線的繪制。 （2）第二行語句的功能是（ ） （A）計(jì)算射程； （B）計(jì)算飛行時(shí)間； （C）計(jì)算航點(diǎn)的橫坐標(biāo)； （D）計(jì)算航點(diǎn)的縱坐標(biāo)。 2關(guān)于“牟合方蓋

30、”的實(shí)驗(yàn)程序如下 h=2*pi/100;t=0:h:2*pi; r=0:0.05:1;x=r'*cos(t);y=r'*sin(t); z=sqrt(1-x.2); %第三行 meshz(x,y,z),axis off colormap(0 0 1) view(-47,56),hold on x1=cos(t);y1=sin(t);z1=abs(sin(t); plot3(x1,y1,z1,'ro')； （1）下面有關(guān)程序的功能的說法確切的是（ ） （A）繪圓柱面 x2 + y2 = 1, x2 + z2 = 1 的交線； （B）繪圓柱面 x2 + y2 = 1

31、, x2 + z2 = 1 所圍區(qū)域的邊界曲面； （C）繪圓柱面 x2 + y2 = 1, x2 + z2 = 1 的交線及所圍區(qū)域的邊界曲面； （D）繪圓柱面 x2 + y2 = 1, x2 + z2 = 1 的交線及所圍區(qū)域的邊界曲面的上半部分。 （2）關(guān)于第三行語句錯(cuò)誤的解釋是（ ） （A）z 是矩形域上曲頂柱面高度值； （B）z 是與 y 同型的矩陣； （C）z 是圓域上曲頂柱面高度值； （D）z 是與 x 同型的矩陣 3中國農(nóng)歷年由天干(10 干)和地支(12 支)相配而成，計(jì)算農(nóng)歷年的 MATLAB 程序如下 year=input('input year:=');

32、第 2 頁 共 5頁 學(xué)院 姓名 學(xué)號 任課老師 選課號 密封線以內(nèi)答題無效 S1='辛壬癸甲乙丙丁戊己庚'； S2='酉戍亥子丑寅卯辰巳午未申' k1=mod(year,10); k2=mod(year,12); if k1=0,k1=10;end %第五行 if k2=0,k2=12;end s1=S1(k1); s2=S2(k2)； strcat(int2str(year),'年是', s1,s2,'年') （1）輸入 2010，實(shí)驗(yàn)程序的結(jié)果將給出（ ） （A）2010 年是丁亥年； （B）2010 年是庚寅年； （C）2

33、010 年是戊子年； （D）2010 年是丙戍年 （2）第五行語句的功能是（ ） （A）當(dāng)年份是 12 的倍數(shù)時(shí)定位為地支 12 （B）當(dāng)年份是 12 的倍數(shù)時(shí)定位為天干 12; （C）當(dāng)年份是 10 的倍數(shù)時(shí)定位為地支 10； （D）當(dāng)年份是 10 的倍數(shù)時(shí)定位為天干 10 4 非負(fù)函數(shù) y =f(x)在有限區(qū)間上的圖形為上半平面的一條曲線，曲線繞 x 軸旋轉(zhuǎn)時(shí)，產(chǎn)生以 x 為對稱軸 2 的旋轉(zhuǎn)曲面，其體積 V = f ( x)2 dx 0 syms a b x f=exp(a*x)*sin(b*x); f1=subs(f,a,-0.2); f2=subs(f1,b,0.5); V=pi*i

34、nt(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(0 0 0) axis off view(-17,54) （1） 關(guān)于程序的功能確切的說法（ ） (A)計(jì)算曲線段 f(x)=exp(a x)sin(b x), 0 x 2 繞 X 軸旋轉(zhuǎn)的旋轉(zhuǎn)曲面體積 (B)計(jì)算曲線段 f(x)=exp(-0.2x)sin(0.5x),

35、 0 x 2 繞 X 軸旋轉(zhuǎn)的旋轉(zhuǎn)曲面體積 (C) 計(jì)算曲線段 f(x)=exp(a x)sin(b x), 0 x 2 繞 X 軸旋轉(zhuǎn)的旋轉(zhuǎn)曲面體積并繪圖 （D）計(jì)算曲線段 f(x)=exp(-0.2x)sin(0.5x), 0 x 2 繞 X 軸旋轉(zhuǎn)的旋轉(zhuǎn)曲面體積 并繪圖 （2） 由第九行至第十一行語句可得旋轉(zhuǎn)曲面的方程（ ） 2 2 2 2 2 (A) f ( x ) = y + z + t （B） f ( x ) = y + z 2 2 2 2 2 （C） y = f ( x + z ) (D) y = f ( x + z + t ) 三、程序填空（32 分） 2 2 2 2 2 2 2 1維維安尼（Viviani）體是圓柱體( x R/2) + y R /4 被球面 x + y + z = R 所割下的立體。下 面的實(shí)驗(yàn)程序功能是取 R=2 求體積上半部分，先利用符號計(jì)算處理重積分并轉(zhuǎn)換為數(shù)值數(shù)據(jù)，再用蒙特 卡羅方法計(jì)算體積做對比。完成下面程序填空 第 3 頁 共 5頁 學(xué)院 姓名 學(xué)號 任課老師 選課號 密封線以內(nèi)答題無效 syms x y; f=sqrt(4-x2-y2);