看个diff函数小问题

liangwy0106 · 发表于 2009-3-23 10:29

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

function [f,df]=detaf(x)
f=x(1)^4+25*x(2)^4+x(1)^2*x(2)^2;
df(1)=4*x(1)^3+2*x(1)*x(2)^2;
df(2)=100*x(2)^3+2*x(1)^2*x(2);

这个是M文件里面的函数我想用diff函数来分别求x(1)和x(2)的一阶导数
df(1)=diff(f,x(1),1)
df(2)=diff(f,x(2),1)
可是总会有错误:@( 能不能帮我看下! 谢谢！！

[ 本帖最后由 ChaChing 于 2009-3-23 13:20 编辑 ]

ChaChing · 发表于 2009-3-23 11:21

LZ再看看diff的帮助文档, 可以df(1)=diff(f,x(1),1)如此用吗?
看看gradient 是否LZ要的

liangwy0106 · 发表于 2009-3-23 11:44

diff 对于一个变量可以这么用但是我不知道两个变量x(1),x(2)行不行所以在这里问一下

ChaChing · 发表于 2009-3-23 13:41

为何需要用matlab算!?
df(1)对x(1)之微分不就是12*x(1)^2+2*x(2)^2;
df(2)对x(2)之微分不就是300*x(2)^2+2*x(1)^2;

liangwy0106 · 发表于 2009-3-23 14:22

我是想实现这么个功能输入一个函数f=x(1)^4+25*x(2)^4+x(1)^2*x(2)^2;
对x(1)求偏导，得到df(1)=4*x(1)^3+2*x(1)*x(2)^2;
对x(2)求偏导，得到df(2)=100*x(2)^3+2*x(1)^2*x(2);

大家现在明白了吗？有可能我刚才没有说清楚

ChaChing · 发表于 2009-3-23 18:46

help jacobian

还有建議樓主看下本版規則並加強發問題方式!

[ 本帖最后由 ChaChing 于 2009-3-23 18:48 编辑 ]

liangwy0106 · 发表于 2009-3-23 22:18

jacobian好像也不行....

ChaChing · 发表于 2009-3-23 23:43

f=x^4+25*y^4+x^2*y^2; v=[x,y]; jacobian(f,v)

liangwy0106 · 发表于 2009-3-24 09:20

谢谢楼上的这位朋友！
当变量为x,y时，你这个方程可以实现。但是当变量为x(1),x(2)时就会出错
能不能试一下然后告诉我为什么呢？谢谢

ydlcsu · 发表于 2009-3-24 10:49

没见过变量名可以为x(1),x(2)，你换成x1，x2试试吧

ChaChing · 发表于 2009-3-24 11:07

matlab仅是一种工具, 建议楼主配合现有的功能想法来实现

账号		自动登录	找回密码
密码			我要加入

[综合讨论] 看个diff函数小问题

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复 9楼 liangwy0106 的帖子

评分