| 查看: 1316 | 回复: 7 | |||
[交流]
请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) 感激不尽！！ |
回复此楼» 猜你喜欢
烦躁
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 为毛每次合成出来的羧酸化合物的氢谱羧酸出峰都看不见啊~~~
已经有7人回复
- 2012考研政治终级笔记-马哲+政经+毛概+邓三【还没有复习好的赶紧看！自创无重复】
已经有135人回复
- 【求助/交流】大家帮忙看看这是什么虫卵（对毛毛虫及虫卵感到恶心的朋友慎入）
已经有28人回复
- 【分享】方舟子最新力作《大象为什么不长毛》【已搜索无重复】
已经有25人回复
- 【求助】怎样将图中毛糙的边改成光滑的边呢？
已经有4人回复
- 【求助】关于羊毛用阻燃剂
已经有4人回复
2楼2012-08-06 21:07:48
3楼2012-08-06 21:11:53
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
polarist: 金币+10 2012-08-06 22:48:29
小木虫: 金币+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朵）
polarist4楼2012-08-06 21:16:56
6楼2012-08-06 22:48:16
7楼2012-08-06 22:59:11
8楼2012-08-06 23:30:06
简单回复
yfliah5楼
2012-08-06 21:27
回复
