matlab计算导数并求值(matlab 计算导数)

用matlab怎么求在给定x为0y为1二阶导数的值

matlab求导数的方法:1、打开matlab软件;2、输入指令“sysm x”定义一个符号变量,并输入指令“f1=atan(x)”定义一个函数;3、输入“df1=diff(f1,x)”指令求出导数;4、输入“subplot”指令输出即可。

matlab中对函数求导并求极值

用matlab对函数求导及其极值 ,可以先用diff()函数求出y\’,然后令y\’=0,用vpasolve()解出x值。实现代码如下:

syms x

z=exp(4.1399-0.0235*x);

y=0.000157*z*(1-z)*(1+z)^3/((1+z)^4+0.00004476*z^2)^(3/2);

dy=diff(y,1)

x=vpasolve(dy==0)

运行结果

x=120.125

matlab数值导数计算

先确定 假设在[-3,3]区间(定义域)内以0.01为步长求数值导数。

第一种:首先用一个4次多项式p(x)拟合y=f(x),并对p(x)求一般意义下的导数dp(x),然后求出dp(x)在假设点的值;

第二种:直接求y=f(x)在假设点的数值导数;

第三种:先手工算出y=f(x)导函数,然后求在假设点的值。

程序如下:

f=inline(\’sqrt(x.^3+2*x.^2-x+12)+7*x.^2+36\’);

x=-3:0.01:3;

p=polyfit(x,f(x),4);

dp=polyder(p);

dpx=polyval(dp,x);

dx=diff(f([x,3.01]))/0.01;

g=inline(\'(3*x.^2+4*x-1)./sqrt(x.^3+2*x.^2-x+12)/2+14*x\’);

gx=g(x);

plot(x,dpx,x,dx,\’.\’,x,gx,\’-\’);

matlab求导后赋值计算

问题描述:

MATLAB 函数先求导再赋值

简单的说就是建立如下m文件

function j=myfun(x)

j=diff(x^2-1)

进行如下调用

s=myfun(2)

得到

s=

[ ]

也就是说在函数myfun调用时,是先赋值后求导的。

而我希望得到先求导后赋值的结果s=4

其实我的实际程序是如下结构的

建立如下m文件

function j=myfun(x,a,b)

syms z1 z2 z3

t=@(z1,z2,z3)z1^3-z2^2+z3;

to=t(a,b,x(1))

tc=t(b,a,x(2))

m=[to;tc]

d=[a,b]

j=jacobian(m,d)

进行调用

s=@(x)myfun(x,2,3)

fminsearch(ss,[1,2])

得到

??? Function \’jacobian\’ is not defined for values of class \’double\’

……

我觉得这个和第一个是同样的错误,因为先赋值了,所以求不出雅克比矩阵了。

请不要提供手算导数再赋值的答案,因为我的真正的程序方程比较复杂,如下

function j=myfun(x,a,b,c)

syms z1 z2 z3

t=@(z1,z2,z3)z1*(z2+z3)/((z3^2+(z1+z2)*z3)*exp((z2+z3)*c)+z1*z2)+exp((z2+z3)*c);

to=mt(1/a,1/b,x(1));

tc=mt(1/b,1/a,x(2));

m=[to;tc];

d=[a,b];

j=det(inv(jacobian(m,d)));

如果用syms是可以用subs部分赋值,但下一步:fminsearch(s,[1,2])就无法进行了。

要如何再转换成函数呢?

1个回答

分类:

综合

2014-09-30

问题解答:

我来补答

1、利用微分diff的零点,驻点值和边界值进行比较;

2、利用fmincon等优化函数

展开全文阅读

怎么用MATLAB求微分并求值?

用MATLAB求微分并求值的方法:

f =x^2+exp(x^2)+sin(x)*cos(2*x+1)

1、diff(S,\’v\’):对符号对象S 中指定的符号变量v,求其1阶导数。

2、diff(S):对符号对象S 中的默认的独立变量求其1阶导数。

3、diff(S,n):对符号对象S中的默认的独立变量求其n阶导数。

4、diff(S,\’v\’,n):对符号对象S 中指定的符号变量v求其n阶导数。

即:

syms x n

y=sin(x)^n*cos(n*x);

Xd=diff(y)

Xd =

sin(x)^n*n*cos(x)/sin(x)*cos(n*x)-sin(x)^n*sin(n*x)*n

Nd=diff(y, n)

Nd =

sin(x)^n*log(sin(x))*cos(n*x)-sin(x)^n*sin(n*x)*x

而微分值则要按照上述方法,然后调用函数subs()完成所求值就能求得了。

MATLAB是一款美国MathWorks公司出品的商业数学软件。MATLAB是矩阵实验室(Matrix Laboratory)的简称。MATLAB用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境,主要包括MATLAB和Simulink两大部分。MATLAB目前只有英语和日语版,不提供中文版。

0

评论0