24小时热门版块排行榜    

查看: 1069  |  回复: 7

polarist

木虫 (著名写手)


[交流] 请Matlab高手帮我看看这个数值微分程序,为毛算不出啊!!!

clear all;
x_num=2;
f=sym('log(besselk(x,y))');
dfdx=diff(f,'x',2)%对f求x偏导
dfdx_num=subs(dfdx,'x',x_num)

感激不尽!!
回复此楼

» 猜你喜欢

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

» 抢金币啦!回帖就可以得到:

查看全部散金贴

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

kezhou1980

铜虫 (正式写手)



小木虫: 金币+0.5, 给个红包,谢谢回帖
这个真不懂,对不住啊,兄弟
2楼2012-08-06 21:07:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dlcheng

金虫 (职业作家)



小木虫: 金币+0.5, 给个红包,谢谢回帖
这个可不懂,对不住啊,兄弟
3楼2012-08-06 21:11:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangming_g

铜虫 (初入文坛)


★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
polarist: 金币+10 2012-08-06 22:48:29
我在matlab2011b上试了一下,好像有结果呀
dfdx =

diff(besselk(x, y), x, x)/besselk(x, y) - diff(besselk(x, y), x)^2/besselk(x, y)^2


dfdx_num =

(besselk(1, y)/y + (2*besselk(2, y))/y^2 + (2*((2*besselk(2, y))/y + besselk(1, y)))/y + besselk(0, y))/besselk(2, y) - (besselk(1, y)/y + (2*besselk(2, y))/y^2 + (2*((2*besselk(2, y))/y + besselk(1, y)))/y + besselk(0, y))^2/besselk(2, y)^2
不知道是不是你想要的?

» 本帖已获得的红花(最新10朵)

4楼2012-08-06 21:16:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

polarist

木虫 (著名写手)


送鲜花一朵
引用回帖:
4楼: Originally posted by yangming_g at 2012-08-06 21:16:56
我在matlab2011b上试了一下,好像有结果呀
dfdx =

diff(besselk(x, y), x, x)/besselk(x, y) - diff(besselk(x, y), x)^2/besselk(x, y)^2


dfdx_num =

(besselk(1, y)/y + (2*besselk(2, y))/y^2 + ...

我的是matlab7
谢谢啊
我下载个2011b测试下!
6楼2012-08-06 22:48:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ytinghost3

铁杆木虫 (著名写手)



小木虫: 金币+0.5, 给个红包,谢谢回帖
dfdx =

diff(besselk(x, y), x, x)/besselk(x, y) - diff(besselk(x, y), x)^2/besselk(x, y)^2


dfdx_num =

diff(besselk(2, y), 2, 2)/besselk(2, y) - diff(besselk(2, y), 2)^2/besselk(2, y)^2
7楼2012-08-06 22:59:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

polarist

木虫 (著名写手)


楼上的运行下这个看有结果吗? 谢谢!
clear all;
x_num=2;
y_num=2;
f=sym('log(besselk(x,y))');
dfdx=diff(f,'x',2);%对f求x偏导
dfdx_num=subs(dfdx,{x,y},{x_num,y_num})
8楼2012-08-06 23:30:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
yfliah5楼
2012-08-06 21:27   回复  
相关版块跳转 我要订阅楼主 polarist 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 焦虑 +7 水冰月月野兔 2026-03-13 9/450 2026-03-16 10:00 by Quakerbird
[考研] 化学调剂0703 +7 啊我我的 2026-03-11 7/350 2026-03-15 23:03 by 凌千颂111
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
[考研] 080500,材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
[基金申请] 现在如何回避去年的某一个专家,不知道名字 +3 zk200107 2026-03-12 6/300 2026-03-14 17:13 by zk200107
[考研] 【0703化学调剂】-一志愿华中师范大学-六级475 +5 Becho359 2026-03-11 5/250 2026-03-14 11:35 by 哦哦123
[考研] 一志愿天津大学,英一数二305分求调剂,四六级已过 +8 小小番的茄 2026-03-09 8/400 2026-03-14 01:53 by JourneyLucky
[考研] 085600求调剂 +3 a邵星池 2026-03-09 3/150 2026-03-14 01:32 by JourneyLucky
[考研] 0703求调剂 +7 jtyq001 2026-03-10 7/350 2026-03-14 01:06 by JourneyLucky
[考研] b区环境工程求调剂 +4 Maps1 2026-03-10 6/300 2026-03-14 00:23 by JourneyLucky
[考研] 311求调剂 +8 zchqwer 2026-03-10 8/400 2026-03-14 00:01 by JourneyLucky
[考研] 279求调剂 +3 Dizzy123@ 2026-03-10 3/150 2026-03-13 23:02 by JourneyLucky
[考研] 考研调剂 +4 芬达46 2026-03-12 4/200 2026-03-13 16:04 by ruiyingmiao
[考研] 材料专硕350 求调剂 +4 王金科 2026-03-12 4/200 2026-03-13 16:02 by ruiyingmiao
[考研] 314求调剂 +7 无懈可击的巨人 2026-03-12 7/350 2026-03-13 15:40 by JourneyLucky
[考研] 274求调剂 +3 S.H1 2026-03-12 3/150 2026-03-13 15:15 by JourneyLucky
[考研] 328化工专硕求调剂 +4 。,。,。,。i 2026-03-12 4/200 2026-03-13 14:44 by JourneyLucky
[考研] 289求调剂 +3 李政莹 2026-03-12 3/150 2026-03-13 11:02 by 求调剂zz
[考研] 268求调剂 +4 好运连绵不绝 2026-03-12 4/200 2026-03-13 10:45 by hyswxzs
[考研] 一志愿河海大学085900土木水利专硕279求调剂不挑专业 +4 SunWwWwWw 2026-03-10 8/400 2026-03-13 02:23 by SunWwWwWw
信息提示
请填处理意见