打印实验3.docVIP

一、实验目的 掌握 MATLAB 的数值运算及其运算中所用到的函数，掌握结构数组和细胞数组的 操作。 二、实验内容： (1) 多项式运算。 (2) 多项式插值和拟合。 (3) 数值微积分。 (4) 结构数组和细胞数组。 三、实验步骤： 1. 多项式运算 (1) 多项式表示。在MATLAB 中，多项式表示成向量的形式。 如：在MATLAB 中表示为 S=[ 1 3 -5 0 9] S=[1 3 -5 0 9] S =1 3 -5 0 9 (2) 多项式的加减法相当于向量的加减法，但须注意阶次要相同。如不同，低阶的要 补0。如多项式 S1=[0 0 2 3 11 ] S2=[1 3 -5 4 7 ] S3=S1+S2 S1=[0 0 2 3 11] S1 = 0 0 2 3 1 S2 = 1 3 -5 4 7 S3=S1+S2 S3 = 1 3 -3 7 18 (3) 多项式的乘、除法分别用函数conv 和deconv 实现 S1=[ 2 3 11 ] S2=[1 3 -5 4 7 ] S3=conv(S1,S2) S4=deconv(S3,S1) S1=[2 3 11] S1 =2 3 11 S2=[1 3 -5 4 7] S2 =1 3 -5 4 7 S3=conv(S1,S2) S3 =2 9 10 26 -29 65 77 S4=deconv(S3,S1) S4 = 1 3 -5 4 7 (4) 多项式求根用函数roots S1=[ 2 4 2 ] roots(S1) S1 = 2 4 2 roots(S1) ans = -1 -1 (5) 多项式求值用函数polyval S1=[ 2 4 1 -3 ] S1 = 2 4 1 -3 polyval(S1,3) %计算x＝3 时多项式的值 x＝1:10 polyval(S1,3) ans = 90 x=1:10 x = 1 2 3 4 5 6 7 8 9 10 y=polyval(S1,x) %计算x 向量对应的值得到y 向量 polyval(S1,x) ans = Columns 1 through 5 4 31 90 193 352 Columns 6 through 10 579 886 1285 1788 2407 2. 多项式插值和拟合 有一组实验数据如附表1-1 所示。 请分别用拟合(二阶至三阶)和插值(线性和三次样条)的方法来估测 X=9.5 时Y 的值。 以下是实现一阶拟合的语句。 x=1:10 x = 1 2 3 4 5 6 7 8 9 10 y=[16 32 70 142 260 436 682 1010 1432 1960] y =Columns 1 through 5 16 32 70 142 260 Columns 6 through 10 436 682 1010 1432 1960 p1=polyfit(x,y,1) %一阶拟合 pl = 204.8000 -522.4000 y1=polyval(p1,9.5) %计算多项式p1 在x＝9.5 时的值 y1 = 1.3208e+003 3. 数值微积分 (1) 差分使用diff 函数实现。 x=1:2:9 x = 1 3 5 7 9 diff(x) ans = 2 2 2 2 (2) 可以用因变量和自变量差分的结果相除得到数