电-Matlab基础及应用3-091130.pptVIP

第3 章 Matlab数值计算 3.1 多项式 1、多项式的建立 Matlab 采用列向量表示多项式，将多项式的系数按降幂次序存放在列向量中。 2、多项式的四则运算 (1).多项式的加减运算： 相同次数的，直接相加减 不同次数的，在较低次幂的多项式系数前补0 (2).多项式乘法运算 函数conv(p1,p2)用于求多项式p1和p2的乘积。这里，p1、p2是两个多项式系数向量。 (3).多项式除法 函数[Q,r]=deconv(p1,p2)用于对多项式p1和p2作除法运算。 其中Q返回多项式p1除以p2的商式，r返回p1除以p2的余式。这里，Q和r仍是多项式系数向量。 deconv是conv的逆函数，即有p1=conv(p2,Q)+r。 3、多项式求值函数：polyval(p,x) 4、多项式求根函数：r=roots(p) 5、多项式的求导函数是p=polyder(a) 3.2 插值和拟合 如果测量值是准确的，没有误差，一般用插值；如果测量值与真实值有误差，一般用曲线拟合。 (1) 多项式的插值和拟合 3.3 数值微积分 数值微分 数值积分 用两种不同的方法求定积分 先建立一个函数文件ex.m： function ex=ex(x) ex=exp(-x.^2); 也可不建立关于被积函数的函数文件，而使用语句函数(内联函数inline)求解，命令如下： g=inline(exp(-x.^2)); %定义一个语句函数I=quadl(g,0,1) %注意函数名不加号 I = 0.7468 3.二重积分数值求解 使用MATLAB提供的dblquad函数就可以直接求出上述二重定积分的数值解。该函数的调用格式为： I=dblquad(f,a,b,c,d,tol,trace) 该函数求f(x,y)在[a,b]×[c,d]区域上的二重定积分。参数tol，trace的用法与函数quad完全相同。 3.5 稀疏矩阵 稀疏矩阵的建立 稀疏矩阵的查看和运算 第4章 结构数组与细胞数组 4.1 结构数组 结构数组的格式: 结构数组名.域名.子域名… 结构数组的创建 赋值法 struct_name(record#).field_name=data 创建1*1结构数组时可省略记录号。 1*2的结构数组，并且当结构数组包含两个以上的元素时，查看结构数组不显示各个元素的值，而是显示数组的结构信息。 利用函数struct创建 函数基本形式为： struct_name = struct(field1,{},field2,{}, ...) struct_name = struct(field1,val1,field2,val2, ...) 注意：续行符“…” 如果域没有值，则一定要赋空值，不能空着。 多个元素域值相同，可以只赋一次值。 结构数组的操作 4.2 细胞数组 细胞数组的创建 通过赋值 cell_name{i,j}={value} 利用函数cell创建细胞数组 cell_name=cell(n) cell_name=cell(m,n)或cell_name=cell([m n]) cell_name=cell(m,n,p,…)或cell_name=cell([m n p…]) cell_name=cell(size(A)) 细胞数组的操作 * 多项式 插值和拟合 数值微积分 线性方程组的数值解 稀疏矩阵 常微分方程的数值解 p：多项式的系数向量 x：为一个向量或者矩阵，要求解的x点 p:多项式的系数向量 r:方程p(x)=0的全部根(重根、复根)；r(1),r(2),…,r(n)分别代表多项式的n个根。 由给定的根创建多项式：p=poly(r) r:方程的根 p:多项式的系数向量 p=polyder(a)：求多项式a的导函数 p=polyder(a,b)：求a·b的导函数 [p,q]=polyder(a,b)：求a/b的导函数，导函数的分子存入p，分母存入q。 上述函数中，参数a,b是多项式的向量表示，结果p,q也是多项式的向量表示。 数值插值的原理：参考数值计算方法的有关内容 1、一维数值插值：被插值函数是一个单变量函数。 Y1=interp1(X,Y,X1,’method’) 根据X,Y的值，计算函数在X1处的值。 X,Y：分别描述采样点