計(jì)算機(jī)仿真復(fù)習(xí)題

.數(shù)字1.5e2,1.5e3中的哪個(gè)與1500相同？數(shù)字相同？［答案1.5e3答案］,.請指出如下5個(gè)變量名中，哪些是合法的？［答案xyz_3ABCDefgh答案］abcd-2xyz_33chana變量ABCDefgh.在Matlab中，比1大的最小數(shù)是多少？［答案］1+eps.指出指令clear,elf,clc各有什么作用？［答案clear:清除Matlab工作空間中的變量Clf:清除圖形窗口中的圖形Clc:清除當(dāng)前指令行窗口中的內(nèi)容.求矩陣行列式的指令是［答案］det.Matlab繪圖中，用于多子圖的指令是？［答案］subplot.指定繪圖使用菱形符藍(lán)色虛劃線的字符號(hào)是？［答案］'d-b'.視覺控制指令View(［az,el］)中，az表示？［答案方位角.用戶想計(jì)算y=2sin(0.3n)/1+5的值，那么應(yīng)依次輸入：2*sin(0.3*pi)/(l+sqrt(5))的值,.創(chuàng)建符號(hào)表達(dá)式f(x)=sinx+x［答案?f=sym('sin(x)+x')］.計(jì)算以上表達(dá)式在x=0,x=n/4,x=2Ji計(jì)算以上表達(dá)式在>>x(3)處的值。ans=0.2［答案答案］?subs(f,x,0)答案?x(［l479］)ans=0ans=00.30.60.8?subs(f,x,pi/4)?x(5:end)ans=1.4925ans=0.40.50.60.70.80.9?subs(f,x,2*pi)?x(find(x.0.4))ans=6.2832ans=0.50.60.70.80.91.032f.設(shè)x為符號(hào)變量，(x)=x4+x2+l,g(x)=x+4x+5x+8,試進(jìn)行如下運(yùn)算。①f(x)+g(x)②f(x)Xg(x)?f=sym('xA4+xA2+1');?g=sym('xA3+4*xA2+5*x+8');?件gans=xA4+5*xA2+9+xA3+5*x?f*gans=(xA4+xA2+l)*(xA3+4*xA2+5*x+8)兀.求cos(4x+3)dx和6cos(4x+3)dx求0>>symsx?f=cos(4*x+3);?int(f)ans=l/4*sin(4*x+3)?int(f,0,pi/6)ans=l/8*3A(l/2)*cos(3)-3/8*sin(3).已知x的閉區(qū)間［0,10*pi］內(nèi),y=sin(x),z=cos(x),寫出繪制如圖示的三維螺旋線的指令》x=0:pi/50:1O*pi?y=sin(x);z=cos(x)ffI?plot3(y,z,x) ?gridon?fitle('el;x')?xlabel('sin(t)')?ylabel('cos(t)')?zlabel()?text(0,0,0,'origin')xx<l16.有一個(gè)函數(shù)y=2xA2-l 1.5x<103x+2 x>=10編寫一個(gè)M文件輸入x,輸出y.functiony=y(x)ifx<lY=x;elseifx>=l&x<10Y=2*xA2-l;elseY=3*x=2end17.對于如下簡單系統(tǒng)進(jìn)行Simulink仿真。其中u(t)為系統(tǒng)輸入，(t)為系統(tǒng)輸出，t為系統(tǒng)時(shí)間。y(t)=2u(t),t>25lOu(t),tW25y很顯然，此系統(tǒng)服從簡單系統(tǒng)的條件，為一簡單系統(tǒng)。系統(tǒng)輸出僅由系統(tǒng)當(dāng)前時(shí)刻的輸入決定。一、填空題Matlab的命令提示符是（?）<,命令行以（回車）結(jié)束。MATLAB的一個(gè)命令行如有多條命令，命令間以（逗號(hào)）分隔。MATLAB的命令行以分號(hào)（；）結(jié)尾表示（不希望顯示命令運(yùn)算結(jié)果）。4、MATLAB的M文件由兩種類型分別是（命令文件或腳本文件）和（函數(shù)文件）。5、命令文件中的變量是（全局變量），在命令文件執(zhí)行后（保存）（填保存/不保存）于工作空間;函數(shù)文件中的變量一般是局部變量，調(diào)用執(zhí)行后（不保存）（填保存/不保存）于工作空間。6,函數(shù)文件保存時(shí)，文件名（必須是函數(shù)定義語句中的函數(shù)名）。（holdon）命令打開圖形保持功能。（%）符號(hào)表示其后為程序注釋。9,已知矩陣A,則矩陣A的轉(zhuǎn)置可表示為（AD,共碗轉(zhuǎn)置可表示為（AD10、Matlab命令窗口中的“>>”符號(hào)表示（）11、在MATLAB中，圖形的繪制必須在（圖形）或（Figure）窗口中進(jìn)行。二、選擇題1、說明下列文件的后綴A、命令文件；B、Simulink仿真文件；C、數(shù)據(jù)文件；D、函數(shù)文件（1）.mdl（2）.m（3）.dat（4）2、MATLAB是基于（函數(shù)）運(yùn)行的軟件，函數(shù)根據(jù)設(shè)計(jì)形式可以分為（內(nèi)建函數(shù)或預(yù)定義函數(shù)）和（用戶自定義函數(shù)）兩類；根據(jù)功能又可以分為（通用函數(shù)）和專用函數(shù)。3、MATLAB的賦值符是（1）,相等符號(hào)是（2）(1)= (2)=(3)=4、A(:J)表示(1)；A(i,:)表示(3)；A(iJ)表示(2)；A(i:i+m,:)表示(4).(1)取A矩陣的第j列全部元素；(2)取A矩陣第i行、第j列的元素。(3)A矩陣第i行的全部元素；(4)取A矩陣第i?i+m行的全部元素；5、說明下列MATLAB數(shù)據(jù)對象的分隔符A、字符串(2)；B、向量(3) C、函數(shù)的輸出變量(1)；D、矩陣(1)；E、矩陣中的行元素之間(4)；F、矩陣中的列之間(5)(1)方括號(hào)口 (2)單引號(hào)''(3)圓括號(hào)()(4)逗號(hào)，或空格(5)分號(hào)；6、命令subplot的用法和作用是(用來同時(shí)畫出數(shù)個(gè)小圖形于同一視窗中)，subplot(2,2,l)在一張圖中同時(shí)顯示(4)個(gè)子圖。7、說明下列特殊變量的含義或取值A(chǔ)、inf()；B、pi() ；C、NAN(不定項(xiàng))；D、ans()；E、：F、()。8、在MATLAB的命令窗口中執(zhí)行(C)命令，將命令窗口的顯示內(nèi)容清空。A.clearB.clcC.echooffD.cd9、下列變量名中A是合法的。A.char_l,i,jB.x*y,a.lC.x\y,al234D.end,lbcxB。10、已知x為一個(gè)向量，計(jì)算ln(x)的MATLAB命令是A.A.ln(x)B.log(x)C.Ln(x)D.lglO(x)11,已知a=0:4,b=l:5,下面說法中錯(cuò)誤的為DA.a+bB.a./bC.a'*bD.a*b。12、計(jì)算三個(gè)多項(xiàng)式sl、s2和s3的乘積，則算式為CA.conv(sI,s2,s3) B.s1*s2*s3 C.conv(conv(sI,s2),s3) D.conv(s1*s2*s3)13、運(yùn)行以下命令：?x=[l23;456];?y=x+x*i?plot(y)則在圖形窗口繪制C條曲線。A.3B.2 C.6 D.414、subplot(2,l,l)是指BA.兩行一列的上圖 B.兩行一列的下圖C.兩列一行的左圖D.兩列一行的右圖15、極坐標(biāo)圖是使用B來繪制的。A.原點(diǎn)和半徑 B.相角和距離 C.縱橫坐標(biāo)值 D.實(shí)部和虛部16、if結(jié)構(gòu)的開始是“if”命令，結(jié)束是BA.EndifB.endC.End一、單項(xiàng)選擇題.符號(hào)計(jì)算與一般數(shù)值計(jì)算有很大區(qū)別，它得到準(zhǔn)確的符號(hào)表達(dá)式。在MATLAB命令窗口中鍵入命令symsx,yl=sqrt(x);y2=xA2；int(yl-y2,x,0,l)?屏幕顯示的結(jié)果是(D)(A)yl=xA(l/2)(B)ans=2/3； (C)y2=xA2； (D)ans=1/3.在MATLAB命令窗口中鍵入命令A(yù)=[l42;312;615];det(A(l:2,2:3).*A(l:2,2:3))o結(jié)果是(B)(A)ans=-143(B)ans=60(C)ans=-16(D)ans=-19.設(shè)n階方陣A的特征值為：Ai(i=l,2,…，n),稱P(A)=max|入i|為矩陣Ai的譜半徑，則下列MATLAB求譜半徑命令是(A)(A)max(abs(eig(A)))；(B)abs(max(eig(A)))；(C)max(norm(eig(A)))； (D)norm(max(eig(A))).MATLAB系統(tǒng)運(yùn)行時(shí)，內(nèi)存中有包括X和Y在內(nèi)的多個(gè)變量(數(shù)據(jù))，要?jiǎng)h除所有變量(數(shù)據(jù))，應(yīng)該使用的命令是(A)(A)clear；(B)clc；(C)home；(D)clearXY.用賦值語句給定x數(shù)據(jù)，計(jì)算(7sin(3+2x)+e21n3嚴(yán)對應(yīng)的MATLAB表達(dá)式是(A)(A)sqrt(7*sin(3+2*x)+exp(2)*log(3))(B)sqrt(7sin(3+2x)+exp(2)log(3))(C)sqr(7*sin(3+2*x)+eA2*log(3))(D)sqr(7sin(3+2x)+eA2log(3)).在MATLAB命令窗口中輸入命令data=[4123131324];y=hist(data,4),hist是用來計(jì)算data數(shù)組的分布情況，也就是分別統(tǒng)計(jì)1,2,3.4每個(gè)數(shù)出現(xiàn)的次數(shù).其中，data中有3個(gè)1,2個(gè)33個(gè)3,2個(gè)4.所以輸出是3232.結(jié)果是(B)(A)y=4123； (B)y=3232； (C)y=1324； (D)y=4211.在MATLAB命令窗口中鍵入A=magic(6);B=A(2:5,1:2:5)將得到矩陣B,B是(C)(A)2行5列矩陣；(B)4行兩列矩陣；(C)4行3列矩陣；(D)4行5列矩陣.MATLAB繪三維曲面需要構(gòu)建網(wǎng)格數(shù)據(jù)，語句[x,y]=meshgrid(-2:2)返回?cái)?shù)據(jù)中：(D)(A)x是行向量，y是列向量；(B)x是列向量，y是行向量；(C)x是行元素相同的矩陣；(D)x是列元素相同矩陣9.下面有關(guān)MATLAB函數(shù)的說法，哪一個(gè)是錯(cuò)誤的(D)(A)函數(shù)文件的第一行必須由function開始，并有返回參數(shù)，函數(shù)名和輸入?yún)?shù)；MATLAB的函數(shù)可以有多個(gè)返回參數(shù)和多個(gè)輸入?yún)?shù)；(C)如果函數(shù)文件內(nèi)有多個(gè)函數(shù)，則只有第一個(gè)函數(shù)可以供外部調(diào)用；(D)在函數(shù)中可以用nargin檢測用戶調(diào)用函數(shù)時(shí)的輸出參數(shù)個(gè)數(shù)10.將帶小數(shù)的實(shí)數(shù)處理為整數(shù)稱為取整，常用四種取整法則是：向正無窮大方向取整、向負(fù)無窮大方向取整、向零方向取整和四舍五入取整。MATLAB提供了如下四個(gè)取整函數(shù)，若a=-1.4,對a取整的結(jié)果是-1,則不應(yīng)該選用下面哪個(gè)函數(shù)。(A)(A)floor；(B)round；(C)ceil；(D)fix一、單項(xiàng)選擇題I.在MATLAB命令窗口中，鍵入命令symsx；int(x*sin(x))。結(jié)果是(A)(A)ans=sin(x)-x*cos(x)；(B)ans=cos(x)+x*sin(x)；ans=sin(x)-cos(x)；(D)ans=-l/2*cos(x)*sin(x)+l/2*x.在MATLAB命令窗口中，鍵入命令symsx,F=1/(2+cos(x));ezplot(diff(F)),結(jié)果是(D)(A)繪出函數(shù)F在［0,2n］的圖形；(B)繪出函數(shù)F在［-2n,2n］的圖形；(C)繪函數(shù)F的導(dǎo)函數(shù)在［0,2n］的圖形；(D)繪函數(shù)F的導(dǎo)函數(shù)在［-2",2n］的圖形.在MATLAB命令窗口中鍵入命令B=［8,l,6;3,5,7;4,9,2］；B*B(:,2)。結(jié)果是(B)(A)ans=(B)ans=(C)ans=(D)ans=916767676791676767679167.MATLAB命令x=3:2:100將創(chuàng)建等差數(shù)列，該數(shù)列是((D))(A)以3為初值的98個(gè)數(shù),；(B)以100為終值的98的個(gè)數(shù)；(C)以99為終值的97個(gè)數(shù)；(D)以3為初值的49個(gè)數(shù)。.MATLAB語句strcat(int2str(2008)；年是s；年)的功能是(C)(A)將數(shù)據(jù)2008轉(zhuǎn)換為符號(hào)；(B)將數(shù)據(jù)2008與符號(hào)變量合并；(C)將幾個(gè)符號(hào)變量合并為一個(gè)；(D)將符號(hào)變量轉(zhuǎn)換為數(shù)值變量；.MATLAB的語句colormap(001)(D)(A)將三維網(wǎng)面圖確定為紅色；(B)將三維網(wǎng)面圖確定為綠色；(C)將三維網(wǎng)面圖確定為藍(lán)色；(D)語句使用格式錯(cuò)誤.設(shè)a,b,c表示三角形的三條邊，表達(dá)式a+b<c|a+c<b|b+c<a, (B)(A)是三條邊構(gòu)成三角形的條件：(B)是三條邊不構(gòu)成三角形的條件；（C）構(gòu)成三角形時(shí)邏輯值為真；（D）不構(gòu)成三角形時(shí)邏輯值為假三、簡答題：I、下面的變量名那些合法那些不合法.為什么？dogl 合法Idog 不合法，應(yīng)以字母開頭Do_you_know_the_way_to_sanJose合法d.help不合法，應(yīng)以字母開頭e.Whafsup?不合法，，和？不是合法字符。2、看下面的數(shù)組回答有關(guān)問題■1.10.02.1-3.56.0-0.01.1-6.62.83.4array\-2.10.10.3-0.41.3-1.45.10.0 1.10.0a.array1的大小是多少？5*4b.arrayl(4,l)的值是多少？-1.4-1.10.0-c.arrayl（:,l:2）的大小和值為多少?大小為4*2,值為《；；-1.45.13、假設(shè)a,b,c和d的定義如下分別指出下列表達(dá)式的運(yùn)算結(jié)果a+b(b)a.*c(c)a*b(d)a*c(e)a+c(f)a+d(g)a.*d(h)a*d答案：(a)這是一個(gè)數(shù)組或矩陣加法:10a=這是一個(gè)數(shù)組乘法： L21.「T0101(c)這是一個(gè)矩陣乘法：「一121f=-25_(d)這是一個(gè)矩陣乘法：°=38一」O(e)操作非法，兩數(shù)組形不同(f)數(shù)組與標(biāo)量的加法：a+d(g)數(shù)組乘法:50105(h)矩陣乘法:105-"Jo4、檢測下面的MATLAB語句，是對是錯(cuò)？正確的，輸出結(jié)果如何，錯(cuò)誤的，錯(cuò)在哪里?ifvolts>125disp('WARNING:Highvoltageonline.');ifvolts<105dispC'WARNING:Lowvoltageonline/);elsedisp(tLinevoltageiswithintolerances.*);end不正確，第二個(gè)if應(yīng)該用elseif語句當(dāng)電壓大于125,警告電壓過高；當(dāng)電壓小于105,警告電壓過低，否則（即電壓大于或等于105且小于或等于125時(shí)，顯示電壓在允許范圍內(nèi)）。.color='yellow1;switch(color);caseTed',disp('Stopnow!*);case'yellow',disp(Treparetostop.*);case'green,,dispCTroceedthroughintersection.*);otherwise,disp(,Illegalcolorencountered.*);end正確，這是十字路口交通信號(hào)燈的交通規(guī)則要求，根據(jù)燈的顏色顯示停車、等待停車或通行命令4、下面的語句用來判斷一個(gè)人的體溫是否處于危險(xiǎn)狀態(tài)(溫度用的是華氏計(jì)量)。這些語句是否正確？如果不正確，指出錯(cuò)在那里？應(yīng)該如何改正？iftemp<97.5dispCTemperaturebelownormal*);elseiftemp>97.5disp(*Temperaturenormal');elseiftemp>99.5disp(,Temperatureslightlyhigh*);elseiftemp>103.0dispCTemperaturedangerouslyhigh*);end正常體溫的數(shù)值范圍不正確，原因是elsfif的設(shè)置條件順序不正確，而且if條件不完備?？梢愿某梢韵滦问剑篿ftemp<97.5disp(*Temperaturebelownormal');elseiftemp>103.0disp(*Temperaturedangerouslyhigh*);elseiftemp>99.5disp(*Temperatureslightlyhigh*);elsedisp('Temperaturenormal*);end6、請你確定下列兩題中函數(shù)的調(diào)用是否正確。如果它是錯(cuò)誤的，指出錯(cuò)誤所在。不正確，test1函數(shù)有兩個(gè)輸入變量，調(diào)用時(shí)，必須給出兩個(gè)實(shí)參out=test1(6);functionres=testl(x,y)res=sqrt(x.A2+y.A2);函數(shù)調(diào)用正確(該例說明：MATLAB函數(shù)允許改變輸入輸出參數(shù)的個(gè)數(shù)。)out=test2(12);functionres=test2(x,y)%nargin函數(shù)返回調(diào)用這個(gè)函數(shù)時(shí)所需要的實(shí)際輸入?yún)?shù)的個(gè)數(shù)error(nargchk(l,2,nargin));%nargout函數(shù)：返回調(diào)用這個(gè)函數(shù)時(shí)所需要的實(shí)際輸出參數(shù)％的個(gè)數(shù)；nargchk函數(shù)：如果一個(gè)函數(shù)被調(diào)用時(shí)參數(shù)過多或過少，那么nargchk函數(shù)將返回一個(gè)標(biāo)準(zhǔn)錯(cuò)誤信息ifnargin=2res=sqrt(x.A2+y.人2);elseres=x; %允許輸入函數(shù)不等于2end7、簡要說明函數(shù)文件與命令文件的區(qū)別(1)函數(shù)文件的第一行必須包含function字符，一般以return結(jié)束，并且第一個(gè)注釋行專用于解釋函數(shù)功能；命令文件無此要求。(2)函數(shù)文件的第一行必須指定函數(shù)名、輸入?yún)?shù)及輸出參數(shù)，命令文件無此要求(即腳本文件沒有輸入?yún)?shù)和輸出參數(shù))。(3)一個(gè)函數(shù)文件可以含0個(gè)、1個(gè)或多個(gè)輸入?yún)?shù)和返回值。(4)函數(shù)文件要在文件的開頭定義函數(shù)名，如function[yl,y2]=fiinc(x,a,b,c),則該函數(shù)文件名必須存為func.m,而命令文件無此要求。(4)命令文件的變量在文件執(zhí)行結(jié)束以后仍然保存在內(nèi)存中而不會(huì)丟失(視為全局變量)，而函數(shù)文件的變量僅在函數(shù)運(yùn)行期間有效(除非用global把變量說明成全局變量，否則函數(shù)文件中的變量均為局部變量)，當(dāng)函數(shù)運(yùn)行完畢后，這些變量也就消失了。8、矩陣運(yùn)算和矩陣元素運(yùn)算有何不同？舉例說明。矩陣運(yùn)算滿足線性代數(shù)的矩陣運(yùn)算法則；矩陣元素運(yùn)算又稱為數(shù)組運(yùn)算，一般是兩個(gè)矢量矩陣之間的對應(yīng)元素進(jìn)行相應(yīng)運(yùn)算。9、ones和eye兩函數(shù)功能分別是什么？Ones::構(gòu)建全部矩陣元素均為1的矩陣。Eye：構(gòu)建單位矩陣，即除主對角線為1夕卜，其他元素均為0的矩陣；r兀2乃1[而而㈤9、用冒號(hào)表達(dá)式建立一個(gè)含有一百個(gè)元素的數(shù)組angles=Angles=(.01:.01:l)*pi10、已知A(:,:)=0.7094 0.6551 0.95970.7547 0.1626 0.34040.2760 0.1190 0.58530.6797 0.4984 0.2238,求下列命令的執(zhí)行結(jié)果?size(A) %獲取數(shù)組尺寸ans=43?numel(A)%獲取數(shù)組元素的個(gè)數(shù)ans=1211、MATLAB中的向量和矩陣有什么不同(1)向量是一維數(shù)組，包括行向量和列向量；只有一個(gè)元素的向量就是標(biāo)量；(2)矩陣包括線性代數(shù)中的二維矩陣和多維矩陣。12、有幾種建立矩陣的方法？各有什么優(yōu)點(diǎn)？？1、直接輸入法，例：a=[379;49101],輸入方法方便簡捷2、通過M文件建立矩陣，適用于建立尺寸較大的矩陣，并且易于修改3、由函數(shù)建立，如y=sin(x),可以用函數(shù)內(nèi)部建立一些特殊的矩陣4、通過數(shù)據(jù)文件建立，該方法可以調(diào)用其他軟件產(chǎn)生數(shù)據(jù)13、在sin(x)運(yùn)算中，x是角度還是弧度？在sin(x)運(yùn)算中，x是弧度，MATLAB規(guī)定所有的三角函數(shù)運(yùn)算都是按弧度進(jìn)行運(yùn)算。14、四、程序分析題寫出以下命令的執(zhí)行結(jié)果函數(shù)印per和lower分別把一個(gè)字符串中所有轉(zhuǎn)化大定和小寫。?result=upper(*Thisistest1!')result=THISISTEST1!?result=lower(*Thisistest2V)result=thisistest2!寫出以下命令的執(zhí)行結(jié)果?a=[l23;456]a=2 34 5 6?b=mat2str(a)b=[123;456]3、下面的語句產(chǎn)生的result的值是多少？(a)x=12+i*5;y=5-i*13;result=x>y;result=1,因?yàn)楸容^的是數(shù)的實(shí)部。(b)x=12+i*5;y=5-i*13;result=abs(x)>abs(y);result=0,因?yàn)閮蓚€(gè)數(shù)的絕對值是確定的。(c)x=12+i*5;y=5-i*13;result=real(x)-imag(y);result=254.說明下列程序段的功能：根據(jù)輸入的數(shù)值分?jǐn)?shù)，輸出等級(jí)分?jǐn)?shù)A級(jí)、B級(jí)和C級(jí)，滿足grade>95A95grade>86B86grade>76C76grade>66D66grade>0Fifgrade>95.0disp(*ThegradeisA.');elseifgrade>86.0disp(*ThegradeisB.');elseifgrade>76.0disp(*ThegradeisCJ);elseifgrade>66.0disp(*ThegradeisD/);elsedispC'ThegradeisF.');endifgrade>95.0disp(*ThegradeisA.*);elseifgrade>86.0disp(fThegradeisB.1);elseifgrade>76.0disp(*ThegradeisCJ);elseifgrade>66.0disp(*ThegradeisDJ);elsedisp('ThegradeisF/);endendendend5、說明下列程序段功能：根據(jù)輸入自變量x的值,計(jì)算分段函數(shù)y的值,滿足：丫4伙，x<。lg(^+vl+x2)/2,x>0x=input(，請輸入x的值ifx<=0y=(x+sqrt(pi))/exp(2);elsey=log(x+sqrt(l+x*x))/2;end6、說明下列程序段功能分析下列程序繪制的曲線。x1=linspace(0,2*pi,l00);x2=linspace(0,3*pi,100);x3=linspace(0,4*pi,100);yl=sin(xl);y2=l+sin(x2);y3=2+sin(x3);x=[xl;x2;x3]1;y=[yi;y2;y3]‘；plot(x,y)在一個(gè)圖形窗口繪制出三條二維曲線，y1=sinx,y2=1+sinx,y3=2+sinx,橫軸x取值范圍都是0~100,步長分別為2冗、3n和4n。分析下列程序繪制的曲線。t=0:pi/100:20*pi;x=sin(t);y=cos(t);z=t.*sin(t).*cos(t);plot3(x,y,z);三維曲線title(*Linein3-DSpace*);xlabel(,X,);ylabel(,Y,);zlabel(,Z,);gridon;在X、Y、Z坐標(biāo)系繪制帶柵格的三維曲線，標(biāo)題為Linein3?DSpace。分析下列程序的功能。[x,y]=meshgrid(0:0.25:4*pi);z=sin(x+sin(y))-x/l0;mesh(x,y,z);三維網(wǎng)格圖axis([04*pi04*pi-2.51]);繪制三維曲面圖z=sin(x+sin(y))-x/10o7、創(chuàng)建符號(hào)變量有幾種方法？MATLAB提供了兩種創(chuàng)建符號(hào)變量和表達(dá)式的函數(shù)：sym和symsosym用于創(chuàng)建一個(gè)符號(hào)變量或表達(dá)式，用法如x=sym(，x，)及f=sym(tx+y+z,),syms用于創(chuàng)建多個(gè)符號(hào)變量，用法如symsxyz。f=sym('x+y+z')相當(dāng)于symsxyzf=x+y+z五、編程題I、骰子模擬。模擬擲骰子的情況在現(xiàn)實(shí)中非常有用。編寫一個(gè)maltab程序模擬擲骰子，每次產(chǎn)生一個(gè)1到6之間的隨機(jī)整數(shù)。(可以編寫命令或函數(shù)文件)。命令文件dice.m函數(shù)文件fdice.m2、某商場對顧客所購買的商品實(shí)行打折銷售，標(biāo)準(zhǔn)如下(商品價(jià)格用price來表示):price<200 沒有折扣200^price<500 3%折扣500^price<1000 5%折扣1000^price<2500 8%折扣2500^price<5000 10%折扣5000^price 14%折扣編程輸入所售商品的價(jià)格，并求其實(shí)際銷售價(jià)格。程序如下：price=input(,W輸入商品價(jià)格')；switchfix(price/100)case{0,1} %價(jià)格小于200rate=0;case{2,3,4}%價(jià)格大于等于200但小于500rate=3/100;casenum2cell(5:9)%價(jià)格大于等于500但%小于1000rate=5/100;casenum2ceH(10:24)%價(jià)格大于等于1000%但小于2500rate=8/100;casenum2cell(25:49)%價(jià)格大于等于2500%但小于5000rate=10/100;otherwise %價(jià)格大于等于5000rate=14/100;endprice=price*(l-rate)%輸出商品實(shí)際銷售價(jià)格3、符號(hào)運(yùn)算題(每題5分，共15分)Calculatethethird-orderderivative(三階導(dǎo)數(shù))offollowingfunctionbysymbolicoperation:sinx+e2xP='sinx+exp(2*x);'DifT(p,'x',3);ans=8*exp(2*x)或p=sym(<sinx4-exp(2*x)>);Solvethefollowingsystemofequationsbysymbolicoperation:x+2xy-y=32x2x-3x+5=0Pl=,xA2+2*x*y-y=3,;P2=,2*xA2-3*x+5=0,;[x,y]=solve(Pl,P2)x=3/4+l/4*i*31A(l/2)3/4-l/4*i*31A(l/2)y=-29/64-19/64*i*31A(1/2)-29/64+19/64*i*31A(1/2)4、編寫MATLAB腳本文件，計(jì)算Fibonacci數(shù)列中的第一個(gè)大于10000的元素及n值。Fibonacci數(shù)列從［1,1］開始，由數(shù)列的最后兩個(gè)元素之和生成新的元素，依次遞推。clear,closeallN=inpu(輸入最大數(shù)值范圍N=，)％等待用戶輸入數(shù)據(jù)HU］；i=l； %變量的初始化whilef⑴+f(i+l)<N%循環(huán)條件檢驗(yàn)Ri+2尸f(i+l)+f(i);i=i+l;%求fibonnaci數(shù)的算式end5、編寫MATLAB程序，以0.2為步長繪制正弦曲線sin(nt)和余弦曲線cos(nt)在［0,2冗］之間的圖形，分別位于圖形窗口的二個(gè)子圖中，正弦曲線以紅色虛線繪制，并且用星號(hào)標(biāo)注每個(gè)節(jié)點(diǎn)，余弦曲線以紅色實(shí)線繪制，用圓標(biāo)注每個(gè)節(jié)點(diǎn)。6、在0WxW2兀區(qū)間內(nèi)，繪制曲線y=2e-0.5xcos(4nx)程序如下：x=0:pi/100:2*pi;y=2*exp(-0.5*x).*cos(4*pi*x);plot(x,y)7、采用圖形保持，在同一坐標(biāo)內(nèi)繪制曲線yl=0.2e?0?5Xcos(4冗x)和y2=2e?°?5Xco§(冗x)。程序如下：x=0:pi/100:2*pi;yl=0.2*exp(-0.5*x).*cos(4*pi*x);plot(x,yl)holdony2=2*exp(-0.5*x).*cos(pi*x);plot(x,y2);holdoff8、在同一坐標(biāo)內(nèi)，分別用不同線型和顏色繪制曲線yl=0.2e-0.5xcos(4冗x)和y2=2e-0.5xcos(nx),標(biāo)記兩曲線交叉點(diǎn)。程序如下：x=linspace(0,2*pi,1000);yl=0.2*exp(-0.5*x).*cos(4*pi*x);y2=2*exp(-0.5*x).*cos(pi*x);k=find(abs(y1-y2)<1e-2); %查找yl與y2相等點(diǎn)(近似相等)的下標(biāo)xl=x(k); %取yl與y2相等點(diǎn)的x坐標(biāo)y3=0.2*exp(-0.5*xl).*cos(4*pi*xl);%求yl與y2值相等點(diǎn)的y坐標(biāo)plot(x,y1,x,y2；k:\xl,y3；bp*);9、在QWxW2兀區(qū)間內(nèi)，繪制曲線yl=2e05x和y2=cos(4nx),并給圖形添加以下標(biāo)注。(1)標(biāo)題為：xfrom0to2Jr；(2)X軸標(biāo)簽為：VariableX,Y軸標(biāo)簽為：VariableY：(3)在坐標(biāo)(0.8,1,5)處添加標(biāo)注曲線yl=2eq5x,在坐標(biāo)(2.5,1.1)處添加標(biāo)注曲線y2Kos(4nx)；(4)添加圖例Yl、Y2.程序如下：x=0:pi/100:2*pi;yl=2*exp(-0.5*x);y2=cos(4*pi*x);plot(x,yl,x,y2)%加圖形標(biāo)題%加%加圖形標(biāo)題%加X軸說明%加Y軸說明%在指定位置添加圖形說明%加圖例xlabel('VariableX');ylabelfVariableY*);text(0.8,l.5；曲線y1=2eA{-0.5x}t);text(2.5,L1；曲線y2=cos(4!\pi}x),);legend(fcyr；y2')10、求下列聯(lián)立方程的解3x+4y-7z—12w=45x-7y+4z+2vv=4x+8-5卬=9-6x+5y-2z+lOvv=4?a=[34-7-12;5-742;180-55-65-210];?y=[4494];?x=inv(a)*y'2.94471.5512-0.60171.270911、（共15分,每小題5分）利用函數(shù)int計(jì)算:f2%—dxJ(x2+l)(x2+x)?int(1/((xA2+1)*(xA2+x)))ans=-1/4*log(xA2+1)-1/2*atan(x)+log(x)-1/2*log(x+1)(2)Jj(x+y)e~xydxdy?symsxy?int(int((x+y)*exp(-x*y)))ans=-l/y*(-l/y*(-l/y*(-exp(-x*y)*x*y-exp(-x*y))+exp(-x*y)/y)-exp(-x*y))?symsxy?f=(x+y)*exp(-x*y);?F=int(int(f,x),y)F=x*Ei(l,x*y)+x*(l/x/y*exp(-x*y)-Ei(l,x*y))+l/x*exp(-x*y)?simple(F)simplify:exp(-x*y)*(x+y)/x/y(3)2=向。+3＞求?symsxy?y=sin(x)A4+cos(x)A4;?subs(y,x,10)ans=0.583312、(共10分)編制m文件，等待鍵盤輸入，輸入密碼555,密碼正確，顯示輸入密碼正確，程序結(jié)束；否則提示，重新輸入。a=input('inputpasswordplease')；while(a?=555)a=input('wrongpassword,inputagain：')；end13、（共14分）編寫m函數(shù)文件求底面半徑為r、高為h的圓柱體的體積。functionttr=input（'輸入底面半徑'）;h=input（'輸入高度'）；v=pi*r*r*h;disp（v）;14、（共14分）從鍵盤輸入若干個(gè)數(shù)，當(dāng)輸入0時(shí)結(jié)束輸入，求這些數(shù)的平均值和立方和。sum=0;sum3=0;n=0;a=input（，請輸入數(shù)，）;while（a?=0）sum=sum+a;sum3=sum34-aA3；n=n+l;a=input（'輸入數(shù)：'）；endifn?=0aver=sum/n;enddisp('平均數(shù)為:');disp(aver);dispC立方和為:9);disp(sum3);15、設(shè)x=rcost+3t,y=rsint+3,分別令r=2,3,4,畫出參數(shù)t=0~10區(qū)間生成的x~y曲線。>t=linspace(0,10);>rl=2;>xl=(rl*cos(t)+3*t);>yl=rl*sin(t)+3;?r2=3;>x2=(r2*cos(t)+3*t);>y2=r2*sin(t)+3;?r3=4;>x3=(r3*cos(t)+3*t);>y3=r3*sin(t)+3;?plot(xl,yl,'r,x2,y2,'b',x3,y3,'m')16、設(shè)f(x)=x"-4x'+3x2-2x+6(1)在x=-2,8］之間取100個(gè)點(diǎn)，畫出曲線，看它有幾個(gè)過零點(diǎn)。(提示:用polyval函數(shù))>>x=linspace(2,8,100);>>y=polyval([10-43-26],x);?plot(x,y,'b',x,0,'y')(2)用roots函數(shù)求此多項(xiàng)式的根。t=[l0-43-26]p=roots(t)17、繪制在定義域x=[-2,2],y=[-2,2]內(nèi)的曲面。(利用meshgrid)x=-2:2;y=x;[X,Y]=meshgrid(x,y);Z=X"2*exp(-(X"2+Y"2));mesh(X,Y,Z)4.求代數(shù)方程3x5+4x"+7x3+2x2+9x+12=0的所有根。(利用roots函數(shù))p=[3472912];roots(p)18、把1開五次方，并求其全部五個(gè)根。(提示：解x‘-1=0)p=[l0000-1];roots(p)(-x2-4x-3)/2 -34xv-l19>y=-x2+l-1<x<1 ,回出在[-3,3]上的曲線。(-x2+4x-3)/21^x<3xl=-3:0.01:-l;yl=(-xl.A2-4.*xl-3)/2;x2=-l:0.01:l;y2=-x2.A2+l;x3=l:0.01:3;y3=(-x3.A2+4.*x3-3)/2;plot(xl,yl);holdon;plot(x2,y2);holdon;plot(x3,y3)20、設(shè)”?！?，求尸x=symCx');y=x*exp(2*x);diff(y,20)21>設(shè)+求小。x=sym(,x');y=(sin(x))"4+(cos(x))"4;diff(y,10)22、設(shè)?=八也&十w石§,求無x=sym(,x');y=sym(,y*);z=y*4*((sin(x))4)+x4*((cos(y))4);int(z,x)23、f—一dxJ(1+1)(/+x)x=sym(，x');int(1/(x2+l)*(x*2+x),x)24、慳聲x=sym(,x);int((sqrt((9-x*2)*3)/x6),x)25、,今x=§ym('x');iiit(x+xA3,x,l,2)26、求級(jí)數(shù)的和：n=l2sym(4n,);symsum(n+2/2An,nJJnf)27、利用函數(shù)int計(jì)算一隹不施積舒’x=sym('x');y=sym('y');z=(x+y)*exp(-x*y);a=int(z,x);int(a,y)28、編制m文件，等待鍵盤輸入，輸入密碼123,密碼正確，顯示輸入密碼正確，程序結(jié)束；否則提示，重新輸入。方法一：functionf=checkpassword=input(*Pleaseinputthepassword:\nf);if(password=123)dispC輸入密碼正確，)elsecheck;end方法二：functionpwn=0;k=0;n=inputC請輸入密碼(整數(shù))\n);while(k=0)ifn=123dispC輸入正確，)k=l;elsek=0;n=input。密碼錯(cuò)誤，請重新輸入密碼(整數(shù))\n);endend29、編制m文件，輸入n(正整數(shù))，顯示所有小于n的質(zhì)數(shù)。%第一m文件(用于判斷某數(shù)是否為素?cái)?shù))：functionleap=panduan(n)leap=l;fori=2:sqrt(n+l)ifmod(nj)=0leap=0;break;endend%第二個(gè)m文件(用于讀進(jìn)n個(gè)數(shù)并調(diào)用panduan.m判斷讀數(shù)是否為素?cái)?shù)，并把是素?cái)?shù)的顯示出來)：functionzhishu(n)form=2:nifpanduan(m)=ldisp(m);endend30、編寫一個(gè)M函數(shù)文件fun_es(x),計(jì)算如下函數(shù)：＞=0.5/"-3皿)，其中參數(shù)可以為標(biāo)量，也可以為向量。在Matlab指令窗口鍵入如下命令檢驗(yàn)該函數(shù):fun_es(3)functionfunesy=0.5.*exp(x./3)-x.2.*sin(x)31、下面程序?qū)崿F(xiàn)什么功能？forn=100:200ifrem(n,77)~=0continueendbreakendn求100到200之間，第一個(gè)77的整倍數(shù)32、編寫m函數(shù)文件求半徑為r的圓的面積和周長。functionf=circle(r)S=pi*J2L=2*pi*r33、數(shù)組［1：99］中素?cái)?shù)(即質(zhì)數(shù))元素的和%第一m文件(用于判斷某數(shù)是否為素?cái)?shù))：functionleap=panduan(n)leap=l;fbri=2:sqrt(n+l)ifmod(n,i)=0leap=0;break;endend%第二個(gè)m文件(用于讀進(jìn)1—99并調(diào)用panduan.m判斷讀數(shù)是否為素?cái)?shù)，并把是素?cái)?shù)的加起來,最后輸出)：functiontotal(n)k=0;form=l:nifpanduan(m)=lk=k+m;endendk34、從鍵盤輸入若干個(gè)數(shù)，當(dāng)輸入0時(shí)結(jié)束輸入，求這些數(shù)的平均值以及平方和。%求n個(gè)數(shù)的平均值及平方和functionpingfang11=0;total1=0;total2=0;a=input('inputa:\n')；while(a?=0)1=1+1;totall=totall+a;tota!2=total2+aA2;a=input(M);endavg=totall/ltotal235、輸入x,y的值，并將它們的值互換后輸出。functionf=change(x,y)x=input(*Inputx:\rf);y=input(*Inputy:\n')；a=x;x=y;y=a;x,y36、利用rand產(chǎn)生10個(gè)隨機(jī)數(shù)，利用for循環(huán)對其進(jìn)行排序(從大到小)。disp(,排列前:')a=rand(10,l)fori=l:10forj=l:10-iifa(j+l,l)>ad，l)b=aG，l);a(j,l)=a(j+l,l);a(j+l,l)=b;endendenddisp。排列后:')37、角度仙[304560],求工的正弦、余弦、正切和余切。?x=[304560];?xl=x/180*pi;?sin(xl)ans=0.5000 0.7071 0.8660?cos(xl)ans=0.8660 0.7071 0.5000?tan(x1)ans=0.5774 1.0000 1.7321?cot(xl)ans=1.7321 1.0000 0.577438、用四舍五入的方法將數(shù)組[2.45686.39823.93758.5042]取整。?b=[2.45686.39823.93758.5042];?round(b)ans=39、基本編程題(每小題10分，共30分)1.Writeaprogramtocomputethefollowingexpressions10Z(4-可⑸-虧)%=1u ,， …,10；10 10 2-利￡(%-%)V*=iHere,wesupposethevariablexhasexistedintheworkspace.fori=l:10forj=l:10xbar(i)=xbar(i)+x(i,j);endxbar(i)=xbar(i)/10;endfori=l:10forj=l:10tl=0;t2=0;t3=0;fork=l:3tl=tl+(x(i,k)-xbar(i))*(x(j,k)-xbar(j));t2=t2+