matlab常用工具[]

《matlab常用工具[]》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、Matlab工具一、多项式1．找出多项式的根，即多项式为零的值。在MATLAB里，多项式由一个行向量表示，它的系数是按降序排列。例如，输入多项式x4－12x3＋0x2＋25x＋116»p=[1-12025116]p=1-12025116注意，必须包括具有零系数的项。给出这种形式，用函数roots找出一个多项式的根。»r=roots(p)%matlab中对字母严格区分大分小写r=11.74732.7028-1.2251+1.4672i-1.2251-1.4672i因为在MATLAB中，无论是一个多项式，还是它的根，都是

2、向量，MATLAB按惯例规定，多项式是行向量，根是列向量。给出一个多项式的根，也可以构造相应的多项式。在MATLAB中，命令poly执行这个任务。»pp=poly(r)pp=1.0e+002*Columns1through40.0100-0.12000.00000.2500Column51.1600+0.0000i»pp=real(pp)%throwawayspuriousimaginarypartpp=1.0000-12.00000.000025.0000116.0000MATLAB当用根重组多项式时，如果一些根有

3、虚部，由于截断误差，则poly的结果有一些小的虚部，这是很普通的。消除虚假的虚部，只要使用函数real抽取实部。2．导数MATLAB为多项式求导提供了函数polyder。»gg=162048697244»h=polyder(g)h=63080144138723．估值根据多项式系数的行向量，可对多项式进行估值。在MATLAB中，这由函数polyval来完成。»x=linspace(-1,3);%choose100datapointsbetween-1and3.»p=[14-7-10];%usespolynomialp(

4、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')%helpplot4．图像的绘制plot在线性坐标系中绘制二维数据plot3绘制三维数据plot函数调用格式如下:plot(Y)plot(X1,Y1,…)plot(X1,Y1,Linespec,…)plot(…,‘PropertName’,PropertyValue,…

5、)h=plot(…)其中,Plot(Y)中如Y为一个矢量,则绘制连接(i,Y(i))的折线;若Y中的1个元素为复数,则绘制(real(Y(i)),imag(Y(i))的折线;若Y为矩阵,则将每一列都视为一个矢量绘制一条折线.。5．说到绘图，只要计算函数在某一区间的值，并且画出结果向量，这样就得到了函数的图形。在大多数情况下，这就足够了。然而，有时一个函数在某一区间是平坦的并且无激励，而在其它区间却失控。在这种情况下，运用传统的绘图方法会导致图形与函数真正的特性相去甚远。MATLAB提供了一个称为fplot的巧妙的绘图

6、函数。该函数细致地计算要绘图的函数，并且确保在输出的图形中表示出所有的奇异点。该函数的输入需要知道以字符串表示的被画函数的名称以及2元素数组表示的绘图区间。>>f='2*exp(-x).*sin(x)';>>fplot(f,[08]);>>title(f),xlabel('x')的曲线6．极小化作图除了提供视觉信息外，还常常需要确定一个函数的其它更多的特殊属性。在许多应用中，特别感兴趣的是确定函数的极值，即最大值（峰值）和最小值（谷值）。数学上，可通过确定函数导数（斜率）为零的点，解析上求出这些极值点。显然，如果定义

7、的函数简单，则这种方法常常奏效。然而，即使很多容易求导的函数，也常常很难找到导数为零的点。在这种情况下，以及很难或不可能解析上求得导数的情况下，必须数值上寻找函数的极值点。MATLAB提供了两个完成此功能的函数fmin和fmins。这两个函数分别寻找一维或n维函数的最小值。这里仅讨论fmin。有关fmins的详细信息，参阅《MATLAB参考指南》。因为f(x)的最大值等于-f(x)的最小值，所以，上述fmin和fmins可用来求最大值和最小值。如果还不清楚，把上述图形倒过来看，在这个状态下，峰值变成了谷值，而谷值则变

8、成了峰值。考虑上述例子，从图中可知，在xmax=0.7附近有一个最大值，并且在xmin=4附近有一个最小值。而这些点的解析值为：和。用文本编辑器编写一个脚本M文件，并用fmin寻出数值上极值点，给出函数主体如下：%ex_fmin.mfn='2*exp(-x)*sin(x)';%definefunctionforminxmin=fmin(fn,2