1(III)仿真软件MATLAB--数值计算和编程.pptVIP

1.4 MATLAB的数值运算 多项式的表示 用行向的一维数组表示多项式，数组元素为多项式的系数，并按照从高阶到低阶的顺序排列 例：p(x)=x3-4x2+5x-2，表示为 p=[1 -4 5 -2] 若某项的系数为0，则必须用0补齐这些系数 例： p(x)=x5-3x3+2x p=[1 0 -3 0 2 0] 创建多项式 系数矢量直接输入法 p=[1 -4 3 2] 特征多项式输入法 a=[1 2 3;2 4 2;2 4 5]; p =poly(a) 根矢量创建多项式 r=[-0.5 -0.3+0.4i -0.3-0.4i]; p =poly(r) 多项式运算 求多项式的值： polyval(p,x);计算多项式在指定点值， polyvalm(p,x);可以接受二维数组形式的输入参数，要求方阵 求多项式的根； r=roots(p) 多项式的乘除运算； c=conv(a,b)， [q,r]=deconv(v,u), q是商多项式,r是余多项式 多项式的微积分； polyint(p),微分见下页 多项式的部分分式展开；见后页 多项式的拟合，见后页 多项式的微分 k=polyder(p),直接计算p的微分多项式k，表示k=p’ k=polyder(a,b),计算conv(a,b)的微分多项式k，表示k=(conv(a,b))’ [q,d]=polyder(a,b),计算分式a/b的微分结果，相当于q/d=(a/b)’ p=[1 2 4 2 3 5]; polyder(p) ans = 5 8 12 4 3 多项式的部分分式展开 p1, p2, …,pn, 称为极点，r1, r2, …rn称为留数，k(x)称为直项。 [r,p,k]=residue(b,a) [b,a]=residue(r,p,k) 多项式的拟合 Polyfit; P=polyfit(x,y,n); [p,S]=polyfit(x,y,n); 其中，x为自变量，y为函数，n是阶次，S是误差估计的结构。 例： x = 0: 0.1: 2.5; y=exp(x); p = polyfit(x,y,6) p = 0.0051 -0.0074 0.0715 0.1386 0.5124 0.9979 1.0000 小结：MATLAB中的多项式函数 1.5 MATLAB的控制语句 基本原则 循环语句 条件转移语句 程序的调试与优化 1.5.1 MATLBA程序的基本设计原则 ％后面的内容是程序的注解，要善于运用注解使程序更具可读性; 养成在主程序开头用clear指令清除变量的习惯，以消除工作空间中其他变量对程序运行的影响。但注意在子程序中不要用clear; 参数值要集中放在程序的开始部分，以便维护。 要充分利用MATLAB工具箱提供的指令来执行所要进行的运算，在语句行之后输入分号使其及中间结果不在屏幕上显示，以提高执行速度 input指令可以用来输入一些临时的数据；而对于大量参数，则通过建立一个存储参数的子程序，在主程序中用子程序的名称来调用; 程序尽量模块化，也就是采用主程序调用子程序的方法，将所有子程序合并在一起来执行全部的操作; 充分利用Debugger来进行程序的调试（设置断点、单步执行、连续执行），并利用其他工具箱或图形用户界面（GUI）的设计技巧，将设计结果集成到一起; 设置好MATLAB的工作路径，以便程序运行 MATLAB程序的基本组成结构 ％说明 清除命令：清除workspace中的变量和图形（clear,close） 定义变量：包括全局变量的声明及参数值的设定 逐行执行命令：指MATLAB提供的运算指令或工具箱 … … … 提供的专用命令 控制循环 ：包含for,if then,switch,while等语句 逐行执行命令 … … … end 绘图命令：将运算结果绘制出来 当然更复杂程序还需要调用子程序，或与simulink以及其他应用程序结合起来。 M文件的编辑及MATLAB工作路径的设置 进入MATLAB的Editor/Debugger窗口来编辑程序; 在编辑环境中，文字的不同颜色显示表明文字的不同属性。 绿色：注解；黑色：程序主体；红色：属性值的设定；蓝色：控制流程; 在运行程序之前，必须设置好MATLAB的工作路径，使得所要运行的程序及运行程序所需要的其他文件处在当前目录之下，否则可能导致无法读取某些系统文件或数据，从而程序无法执行 通过cd指令在命令窗口中可以更改、显示当前工作路径。 通过路径浏览器（path