资源描述：
《matlab 1(III) 数值计算编程和绘图ppt课件.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、1.4MATLAB的数值运算多项式的表示用行向的一维数组表示多项式，数组元素为多项式的系数，并按照从高阶到低阶的顺序排列例：p(x)=x3-4x2+5x-2，表示为p=[1-45-2]若某项的系数为0，则必须用0补齐这些系数例：p(x)=x5-3x3+2xp=[10-3020]创建多项式系数矢量直接输入法p=[1-432]特征多项式输入法a=[123;242;245];p=poly(a)根矢量创建多项式r=[-0.5-0.3+0.4i-0.3-0.4i];p=poly(r)多项式运算加减＋－：多项式的乘除运算；乘法：c=conv(a,b)除法：[q,r]=deconv(v,u)
2、,q是商多项式,r是余多项式求多项式的值：polyval(p,x);计算多项式在指定点值，polyvalm(p,x);多项式的系数为矩阵的求值，要求p为方阵求多项式的根；r=roots(p)多项式的微分k=polyder(p),直接计算p的微分多项式kk=polyder(a,b),计算conv(a,b)的微分多项式kp=[124235];>>polyder(p)ans=581243多项式的部分分式展开p1,p2,…,pn,称为极点，r1,r2,…rn称为留数，k(x)称为直项。[r,p,k]=residue(b,a)[b,a]=residue(r,p,k)多项式的拟合p=pol
3、yfit(x,y,n);x，y为利用最小二乘法进行拟合的数据，n为拟合的次数[p,s]=polyfit(x,y,n);其中，x为自变量，y为函数，n是阶次，S是误差估计的结构。例：用6阶多项式对sint函数在[0,2*pi]进行拟合x=0:2*pi/80:2*pi;y=sin(x);p=polyfit(x,y,6);y1=polyval(y,x);plot(x,y,’r’,x,y1)小结：MATLAB中的多项式函数r=roots(p)计算多项式p的根k=polyder(p)计算多项式的微分p=poly(r),p=poly(A)由根r创建多项式p，或计算A的特征多项式pk=pol
4、yder(a,b)计算两个多项式乘积的微分v=polyval(p,x),v=polyvalm(p,A)多项式求值[q,d]=polyder(a,b)计算两个多项式相除的微分p=conv(a,b)多项式相乘[r,p,k]=residue(a,b)计算两个一般形式的多项式相除计算部分分式展开式[p,r]=deconv(a,b)多项式相除[a,b]=residue(r,p,k)从部分分式展开式形式计算两个相除多项式的一般形式r=polyfit(x,y,n)将数据点x，y拟合成n次多项式1.5MATLAB的控制语句基本原则循环语句条件转移语句程序的调试与优化1.5.1MATLBA程序的
5、基本设计原则％后面的内容是程序的注解，要善于运用注解使程序更具可读性;养成在主程序开头用clear指令清除变量的习惯，以消除工作空间中其他变量对程序运行的影响。但注意在子程序中不要用clear;参数值要集中放在程序的开始部分，以便维护。要充分利用MATLAB工具箱提供的指令来执行所要进行的运算，在语句行之后输入分号使其及中间结果不在屏幕上显示，以提高执行速度input指令可以用来输入一些临时的数据；而对于大量参数，则通过建立一个存储参数的子程序，在主程序中用子程序的名称来调用;程序尽量模块化，也就是采用主程序调用子程序的方法，将所有子程序合并在一起来执行全部的操作;充分利用De
6、bugger来进行程序的调试（设置断点、单步执行、连续执行），并利用其他工具箱或图形用户界面（GUI）的设计技巧，将设计结果集成到一起;设置好MATLAB的工作路径，以便程序运行M文件的编辑及MATLAB工作路径的设置进入MATLAB的Editor/Debugger窗口来编辑程序;在编辑环境中，文字的不同颜色显示表明文字的不同属性。绿色：注解；黑色：程序主体；红色：属性值的设定；蓝色：控制流程;在运行程序之前，必须设置好MATLAB的工作路径，使得所要运行的程序及运行程序所需要的其他文件处在当前目录之下，否则可能导致无法读取某些系统文件或数据，从而程序无法执行1.5.2循环语句
7、--for-end基本格式for循环变量＝起始值：步长：终止值循环体end步长缺省值为1，可以在正实数或负实数范围内任意指定。对于正数，循环变量的值大于终止值时，循环结束；对于负数，循环变量的值小于终止值时，循环结束。循环结构可以嵌套使用。书写格式不必太过于拘泥，在Editor编辑器中会自动进行处理e1_3for循环的嵌套e1_4e1_5While循环基本格式while表达式循环体end若表达式为真，则执行循环体的内容，执行后再判断表达式是否为真，若不为真，则跳出循环体，向下继续执行Whil
