matlab實(shí)驗(yàn)(西北農(nóng)林科技大學(xué))

水利與建筑工程學(xué)院課程實(shí)驗(yàn)報(bào)告課程名稱MATLAB與科學(xué)計(jì)算 實(shí)驗(yàn)室測(cè)量測(cè)繪實(shí)驗(yàn)中心計(jì)算機(jī)室 學(xué)生姓名 學(xué)生學(xué)號(hào)20130 提交時(shí)間2016-06-02 成績(jī)王增紅1(1)熟練掌握MATLAB的啟動(dòng)和退出、MATLAB的命令窗口；常用選單和工具欄，通過(guò)實(shí)例初步(2)熟練掌握MATLAB變量的使用；(3)熟練掌握矩陣的創(chuàng)建、掌握MATLAB的矩陣和數(shù)組的運(yùn)算；(4)熟悉MATLAB多項(xiàng)式的運(yùn)用，掌握MATLAB的擬合和插值方法、積分運(yùn)算和統(tǒng)計(jì)運(yùn)算。二。 (1)用MALTAB命令和M程序兩種方式求解三元一次方程組。答：以求解三元一次方程組{為例：a=[2-13;31-5;4-11];b=[5;5;9];x=a\bx=2-10>>a=[2-13;31-5;4-11];>>b=[5;5;9];>x=a\bx=2-10>>a=fix(20*rand(3,3))a=fix(20*rand(3,3))a=b=magic(3);16185A=a+b*i;181210r=abs(A)2119theta=angle(A)*180/pi>>b=magic(3);運(yùn)行結(jié)果：>>A=a+b*I;a=>>r=abs(A)15013r=19161517.888518.02787.810213181418.248313.000012.2066r=4.47219.055419.105017.00001.000014.3178>>theta=angle(A)*180/pi19.235416.763116.5529theta=13.601520.124614.142126.56513.179850.1944theta=9.462322.619934.992028.072590.000024.775163.434983.65986.00908.972617.354025.016917.102726.56518.13012a(3)創(chuàng)建一個(gè)4╳4矩陣，求逆矩陣、產(chǎn)生對(duì)角矩陣、求矩陣的秩等。a>>>>a=fix(20*rand(4,4))a=2101190>>ni=inv(a)ni=-0.00540.0664-0.0079-0.0536-0.1375-0.07220.18460.04746360.02510.0863-0.11740.050735-0.02940.00740.1740-0.1029>>duijiao=diag(a)duijiao=15>>zhi=rank(a)zhi=4(4)創(chuàng)建兩個(gè)4╳4矩陣進(jìn)行加、減、乘法和左除運(yùn)算。>>>>a=fix(10*rand(4,4)+1)=289>>b=fix(10*rand(4,4)+1)b=846abans=-160-1a+bans=76a*bans=a\bans=-2.1370110-11.2877920.0000-3.08226.2466-14.972686-5.0000-0.931552128.0000-3.73979.2192-22.7534-4-80-55-393936384237549625079362363((5)創(chuàng)建一個(gè)一元四次多項(xiàng)式，進(jìn)行多項(xiàng)式求值，求根計(jì)算。>>p=[5,4,3,2,1]p=54321>>p=[5,4,3,2,1];ans=>>roots(p)%計(jì)算多項(xiàng)式的根ans=0.1378+0.6782i0.1378-0.6782i-0.5378+0.3583i-0.5378-0.3583i3(6)創(chuàng)建一個(gè)一元三次多項(xiàng)式，取十個(gè)不同自變量值，計(jì)算多項(xiàng)式的值，并分別進(jìn)行一階、二階擬合，圖示原始數(shù)據(jù)和擬合結(jié)果>>>>p=fix(10*rand(1,4)+1)p=7785>>x=1:2:20;>>y=polyval(p,x);>>p1=polyfit(x,y,1);p2=polyfit(x,y,2);>>y1=polyval(p1,x);y2=polyval(p2,x);>>plot(x,y1,'r')>>holdon>>plot(x,y2)(7)自建一個(gè)一元三次多項(xiàng)式，取十個(gè)不同自變量值，計(jì)算多項(xiàng)式的值，并用線性、相鄰點(diǎn)和三次插值，圖示其插值結(jié)果。>>>>p=fix(10*rand(1,4)+1)p=7136>>x=1:3:30;>>y=polyval(p,x);>>y01=interp1(x,y,12)%線性插值y01=>>y02=interp1(x,y,12,'nearest')%相鄰點(diǎn)插值y02=>>y03=interp1(x,y,12,'cubic')%三次插值y03=e004>>plot(x,y)>>holdon>>plot(12,y01,'or')>>plot(12,y02,'*g')>>plot(12,y03,'+')(8)給出多組溫度的測(cè)定值(內(nèi)含平均溫度、最低、最高溫度)，并求其平均溫度、最低、最高溫度的平均值、標(biāo)準(zhǔn)差等。>>a=fix(40*rand(6,3)+1)4a=a393228203931332730621617342737387>>mean(a)%各列的平均值ans=25.333328.666723.1667>>std(a)%各列的標(biāo)準(zhǔn)差ans=13.063913.76479.5795水利與建筑工程學(xué)院課程實(shí)驗(yàn)報(bào)告學(xué)年學(xué)期2015-2016學(xué)年第2學(xué)期 課程名稱MATLAB與科學(xué)計(jì)算實(shí)驗(yàn)室測(cè)量測(cè)繪實(shí)驗(yàn)中心計(jì)算機(jī)室 專業(yè)年級(jí)土木13級(jí) 學(xué)生姓名學(xué)生學(xué)號(hào)2016-06-02成績(jī)?nèi)握n教師王增紅 1ans=xans=xx+x+10)*(5*x^3+3*x^2+x+20)(1)熟練掌握MATLAB符號(hào)表達(dá)式的創(chuàng)建、代數(shù)運(yùn)算及化簡(jiǎn)。(2)熟悉符號(hào)方程的求解。(3)熟練掌握MATLAB二維曲線、三維圖形的繪制。(4)熟練掌握各種特殊圖形的繪制。四。(1)使用sym和syms命令創(chuàng)建符號(hào)變量和符號(hào)表達(dá)式；>>sym('x','real')ans=x>>sym('y=3*x^8+7*x^6+x+1')ans=y=3*x^8+7*x^6+x+1>>symsxyreal>>yy=x+i*y^3yy=y^3*i+x(2)用常數(shù)替換符號(hào)變量。>>z=x^6+y^5+x*yx^6+x*y+y^5ans=y^5+5*y+15625(3)尋找兩個(gè)一元三次符號(hào)表達(dá)式，分別進(jìn)行其符號(hào)表達(dá)式的加、減、乘等運(yùn)算；兩兩個(gè)一元三次符號(hào)表達(dá)式>>y1=sym('3*x^3+2*x^2+x+10')y1=3*x^3+2*x^2+x+10>>y2=sym('5*x^3+3*x^2+x+20')y2=5*x^3+3*x^2+x+20加、減、乘運(yùn)算>>y1+y2%符號(hào)表達(dá)式相加ans=8*x^3+5*x^2+2*x+30>>y1-y2%符號(hào)表達(dá)式相減ans=-2*x^3-x^2-10>>y1*y2%符號(hào)表達(dá)式相乘(4)尋找一個(gè)可化簡(jiǎn)的符號(hào)表達(dá)式，并化簡(jiǎn)；y=yx^7+x^5+1>>factor(y)ans=(x^2+x+1)*(x^5-x^4+x^3-x+1)2x=-2*pi:0.01:2*pi;subplot(221)plot(x,sin(x))title('y=sin(x)')xlabel('x')ylabel('y')subplot(222)plot(x,cos(x))title('y=cos(x)')xlabel('x')ylabel('y')subplot(223)plot(x,x.^2)title('y=x^2')xlabel('x')ylabel('y')subplot(224)plot(x,exp(x))axis([-5,5,0,12])title('y=e^x')xlabel('x')ylabel('y')(6)繪制特殊圖形。參照教材中相關(guān)內(nèi)容，選取不同的數(shù)據(jù)繪制三維條形圖；yy=fix(50*rand(4,3)+1)x=2:5bar3(x,y)3(7)隨機(jī)產(chǎn)生正態(tài)分布的數(shù)據(jù)100個(gè)，給出統(tǒng)計(jì)數(shù)據(jù)，并繪制餅狀圖。a=fix(100*randn(10,10))s1=0;s2=0;s3=0;s4=0;forn=1:100ifa(n)<=-100s1=s1+1;elseifa(n)<=0s2=s2+1;elseifa(n)<100s3=s3+1;elses4=s4+1;y=[s1,s2,s3,s4]explode=[1111];pie(y,explode,{'(<=-100)','(-99~0)','(1~99)','(>=100)’})(7)題要求用M文件表達(dá)。水利與建筑工程學(xué)院課程實(shí)驗(yàn)報(bào)告課程名稱MATLAB與科學(xué)計(jì)算 B實(shí)驗(yàn)室測(cè)量測(cè)繪實(shí)驗(yàn)中心計(jì)算機(jī)室 級(jí)學(xué)生姓名 學(xué)生學(xué)號(hào)1 提交時(shí)間2016-06-02 成績(jī)王增紅1(1)熟練掌握MATLAB的程序流程控制結(jié)構(gòu)。(2)熟練掌握M文件的結(jié)構(gòu)。(3)熟練掌握函數(shù)調(diào)用和參數(shù)傳遞。functioncos234(x,n)f=1;sum=1;i=1;whilei<=nf=f*(-1)*x^2/(2*i-1)/(2*i);sum=sum+f;functionz=cos123(x,n)yform=1:ny=y+(-1)^m*x^(2*m)/jiecheng(2*m);yfunctions=jiecheng(n)fori=1:nFor語(yǔ)句functionz=cos222(x,n)yform=1:nfori=1:2*my=y+(-1)^m.*x^(2*m)/s;y>>cos234(pi,30)-1>>cos123(pi,32)y=-1.0000>>cos222(pi,34)y=-1.00002(2)ln(1+x)=x++…+(1)n+…(-1＜x≤1)n≥2000即可functionzky(x,n)fori=1:ns=s+(-1)^i*x^(i+1)/(i+1);sfunctionzky2(x,n)sxjwhilej<=ns=s+(-1)^j*x^(j+1)/(j+1);sfunctionzky3(x,n)fori=1:ns=s+xx(x,i);sfunctio