简单的说就是建立如下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])就无法进行了。
要如何再转换成函数呢？
|
