24小时热门版块排行榜    

Znn3bq.jpeg
查看: 3516  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 294求调剂 +6 淡然654321 2026-04-17 7/350 2026-04-17 13:32 by 学员JpLReM
[有机交流] 二苯甲酮酸类衍生物 50+3 小白爱主人 2026-04-17 4/200 2026-04-17 11:18 by czyzsu
[考研] 一志愿中科大材料与化工,353分还有调剂学校吗 +9 否极泰来2026 2026-04-15 11/550 2026-04-17 11:08 by liuzhexuan12
[考研] 335求调剂 +20 想上岸呀!! 2026-04-12 23/1150 2026-04-17 10:50 by cuisz
[考研] 304求调剂 +6 castLight 2026-04-16 6/300 2026-04-17 10:34 by 一切OK
[考研] 294求调剂 +14 淡然654321 2026-04-15 14/700 2026-04-16 21:01 by lpl364211
[考研] 297,工科调剂?河南农业大学本科 +14 河南农业大学-能 2026-04-14 14/700 2026-04-16 14:41 by dingyanbo1
[考博] 申博自荐 +3 Linxia林夏 2026-04-13 3/150 2026-04-16 12:55 by 墨荷之露
[考研] 279学硕食品专业求调剂院校 20+7 孤独的狼爱吃羊 2026-04-12 29/1450 2026-04-16 09:00 by screening
[考研] 复试调剂 +21 积极向上; 2026-04-10 23/1150 2026-04-15 12:50 by 西北望—风沙
[考研] 化工学硕294分,求导师收留 +32 yzyzx 2026-04-12 36/1800 2026-04-14 17:45 by lhj2009
[考研] 各位老师好,求调剂,本科211,一志愿天津大学生物与医药学硕,差两名录取。 +11 路六六jjj 2026-04-13 11/550 2026-04-14 16:01 by zs92450
[考研] 求调剂 +12 何气正 2026-04-13 13/650 2026-04-14 14:47 by zs92450
[考研] 食品与营养(0955)271求调剂 +15 升格阿达 2026-04-12 16/800 2026-04-14 13:18 by 浮若_安生
[考研] 一志愿中南大学 0855 机械 286 求调剂 +11 不会吃肉 2026-04-12 11/550 2026-04-13 21:59 by bljnqdcc
[考研] 求调剂 +3 胃痉挛累了 2026-04-11 5/250 2026-04-11 14:13 by luhong1990
[考研] 0854调剂 +5 音像店听花鼓戏 2026-04-10 5/250 2026-04-11 10:49 by qingpingzhu
[考研] 一志愿985机械学硕380求调剂 +5 关关雎鸠10 2026-04-11 5/250 2026-04-11 10:10 by 知念。A
[考研] 263能源动力专硕求调剂 +3 加大号饭盒袋 2026-04-10 3/150 2026-04-10 22:23 by 286640313
[考研] 368求调剂 +3 17385968172 2026-04-10 3/150 2026-04-10 20:12 by 电子奥德彪
信息提示
请填处理意见