应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 程序代码 » 怎么使用oringin或spss自定义函数实现非线性拟合y=b+A*cos(c*t−d)
61/1返回列表
查看: 1451  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

钟声入翠微

铜虫 (小有名气)

[求助] 怎么使用oringin或spss自定义函数实现非线性拟合y=b+A*cos(c*t−d) 已有1人参与

请问怎么使用oringin或spss自定义函数实现非线性拟合
公式：y=b+A*cos(c*t−d)，其中c=2π/365

数据拟合，根据残差分析剔除绝对值大于2的数据点，然后重新拟合，直到无异常点
本人最近一直在尝试，但是编程基础太差，自定义函数总是不成功，希望大家可以帮帮忙。
数据
1987-05-15        -3.64
1987-06-15        -8.6
1987-07-15        -3.91
1987-08-15        -5.53
1987-09-15        -8.64
1987-10-15        -13.88
1987-11-15        -4.63
1987-12-15       
1988-01-15        -21.82
1988-02-15        -13.01
1988-03-15        -15.37
1988-04-15       
1988-05-15        -8.22
1988-06-15        -7.36
1988-07-15        -10.09
1988-08-15        -9.67
1988-09-15        -6.98
1988-10-15        -7.97
1988-11-15       
1988-12-15       
1989-01-15       
1989-02-15        -6.36
1989-03-15        -5.3
1989-04-15        -4.16
1989-05-15        -2.59
1989-06-15        -5.99
1989-07-15        -7.58
1989-08-15        -9.1
1989-09-15        -7.61
1989-10-15        -8.09
1989-11-15        -10.9
1989-12-15       
1990-01-15        -10.96
1990-02-15        -12.61
1990-03-15        -14.3
1990-04-15        -6.85
1990-05-15        0.73
1990-06-15        -5.94
1990-07-15        -7.67
1990-08-15        -6.81
1990-09-15        -4.96
1990-10-15        -9.61
1990-11-15        -13.61
1990-12-15       
1991-01-15       
1991-02-15       
1991-03-15        -11.79
1991-04-15        -7.43
1991-05-15        -4.79
1991-06-15        -4.16
1991-07-15        -5.86
1991-08-15        -7.83
1991-09-15        -5.3
1991-10-15        -4.8@月只蓝
回复此楼

» 猜你喜欢
1楼 2020-06-08 18:49:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)
你这数据输入格式就不对，t应该是时间序列如1.2.3... 而不是具体年月。
剔除的部分对应时间也不要输入了
发自小木虫Android客户端
赞一下(1人) 回复此楼
数值计算
2楼2020-06-08 19:53:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

自为胜利而来

至尊木虫 (著名写手)
百度一下吧，origin自定义拟合函数。照着做一遍

发自小木虫Android客户端
赞一下 回复此楼
道系青年
3楼2020-06-08 19:56:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

钟声入翠微

铜虫 (小有名气)
引用回帖:
2楼: Originally posted by 独孤神宇 at 2020-06-08 19:53:45
你这数据输入格式就不对，t应该是时间序列如1.2.3... 而不是具体年月。
剔除的部分对应时间也不要输入了

我给出的是原始的下载下来的数据，我知道t是时间序列，但是我主要是对于自定义函数编程方面无能为力。还需要大佬帮帮忙，谢谢。

发自小木虫Android客户端
赞一下 回复此楼
4楼2020-06-08 20:38:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hzlhm

至尊木虫 (著名写手)

【答案】应助回帖

感谢参与，应助指数 +1
可以用插值的方法补齐空缺的数据，再用拟合的方法求出方程的系数。或先拟合系数，再补全空缺的数据。

发自小木虫Android客户端
赞一下 回复此楼
QQ：2120156492
5楼2020-06-10 07:00:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

钟声入翠微

铜虫 (小有名气)
引用回帖:
5楼: Originally posted by hzlhm at 2020-06-10 07:00:22
可以用插值的方法补齐空缺的数据，再用拟合的方法求出方程的系数。或先拟合系数，再补全空缺的数据。

请问如何拟合系数呢？

发自小木虫Android客户端
赞一下 回复此楼
6楼2020-06-10 09:03:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 钟声入翠微 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 17:01 by BCdvBAEpt3zm
[博后之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 16:58 by BCdvBAEpt3zm
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 18/900 2026-11-01 16:48 by BCdvBAEpt3zm
[基金申请] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +6 X5bHxOULxFqb 2026-10-30 22/1100 2026-11-01 16:35 by BCdvBAEpt3zm
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +5 X5bHxOULxFqb 2026-10-30 22/1100 2026-11-01 16:35 by BCdvBAEpt3zm
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 hKe7U1kad5Sg 2026-10-31 7/350 2026-11-01 15:08 by BCdvBAEpt3zm
[基金申请] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 7/350 2026-11-01 15:04 by BCdvBAEpt3zm
[考研] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 xvJS1BYDxvBS 2026-10-31 10/500 2026-11-01 14:50 by BCdvBAEpt3zm
[论文投稿] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 xvJS1BYDxvBS 2026-10-31 7/350 2026-11-01 14:50 by BCdvBAEpt3zm
[找工作] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 14:48 by BCdvBAEpt3zm
[硕博家园] 越来越担心硕士毕业后没有竞争力了 +9 zuesli 2026-10-27 15/750 2026-11-01 12:59 by Equinoxhua
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 3/150 2026-11-01 12:53 by BCdvBAEpt3zm
[考博] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 6/300 2026-11-01 12:53 by BCdvBAEpt3zm
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 xvJS1BYDxvBS 2026-10-31 5/250 2026-11-01 12:51 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 X5bHxOULxFqb 2026-10-30 18/900 2026-11-01 12:11 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 xvJS1BYDxvBS 2026-10-31 9/450 2026-11-01 10:53 by BCdvBAEpt3zm
[基金申请] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 14/700 2026-11-01 10:43 by BCdvBAEpt3zm
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-01 10:27 by BCdvBAEpt3zm
[博后之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-01 05:22 by YXM1FDCSE8Bn
[考博] 2027博士申请-临床检验诊断学、病原生物学、微生物方向 50+3 yyh@ 2026-10-28 4/200 2026-10-31 12:47 by 可淡不可忘
信息提示
关闭
请填处理意见
关闭