matlab求导赋值,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等优化函数

展开全文阅读