24小时热门版块排行榜    

查看: 3429  |  回复: 13

smart824

新虫 (初入文坛)

[求助] 离散点求导数?

这样,我有两组离散的数据,一组是自变量,另一组是应变量,现在我想求他们的导数,我自己已经用差分法编了程序算,但是有几个点总是不太好,所以想请问一下大家能不能给我提供另外一种好的方法,万分感谢!
回复此楼

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

yalefield

金虫 (文坛精英)

老汉一枚

公式有小错啦!
CODE:
y[i] = (22 * (y[i - 3] - y[i+3]) - 67 * (y[i - 2] - y[i+2]) - 58 * (y[i - 1] - y[i + 1])) / 252;

3楼2012-03-15 11:38:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

smart824

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by yalefield at 2012-03-15 11:38:08:
公式有小错啦!

y = (22 * (y - y) - 67 * (y - y) - 58 * (y - y)) / 252;

谢谢啊,我先好好看看你给我说的!
4楼2012-03-15 13:04:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

锐利的碎片

木虫 (正式写手)

star watcher

【答案】应助回帖

感谢参与,应助指数 +1
样条拟合也可以,你看下附件里别人写的程序。

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : diff.rar
  • 2012-03-16 14:11:16, 2.59 K
11楼2012-03-16 14:11:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

yalefield

金虫 (文坛精英)

老汉一枚

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
dbb627: 金币+2, 感谢应助 2012-03-15 14:10:33
建议先插值、平滑。
然后,用Savitzky-Golay求导
【要求自变量是等距离变化的!所以,如果原来的自变量不是均匀的,就要插值咯】。
满足自变量等距的条件后,用Savitzky-Golay多项式(如7点3次)求导【下面的这些y应该是等距的x所对应的y哟】:

一阶导数【注意下标范围】:
CODE:
y[i] = (22 * (y[i - 3] - y[i+3] - 67 * (y[i - 2] - y[i+2]) - 58 * (y[i - 1] - y[i + 1])) / 252;

二阶导数:
您自己推导?还是另外开贴?
2楼2012-03-15 11:37:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

smart824

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by yalefield at 2012-03-15 11:37:15:
建议先插值、平滑。
然后,用Savitzky-Golay求导
【要求自变量是等距离变化的!所以,如果原来的自变量不是均匀的,就要插值咯】。
满足自变量等距的条件后,用Savitzky-Golay多项式(如7点3次)求导【下面的这 ...

请问:什么是Savitzky-Golay多项式?
5楼2012-03-15 14:14:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yalefield

金虫 (文坛精英)

老汉一枚

【答案】应助回帖

★ ★
cenwanglai: 金币+2, ~~ 2012-03-15 18:30:03
http://wenku.baidu.com/view/9181c31f59eef8c75fbfb3e3.html
哦,这是解释用于平滑的。算导数也用这个多项式。
6楼2012-03-15 16:22:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

smart824

新虫 (初入文坛)

引用回帖:
6楼: Originally posted by yalefield at 2012-03-15 16:22:53:
http://wenku.baidu.com/view/9181c31f59eef8c75fbfb3e3.html
哦,这是解释用于平滑的。算导数也用这个多项式。

这篇文章我看到啦,但是讲的挺复杂的,我看你的公式很简单啊
7楼2012-03-15 18:04:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yalefield

金虫 (文坛精英)

老汉一枚


cenwanglai: 金币+1, 随礼~~ 2012-03-15 19:15:09
引用回帖:
7楼: Originally posted by smart824 at 2012-03-15 18:04:03:
这篇文章我看到啦,但是讲的挺复杂的,我看你的公式很简单啊

你是说俺?
飞花摘叶,吐故纳新,
承前启后,治病救人,
简单公式,复杂论文,
人生指南,地方天浑。
8楼2012-03-15 19:12:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

smart824

新虫 (初入文坛)

引用回帖:
8楼: Originally posted by yalefield at 2012-03-15 19:12:31:
你是说俺?
飞花摘叶,吐故纳新,
承前启后,治病救人,
简单公式,复杂论文,
人生指南,地方天浑。

那你的式子中的那些常数之类的是什么意思啊?我不太理解啊
9楼2012-03-15 19:26:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yalefield

金虫 (文坛精英)

老汉一枚

引用回帖:
9楼: Originally posted by smart824 at 2012-03-15 19:26:41:
那你的式子中的那些常数之类的是什么意思啊?我不太理解啊

俺隐退了。
10楼2012-03-15 21:11:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 smart824 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 263求调剂 +6 yqdszhdap- 2026-03-22 9/450 2026-03-23 12:57 by yqdszhdap-
[考研] 接收2026硕士调剂(学硕+专硕) +3 allen-yin 2026-03-23 4/200 2026-03-23 12:27 by allen-yin
[考研] 求调剂材料学硕080500,总分289分 5+3 @taotao 2026-03-19 21/1050 2026-03-23 10:17 by 冠c哥
[考研] 一志愿西安交通大学材料工程专业 282分求调剂 +11 枫桥ZL 2026-03-18 13/650 2026-03-22 20:26 by edmund7
[考研] 求调剂 +5 Zhangbod 2026-03-21 7/350 2026-03-22 13:13 by Zhangbod
[考研] 生物学一志愿985,分数349求调剂 +4 zxts12 2026-03-21 7/350 2026-03-22 09:57 by zxts12
[考研] 一志愿南大,0703化学,分数336,求调剂 +3 收到VS 2026-03-21 3/150 2026-03-21 18:42 by 学员8dgXkO
[基金申请] 学校已经提交到NSFC,还能修改吗? 40+4 babangida 2026-03-19 9/450 2026-03-21 16:12 by babangida
[考研] 265求调剂 +12 梁梁校校 2026-03-19 14/700 2026-03-21 13:38 by lature00
[考研] 265求调剂 +9 梁梁校校 2026-03-17 9/450 2026-03-21 02:17 by JourneyLucky
[考研] 085700资源与环境308求调剂 +12 墨墨漠 2026-03-18 13/650 2026-03-21 01:42 by JourneyLucky
[考研] 材料 336 求调剂 +3 An@. 2026-03-18 4/200 2026-03-21 01:39 by JourneyLucky
[考研] 一志愿重庆大学085700资源与环境专硕,总分308求调剂 +3 墨墨漠 2026-03-18 3/150 2026-03-21 00:39 by JourneyLucky
[考研] 296求调剂 +6 www_q 2026-03-18 10/500 2026-03-20 23:56 by JourneyLucky
[考研] 274求调剂 +10 S.H1 2026-03-18 10/500 2026-03-20 23:51 by JourneyLucky
[考研] 材料与化工 322求调剂 +4 然11 2026-03-19 4/200 2026-03-20 22:12 by luoyongfeng
[考研] 一志愿南理工085701环境302求调剂院校 +3 葵梓卫队 2026-03-20 3/150 2026-03-20 19:28 by zhukairuo
[考研] 298-一志愿中国农业大学-求调剂 +9 手机用户 2026-03-17 9/450 2026-03-20 14:24 by 无懈可击111
[论文投稿] 申请回稿延期一个月,编辑同意了。但系统上的时间没变,给编辑又写邮件了,没回复 10+3 wangf9518 2026-03-17 4/200 2026-03-19 23:55 by babero
[考研] 085600材料与化工求调剂 +6 绪幸与子 2026-03-17 6/300 2026-03-19 13:27 by houyaoxu
信息提示
请填处理意见