資源描述:
《matlab結(jié)課論文.doc》由會員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在學(xué)術(shù)論文-天天文庫。
1、北京林業(yè)大學(xué)全校公共選修課結(jié)課作業(yè)MATLAB程序應(yīng)用姓名:課程序號:51班級:電氣10-1學(xué)號:1010541182011年12月2011-2012學(xué)年MATLAB程序應(yīng)用結(jié)課作業(yè)-51-喬林翰1、(指選題1)有一天小猴摘下了若干個桃子,當(dāng)即吃掉了一半,還覺得不過隱,又多吃了2個。第二天接著吃下剩下的一半,又多吃了2個。以后每天都是吃掉尚存的一半零兩個。到第十天早上,小猴準(zhǔn)備吃桃子時,看到只剩下n個桃子了。問小猴第一天共摘下多少個桃子?1)寫出求解小猴第一天共摘下多少個桃子的數(shù)學(xué)表達(dá)式2)根據(jù)確定的數(shù)學(xué)表達(dá)式,利用mat
2、lab編寫出一個函數(shù)文件用于求解3)編寫一個命令文件,完成輸入一個n值,然后調(diào)用上述函數(shù)求解,并將結(jié)果輸出顯示出來【思路和方法】假設(shè)猴子第一天吃掉……因此本題可以轉(zhuǎn)化利用matlab求解級數(shù)的問題……(正文格式:中文宋體,西文timesNewRoman,1.25倍行距)【源程序】inputn%輸入n【運(yùn)行結(jié)果】>>display('Pleaseinputthevalueofn');n=input('n=');fori=1:9;n=(n+2)*2;endm=nPleaseinputthevalueofnn=112011-20
3、12學(xué)年MATLAB程序應(yīng)用結(jié)課作業(yè)-51-喬林翰2、(指選題2)1、分別用矩陣求逆、矩陣除法、矩陣分解求下面線性方程組的解;矩陣除法:A=1-23-401-1113010-731>>b=[4,-3,1,-3]b=4-31-3>>x=A/bx=0.7143-0.2000-0.22860.6000矩陣求逆A=[1,-2,3,-4;0,1,-1,1;1,3,0,1;0,-7,3,1]A=1-23-401-11112011-2012學(xué)年MATLAB程序應(yīng)用結(jié)課作業(yè)-51-喬林翰13010-731>>b=[4,-3,1,-3]b=
4、4-31-3x=inv(A)*bx=0.7143-0.2000-0.22860.6000矩陣分解法(LU分解):[L,U]=lu(A);>>x=U/(L/b)x=0.7143-0.2000-0.22860.60004、某鄉(xiāng)鎮(zhèn)企業(yè)某產(chǎn)品年銷售額與所獲利潤從2000年到2010年的數(shù)據(jù)如下:年度20002001200220032004200520062007200820092010銷售額x6.17.59.410.714.617.421.124.429.832.934.3純利潤4.56.48.38.49.711.513.715.
5、417.720.522.3試求對的經(jīng)驗(yàn)回歸直線方程,并做回歸分析,對回歸的參數(shù)進(jìn)行區(qū)間估計(jì);x=[6.1,7.5,9.4,10.7,14.6,17.4,21.1,24.4,29.8,32.9,34.3];>>y=[4.5,6.4,8.3,8.4,9.7,11.5,13.7,15.4,17.7,20.5,22.3];>>p=polyfit(x,y,1)112011-2012學(xué)年MATLAB程序應(yīng)用結(jié)課作業(yè)-51-喬林翰p=0.56391.9096xi=6.1:34.3;>>yi=polyval(p,xi)yi=Columns
6、1through85.34915.91306.47687.04077.60458.16848.73229.2961Columns9through169.860010.423810.987711.551512.115412.679213.243113.8069Columns17through2414.370814.934615.498516.062316.626217.190017.753918.3177Columns25through2918.881619.445420.009320.573121.1370plot(x,y
7、,':o',xi,yi,'-*')112011-2012學(xué)年MATLAB程序應(yīng)用結(jié)課作業(yè)-51-喬林翰112011-2012學(xué)年MATLAB程序應(yīng)用結(jié)課作業(yè)-51-喬林翰3、(自選題)>>surf(x,y,z)title('surf(x,y,z)')t=0:pi/20:2*pi;[x,y,z]=cylinder(2+sin(t),30);subplot(2,1,2);surf(x,y,z);112011-2012學(xué)年MATLAB程序應(yīng)用結(jié)課作業(yè)-51-喬林翰4、(自選題)我們是電氣工程及其自動化專業(yè),對電路分析有時會需要進(jìn)
8、行傅里葉變換和拉普拉斯變換,通過學(xué)習(xí)這門課,我們發(fā)現(xiàn)matlab可以實(shí)現(xiàn)對一個函數(shù)進(jìn)行快速傅氏變換(fft),以下為具體操作:對函數(shù)x(t)=12sin(10πt+π/2)+10(cos20πt)進(jìn)行fft:N=128N=128>>T=1T=1>>t=linspace(0,T,N);>>x=12*sin