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

1、學(xué)年學(xué)期課程名稱實(shí)驗(yàn)名稱實(shí)驗(yàn)室專業(yè)年級(jí)學(xué)生姓名學(xué)生學(xué)號(hào)提交時(shí)間成績(jī)?nèi)握n教師課程實(shí)驗(yàn)報(bào)告2015-2016學(xué)年第2學(xué)期MATLAB與科學(xué)計(jì)算實(shí)驗(yàn)一MATLAB環(huán)境與數(shù)值計(jì)算測(cè)量測(cè)繪實(shí)驗(yàn)中心計(jì)算機(jī)室土木13級(jí)201302016-06-02王增紅水利與建筑工程學(xué)院實(shí)驗(yàn)四Simulink仿真1、目的和要求掌握Simulink的模型窗口、創(chuàng)立過(guò)程，會(huì)設(shè)置模塊和信號(hào)線參數(shù)。2、內(nèi)容和步驟參見(jiàn)教材7.5.2小節(jié)內(nèi)容。3、實(shí)驗(yàn)報(bào)告提交要求(1) 使用Simulink模型窗口用兩種方法建立二階系統(tǒng)的仿真模型，并寫(xiě)出具體步驟。法一： 翻開(kāi)simulink窗口，新建一個(gè)模型。 在MSources模塊庫(kù)選擇Step”

2、模塊，在Conyinuous”模塊庫(kù)選擇TransferFen”模塊，在MathOperations*模塊庫(kù)選擇Sum”模塊，在“Sinks”模塊庫(kù)選擇Scope模塊。連接各模塊，從信號(hào)線引出分支點(diǎn)，構(gòu)成環(huán)閉系統(tǒng)。 設(shè)置.模塊參數(shù)，翻開(kāi)“Sum”模塊“參數(shù)設(shè)置”對(duì)話框，將“Iconshape設(shè)置為“rectangular”,將“Lisiofsigns*設(shè)置為“I。在TransferFcnw模塊的“參數(shù)設(shè)置”對(duì)話框中，將墳?zāi)苟囗?xiàng)式“Denominator設(shè)置為“10.60。在“Step”模塊的“參數(shù)設(shè)置”對(duì)話框中，將“Steptime”修改為0。 添加信號(hào)線文本注釋。 方針并分析。在Simuli

3、nk模型窗口,選擇菜單Simulation-*Simulationparameters,w命令，在“Solver*頁(yè)將“Stoptime”設(shè)置為15,單機(jī)“Startsimulation*按鈕，示波器顯示時(shí)間為15s。翻開(kāi)示波器的“y坐標(biāo)設(shè)置對(duì)話框，將y坐標(biāo)的“Y-min改為0,“Y-max改為2。運(yùn)行后得到以下圖。法二：將法一中系統(tǒng)的構(gòu)造修改成使用積分模塊(integrator)和零極點(diǎn)模塊(zero-pole)串聯(lián)，反響使用“mathoperations”模塊庫(kù)中的“gain”模塊構(gòu)造成反響環(huán)的增益為T(mén),二階系統(tǒng)構(gòu)造如圖所示，因此運(yùn)行后的波形備戰(zhàn)考試i千錘百練備戰(zhàn)考試i千錘百練圖也與法一所

4、得到的波形圖一樣。反，斗.目回區(qū)!圖二并進(jìn)展仿真。使用Simulink模型窗口創(chuàng)立如下所示的模型,SineWavelIntegrator步驟： 選擇“sources”模塊庫(kù)中選擇兩個(gè)sinewave”模塊，在MMathOperations”中選擇sum”模塊，在sinks”模塊中選擇scope”模塊，在“continuous模塊庫(kù)中選擇inlegralor模塊。 然后用信號(hào)線連接個(gè)模塊，構(gòu)成閉環(huán)系統(tǒng)圖二圖二 仿真并分析，單擊工具欄中的startsimulation*按鈕開(kāi)場(chǎng)仿真，在示波器上顯示仿真波形。實(shí)驗(yàn)一MATLAB環(huán)境與數(shù)值計(jì)算1、目的和要求(1) 熟練掌握MATLAB的啟動(dòng)和退出、MA

5、TLAB的命令窗口：常用選單和工具欄，通過(guò)實(shí)例初步認(rèn)識(shí)對(duì)MAT1.AB的命令和程序運(yùn)行等。(2) 熟練掌握MATLAB變量的使用；(3) 熟練掌握矩陣的創(chuàng)立、掌握MATLAB的矩陣和數(shù)組的運(yùn)算：(4) 熟悉MATLAB多項(xiàng)式的運(yùn)用，掌握MATLAB的擬合和插值方法、枳分運(yùn)算和統(tǒng)計(jì)運(yùn)算。2、內(nèi)容和步驟參見(jiàn)教材實(shí)驗(yàn)一、二。3、實(shí)驗(yàn)報(bào)告提交要求(1)用MALTAB命令和M程序兩種方式求解三元一次方程組。(2a-b+3c=5答：以求解三元一次方程組3a+b-5c=5為例:(4ab+c=9(2)用MALTAB命令和M程序兩種方式求解3X3復(fù)數(shù)矩陣的模和相角。MALTAB命令M程序a=2-13;31-5:

6、4-11：a=2-13;31-5;4-11;b=5;5;9;b=5;5;9;x=abx=abX=運(yùn)行結(jié)果：2X=-120-10答:MALTAB命令M程序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/pib=magic(3)運(yùn)行結(jié)果：A=a+b+l;a=r=abs(A)15013r=19161517.888518.02787.810213181418.248313.000012.2066r=4.47219.055419.105017.00001.

7、000014.3178theta=ang1e(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.1301(3)創(chuàng)立一個(gè)4X4矩陣,求逆矩陣、產(chǎn)生對(duì)斤片矩陣、求矩陣的秩等。a=fix(20*rand(4,4)duijiao=diag(a)a=duijiao=12610110113311159

8、151215180655ni=inv(a)zhi=rank(a)ni=zhi=-0.00540.0251-0.029440.0664-0.00790.08630.0074-0.0536-0.1375-0.11740.1740-0.07220.18460.0507-0.10290.0474(4)創(chuàng)立兩個(gè)4X4矩陣進(jìn)展加、減、乘法和左除運(yùn)算。a=fix(10*rand(4,4)+1)a+ba=ans=372371281189291013121331033619696691013121711b=fix(10*rand(4,4)+1)a*hb=ans=4568537911867241041191482

9、16121393662100151857681133195203136ababans=ans=-12-4-514.0000-5.000020.000028.000065-85-2.13701.6438-3.0822-3.7397010-35.4110-0.93156.24669.2192-1019-11.28773.4521-14.9726-22.7534(5)創(chuàng)立一個(gè)一元四次多項(xiàng)式，進(jìn)展彳$項(xiàng)式求值,求根計(jì)算。解：以函數(shù)y=5x+4x3+3x2+2x+l為例:P=5.4,3,2,1roots(p)%計(jì)算多項(xiàng)式的根P=ans=543210.1378+0.6782ip=5,4,3,2,1;0.1

10、378-0.6782ipolyval(p,3)%計(jì)算x=3時(shí)多項(xiàng)式的值-0.5378+0.3583ians=-0.5378-0.3583i547(6) 創(chuàng)立一個(gè)一元三次多項(xiàng)式，取十個(gè)不同自變?cè)牢唬?jì)算多項(xiàng)式的位，并分別進(jìn)展一階、二階擬合，圖示原始數(shù)據(jù)和擬合結(jié)果解：p=fix(10*rand(l,4)+l)P=7785x=l:2:20;y=polyval(p,x):pl=polyfit(x,y,1):p2=polyfit(x,y,2);yl=polyval(pl,x):y2=polyval(p2,x);plot(x,yl,r)holdonplot(x,y2)(7) 自建一個(gè)一元三次多項(xiàng)式，取十個(gè)

11、不同自變量值，計(jì)算多項(xiàng)式的值，并用線性、相鄰點(diǎn)和三次插恒，圖示其插值結(jié)果。p=fix(10*rand(l,4)+1)P=7136x=l:3:30;y=polyval(p,x):y01=interpl(x,y,12)%線性插值yOl=12774y02=interpl(x,y,12,nearest*)%相鄰點(diǎn)插值y02=15593y03=interpl(x,y,12,cubic*)%三次插值y03=1.2310e+004plot(x,y)holdonplot(12,yOl,or)plot(12,y02,*g)plot(12,y03,+)(8) 給出多組溫度的測(cè)定值(內(nèi)含平均溫度、最低、最高溫度)，

12、并求其平均溫度、最低、最高溫度的平均值、標(biāo)準(zhǔn)差等。解：a=fix(40*rand(6,3)+1)393228203931332730621617342737387mean(a)%各列的平均值ans=25.333328.666723.1667std(a)餡各列的標(biāo)準(zhǔn)差ans=13.063913.76479.5795學(xué)年學(xué)期課程名稱實(shí)驗(yàn)名稱實(shí)驗(yàn)室專業(yè)年級(jí)學(xué)生姓名學(xué)生學(xué)號(hào)提交時(shí)間成績(jī)?nèi)握n教師課程實(shí)驗(yàn)報(bào)告2015-2016學(xué)年第2學(xué)期MATLAB與科學(xué)計(jì)算實(shí)驗(yàn)二MATLAB的符號(hào)計(jì)算與可視化測(cè)量測(cè)繪實(shí)驗(yàn)中心計(jì)算機(jī)室土木13級(jí)2016-06-02王增紅水利與建筑工程學(xué)院實(shí)驗(yàn)二MATLAB的符號(hào)計(jì)算與可視

13、化1、目的和要求(1) 熟練掌握MATLAB符號(hào)表達(dá)式的創(chuàng)立、代數(shù)運(yùn)算及化簡(jiǎn)。(2) 熟悉符號(hào)方程的求解。(3) 熟練掌握MATLAB二維曲線、三維圖形的繪制。(4) 熟練掌握各種特殊圖形的繪制。2、內(nèi)容和步驟參見(jiàn)教材實(shí)驗(yàn)三、四。3、實(shí)驗(yàn)報(bào)告提交要求、1)使用sym和syms命令創(chuàng).符號(hào)變量和符號(hào)表達(dá)式；sym命令syms命令symCx,real)symsxyrealans=yy=x+i*y3Xyy=symCy=3*x8+7*x*6+x+r)ans=y=3*x8+7*x6+x+1y3*i+x(2) 用常數(shù)替換符號(hào)變量。解：symsxyz=x*6+y*5+x*yz=x6+x*y+y*5subs(

14、z,x,5)驗(yàn)用5昔換了xans=y*5+5*y+15625(3) 尋找兩個(gè)一元三次符號(hào)表達(dá)式，分別進(jìn)展其符號(hào)表達(dá)式的加、減、乘等運(yùn)算:兩個(gè)一元三次符號(hào)表達(dá)式加、減、乘運(yùn)算yl=symC3*x3+2*x2+x+10,)yl+y2%符號(hào)表達(dá)式相加yi=ans=3*x3+2*x2+x+108*x3+5*x2+2*x+30y2=symC5*x*3+3*x2+x+20*)yl-y2%符號(hào)表達(dá)式相減y2=ans=5*x3+3*x2+x+20-2*x3-x2-10yl*y2%符號(hào)表達(dá)式相乘ans=(3*x3+2*x2+x+10)*(5*x3+3*x2+x+20)(4) 尋找一個(gè)可化簡(jiǎn)的符號(hào)表達(dá)式，并化簡(jiǎn)：

15、解：y=sym(x7+x5+l)y=x7+x5+1factor(y)ans=(x2+x+l)*(x”5-x4+x3-x+1)繪制二維曲線，尋找四種二維的曲線，把圖形窗口分割為2列2行，并分別標(biāo)明圖名、坐標(biāo)等。解:M文件x=-2*pi:0.01:2*pi;subplot(221)plot(x,sin(x)title(*y=sin(x)1)xlabel(x)ylabel(y*)subplot(222)plot(x,cos(x)title(y=cos(x)xlabel(x)ylabel(y*)subplot(223)plot(x,x.A2)title(*y=xA2)xlabel(x)ylabel(y

16、*)subplot(224)plot(x,exp(x)axis(-5,5,0,12)title(，y=e*)xlabel(x*)ylabel(y*)繪制特殊圖形。參照教材中相關(guān)內(nèi)容，選取不同的數(shù)據(jù)繪制三維條形圖:M文件y-fix(50*rand(4z3)+1)x=2:5bar3(xzy)隨機(jī)產(chǎn)生正態(tài)分布的數(shù)據(jù)100個(gè)，給出統(tǒng)計(jì)數(shù)據(jù)，并繪制餅狀圖。解:M文件a=fix(100*randn(10,10)sl=0;s2=0;s3=0;s4=0;forn=l:100ifa(n)=-100sl=sl+l;elseifa(n)=0s2=s2+l;elseifa(n)100s3=s3+l;elses4=s4

17、+l;endendy=slzs2,s3zs4explode=(l111);pie(yzexplode,(1(=100)注:所有題目要求不同的學(xué)生采用不同的符號(hào)數(shù)字，報(bào)告中應(yīng)列明命令或M文件及結(jié)果，(5)(7)題要求用M文件表達(dá)。課程實(shí)驗(yàn)報(bào)告學(xué)年學(xué)期2015-2016學(xué)年第2學(xué)期課程名稱MATLAB與科學(xué)計(jì)算實(shí)驗(yàn)名稱實(shí)驗(yàn)三MATLAB程序設(shè)計(jì)實(shí)驗(yàn)室測(cè)量測(cè)繪實(shí)驗(yàn)中心計(jì)算機(jī)室專業(yè)年級(jí)級(jí)學(xué)生姓名學(xué)生學(xué)號(hào)1提交時(shí)間2016-06-02成績(jī)?nèi)握n教師王增紅水利與建筑工程學(xué)院備戰(zhàn)考試千錘百練實(shí)驗(yàn)三MATLAB程序設(shè)計(jì)1、目的和要求(1) 熟練掌握MATLAB的程序流程控制構(gòu)造。(2) 熟練掌握M文件的構(gòu)造。(

18、3) 熟練掌握函數(shù)調(diào)用和參數(shù)傳遞。2、內(nèi)容和步驟參見(jiàn)教材實(shí)驗(yàn)五。3、實(shí)驗(yàn)報(bào)告提交要求用for語(yǔ)句、while和函數(shù)調(diào)用三種方式分別對(duì)卜述兩題進(jìn)展編程，并上機(jī)進(jìn)展操作、運(yùn)行出結(jié)果。(1)cosx=11+(1)”n，30即可2!4!(2)！While語(yǔ)句：functioncos234(x,n)f=l：sum=l;i=l;whilei=nf=f*(-l)*x2/(2*i-l)/(2*i);sum=sum+f;i=i+l;endsum結(jié)果：cos234(pi,30)sum=-1IcomrrondWindowHcos234(pi,30)SUB=-1函數(shù)調(diào)用functionz=cosl23(x,n)y=l;form=l:ny=y+(-l)m*x(2*m)/jiecheng(2*m);endyfunctions=jiecheng(n)s=l;fori=l:ns=i*s;end結(jié)果：cosl23(pi,32)y=-1.0000cosl23(pi,32)y=-1.0000lFor語(yǔ)句functionz=cos222(x,n)y=l;form=l:ns=l;fori=l:2*ms=s*i;end