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轴，选择方便的原点，并将从最西边点到最东边界点在*轴的区间适当划分假设干级，在每个分点的Y方向测出南边界点和北边界点的Y坐标Y1和Y2,这样就得到了下表，根据地图比例知道18mm相当于40km，试有测量数据计算瑞士国土近似面积，与其准确值41228km^2比拟。*710.51317.53440.544.548566168.576.580.591Y14445475050383030343634414546Y24459707293100110110110117118116118118*96101104106.5111.5118123.5136.5142146150157158Y143373328326555545250666668Y2121124121121121116122838182868568构造数组与细胞数组构造数组的创立构造数组的操作练习：创立一构造数组stusorce，其域为：No,Name,English，Math，Chinese，Total，Average。构造数组的大小为2*2。细胞数组的创立细胞数组的操作练习：创立一大小为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