6MATLAB数值计算.pptVIP

5. 多项式的导函数 对多项式求导数的函数是： p=polyder(P) 求多项式P的导函数 p=polyder(P,Q) 求P*Q的导函数 [p,q]=polyder(P,Q) 求P/Q的导函数，导函数的分子存入p，分母存入q。 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 例6.11 求有理分式的导数。 命令如下： P=[3,5,0,-8,1,-5]; Q=[10,5,0,0,6,0,0,7,-1,0,-100]; [p,q]=polyder(P,Q) Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 6.1.5 函数的最大值与最小值 MATLAB中用于求最小值的函数是： fmin(f,a,b) 求单变量函数f(x)在区间(a,b)上的最小值点。 fmins(F,X0) 求多变量函数F(x)在估计值X0附近的最小值点。 MATLAB没有专门提供求函数最大值点的函数，但只 要注意到-f(x)在区间(a,b)上的最小值点就是f(x)在(a,b) 的最大值点，所以fmin(-f,a,b)返回函数f(x)在区间(a,b) 上的最大值。 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 例6.12 求函数f(x)在区间(-10，1)和(1，10)上的最小值点。 首先建立函数文件fx.m： function f=f(x) f=x-1/x+5; return 上述函数文件也可用一个语句函数代替： Ff=inline( x-1/x+5 ) 再在MATLAB命令窗口，输入命令： fmin(fx,-10,-1) %求函数在区间(-10，-1)内的最小值点 fmin(ff,1,10) %求函数在区间(1，10)内的最小值点。 注意函数名f不用加 有时为了描述某个数学函数的方便，可以用inline()函数来直接编写该函数，形式相当于M-函数，但无编写一个真正的MATLAB文件，就可以描述出某种数学关系。其调用格式为fun=inline(‘函数内容’,自变量列表) Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 例6.13 设有函数f(x,y,z)，求函数f在(0.5,0.5,0.5)附近的最小值。 建立函数文件fxyz.m： function f=f(u) x=u(1);y=u(2);z=u(3); f=x+y.^2./x/4+z.^2./y+2./z; return 在MALAB命令窗口，输入命令： U=fmins(fxyz,[0.5,0.5,0.5]) %求函数的最小值点 fxyz(U) %求函数的最小值 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 6.2 傅立叶分析 MATLAB中，提供了对向量(或直接对矩阵的行或列)进行离散傅立叶变换的函数， 其调用格式是： Y=fft(X,n,dim) (1)当X是一个向量时，返回对X的离散傅立叶变换。 (2)当X是一个矩阵时，返回一个矩阵并送Y，其列(行)是对X的列(行)的离散傅立叶变换。 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 例6.14 求X=(1，0，-3，5，2)的离散傅立叶变换。 在MATLAB命令窗口，输入命令： X=[1,0,-3,5,2]; Y=fft(X) %对X进行变换 3. 离散傅立叶变换的逆变换 MATLAB中，对向量(或直接对矩阵的行或列)进