matlab工具.doc.docVIP

Matlab工具 一、多项式 1．找出多项式的根，即多项式为零的值。在MATLAB里，多项式由一个行向量表示，它的系数是按降序排列。例如，输入多项式x4－12x3＋0x2＋25x＋116 ? p=[1 -12 0 25 116] p = 1 -12 0 25 116 注意，必须包括具有零系数的项。给出这种形式，用函数roots找出一个多项式的根。 ? r=roots(p)matlab中对字母严格区分大分小写 r = 11.7473 2.7028 -1.2251 + 1.4672i -1.2251 - 1.4672i 因为在MATLAB中，无论是一个多项式，还是它的根，都是向量，MATLAB按惯例规定，多项式是行向量，根是列向量。给出一个多项式的根，也可以构造相应的多项式。在MATLAB中，命令poly执行这个任务。 ? pp=poly(r) pp = 1.0e+002 * Columns 1 through 4 0.0100 -0.1200 0.0000 0.2500 Column 5 1.1600 + 0.0000i ? pp=real(pp) %throw away spurious imaginary part pp = 1.0000 -12.0000 0.0000 25.0000 116.0000 MATLAB当用根重组多项式时，如果一些根有虚部，由于截断误差，则poly的结果有一些小的虚部，这是很普通的。消除虚假的虚部，只要使用函数real抽取实部。 2．导数 MATLAB为多项式求导提供了函数polyder。 ? g g = 1 6 20 48 69 72 44 ? h=polyder(g) h = 6 30 80 144 138 72 估值 根据多项式系数的行向量，可对多项式进行估值。在MATLAB中，这由函数polyval来完成。 ? x=linspace(-1, 3) ; % choose 100 data points between -1and 3. ? p=[1 4 -7 -10] ; % uses polynomial p(x) = x3＋4x2－7x－10 ? v=polyval(p , x) ; 计算x值上的p(x)，把结果存在v里。然后用函数plot绘出结果。 ? plot(x , v),title( x^3+4x^2-7x-10 ), xlabel( x ) %help plot 图像的绘制 plot 在线性坐标系中绘制二维数据 plot3 绘制三维数据 plot 函数调用格式如下: plot(Y) plot(X1,Y1,…) plot(X1,Y1, Linespec, …) plot(…, ‘PropertName’, Property Value, …) h=plot(…) 其中, Plot(Y)中如Y为一个矢量,则绘制连接(i, Y(i))的折线; 若Y中的1个元素为复数,则绘制(real(Y(i)), imag(Y(i))的折线; 若Y为矩阵,则将每一列都视为一个矢量绘制一条折线.。 说到绘图，只要计算函数在某一区间的值，并且画出结果向量，这样就得到了函数的图形。在大多数情况下，这就足够了。然而，有时一个函数在某一区间是平坦的并且无激励，而在其它区间却失控。在这种情况下，运用传统的绘图方法会导致图形与函数真正的特性相去甚远。MATLAB提供了一个称为fplot的巧妙的绘图函数。该函数细致地计算要绘图的函数，并且确保在输出的图形中表示出所有的奇异点。该函数的输入需要知道以字符串表示的被画函数的名称以及2元素数组表示的绘图区间。 f=2*exp(-x) .* sin(x); fplot(f , [0 8]); title(f) , xlabel(x) 的曲线 极小化 作图除了提供视觉信息外，还常常需要确定一个函数的其它更多的特殊属性。在许多应用中，特别感兴趣的是确定函数的极值，即最大值（峰值）和最小值（谷值）。数学上，可通过确定函数导数（斜率）为零的点，解析上求出这些极值点。显然，如果定义的函数简单，则这种方法常常奏效。然而，即使很多容易求导的函数，也常常很难找到导数为零的点。在这种情