MATLAB數(shù)值運(yùn)算實(shí)驗(yàn)報(bào)告

-.z實(shí)驗(yàn)報(bào)告實(shí)驗(yàn)報(bào)告系(部)：信息工程班級(jí)：**：**：課程：MATLAB實(shí)驗(yàn)名稱：Matlab數(shù)值運(yùn)算-.z目錄TOC\o"1-3"\h\u一.實(shí)驗(yàn)?zāi)康?二.實(shí)驗(yàn)容2三.實(shí)驗(yàn)步驟2四.實(shí)驗(yàn)具體過(guò)程及數(shù)據(jù)分析3五.實(shí)驗(yàn)原始記錄3六.實(shí)驗(yàn)心得、體會(huì)及思考3一.實(shí)驗(yàn)?zāi)康恼莆誐ATLAB的數(shù)值運(yùn)算及其運(yùn)算中所用到的函數(shù)，掌握構(gòu)造數(shù)組和細(xì)胞數(shù)組的操作。二.實(shí)驗(yàn)容多項(xiàng)式運(yùn)算。多項(xiàng)式插值和擬合。數(shù)值微積分。構(gòu)造數(shù)組和細(xì)胞數(shù)組。三.實(shí)驗(yàn)步驟多項(xiàng)式運(yùn)算(1)多項(xiàng)式表示。在MATLAB中，多項(xiàng)式表示成向量形式。如：s^4+3s*s^3-5*s^2+9>>S=[13-509](2)多項(xiàng)式的加減法相當(dāng)于向量的加減法，但須注意階次要一樣。如不同，低階要補(bǔ)0。如多項(xiàng)式2*s^2+3*s+9與多項(xiàng)式s^4+3*s^3-5*s^2+4s+7相加。(3)多項(xiàng)式的乘、除法分別用函數(shù)conv和deconv實(shí)現(xiàn)。(4)多項(xiàng)式求根用函數(shù)roots(5)多項(xiàng)式求值用函數(shù)polyval練習(xí)1：求(s^2+1)(s+3)(s+1)/(s^3+2*s+1)的"商〞及"余〞多項(xiàng)式多項(xiàng)式插值和擬合有一組實(shí)驗(yàn)數(shù)據(jù)如表所示*12345678910Y163270142260436682101014321960請(qǐng)分別用擬合〔二階至三階〕和插值〔線性和三次樣條〕的方法來(lái)估測(cè)*=9.5時(shí)Y的值。數(shù)值微積分差分使用diff函數(shù)的實(shí)現(xiàn)可以用因變量和自變量差分的結(jié)果相除得到數(shù)值微分Cumsum函數(shù)求累計(jì)積分，trapz函數(shù)用梯形法求定積分，即曲線的面積練習(xí)：如圖瑞士地圖，為了算出其國(guó)土面積，首先對(duì)地圖作如下測(cè)量：以由西向向?yàn)?軸，由南到北方為Y軸，選擇方便的原點(diǎn)，并將從最西邊點(diǎn)到最東邊界點(diǎn)在*軸的區(qū)間適當(dāng)劃分假設(shè)干級(jí)，在每個(gè)分點(diǎn)的Y方向測(cè)出南邊界點(diǎn)和北邊界點(diǎn)的Y坐標(biāo)Y1和Y2,這樣就得到了下表，根據(jù)地圖比例知道18mm相當(dāng)于40km，試有測(cè)量數(shù)據(jù)計(jì)算瑞士國(guó)土近似面積，與其準(zhǔn)確值41228km^2比擬。*710.51317.53440.544.548566168.576.580.591Y14445475050383030343634414546Y24459707293100110110110117118116118118*96101104106.5111.5118123.5136.5142146150157158Y143373328326555545250666668Y2121124121121121116122838182868568構(gòu)造數(shù)組與細(xì)胞數(shù)組構(gòu)造數(shù)組的創(chuàng)立構(gòu)造數(shù)組的操作練習(xí)：創(chuàng)立一構(gòu)造數(shù)組stusorce，其域?yàn)椋篘o,Name,English，Math，Chinese，Total，Average。構(gòu)造數(shù)組的大小為2*2。細(xì)胞數(shù)組的創(chuàng)立細(xì)胞數(shù)組的操作練習(xí)：創(chuàng)立一大小為2*2細(xì)胞數(shù)組stucell，其元素的類型分別為：構(gòu)造類型、字符串、矩陣和細(xì)胞類型。四.實(shí)驗(yàn)具體過(guò)程及數(shù)據(jù)分析1.>>S1=[242]S1=242>>roots(S1)ans=-1-1>>S=[13-509]S=13-509>>S1=[2311]S1=2311>>S2=[13-547]S2=13-547>>S3=conv(S1,S2)S3=291026-296577>>S4=deconv(S3,S1)S4=13-547>>S1=[242]S1=242>>polyval(S1,3)ans=32>>*=1:10*=12345678910>>y=polyval(S1,*)y=81832507298128162200242練習(xí)1：>>clearall>>s1=[101]s2=[13]s3=[11]s4=conv(s1,s2)Y=conv(s4,s3)*=[1021][Q,R]=deconv(Y,*)s1=101s2=13s3=11s4=1313Y=14443*=1021Q=14R=002-5-1>>poly2sym(Q)ans=*+4>>poly2sym(R)ans=2**^2-5**-1>>conv(Q,*)+R-Yans=00000>>s1=[101]s2=[13]s3=[11]s4=conv(s1,s2)Y=conv(s4,s3)*=[1021][Q,R]=deconv(Y,*)s1=101s2=13s3=11s4=1313Y=14443*=1021Q=14R=002-5-1poly2sym(Q)ans=*+4poly2sym(R)ans=2**^2-5**-1conv(Q,*)+R-Yans=000002.*=1:10y=[163270142260436682101014321960]p1=polyfit(*,y,1)y1=polyval(p1,9.5)3.*=1:2:9diff(*)*=linspace(0,2*pi,100);y=sin(*);plot(*,y)y1=diff(y)./diff(*);plot(*(1:end-1),y1)*=ones(1,10)cumsum(*)*=linspace(0,pi,100);y=sin(*);s=trapz(*,y)練習(xí)2：*=[710.51317.53440.544.548566168.576.580.59196101104106.5111.5118123.5136.5142146150157158];y1=[444547505038303034363441454643373328326555545250666668];y2=[4459707293100110110110117118116118118121124121121121116122838182868568];*=*./18*40;Y1=y1./18*40;Y2=y2./18*40;t1=trapz(*,Y1),t2=trapz(*,Y2),t=t2-t14.〔1〕>>student.number='20050731001';>>='jack';>>student(2).number='20050731002';>>student(2).name='lucy'；〔2〕>>student(1).subject=[]>>student(1).sorce=[]>>student>>fieldnames(student)>>getfield(student,{2},'name')>>student=rmfield(student,'subject')>>student=setfield(student,{1},'sorce',90);>>stuent(2).sorce=88;〔3〕>>A={'Howareyou!',ones(3);[12;34],{'cell'}};>>B(1,1)={'Helloworld'};>>B(2,1)={[1234]};〔4〕>>ans1=A(1,1