MATLAB程序運行結(jié)果

closeall%關(guān)閉打開了的所有圖形窗口clc %清屏命令clear %清除工作空間中所有變量%定義時間范圍t=[0:pi/10:8*pi];y=sin(t);plot(t,y,'b:square')closeallclcclear%定義時間范圍t=[0:pi/20:9*pi];gridonholdon %允許在同一坐標(biāo)系下繪制不同的圖形plot(t,sin(t),'r:*')plot(t,cos(t))plot(t,-cos(t),'k')%gridon %在所畫出的圖形坐標(biāo)中添加?xùn)鸥?，注意用在plot之后4-1:closeallclcclear%定義時間范圍t=[0:pi/20:9*pi];holdon %允許在同一坐標(biāo)系下繪制不同的圖形plot(t,sin(t),'r:*')plot(t,cos(t))plot(t,-cos(t),'k')gridon %在所畫出的圖形坐標(biāo)中添加?xùn)鸥?，注意用在plot之后holdoff %覆蓋舊圖，自動把柵格去掉，且若要在加入柵格就%必須把gridon加在plot后面plot(t,-sin(t))gridon%主程序exp2_10.mglobala%聲明變量a為全局變量x=1:100;a=3;c=prods(x)%調(diào)用子程序prods.m%子程序prods.m%functionresult=prods(x)%globala%result=a*sum(x);%聲明了與主程序一樣的全局變量a,以便在子程序中可以%使用主程序中定義的變量答案：15150exmdl2_1.mclearcloseallclc

t=[0:pi/20:5*pi];figure(1)plot(t,out)gridonxlabel('time')ylabel('magnitude')clc%清屏clear%從內(nèi)存中清除變量和函數(shù)moreonechoon%求矩陣與矩陣的乘積，矩陣與向量的乘積A=[567;946;436]B=[345;579;731]X=[5;7;8]C=A*BY=A*Xmoreoffechooff答案：%求矩陣與矩陣的乘積，矩陣與向量的乘積A=[567;946;436]A=567946436B=[345;579;731]B=345579731X=[5;7;8]X=578C=A*BC=948386898287695553Y=A*X12312189moreoffechooffexp2_2.mclcclearmoreonechoon%為便于理解，在程序等執(zhí)行過程中顯示程序的表達式a=16;b=12;c=3;d=4;e=a+b-c*df=e/2k=e\2h=cA3g=e+f+...2+1-9aa=sin(g)abs(aa)bb=2+3jcc=conj(bb)rbb=real(bb)log(rbb)sqrt(rbb)exp(rbb)echooffmoreoffa=16;b=12;c=3;d=4;e=a+b-c*de=16f=e/2f=8k=e\2k=0.1250h=cA3h=27g=e+f+...2+1-918aa=sin(g)aa=-0.7510abs(aa)ans=0.7510bb=2+3jbb=2.0000+3.0000icc=conj(bb)cc=2.0000-3.0000irbb=real(bb)rbb=2log(rbb)ans=0.6931sqrt(rbb)ans=1.4142exp(rbb)ans=7.3891echooffans=7.3891%繪制單位圓clearcloseallclc%定義時間范圍t=[0:0.01:2*pi];x=sin(t);y=cos(t);plot(x,y)axis([-1.51.5-1.51.5])%限定x軸和y軸的顯示范圍gridonaxis('equal')%axis([xmin,xmax,ymin,ymax])函數(shù)來調(diào)整圖軸的范圍:答案：

exp2_5_.mclearcloseallclct=[0:pi/20:5*pi];plot(t,sin(t),'r:*')axis([05*pi-1.51.5])%給x軸和y軸命名xlabel('t(deg)')ylabel('magnitude')%給圖形加標(biāo)題title('sinewavefromzeroto5\pi')%在指定位置創(chuàng)建說明性文字text(pi/2,sin(pi/2),'\bullet\leftarrowThesin(t)att=2')EileEditViewInsertloolsIlesktupWindowHelp□qa曉食q鋼?聶口屈石口CLIpnl'EHEUU%圖形文字標(biāo)示命令的使用clearcloseallclct=[0:pi/20:5*pi];plot(t,sin(t),'r:*')axis([05*pi-1.51.5])%給x軸和y軸命名xlabel('t(deg)')ylabel('magnitude')%給圖形加標(biāo)題title('sinewavefromzeroto5\pi')%在指定位置創(chuàng)建說明性文字text(pi/2,sin(pi/2),'\bullet\leftarrowThesin(t)att=2')%輸入特定的字符%\pi%\alpha%\beta%\leftarrow%\rightarrow%\bullet(點號)holdonplot(t,cos(t))%區(qū)分圖形上不同的曲線legend('sin(t)','cos(t)')%用鼠標(biāo)在特定位置輸入文字gtext('文字標(biāo)示命令舉例')holdoffFigure1口回岡EileEdit里iEwInsertI.uolsHeEktopWindowHelp□qalteqn?l^s|□[b]|a□sinewavefromzeroto5nTOC\o"1-5"\h\z1.5i 1 1 ■■文字標(biāo)示命令舉例|…爭一而cos^t)mpnlcEiELUcos^t)mpnlcEiELU-1.51 1 1 L0 5 10 15t(dmg)exp2_6.m%圖形分割命令的使用clearcloseallclct=[0:pi/20:5*pi];figure(1)subplot(321)plot(t,sin(t))axis([016-1.51.5])xlabel('t(deg)')ylabel('magnitude')gridontitle('sin(t)')subplot(322)plot(t,-sin(t))axis([016-1.51.5])xlabel('t(deg)')ylabel('magnitude')gridontitle('-sin(t)')subplot(323)plot(t,cos(t))axis([016-1.51.5])xlabel('t(deg)')ylabel('magnitude')gridontitle('cos(t)')subplot(324)plot(t,-cos(t))axis([016-1.51.5])xlabel('t(deg)')ylabel('magnitude')gridontitle('-cos(t)')subplot(325)subplot(326)exp2_7.mclcclear%繪制對應(yīng)于每個輸入X的輸出y的高度條形圖subplot(221)x=[12345678910];y=[56348110356];bar(x,y)%繪制xl在以y1為中心的區(qū)間中分布的個數(shù)條形圖subplot(222)x1=randn(1,1000);%生成1000個各隨機數(shù)yl=-3:0.l:3;hist(x1,y1)%繪制y2對應(yīng)于x2的梯形圖subplot(223)x2=0:0.1:10;y2=1./(x2.人3-2.*x2+4);stairs(x2,y2)

%繪制y3對應(yīng)于x3的散點圖subplot(224)x3=0:0.1:10;y3=l./(x2.人3-2.*x2+4);stem(x3,y3)exp2_8.mechooff%不顯示程序內(nèi)容%clearallclearclca=4;b=6disp（'暫停，請按任意鍵繼續(xù)'）%disp指令可以用來顯示字符pause%暫停，直到用戶按任意鍵echoon%顯示程序內(nèi)容，注意matlab默認是不顯示c=a+b%暫時把控制權(quán)交給鍵盤(在命令窗口中出現(xiàn)k提示符),%輸入return，回車后退出，繼續(xù)執(zhí)行下面的語句。keyboardd=input('Enterthevalueofd=');%提示用戶輸入d的值e=a*24+b*8+d*6-cf=e/ddisp('workpause30second')%disp指令可以用來顯示字符pause(10)%暫停10秒disp(f)%顯示指定變量的結(jié)果，而不顯示變量名稱echooff答案：b=6暫停，請按任意鍵繼續(xù)%顯示程序內(nèi)容，注意matlab默認是不顯示c=a+bc=10%暫時把控制權(quán)交給鍵盤(在命令窗口中出現(xiàn)k提示符),%輸入return，回車后退出，繼續(xù)執(zhí)行下面的語句。keyboardK>>returnd=input('Enterthevalueofd=');%提示用戶輸入d的值Enterthevalueofd=10e=a*24+b*8+d*6-ce=194f=e/df=19.4000disp('workpause30second')%disp指令可以用來顯示字符workpause30secondpause(10)%暫停10秒disp(f)%顯示指定變量的結(jié)果，而不顯示變量名稱19.4000echooffexp2_9.msumhe.m函數(shù)%sumhe.mfunctionresult=sumhe(a,b)%sumhe(a,b)sumtheserialofnumbersfromatobresult=sum(a:b);主程序exp2_9.ma=input('pleaseinputvalueofa=');b=input('pleaseinputvalueofb=');c=sumhe(a,b)%在此例中，主程序exp2_9.m設(shè)定了兩個參數(shù)a和b來調(diào)用%子程序sumhe.m(為一個函數(shù)文件)，%所以執(zhí)行主程序時，子程序就會調(diào)用主程序傳遞過來的%a和b值進行計算，再將結(jié)果返回主程序。結(jié)果：pleaseinputvalueofa=100pleaseinputvalueofb=200c=15150%主程序exp2_10.mglobala%聲明變量a為全局變量x=1:100;a=3;c=prods(x)%調(diào)用子程序prods.m%子程序prods.mfunctionresult=prods(x)globalaresult=a*sum(x);%聲明了與主程序一樣的全局變量a,以便在子程序中可以%使用主程序中定義的變量結(jié)果：c=151509主函數(shù)：x=[123;456]mean(x)自函數(shù)：%meanfunctiony=mean(x,dim)ifnargin==1,%nargin指函數(shù)的輸入變量數(shù)；nargout指函數(shù)的輸出變量數(shù)%DeterminewhichdimensionSUMwillusedim=min(find(size(x)?=l));%size求矩陣的行數(shù)與列數(shù)，find找出非零元素的索引號ifisempty(dim),dim=1;endy=sum(x)/size(x,dim);elsey=sum(x,dim)/size(x,dim);end結(jié)果：x=l2 3456ans=2.50003.50004.5000l2 3456ans=x=123456ans=2.50003.50004.500010自函數(shù)%MYHILB生成一個Hilbert矩陣function[A,B]=myhilb(n,m)ifnargin==1,m=n;endfori=1:nforj=1:mA(i,j)=1/(i+j-1);endendnargoutif(nargout==2),B=A'*A;end主函數(shù)a=4b=4myhilb(a)結(jié)果：

ans=0ans=1.00000.50000.33330.25000.50000.33330.25000.20000.33330.25000.20000.16670.25000.20000.16670.1429exp2_13.mclcclearforj=1:10a(j)=j;ifj>5a(j)=10-j;endendj=1:10;b=a(j)結(jié)果：b=432101234543210exp2_13_.m%寫賦值程序：10%寫賦值程序：10x>=1%y=0-1<x<1%-10x<=1x=input('pleaseinputx=');ifx>=1y=10elseifx>-1&x<1y=0elsey=-10endexp2_14.mclearclccloset=-pi:0.1:pi;trigname=input('Inputtrigfunctionsname:')switchtrignamecase'sin'plot(t,sin(t))case'cos'plot(t,cos(t))otherwisebreak%終止，跳出switch語句endexp2_15.m%curvefittingofsinwaveclcclearx=0:0.1:2*pi;%生成樣本點xy=sin(x)+0.5*rand(size(x));%生成樣本