版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

yg5189

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 2178.9
帖子: 203
在线: 324.8小时
虫号: 1649091
注册: 2012-02-28
专业: 运筹学

[求助] matlab符号函数求导及梯度的问题

比如现有函数
function f=fun(x)
f=2*x(1)*x(1)-x(2)-1;

1.求f对x(1)的导数在[1,2]这点的值
2.求f对x的梯度（也就是分别对x(1),x(2)求导的一个向量）在[1,2]这点的向量（如果是重复使用1中的函数的话那这题可以不用写）

*注意：不是将x(1),x(2)变为x1,x2，因为我实际求的函数有n维，我想用for循环！谢谢

回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

求matlab画一个函数的图已经有6人回复
MATLAB函数速查手册已经有211人回复
关于Matlab算符号积分已经有6人回复
matlab的函数调用问题已经有3人回复
MATLAB中产生高斯白噪声，涉及到awgn和wgn函数已经有3人回复
已知一随机信号的功率谱密度函数，如何利用Matlab产生其单次时域样本已经有6人回复
问个关于matlab里find函数的问题！已经有4人回复
关于数学期望与取符号函数的运算次序已经有8人回复
【求助】用matlab对多项式函数求导已经有3人回复
matlab做函数图象的问题已经有6人回复
【求助】请教一个隐函数的求导问题已经有10人回复

1楼2012-06-30 21:48:15

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

yg5189

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 2178.9
帖子: 203
在线: 324.8小时
虫号: 1649091
注册: 2012-02-28
专业: 运筹学

引用回帖:

2楼: Originally posted by libralibra at 2012-06-30 23:10:24
matlab求导命令diff调用格式:

diff(函数) ，求的一阶导数；

diff(函数， n) ，求的n阶导数(n是具体整数)；

diff(函数，变量名)，求对的偏导数；

谢谢你的回复！
不过我试过diff，好像不行，是不是我用的不对，能不能发一个你写的例子？

赞一下

回复此楼

3楼2012-07-01 08:57:25

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 6 个回答

libralibra

至尊木虫 (著名写手)

骠骑将军

程序强帖: 40
应助: 817 (博后)
金币: 12914.1
红花: 64
帖子: 2238
在线: 287.3小时
虫号: 696514
注册: 2009-02-05
专业: 计算机软件

【答案】应助回帖

感谢参与，应助指数 +1

matlab求导命令diff调用格式:

diff(函数) ，求的一阶导数；

diff(函数， n) ，求的n阶导数(n是具体整数)；

diff(函数，变量名)，求对的偏导数；

赞一下

回复此楼

matlab/VB/python/c++/Java写程序请发QQ邮件:790404545@qq.com

2楼2012-06-30 23:10:24

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

libralibra

至尊木虫 (著名写手)

骠骑将军

程序强帖: 40
应助: 817 (博后)
金币: 12914.1
红花: 64
帖子: 2238
在线: 287.3小时
虫号: 696514
注册: 2009-02-05
专业: 计算机软件

【答案】应助回帖

你的问题貌似没啥其他办法

CODE:

 clc;

syms x y;

f = 2*x*x-y-1;

dx = diff(f,x);

dy = diff(f,y);

disp('df/dx = ');

disp(dx);

disp('df/dy = ');

disp(dy);

grad = [dx,dy]; 

dx12 = subs(dx,x,1);

dy12 = subs(dy,y,2);

disp(' ');

disp(['df/dx @ [1,2] = ',num2str(dx12)]);

disp(['gradient(f,x,y) @ [1,2] = ',num2str(dy12)]);

结果:

CODE:

df/dx = 

4*x

df/dy = 

-1

df/dx @ [1,2] = 4

gradient(f,x,y) @ [1,2] = -1

>>

赞一下

回复此楼

matlab/VB/python/c++/Java写程序请发QQ邮件:790404545@qq.com

4楼2012-07-01 18:29:19

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

libralibra

至尊木虫 (著名写手)

骠骑将军

程序强帖: 40
应助: 817 (博后)
金币: 12914.1
红花: 64
帖子: 2238
在线: 287.3小时
虫号: 696514
注册: 2009-02-05
专业: 计算机软件

【答案】应助回帖

★ ★ ★ ★ ★ ★
yg5189: 金币+6, ★★★很有帮助 2012-07-02 09:35:51

抱歉,最后一句写错了

CODE:

 grad12 = subs(grad,{x,y},{1,2});

disp(['gradient(f,x,y) @ [1,2] = [',num2str(grad12),']']);

结果

CODE:

 df/dx = 

4*x

df/dy = 

-1

df/dx @ [1,2] = 4

gradient(f,x,y) @ [1,2] = [4 -1]

赞一下

回复此楼

matlab/VB/python/c++/Java写程序请发QQ邮件:790404545@qq.com

5楼2012-07-01 18:32:00

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 6 个回答