24小时热门版块排行榜

>论坛更新日志 (3753)
>虫友互识 (401)
>文献求助 (385)
>导师招生 (268)
>考博 (133)
>硕博家园 (107)
>休闲灌水 (101)
>博后之家 (98)
>论文投稿 (85)
>基金申请 (78)
>公派出国 (74)
>找工作 (46)
>教师之家 (45)
>考研 (39)
>有机交流 (31)
>招聘信息布告栏 (26)

返回列表

xinyuan_x

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 129.5
散金: 2
红花: 1
帖子: 61
在线: 7小时
虫号: 760492
注册: 2009-04-30
专业: 近场光学

[交流] 【求助】急求：一个求数值偏导数的错误

各位虫友：
在数值分析偏导数时，无法获得结果，程序代码：
syms x y
z=x.*y.*sin(x.*y);
f1=diff(z,x);
f2=diff(z,y)
[x,y]=meshgrid(-500:5:500);
i=f1;
mesh(x,y,i);
xlabel('x(nm)');
ylabel('y(nm)');
zlabel('一阶偏导数fx');
结果出错，原因为f1 =y*sin(x*y)+x*y^2*cos(x*y)是一个表达式，而不是数值，如何获得f1的数值？

回复此楼

» 猜你喜欢

1楼 2010-09-11 21:16:55

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

lijinfeng042

木虫 (小有名气)

Matlab

仿真EPI: 2
应助: 1 (幼儿园)
金币: 2156.1
散金: 115
帖子: 291
在线: 31.5小时
虫号: 1019062
注册: 2010-05-15
性别: GG
专业: 通信理论与系统

xinyuan_x(金币+2):我用的是subs，还是感谢你。 2010-09-12 10:50:26

引用回帖:

Originally posted by xinyuan_x at 2010-09-11 21:16:55:
各位虫友：
在数值分析偏导数时，无法获得结果，程序代码：
syms x y
z=x.*y.*sin(x.*y);
f1=diff(z,x);
f2=diff(z,y)
[x,y]=meshgrid(-500:5:500);
i=f1;
mesh(x,y,i);
xlabel('x(nm)');
ylabel(' ...

要转换为数值用eval即可根据我的理解,你是使用上一步的x,y值吧? 这样该
i=eval(f1);
这是你要的图?

赞一下

回复此楼

工作了,偶尔会上来~可以关注新浪微博 @云是风的梦_Matlab

2楼2010-09-12 09:05:39

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

相关版块跳转我要订阅楼主 xinyuan_x 的主题更新

返回列表