实验四 数值运算.pptxVIP

实验四 数值运算;一、实验目的 掌握 MATLAB 的数值运算及其运算中所用到的函数。 二、实验内容： (1) 多项式运算。 (2) 多项式插值和拟合。 (3) 数值微积分。 ;?;(3) 多项式的乘、除法分别用函数conv 和deconv 实现 S1=[ 2 3 11 ] S2=[1 3 -5 4 7 ] S3=conv(S1,S2) S4=deconv(S3,S1) (4) 多项式求根用函数roots S1=[ 2 4 2 ] roots(S1) (5) 多项式求值用函数polyval S1=[ 2 4 1 -3 ] polyval(S1,3) %计算x＝3 时多项式的值 x＝1:10 y=ployval(S1,x) %计算x 向量对应的值得到y 向量;2. 多项式插值和拟合 有一组实验数据如下表所示。 请分别用拟合(二阶至三阶)和插值(线性和三次样条)的方法来估测 X=9.5 时Y 的值。 以下是实现一阶拟合的语句。 x=1:10 y=[16 32 70 142 260 436 682 1010 1432 1960] p1=ployfit(x,y,1) %一阶拟合 y1=ployval(p1,9.5) %计算多项式p1 在x＝9.5 时的值;例题;格式：zi=interp2(x,y,z,xi,yi,’method’) x,y为已知的点坐标向量，z为矩阵（x,y对应点的值）xi，yi 为插入点的 X,Y 坐标向量，‘method’：同上 ;例： 假设有一个汽车发动机在转速为 2000r/min 时，温度(单位为℃)与时间(单位为s)的5 个测量值如表所示。;例4-11 某观测站测得某日6:00时至18:00时之间每隔2小时的室内外温度(℃)，用3次样条插值分别求得该日室内外6:30至17:30时之间每隔2小时各点的近似温度(℃)。 设时间变量h为一行向量，温度变量t为一个两列矩阵，其中第一列存放室内温度，第二列储存室外温度。命令如下： h =6:2:18; t=[18,20,22,25,30,28,24;15,19,24,28,34,32,30]; XI =6.5:2:17.5 YI=interp1(h,t,XI,‘spline’) %用3次样条插值计算 ;3. 数值微分 (1) 差分使用diff 函数实现。 x=1:2:9 diff(x) (2) 可以用因变量和自变量差分的结果相除得到数值微分。 x=linspace(0,2*pi,100); y=sin(x); plot(x,y) y1=diff(y)./diff(x); plot(x(1:end-1),y1);4. 数值积分;4. 数值积分;练习;3. 大气压强p 随高度x 变化的理论公式为p =1.0332e?( x+500) / 7756，为验证这一公式，测得某地大气压强随高度变化的一组数据如表所示，试用插值法和拟合法进行计算并绘图，看哪种方法较为合理，且总误差最小;4，利用梯形法和辛普森法求定积分 的值，并对结果进行比较。如果积分区间改为-5～5 结果有何不同？梯形积分中改变自变量x 的维数，结果有何不同？;练习5：图A1 是瑞士地图，为了算出其国土面积，首先对地图作如下测量：以由西向东方向为X 轴，由南到北方向为Y 轴，选择方便的原点，并将从最西边界点到最东边界点在X 轴上的区间适当划分为若干段，在每个分点的Y 方向测出南边界点和北边界点的Y 坐标Y1 和Y2，这样就得到了表1，根据地图比例尺知道18mm 相当于40km，试由测量数据计算瑞士国土近似面积，与其精确值41228km2 比较。地图的数据见下表1-2。