实验三MATLAB数值运算.docVIP

实验三 MATLAB数值运算 一、实验目的 掌握 MATLAB 的数值运算及其运算中所用到的函数，掌握结构数组和细胞数组的操作。 二、实验内容： 1 多项式运算。 2 多项式插值和拟合。 3 数值微积分。 三、实验步骤： 1. 多项式运算 1 多项式表示。在 MATLAB 中，多项式表示成向量的形式。 如： 在 MATLAB 中表示为 S [ 1 3 -5 0 9] 2 多项式的加减法相当于向量的加减法，但须注意阶次要相同。如不同，低阶的要补 0。如多项式与多项式相加。 S1 [0 0 2 3 11 ] S2 [1 3 -5 4 7 ] S3 S1+S2 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＝1:10 y polyval S1,x 2. 多项式插值和拟合 有一组实验数据如附表 1-1 所示。 X 1 2 3 4 5 6 7 8 9 10 Y 16 32 70 142 260 436 682 1010 1432 1960 请分别用拟合 二阶至三阶 和插值 线性和三次样条 的方法来估测 X 9.5 时 Y 的值。以下是实现一阶拟合的语句。 x 1:10 y [16 32 70 142 260 436 682 1010 1432 1960] p1 polyfit x,y,1 y1 polyval p1,9.5 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 3 cumsum函数求累计积分，trapz 函数用梯形法求定积分，即曲线的面积。 x ones 1,10 cumsum x x linspace 0, pi,100 ; y sin x ; S trapz y,x 四、实验报告 1 写出实验内容语句的注释及运算结果； 2 写出实验体会