MATLAB軟件入門 第一次數(shù)學(xué)實(shí)驗(yàn)

重慶大學(xué)學(xué)生實(shí)驗(yàn)報(bào)告實(shí)驗(yàn)課程名稱數(shù)學(xué)實(shí)驗(yàn)開課實(shí)驗(yàn)室學(xué)院年級專業(yè)班學(xué)生姓名學(xué)號開課時(shí)間2011至2012學(xué)年第2學(xué)期總成績教師簽名數(shù)學(xué)與統(tǒng)計(jì)學(xué)院制開課學(xué)院、實(shí)驗(yàn)室：數(shù)學(xué)與統(tǒng)計(jì)DS1422實(shí)驗(yàn)時(shí)間：2012年3月29日課程名稱數(shù)學(xué)實(shí)驗(yàn)實(shí)驗(yàn)項(xiàng)目名稱MATLAB軟件入門實(shí)驗(yàn)項(xiàng)目類型驗(yàn)證演示綜合設(shè)計(jì)其他指導(dǎo)教師肖劍成績實(shí)驗(yàn)?zāi)康腫1]熟悉MATLAB軟件的用戶環(huán)境；[2]了解MATLAB軟件的一般目的命令；[3]掌握MATLAB數(shù)組操作與運(yùn)算函數(shù)；[4]掌握MATLAB軟件的基本繪圖命令；[5]掌握MATLAB語言的幾種循環(huán)、條件和開關(guān)選擇結(jié)構(gòu)?；A(chǔ)實(shí)驗(yàn)一、實(shí)驗(yàn)內(nèi)容1．MATLAB軟件的數(shù)組操作及運(yùn)算練習(xí)；2．直接使用MATLAB軟件進(jìn)行作圖練習(xí)；3．用MATLAB語言編寫命令M-文件和函數(shù)M-文件。二、實(shí)驗(yàn)過程（一般應(yīng)包括實(shí)驗(yàn)原理或問題分析，算法設(shè)計(jì)、程序、計(jì)算、圖表等，實(shí)驗(yàn)結(jié)果及分析）1．設(shè)有分塊矩陣，其中E,R,O,S分別為單位陣、隨機(jī)陣、零陣和對角陣，試通過數(shù)值計(jì)算驗(yàn)證。>>E=eye(3);>>R=rand(3,2);>>O=zeros(2,3);>>a=[23];>>S=diag(a);>>A=[ER;OS];>>B=[ER+R*S;OS^2];>>A^2-B結(jié)果：ans=00000000000000000000000002．某零售店有9種商品的單件進(jìn)價(jià)（元）、售價(jià)（元）及一周的銷量如表1.1，問哪種商品的利潤最大，哪種商品的利潤最??；按收入由小到大，列出所有商品及其收入；求這一周該10種商品的總收入和總利潤。表1.1貨號123456789單件進(jìn)價(jià)7.158.253.2010.306.6812.0316.8517.519.30單件售價(jià)11.1015.006.0016.259.9018.2520.8024.1515.50銷量568120575358039521041538810694M文件程序：function[zs,zl]=ht(t)a=[7.158.253.2010.306.6812.0316.8517.519.30];b=[11.1015.006.0016.259.9018.2520.8024.1515.50];c=[568120575358039521041538810694];[m,mi]=max(b-a)[n,ni]=min(b-a)[l,li]=sort(b.*c)zs=sum(b.*c)zl=sum((b-a).*c)結(jié)果：htm=6.7500mi=2n=2.8000ni=3l=1.0e+004*Columns1through80.39110.45180.63050.94251.07571.80751.95623.1990Column93.8398li=531492876zs=1.4294e+005zl=4.6052e+004ans=1.4294e+0053.用兩種方法在同一個(gè)坐標(biāo)下作出y1=x2,y2=x3,y3=x4y4=x5這四條曲線的圖形，并要求用兩種方法在圖上加各種標(biāo)注。（1）M文件程序：function[x,y]=ht0(x)x=linspace(-1,1,20);y1=x.^2;y2=x.^3;y3=x.^4;y4=x.^5;plot(x,y1)holdon,plot(x,y2);plot(x,y3);plot(x,y4);holdofftext(-0.6,0.4,'\leftarrowy1=x^2')text(-0.9,0.2,'y2=x^3\rightarrow')text(-0.9,-0.1,'y3=x^4\rightarrow')text(-0.8,-0.6,'\leftarrowy4=x^5')xlabel('x')ylabel('y')title('curves')結(jié)果：>>ht0ans=Columns1through8-1.0000-0.8947-0.7895-0.6842-0.5789-0.4737-0.3684-0.2632Columns9through16-0.1579-0.05260.05260.15790.26320.36840.47370.5789Columns17through200.68420.78950.89471.0000圖片：（2）M文件程序：function[x,y]=ht2(x)x=linspace(-1,1,20);y1=x.^2;y2=x.^3;y3=x.^4;y4=x.^5;plot(x,y1,x,y2,x,y3,x,y4)gtext('y1=x^2');gtext('y2=x^3');gtext('y3=x^4');gtext('y4=x^5')xlabel('x')ylabel('y')title('curves')圖片：4．作出下列曲面的3維圖形，；M文件程序：function[Z]=ht3(x,y)x=-1.41:0.5:1.41;y=x;[X,Y]=meshgrid(x,y);R=sqrt(X.^2+Y.^2)+eps;Z=sin(R*pi);mesh(X,Y,Z);結(jié)果：>>ht3ans=-0.0187-0.8474-0.9951-0.9628-0.9960-0.6321-0.8474-0.78420.00600.2656-0.2624-0.9685-0.99510.00600.96870.96840.7736-0.4943-0.96280.26560.96840.38930.9541-0.2902-0.9960-0.26240.77360.95410.4971-0.6833-0.6321-0.9685-0.4943-0.2902-0.6833-0.9915圖片：5．建立一個(gè)命令M-文件：求所有的“水仙花數(shù)”，所謂“水仙花數(shù)”是指一個(gè)三位數(shù)，其各位數(shù)字的立方和等于該數(shù)本身。例如，153是一個(gè)水仙花數(shù)，因?yàn)?53=13+53+33。M文件程序：functionf=fun(t)fort=100:999x=fix(t/100);y=fix((t-x*100)/10);z=t-x*100-y*10;if(x^3+y^3+z^3)==tA=tendend結(jié)果:>>funA=153A=370A=371A=4076．編寫函數(shù)M-文件sq.m：用迭代法求的值。求平方根的迭代公式為迭代的終止條件為前后兩次求出的x的差的絕對值小于105。M文件程序：functionx1=sq(a)x0=eps;x1=(1/2)*(x0+a/x0);while(abs(x1-x0)>=10^(-5))x0=x1;x1=(1/2)*(x0+a/x0);end結(jié)果:>>sq(4)ans=2.0000總結(jié)與體會通過這次的實(shí)驗(yàn)操作讓自己對MATLAB的工作環(huán)境有了一個(gè)大概了解，并通過一定的練習(xí)和看書，是自己對MATLAB的基本函數(shù)有了了解，并能進(jìn)行使用。但在練習(xí)與操作過程中與暴露自己的不足，這表明自己在有些方面還是存在缺陷的，自己要注