數(shù)學(xué)軟件Matlab實(shí)驗(yàn)報(bào)告

PAGE計(jì)算機(jī)與信息學(xué)院（數(shù)學(xué)類課程）實(shí)驗(yàn)報(bào)告課程名稱：數(shù)學(xué)軟件姓名：系：應(yīng)用數(shù)學(xué)系專業(yè)：數(shù)學(xué)與應(yīng)用數(shù)學(xué)年級(jí)：2010級(jí)學(xué)號(hào)：指導(dǎo)教師：職稱：講師201

實(shí)驗(yàn)項(xiàng)目列表序號(hào)實(shí)驗(yàn)項(xiàng)目名稱成績(jī)指導(dǎo)教師1MATLAB運(yùn)算基礎(chǔ)林運(yùn)國(guó)2MATLAB矩陣分析與處理林運(yùn)國(guó)3選擇結(jié)構(gòu)程序設(shè)計(jì)，循環(huán)結(jié)構(gòu)程序設(shè)計(jì)（1）林運(yùn)國(guó)4循環(huán)結(jié)構(gòu)程序設(shè)計(jì)（2），函數(shù)文件林運(yùn)國(guó)5高層繪圖操作林運(yùn)國(guó)67891011121314151617181920PAGE15計(jì)算機(jī)與信息學(xué)院數(shù)學(xué)類實(shí)驗(yàn)報(bào)告（一）系：應(yīng)用數(shù)學(xué)專業(yè)：數(shù)學(xué)與應(yīng)用數(shù)學(xué)年級(jí)：2010級(jí)姓名：***學(xué)號(hào)：102260002051實(shí)驗(yàn)課程：數(shù)學(xué)軟件實(shí)驗(yàn)室號(hào)：_田C-513實(shí)驗(yàn)設(shè)備號(hào)：107實(shí)驗(yàn)時(shí)間：2011/09/25指導(dǎo)教師簽字：成績(jī)：1．實(shí)驗(yàn)項(xiàng)目名稱：MATLAB運(yùn)算基礎(chǔ)2．實(shí)驗(yàn)?zāi)康暮鸵螅?.熟悉啟動(dòng)和退出MTALAB的方法。2.熟悉MATLAB命令窗口的組成。3.掌握建立矩陣的方法。4.掌握MATLAB各種表達(dá)式的書寫規(guī)則以及常用函數(shù)的使用。3．實(shí)驗(yàn)使用的主要儀器設(shè)備和軟件：聯(lián)想微機(jī)；MATLAB6.54．實(shí)驗(yàn)的基本理論和方法:冒號(hào)表達(dá)式，邏輯表達(dá)式，whos，find，length，空矩陣5．實(shí)驗(yàn)內(nèi)容與步驟：1(1)z1=2*sin(pi*85/180)/(1+exp(2))(2)x=[2,1+2*i;-0.455];z2=(1/2)*log(x+sqrt(1+x^2))(3)a=-3.0:0.1:3.0;z3=(exp(0.3*a)-exp(-0.3*a))/2.*sin(a+0.3)+log((0.3+a)/2)(4)t=0:0.5:2.5;z4=t.^2.*(t>=0&t<1)+(t.^2-1).*(t>=1&t<2)+(t.^2-2*t+1).*(t>=2&t<3)2A=[12,34,-4;34,7,87;3,65,7];B=[1,3,-1;2,0,3;3,-2,7];A+6*BB+eye(3)(2)A*BA.*B(3)A^3A.^3(4)A/BB\A(5)[A,B][A([1,3],:);B^2]3A=[1,2,3,4,5;6,7,8,9,10;11,12,13,14,15;16,17,18,19,20;21,22,23,24,25];B=[3,0,16;17,-6,9;0,23,-4;9,7,0;4,13,11];C=A*BD=C(3:5,2:3)Whos4(1)A=100:999;B=rem(A,21);C=length(find(B==0))A='IlOveMatLAb2051';a=find(A>='A'&A<='Z');A(a)=[]6．實(shí)驗(yàn)心得（質(zhì)疑、建議）：學(xué)會(huì)了一些Matlab的基礎(chǔ)，并利用這些知識(shí)進(jìn)行基本的數(shù)值處理。計(jì)算機(jī)與信息學(xué)院數(shù)學(xué)類實(shí)驗(yàn)報(bào)告（二）系：應(yīng)用數(shù)學(xué)專業(yè)：數(shù)學(xué)與應(yīng)用數(shù)學(xué)年級(jí)：2010級(jí)姓名：***學(xué)號(hào)：102260002051實(shí)驗(yàn)課程：數(shù)學(xué)軟件實(shí)驗(yàn)室號(hào)：_田C-513實(shí)驗(yàn)設(shè)備號(hào)：107實(shí)驗(yàn)時(shí)間：2011/10/17指導(dǎo)教師簽字：成績(jī)：1．實(shí)驗(yàn)項(xiàng)目名稱：MATLAB矩陣分析與處理2．實(shí)驗(yàn)?zāi)康暮鸵螅?.掌握生成特殊矩陣的方法。2.掌握矩陣分析的方法。3.用矩陣求逆法解線性方程組。3．實(shí)驗(yàn)使用的主要儀器設(shè)備和軟件：聯(lián)想微機(jī)；MATLAB6.54．實(shí)驗(yàn)的基本理論和方法:單位矩陣eye，隨機(jī)矩陣rand，零矩陣zeros，對(duì)角陣diag，條件數(shù)cond，行列式det，跡rank，秩trace，范數(shù)norm，逆inv5．實(shí)驗(yàn)內(nèi)容與步驟：1E=eye(3);R=rand(3,2);O=zeros(2,3);S=diag([2,3]);A=[E,R;O,S];A^2B=[E,(R+R*S);O,S^2]2H=hilb(5);P=pascal(5);Hh=det(H)Hp=det(P)Th=cond(H)Tp=cond(P)3A=[1,2,3,4,5;6,7,8,9,10;11,12,13,34,53;6,3,23,7,15;6,26,9,41,28];det(A)rank(A)trace(A)norm(A)4A=[-29,6,18;20,5,12;-8,8,5][V,D]=eig(A)5(1)A=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6];b=[0.95；0.67；0.52];x=inv(A)*bB1=[0.95;0.67;0.53];X1=inv(A)*B1N=cond(B)N1=cond(B1)6A=fix(10*rand(5))B=sqrtm(A)C=sqrt(A)6．實(shí)驗(yàn)心得（質(zhì)疑、建議）：學(xué)會(huì)了利用本實(shí)驗(yàn)涉及到的代碼對(duì)矩陣進(jìn)行各種求解，對(duì)計(jì)算和解答方便了非常多。只是求什么量要用哪個(gè)代碼還不熟練，要查表。因?yàn)闆](méi)發(fā)現(xiàn)規(guī)律，又太多，易搞混。在以后的實(shí)驗(yàn)操作中要多多練習(xí)，掌握熟練。計(jì)算機(jī)與信息學(xué)院數(shù)學(xué)類實(shí)驗(yàn)報(bào)告（三）系：應(yīng)用數(shù)學(xué)專業(yè)：數(shù)學(xué)與應(yīng)用數(shù)學(xué)年級(jí)：2010級(jí)姓名：***學(xué)號(hào)：102260002051實(shí)驗(yàn)課程：數(shù)學(xué)軟件實(shí)驗(yàn)室號(hào)：_田C-513實(shí)驗(yàn)設(shè)備號(hào)：107實(shí)驗(yàn)時(shí)間：2011/10/29指導(dǎo)教師簽字：成績(jī)：1．實(shí)驗(yàn)項(xiàng)目名稱：選擇結(jié)構(gòu)程序設(shè)計(jì)，循環(huán)結(jié)構(gòu)程序設(shè)計(jì)（1）2．實(shí)驗(yàn)?zāi)康暮鸵螅?.掌握建立和執(zhí)行M文件的方法。2.掌握利用if語(yǔ)句實(shí)現(xiàn)選擇結(jié)構(gòu)的方法。3.掌握利用switch語(yǔ)句實(shí)現(xiàn)多分支選擇結(jié)構(gòu)的方法。4．掌握try語(yǔ)句的使用。5.掌握利用for語(yǔ)句實(shí)現(xiàn)循環(huán)結(jié)構(gòu)的方法。6.掌握利用while語(yǔ)句實(shí)現(xiàn)循環(huán)結(jié)構(gòu)的方法。7.熟悉利用向量運(yùn)算來(lái)代替循環(huán)操作的方法。3．實(shí)驗(yàn)使用的主要儀器設(shè)備和軟件：聯(lián)想微機(jī)；MATLAB6.54．實(shí)驗(yàn)的基本理論和方法：If語(yǔ)句，switch語(yǔ)句，try語(yǔ)句，循環(huán)結(jié)構(gòu)：for語(yǔ)句、while語(yǔ)句，向量運(yùn)算代替循環(huán)操作的方法。5．實(shí)驗(yàn)內(nèi)容與步驟：3.1x=[-5.0,-3.0,1.0,2.0,2.5,3.0,5.0];y=[];forx0=xifx0<0&x0~=-3y=[y,x0*x0+x0-6];elseifx0>=0&x0<5&x0~=2&x0~=3y=[y,x0*x0-5*x0+6];elsey=[y,x0*x0-x0-1];endendxy3.2if語(yǔ)句:x=input('pleaseenteranumberbetween0to100：');ifx>100|x<0disp('Error!!thenumberyouinputiswrong!!');elseifx<=100&x>=90disp('A');elseifx<=89&x>=80disp('B');elseifx<=79&x>=70disp('C');elseifx<=69&x>60disp('D');elsedisp('E');endendswitch語(yǔ)句:x=input('pleaseenteranumberbetween0to100：');while1ifx<0|x>100disp('Error!!thenumberyouinputiswrong!!')x=input('pleaseinputanumberbetween0to100：');elsebreak;endendswitchfix(x/10)case{9,10}disp('A')case8disp('B')case7disp('C')case6disp('D')otherwisedisp('E')end3.3n=input('inputtheworker’sID：');t=input('inputtheworker’sworkingtime：');ift>120x=(t-120)*84*(1+0.15)+120*84;elseift<60x=t*84-700;elsex=t*84;enddisp([num2str(n),'worker’ssalaryis',num2str(x)]);3.4a=fix(10+90*rand(1,2))x=a(1);y=a(2);t=input('pleaseinputtheyunsuanfuhao：','s');ift=='+'z=x+y;elseift=='-'z=x-y;elseift=='*'z=x*y;elseift=='/'z=x/y;enddisp([num2str(x),t,num2str(y),'=',num2str(z)])3.5a=rand(5,6)n=input('whilelinedoyouwant?：');ifn>5disp('overtherange，thelastlineis：')a(5,:)elsedisp(['No.',num2str(n),'''slineis：'])a(n,:)end4.1循環(huán)結(jié)構(gòu)pi=0;n=input('n=');fori=1:npi=sqrt(pi^2+6*(1/i/i));endpi向量運(yùn)算n=input('n=');i=1:n;f=6./i.^2;pi=sqrt(sum(f))4.2y=0;n=1;while(y<3)y=y+1/(2*n-1);n=n+1;endy=y-1/(2*(n-1)-1)n=n-24.3a=input('a=');b=input('b=');Xn=1;Xn1=a/(b+Xn);n=0;whileabs(Xn1-Xn)>1e-5Xn=Xn1;Xn1=a/(b+Xn);n=n+1;ifn==500break;endendnXn1r1=(-b+sqrt(b*b+4*a))/2r2=(-b-sqrt(b*b+4*a))/24.4fori=1:100ifi==1f(i)=1;elseifi==2f(i)=0;elseifi==3f(i)=1;elsef(i)=f(i-1)-2*f(i-2)+f(i-3);endendmax(f)min(f)sum(f)length(find(f>0))length(find(f==0))length(find(f<0))4.5s=0;n=0;fori=2:49b=i*(i+1)-1;m=fix(sqrt(b));forj=2:mifrem(b,j)==0breakendendifj==mn=n+1;s=s+b;endendns6．實(shí)驗(yàn)心得（質(zhì)疑、建議）：本次實(shí)驗(yàn)內(nèi)容多，而且也比較復(fù)雜，但是上機(jī)時(shí)間太少，個(gè)人覺得完成有點(diǎn)倉(cāng)促。建議老師能多安排上機(jī)時(shí)間，畢竟這也是一門操作性比較強(qiáng)的學(xué)科。計(jì)算機(jī)與信息學(xué)院數(shù)學(xué)類實(shí)驗(yàn)報(bào)告（四）系：應(yīng)用數(shù)學(xué)專業(yè)：數(shù)學(xué)與應(yīng)用數(shù)學(xué)年級(jí)：2010級(jí)姓名：***學(xué)號(hào)：102260002051實(shí)驗(yàn)課程：數(shù)學(xué)軟件實(shí)驗(yàn)室號(hào)：_田C-513實(shí)驗(yàn)設(shè)備號(hào)：107實(shí)驗(yàn)時(shí)間：2011/11/14指導(dǎo)教師簽字：成績(jī)：1．實(shí)驗(yàn)項(xiàng)目名稱：循環(huán)結(jié)構(gòu)程序設(shè)計(jì)（2），函數(shù)文件2．實(shí)驗(yàn)?zāi)康暮鸵螅?.理解函數(shù)文件的概念。2.掌握定義和調(diào)用MATLAB函數(shù)的方法。3．實(shí)驗(yàn)使用的主要儀器設(shè)備和軟件：聯(lián)想微機(jī)；MATLAB6.54．實(shí)驗(yàn)的基本理論和方法：建立M文件，并在命令窗口調(diào)用。5．實(shí)驗(yàn)內(nèi)容與步驟：1建立M文件mat1.mfunctiony=mat1(x)y=[exp(x),log(x),sin(x),cos(x)];在命令窗口調(diào)用mat1.m文件y=mat(1+2*i)%假設(shè)給定的復(fù)數(shù)為1+2*i2建立M文件dierti.mfunction[a,b,N,M]=dierti(m,n,t)A=[m*cos(t*pi/180),-m,-sin(t*pi/180),0;m*sin(t*pi/180),0,cos(t*pi/180),0;0,n,-sin(t*pi/180),0;0,0,-cos(t*pi/180),1];B=[0,9.8*m,0,9.8*n];C=inv(A)*B';a=C(1);b=C(2);N=C(3);M=C(4);在命令窗口調(diào)用M文件m1=input('m1=');m2=input('m2=');theta=input('theta=');[a1,a2,N1,N2]=dierti(m1,m2,theta)建立M文件mat2.mfunctionX=mat2(m1,m2,t)g=9.8;A=[m1*cos(t*pi/180),-m1,-sin(t*pi/180),0;m1*sin(t*pi/180),0,cos(t*pi/180),0;0,m2,-sin(t*pi/180),0;0,0,-cos(t*pi/180),1];B=[0;m1*g;0;m2*g];X=inv(A)*B;在命令窗口調(diào)用M文件X=mat2(1,1,60)3建立M文件mat3.mfunctionflag=mat3(x)flag=1;fori=2:sqrt(x)ifrem(x,i)==0flag=0;break;endend在命令窗口調(diào)用M文件fori=10:99j=10*rem(i,10)+fix(i/10);ifmat3(i)&mat3(j)disp(i)endend4建立M文件fx.mfunctiony=fx(x)y=1./((x-2).^2+0.1)+1./((x-3).^4+0.01);在命令窗口調(diào)用M文件y=fx(3)a=[1,2,3;4,5,6;7,8,9];y=fx(a)5(1)建立M文件mat4.mfunctionf1=mat4(n)f1=n+10*log(n*n+5);在命令窗口調(diào)用M文件y=mat4(40)/(mat4(30)+mat4(20))(2)建立M文件mat5.mfunctionf2=mat5(n)f2=0;fori=1:nf2=f2+i*(i+1);end在命令窗口調(diào)用M文件y=mat5(40)/(mat5(30)+mat5(20))6．實(shí)驗(yàn)心得（質(zhì)疑、建議）：難度有所提高，但實(shí)用性也是非常強(qiáng)大的！通過(guò)自己寫函數(shù)來(lái)完成更復(fù)雜的題目，完成之后很有成就感。但實(shí)驗(yàn)安排的時(shí)間還是不夠，一次實(shí)驗(yàn)只做了4道題目，課本的知識(shí)還掌握不牢固，會(huì)找個(gè)時(shí)間好好鞏固下的。計(jì)算機(jī)與信息學(xué)院數(shù)學(xué)類實(shí)驗(yàn)報(bào)告（五）系：應(yīng)用數(shù)學(xué)專業(yè)：數(shù)學(xué)與應(yīng)用數(shù)學(xué)年級(jí)：2010級(jí)姓名：***學(xué)號(hào)：102260002051實(shí)驗(yàn)課程：數(shù)學(xué)軟件實(shí)驗(yàn)室號(hào)：_田C-513實(shí)驗(yàn)設(shè)備號(hào)：107實(shí)驗(yàn)時(shí)間：2011/11/24指導(dǎo)教師簽字：成績(jī)：1．實(shí)驗(yàn)項(xiàng)目名稱：高層繪圖操作2．實(shí)驗(yàn)?zāi)康暮鸵螅?.掌握繪制二維圖形的常用函數(shù)。2.掌握繪制三維圖形的常用函數(shù)。3.掌握繪制圖形的輔助操作。3．實(shí)驗(yàn)使用的主要儀器設(shè)備和軟件：聯(lián)想微機(jī)；MATLAB6.54．實(shí)驗(yàn)的基本理論和方法：Linespace，plot，subplot，title，steam，bar，stairs，fill（要把顏色代碼填進(jìn)去），polar等。5．實(shí)驗(yàn)內(nèi)容與步驟：1x=linspace(0,2*pi,101);y=(0.5+3*sin(x)./(1+x.^2)).*cos(x);plot(x,y)2(1)x=linspace(-2*pi,2*pi,100);y1=x.^2;y2=cos(2*x);y3=y1.*y2;plot(x,y1,'r-',x,y2,'b:',x,y3,'g--');text(-4,16,'\leftarrowy1=x^2');text(3*pi/2,-0.5,'\uparrowy2=cos(2*x)');text(-1.5*pi,-2.25*pi*pi,'\leftarrowy3=y1*y2');(2)x=linspace(-2*pi,2*pi,100);y1=x.^2;y2=cos(2*x);y3=y1.*y2;subplot(2,2,1);plot(x,y1);title('y1=x^2');subplot(2,2,2);plot(x,y2);title('y2=cos(2*x)');subplot(2,2,3);plot(x,y3);title('y3=x^2*cos(2*x)');（3）x=linspace(-2*pi,2*pi,50);y1=x.^2;subplot(2,2,1);bar(x,y1);title('y1=x^2的條形圖');subplot(2,2,2);stairs(x,y1);title('y1=x^2的階梯圖');subplot(2,2,3);stem(x,y1);title('y1=x^2的桿圖');subplot(2,2,4);fill(x,y1,'b');title('y1=x^2的填充圖');x=linspace(-2*pi,2*pi,50);y2=cos(2*x);subplot(2,2,1);bar(x,y2);title('y2=cos(2*x)的條形圖');subplot(2,2,2);stairs(x,y2);title('y2=co