Matlab教程Ch6-(多项式及其运算).pptVIP

* * 2、for循环结构 使用for语句循环需要预先知道循环体执行的次数 格式 for index=start:increment:end MATLAB语句 end index的取值取决于start和end的值 在for循环语句中，不仅可以使用行向量进行循环迭代的处理，也可以使用矩阵作为循环次数的控制变量 循环的索引值将直接使用矩阵的每一列，循环的次数为矩阵的列数 * * 例：根据麦克劳林公式可以得到 n=10; p=1; s=1; for i=1:n p=p*i; s=s+1/p; fprintf(i=%.0f,s=%.8f\n,i,s) end 试求e的近似值 i=1,s=2i=2,s=2i=3,s=2i=4,s=2i=5,s=2i=6,s=2i=7,s=2i=8,s=2i=9,s=2i=10,s=2* * break语句和continue语句 break语句和continue语句的基本功能 当break语句使用在循环体中的时候，其作用是能够在执行循环体的时候强迫终止循环，即控制程序的流程，使其提前退出循环 使用方法： break; pause语句使程序运行暂停，等待用户按任意键继续。pause语句在程序调试或查看中间结果时经常使用。 使用方法 1.直接使用：pause 2.使用按时间暂停形式： pause(n) * * 例：动态显示数列极限An=(1+1/n)^n →e的逼近过程 hold on axis([0,150,2,2.8]); grid for n=2:2:150 An=(1+1/n)^n; plot(n,An,r.) pause (0.1) fprintf(n=%d,s=%.4f\n,n,An) end * * * 第六章 多项式及其运算 多项式的表示方式 * p=[1 5 0 -3 1]; f=poly2str(p,’x’) sym_p = poly2sym(p) f = x^4 + 5 x^3 - 3 x + 1 sym_p = x^4+5*x^3-3*x+1 * r = [0, 1, -1]; poly(r) %生成多项式 根为 0,1,-1 a = [1 2; 3 4]; poly(a) % 生成矩阵的特征多项式 ans = 1.0000 -5.0000 -2.0000 p = [1 0 -2 -5] %表示多项式 ans = 1 0 -1 0 生成多项式 * p=[1 -6 -72 -27]; r=roots(p) 运行结果： r = 12.1229 -5.7345 -0.3884 r=roots(p) 求多项式p（用系数行向量）表示的根. 多项式求根 * 多项式求值 polyval(p,5) %按数组规则计算5处多项式的值 ans = 110 subs(sym_p,5) % 符号变量代值 ans = 110 Y = 377 179 439 111 81 136 490 253 639 X = [2 4 5; -1 0 3; 7 1 5]; Y = polyvalm(p, X) %按矩阵规则计算X处多项式的值 y=polyval(p,X) X可以是向量也可以是矩阵. * 多项式的算术运算 加法 p = [1 0 -2 -5]; p2 = [0 2 -1 3]; add_p = p+p2 %维数必须相同 add_p = 1 2 -3 -2 减法运算类似 * 多项式的算术运算 conv(P1,P2) 多项式P1和P2的乘积 a = [1 2]; b = [2 0 -1]; c = conv(a, b); poly2sym(c) % 计算 2*x^3+4*x^2-x-2 [q,r]=deconv(P1,P2) 多项式P1和P2作除法运算, q为商，r为余数 [q, r] = deconv(c, a) %多项式c和a作除法运算 q = 2 0 -1 r = 0 0 0 0