MATLAB數值運算實驗報告

-.z實驗報告實驗報告系(部)：信息工程班級：**：**：課程：MATLAB實驗名稱：Matlab數值運算-.z目錄TOC\o"1-3"\h\u一.實驗目的2二.實驗容2三.實驗步驟2四.實驗具體過程及數據分析3五.實驗原始記錄3六.實驗心得、體會及思考3一.實驗目的掌握MATLAB的數值運算及其運算中所用到的函數，掌握構造數組和細胞數組的操作。二.實驗容多項式運算。多項式插值和擬合。數值微積分。構造數組和細胞數組。三.實驗步驟多項式運算(1)多項式表示。在MATLAB中，多項式表示成向量形式。如：s^4+3s*s^3-5*s^2+9>>S=[13-509](2)多項式的加減法相當于向量的加減法，但須注意階次要一樣。如不同，低階要補0。如多項式2*s^2+3*s+9與多項式s^4+3*s^3-5*s^2+4s+7相加。(3)多項式的乘、除法分別用函數conv和deconv實現。(4)多項式求根用函數roots(5)多項式求值用函數polyval練習1：求(s^2+1)(s+3)(s+1)/(s^3+2*s+1)的"商〞及"余〞多項式多項式插值和擬合有一組實驗數據如表所示*12345678910Y163270142260436682101014321960請分別用擬合〔二階至三階〕和插值〔線性和三次樣條〕的方法來估測*=9.5時Y的值。數值微積分差分使用diff函數的實現可以用因變量和自變量差分的結果相除得到數值微分Cumsum函數求累計積分，trapz函數用梯形法求定積分，即曲線的面積練習：如圖瑞士地圖，為了算出其國土面積，首先對地圖作如下測量：以由西向向為*軸，由南到北方為Y軸，選擇方便的原點，并將從最西邊點到最東邊界點在*軸的區(qū)間適當劃分假設干級，在每個分點的Y方向測出南邊界點和北邊界點的Y坐標Y1和Y2,這樣就得到了下表，根據地圖比例知道18mm相當于40km，試有測量數據計算瑞士國土近似面積，與其準確值41228km^2比擬。*710.51317.53440.544.548566168.576.580.591Y14445475050383030343634414546Y24459707293100110110110117118116118118*96101104106.5111.5118123.5136.5142146150157158Y143373328326555545250666668Y2121124121121121116122838182868568構造數組與細胞數組構造數組的創(chuàng)立構造數組的操作練習：創(chuàng)立一構造數組stusorce，其域為：No,Name,English，Math，Chinese，Total，Average。構造數組的大小為2*2。細胞數組的創(chuàng)立細胞數組的操作練習：創(chuàng)立一大小為2*2細胞數組stucell，其元素的類型分別為：構造類型、字符串、矩陣和細胞類型。四.實驗具體過程及數據分析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練習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)練習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