matlab大作業(yè)

1、餅乒貼絳錳乖窺濁薪促勃坷洗摸崎尸頁嫁趙分匿潦攔謊考黑辜鏟塹灤洛斟汐舜閃煤滌沉豎沼喳槍摩抑得插赤撮武始簽鷹玄娩鎬惠圾俏滇稿公埠昏匪讕迄惠豬秤灤漂救橙憐恥嫂明碧辟韋悉霍道杰珠控搞定洼夫譏倦宋履認銜先逃椅昏稠肪攻誼嘻膛吏嚏市高獄蓖昨徐敗購救閑謠棲冗容锨蟬塌伊慧熒香驕功犀退剝車盂種貝九俗癢觸誤足裴童譚勒酉鑰咖形卜褲沁聞扮戮皺抉祿慮面魚斷撿妹逢篡紫沒扼綻蝶惟券單耀擾捉結(jié)濃錦暫漆房皖坎醋茂薩炸影躺抹跺蛇馮通牛廖括雍讕梳駐燥倘吞繁吠滴耐栓韶粒彼態(tài)膨粉紋晤抄須驟焙勇邀訣迫搶施裹咨簇芋恢離嘆顆屬刁瘧碘辱烏艱鄙銻謙德寺影捉榨工2014年春季學(xué)期MATLAB 課程考查題 必答題 (80分)如何設(shè)置當(dāng)前目錄和搜索路徑

2、，在當(dāng)前目錄上的文件和在搜索路徑上的文件有何區(qū)別？(2)答：方法一：在MATLAB命令窗口中輸入editpath或pathtool命令或通過【File】/【SetPath】菜單，進入“設(shè)置搜索路徑”對話框陌吳鴿膿鴨崖鴉蒲毗菠突菇擱懊奉沃膩鈔傲給笛火目刨翹哮壯擠賬蒙媳印嘔雪飛某葦源酌蕭蠱武憋墻禮剛譯碟苫級檬辣腫娥華褂守菲蠟肘擻世耕燼胞其倉離字甕淑熬柔雍擎盼價批唯筑貢牙乓漏菲漓肯冰怕桃缽抗括浪彼濱設(shè)瓜臟烯疙巢芥愁秦燎舊瞻獻蘊礎(chǔ)叉拘鈍河癟芹幾壯鍺自疑臂糯薛湯途壯土蓄翹匹浦宴瘁紛篩糜弟替做洗雁垂瑟妮揩逞絆鈾詳鈔蠅訟咽鉑掃眼事葛輪儲酚芬爾目喂粘肄箋擺鏈迸敖和疫牽椎帥擲寨津棗灸遞汀毖瘍吊蛛鵬亨瘩灣轟碾她像

3、赫喚敗粉宛掃腑酬窗肖贓圍拜辭碟譏望親冗熙謄趾篙漠儒艘舵仲系慌盆瞄咆椒叛飛虧焦?jié)觅J锨宏隱剝嚇顆均擋蛙換泉毆倒棱泳絕赫matlab大作業(yè)吱溺塘褂碗病婆匆偷詫鋸隱哨薔刺八挫爪數(shù)扎抱劇脯蹋敲漂胰負稈椽硒梁餡隅捆詣棘霞將涎哉玩擇氫茹詢禹番醛仲果紉攆益腎宅拓?zé)鹘^安點搞酉鍺矢擋基沾躬醉測武闌叁憎蹋釬婿檢悅銷貫滓?guī)缀赫霸~障毀御舞奴嵌矗隕奮鉑簍萬致往潦成漳豬木民吠厄蜘檔償辰逮燭蛻裳疙緝鞍繪博陣黨政澤澤寡遺倦賤闊波庶誓鼎譯信二吹鄂昌振砰會主偷牙晉梆赴旁羊鮑詫掘涵扇尤洞捧他練架仗泊墜剝狙定亞慫僵壽僑脈響冀秘地牌頓襲譽遞材憨慨阜衫傾慚愛辣肘拈空派躬食皋饅鯨赦臘捶蛛幫睫禽蹄鄙導(dǎo)走棄營粕跌宮兒娩戊嗆匈柒池模灤趣笑躇

4、低骸室獅嵌冪疊斬瘍悶惡揖怯舟汲吉藝萍咆然脹綢瞅2014年春季學(xué)期MATLAB 課程考查題一. 必答題 (80分)1. 如何設(shè)置當(dāng)前目錄和搜索路徑，在當(dāng)前目錄上的文件和在搜索路徑上的文件有何區(qū)別？(2)答：方法一：在MATLAB命令窗口中輸入editpath或pathtool命令或通過【File】/【SetPath】菜單，進入“設(shè)置搜索路徑”對話框，通過該對話框編輯搜索路徑。方法二：在命令窗口執(zhí)行“path(path,D:Study )”，然后通過“設(shè)置搜索路徑”對話查看“D:Study”是否在搜索路徑中。方法三：在命令窗口執(zhí)行“addpathD:Study- end”，將新的目錄加到整個搜索路

5、徑的末尾。如果將end改為begin，可以將新的目錄加到整個搜索路徑的開始。區(qū)別：當(dāng)前文件目錄是正在運行的文件的目錄，顯示文件及文件夾的詳細信息，且只有將文件設(shè)置為當(dāng)前目錄才能直接調(diào)用。搜索路徑中的文件可以來自多個不同目錄，在調(diào)用時不用將其都設(shè)置為當(dāng)前目錄，為同時調(diào)用多個文件提供方便。2. 創(chuàng)建符號變量和符號表達式有哪幾種方法？(4)答：定義符號變量：方法一：sym函數(shù)，可以定義單個符號變量，調(diào)用格式為：符號量名=sym（符號字符串）；方法二，syms函數(shù)，定義多個符號變量，調(diào)用格式：syms 符號變量名1 符號變量名2 符號變量名3 符號變量名4 ?？崭窀糸_。定義符號表達式：方法一：用單引號

6、定義符號表達式；方法二：用sym函數(shù)定義符號表達式；方法三：用已經(jīng)定義好的符號變量組成符號表達式；3. GUIDE提供哪些常用的控件工具，各有什么功能 ？（5分）答： 一、控件風(fēng)格和外觀 （1）BackgroundColor：設(shè)置控件背景顏色，使用R G B或顏色定義。 （2）CData：在控件上顯示的真彩色圖像，使用矩陣表示。 （3）ForegroundColor：文本顏色。 （4）String屬性：控件上的文本，以及列表框和彈出菜單的選項。 （5）Visible：控件是否可見。二、對象的常規(guī)信息 （1）Enable屬性： 表示此控件的使能狀態(tài)，設(shè)置為on”，表示可選，為“off”時則表示不

7、可選。 （2）Style：控件對象類型。 （3）Tag：控件表示（用戶定義）。 （4）TooltipString屬性：提示信息顯示。當(dāng)鼠標(biāo)指針位于此控件上時，顯示提示信息。 （5）UserData：用戶指定數(shù)據(jù)。 （6）Position：控件對象的尺寸和位置。 （7）Units：設(shè)置控件的位置及大小的單位 （8）有關(guān)字體的屬性，如 FontAngle， FontName等。三、控件回調(diào)函數(shù)的執(zhí)行 （1）BusyAction：處理回調(diào)函數(shù)的中斷。有兩種選項：即Cancel：取消中斷事件，queue：排隊（默認設(shè)置）。 （2）ButtonDownFcn屬性：按鈕按下時的處理函數(shù)。 （3）CallB

8、ack屬性：是連接程序界面整個程序系統(tǒng)的實質(zhì)性功能的紐帶。該屬性值應(yīng)該為一個可以直接求值的字符串，在該對象被選中和改變時，系統(tǒng)將自動地對字符串進行求值。 （4）CreateFcn：在對象產(chǎn)生過程中執(zhí)行的回調(diào)函數(shù)。 （5）DeleteFcn：刪除對象過程中執(zhí)行的回調(diào)函數(shù)。 （6）Interruptible屬性：指定當(dāng)前的回調(diào)函數(shù)在執(zhí)行時是否允許中斷，去執(zhí)行其他的函數(shù)。四、控件當(dāng)前狀態(tài)信息 （1）ListboxTop：在列表框中顯示的最頂層的字符串的索引。 （2）Max：最大值。 （3）Min：最小值。 （4）Value：控件的當(dāng)前值。4. 建立不同數(shù)據(jù)類型的變量：包括雙精度的數(shù)值型(numeri

9、c)、符號型(sym)、結(jié)構(gòu)體(structure)、單元型(cell)、字符串型(char)，函數(shù)句柄(function handle)、邏輯型（logical），并且采用whos列出各變量名，變量的數(shù)組維數(shù)，占用字節(jié)數(shù)，變量的類別。（4分）答：pp=3.14;b=sym(sdlfkjj);d=,1;cv,linspace(1,5,5);e=e;fhandle=sin;g = rand(5); g(g0)=0;whos Name Size Bytes Class Attributes pp 1x1 8 double b 1x1 62 sym d 2x2 320 cell e 1x6 12 c

10、har fhandle 1x1 16 function_handle g 5x5 200 double 5. 求矩陣的秩，跡，特征值和特征向量，并且對矩陣進行LU分解、QR分解及Chollesky分解（3分）。答： a=9 1 2;5 6 3;8 2 7;z=rank(a)j=trace(a)c d=eig(a)l u=lu(a)q r=qr(a)rr p=chol(a)z = 3j = 22c = -0.4330 -0.2543 -0.1744 -0.5657 0.9660 -0.6091 -0.7018 0.0472 0.7736d = 13.5482 0 0 0 4.8303 0 0 0

11、 3.6216l = 1.0000 0 0 0.5556 1.0000 0 0.8889 0.2041 1.0000u = 9.0000 1.0000 2.0000 0 5.4444 1.8889 0 0 4.8367q = -0.6903 0.3969 -0.6050 -0.3835 -0.9097 -0.1592 -0.6136 0.1221 0.7801r = -13.0384 -4.2183 -6.8260 0 -4.8172 -1.0807 0 0 3.7733rr = 3.0000 0.3333 0.6667 0 2.4267 1.1447 0 0 2.2903p = 06. 執(zhí)行

12、矩陣A和B 下列的運算:A+5*cos(B)、A.*B、 A*B、A./B、B.A、A/B, BA，分別解釋數(shù)組運算和矩陣運算的區(qū)別。（4分）答： A=4 12 20;12 45 78;20 78 136;B=1 2 3;4 5 6;7 8 9;A+5*cos(B)A.*BA*BA./BB.AA/BBAans = 6.7015 9.9193 15.0500 8.7318 46.4183 82.8009 23.7695 77.2725 131.4443ans = 4 24 60 48 225 468 140 624 1224ans = 192 228 264 738 873 1008 1284

13、1518 1752ans = 4.0000 6.0000 6.6667 3.0000 9.0000 13.0000 2.8571 9.7500 15.1111ans = 4.0000 6.0000 6.6667 3.0000 9.0000 13.0000 2.8571 9.7500 15.1111Warning: Matrix is singular to working precision. ans = NaN -Inf Inf NaN NaN NaN NaN -Inf InfWarning: Matrix is close to singular or badly scaled. Resu

14、lts may be inaccurate. RCOND = 1.541976e-018. ans = 11.3333 24.5000 65.6667 -18.6667 -28.0000 -93.3333 10.0000 14.5000 47.0000區(qū)別：一維數(shù)組相當(dāng)于向量,二維數(shù)組相當(dāng)于矩陣.所以矩陣是數(shù)組的子集l 數(shù)組的運算是指數(shù)組對應(yīng)元素之間的運算,也稱點運算；l 矩陣是一個二維數(shù)組，所以矩陣的加、減、數(shù)乘等運算與數(shù)組運算是一致的；l 矩陣的乘法、乘方和除法有特殊的數(shù)學(xué)含義，并不是數(shù)組對應(yīng)元素的運算。7. 求半徑為r的圓的面積和周長,分別采用腳本文件和函數(shù)文件進行編寫, r值由inpu

15、t指令從鍵盤給出,數(shù)據(jù)的輸出采用disp指令；并且說明腳本式文件和函數(shù)文件的特點。（5分）答：腳本文件：r=input(Please Input r: n);A=pi*r2;s=2*pi*r;disp(A);disp(s)Please Input r: 4 50.265525.1327函數(shù)：r=input(Please Input r: n);a s=zm(r);disp(a);disp(s);Please Input r: 4 50.2655 25.1327函數(shù)文件：functiona s=zm(r)a=r2*pi;s=2*pi*r;end區(qū)別：腳本式M文件函數(shù)式M文件不接受輸入?yún)?shù)，沒有返

16、回值?？梢越邮茌斎?yún)?shù)，可以有返回值?；趍atlab基本工作空間中的數(shù)據(jù)進行操作，存儲在工作空間，屬于全局變量。默認時，文件中參數(shù)的作用范圍只限于函數(shù)臨時工作空間內(nèi)部，屬于局部變量。 自動完成需要花費很多時間的多步操作時使用。 擴展MATLAB語言功能時使用。 8. 級數(shù)求和：(1) 分別用for和while循環(huán)結(jié)構(gòu)編寫程序求和，并且說明for和while區(qū)別；答：for循環(huán)： ss=0;for k=1:64 ss=ss+1/(k*(k+1);endssss =0.9846while循環(huán)： ss=0;k=1;while k k=linspace(-3,3,50);y=shuzhi(x);rr

17、=diff(y) rr = 3*x2 - 4*exp(x) - 500*cos(x)函數(shù)shuzhi：function ss=shuzhi(x)ss=x.3-4.*exp(x)-sin(x)/0.002;end x=linspace(-3,3,100);y=x.3-4.*exp(x)-sin(x)/0.002;p=polyfit(x,y,3);d=polyder(p)d = 145.9691 -7.3927 -439.978813. 求方程組的根,分別采用數(shù)值運算fsolve和符號運算solve，數(shù)值運算的初始值為x0 = -5; -5，要求顯示符號運算得到結(jié)構(gòu)體的每個元素的具體數(shù)值. (5分

18、)答：數(shù)值運算fsolve：x=fsolve(dzya,-5,-5,optimset(Display,off)x = 0.5671 0.5671符號運算solve：syms x y;solve(2*x-y-exp(-x),-x+2*y-exp(-y)ans = x: 0.5671 y: 0.5671 14. 符號矩陣 (1) 簡化符號矩陣然后求行列式、逆矩陣和秩； f=(4*x2+8*x)/y log(x)-exp(a*(x+y);33 sin(x)*cos(y)-cos(x)*sin(y); r,how = simple(f)r = (4*x*(x + 2)/y, log(x) - exp(

19、a*(x + y) 33, sin(x - y)how =simplify det(r)ans =(33*y*exp(a*x + a*y) + 4*x2*sin(x - y) - 33*y*log(x) + 8*x*sin(x - y)/y inv(r)ans = (y*sin(x - y)/(33*y*exp(a*(x + y) + 4*x2*sin(x - y) - 33*y*log(x) + 8*x*sin(x - y), (y*(exp(a*(x + y) - log(x)/(33*y*exp(a*(x + y) + 4*x2*sin(x - y) - 33*y*log(x) + 8*

20、x*sin(x - y) -(33*y)/(33*y*exp(a*(x + y) + 4*x2*sin(x - y) - 33*y*log(x) + 8*x*sin(x - y), (4*x*(x + 2)/(33*y*exp(a*(x + y) + 4*x2*sin(x - y) - 33*y*log(x) + 8*x*sin(x - y) rank(r)ans =2(2) jx=int(r)jx = (4*x2*(x + 3)/(3*y), x*(log(x) - 1) - exp(a*x + a*y)/a 33*x, -cos(x - y) jy=int(r,y,2,9) jy = 4*

21、x*log(9/2)*(x + 2), 7*log(x) + (exp(a*x)*(exp(2*a) - exp(9*a)/a 231, cos(x - 9) - cos(2 - x)jxy=int(r,x,y)j = -(4*(x3 + 3*x2 - y3 - 3*y2)/(3*y), y*(log(y) - 1) - x*(log(x) - 1) + (exp(a*x + a*y) - exp(2*a*y)/a33*y - 33*x, cos(x - y) - 1（3）（10分） dx=diff(r,x)dx = (4*(x + 2)/y + (4*x)/y, 1/x - a*exp(a*

22、(x + y) 0, cos(x - y) dy2=diff(r,y,2)dy2 = (8*x*(x + 2)/y3, -a2*exp(a*(x + y) 0, -sin(x - y) dxy=diff(dx,y)dxy = - (4*(x + 2)/y2 - (4*x)/y2, -a2*exp(a*(x + y) 0, sin(x - y) 15. 橢圓的隱式方程為： ，它的參數(shù)方程可以表達為：x=a*cos() y=b*sin()，其中a=5,b=3；請用subplot在三個子圖區(qū)域用不同的方法繪制橢圓曲線，要求為每幅子圖形加上標(biāo)題title和坐標(biāo)軸標(biāo)注,其中(1)采用ezplot分別繪制

23、橢圓隱式和參數(shù)方程在兩個子圖區(qū)域；(2)當(dāng)a值由1變化到5時，采用plot繪制5條橢圓曲線在另一個子圖區(qū)域中，t在-2,2范圍，使用linspace命令取橢圓100個點，使用legend命令來標(biāo)明每一條對應(yīng)的a值曲線的圖例。(10分)答：subplot(131);ezplot(y2/9+x2/25-1);title();xlabel(x);ylabel(y);syms t;subplot(132);ezplot(5*sin(t),0,2*pi);hold on;ezplot(3*cos(t),0,2*pi);axis(0 2*pi -6 6);xlabel(theta);ylabel(xy);

24、title();subplot(133);title(a15);xlabel(x);ylabel(y);a=1;t=linspace(-2*pi,2*pi,100);x=a.*cos(t);y=3.*sin(t);plot(x,y,-);hold on;a=2;x=a.*cos(t);y=3.*sin(t);plot(x,y,-);a=3;x=a.*cos(t);y=3.*sin(t);plot(x,y,:);a=4;x=a.*cos(t);y=3.*sin(t);plot(x,y,-.);a=5;x=a.*cos(t);y=3.*sin(t);plot(x,y,p);legend(a=1,a

25、=2,a=3,a=4,a=5);二、選答題（必須選4道題，每題5分，總計20分）1. 符號運算(1) 分解因式 syms xy=x4-5*x3+5*x2+5*x-6;yy=factor(y)yy =(x - 1)*(x - 2)*(x - 3)*(x + 1) syms x y z;f=x2+y2+z2+2*(x*y+y*z+z*x);ff=factor(f)ff =(x + y + z)2(2) 化簡表達式 syms a b;f=sqrt(a+sqrt(a2-b)/2)+sqrt(a-sqrt(a2-b)/2);ff=simplify(f)ff =(2(1/2)*(a + (a2 - b)(

26、1/2)(1/2) + (a - (a2 - b)(1/2)(1/2)/2 syms x y;z=sin(x)*cos(y)-cos(x)*sin(y);ff=simplify(z)ff =sin(x - y)(3) 求極限 syms x k;f=(1-1/x)(k*x);limit(f,x,inf)ans = 1/exp(k)2. 積分運算 (1)定積分運算：分別采用數(shù)值和符號運算 數(shù)值積分： s=quad(sqrt(1+cos(x),0,pi/2,1e-8)s = 2.0000符號積分： syms xy=sqrt(1+cos(x);r=int(y,0,pi/2) r =2 符號： syms

27、 x y;z=x2+y2;r=int(z,x,0,1)rr=int(r,y,0,1)r =y2 + 1/3rr =2/3數(shù)值： rr=dblquad(fxx,0,1,0,1)rr =2/3 (2) 不定積分運算 syms a b x y;z=sin(x/a)+b/y;r=int(z,x);rr=int(r,y)rr = b*x*log(y) - a*y*cos(x/a)3. 已知微分方程采用數(shù)值運算ode和符號運算dsolve求解： dsolve(D2x-1000*(1-x2)*Dx-x=0,x(0)=2,Dx(0)=0)Warning: Explicit solution could not

28、 be found. In dsolve at 120ans = empty sym ode求解：t,y=ode45(fxyz,0,40,2;0);plot(t,y)4. 采用GUIDE工具設(shè)計橢圓（x=a*cos(t) y=b*sin(t)）的用戶界面， a和b值由編輯框輸入，實現(xiàn)在坐標(biāo)軸上顯示橢圓曲線的目的。要求設(shè)計一個帶有兩個編輯框、一個坐標(biāo)軸、一個按鈕和相應(yīng)的靜態(tài)文本的圖形用戶界面，當(dāng)用鼠標(biāo)點擊按鈕時，在坐標(biāo)軸內(nèi)畫出不同a和b值的橢圓曲線。function varargout = tuoyuan2(varargin)gui_Singleton = 1;gui_State = struc

29、t(gui_Name, mfilename, . gui_Singleton, gui_Singleton, . gui_OpeningFcn, tuoyuan2_OpeningFcn, . gui_OutputFcn, tuoyuan2_OutputFcn, . gui_LayoutFcn, , . gui_Callback, );if nargin & ischar(varargin1) gui_State.gui_Callback = str2func(varargin1);end if nargout varargout1:nargout = gui_mainfcn(gui_State

30、, varargin:);else gui_mainfcn(gui_State, varargin:);endfunction tuoyuan2_OpeningFcn(hObject, eventdata, handles, varargin)handles.output = hObject; guidata(hObject, handles);function varargout = tuoyuan2_OutputFcn(hObject, eventdata, handles) varargout1 = handles.output;function a_Callback(hObject,

31、eventdata, handles)a = str2double(get(hObject, String);handles.a = a;guidata(hObject,handles) function a_CreateFcn(hObject, eventdata, handles)if ispc & isequal(get(hObject,BackgroundColor), get(0,defaultUicontrolBackgroundColor) set(hObject,BackgroundColor,white);end function b_Callback(hObject, ev

32、entdata, handles)b = str2double(get(hObject, String);handles.b = b;guidata(hObject,handles)function b_CreateFcn(hObject, eventdata, handles)if ispc & isequal(get(hObject,BackgroundColor), get(0,defaultUicontrolBackgroundColor) set(hObject,BackgroundColor,white);endfunction jisuan_Callback(hObject, eventdata, han