【答案】应助回帖

» 猜你喜欢

博士读完未来一定会好吗已经有6人回复
小论文投稿已经有3人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有9人回复
心脉受损已经有3人回复
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有8人回复
申请2026年博士已经有6人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有5人回复

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

1楼2014-12-15 22:20:01

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询

数学EPI: 6
应助: 413 (硕士)
贵宾: 0.002
金币: 7116.5
散金: 15
红花: 63
帖子: 1639
在线: 798.9小时
虫号: 2469437
注册: 2013-05-14
专业: 计算数学与科学工程计算

【答案】应助回帖

感谢参与，应助指数 +1

这恐怕和你数据有关。如果二次项系数确实是正的，那也没有办法。实在要坚持正的二次项系数，那就把数据都反一个符号。

2楼2014-12-15 23:29:03

不倦

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 1663.2
帖子: 89
在线: 214.6小时
虫号: 1626308
注册: 2012-02-19
专业: 天然有机化学

引用回帖:

2楼: Originally posted by pippi6 at 2014-12-15 23:29:03
这恐怕和你数据有关。如果二次项系数确实是正的，那也没有办法。实在要坚持正的二次项系数，那就把数据都反一个符号。

主要是拟合得到的公式还要用于后面求导，如果二次项不为正的话，求导后的单调性不能保证与其他二次项为正的情况不一致。不好说明数据规律。有没有其他软件能做的？

3楼2014-12-15 23:58:14

feixiaolin

荣誉版主 (文坛精英)

专家经验: +518
应助: 942 (博后)
贵宾: 1.275
金币: 2930
散金: 58785
红花: 532
沙发: 11
帖子: 24215
在线: 2601.8小时
虫号: 2139575
注册: 2012-11-21
专业: 光学信息获取与处理
管辖: 数学

可以尝试。
y=k*(x-a)^2+b;
加约束条件 a>0.0035, k>0

4楼2014-12-16 15:12:38

feixiaolin

荣誉版主 (文坛精英)

专家经验: +518
应助: 942 (博后)
贵宾: 1.275
金币: 2930
散金: 58785
红花: 532
沙发: 11
帖子: 24215
在线: 2601.8小时
虫号: 2139575
注册: 2012-11-21
专业: 光学信息获取与处理
管辖: 数学

把数据挂上来，我晚上帮你看一下。

5楼2014-12-16 18:32:44

不倦

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 1663.2
帖子: 89
在线: 214.6小时
虫号: 1626308
注册: 2012-02-19
专业: 天然有机化学

引用回帖:

4楼: Originally posted by feixiaolin at 2014-12-16 15:12:38
可以尝试。
y=k*(x-a)^2+b;
加约束条件 a>0.0035, k>0

就是这个意思，k>0，二次方程开口向上，就是不知道怎么加这个限制条件。

6楼2014-12-16 23:58:45

不倦

木虫 (小有名气)

feixiaolin: 屏蔽内容, 重复发帖 2014-12-17 04:06:14

本帖内容被屏蔽

7楼2014-12-17 00:01:03

不倦

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 1663.2
帖子: 89
在线: 214.6小时
虫号: 1626308
注册: 2012-02-19
专业: 天然有机化学

引用回帖:

5楼: Originally posted by feixiaolin at 2014-12-16 18:32:44
把数据挂上来，我晚上帮你看一下。

0.003469813 2.356978577
0.003410641 2.47182141
0.003353454 2.613079964
0.003298153 2.69577517
0.003244646 2.790551423
0.003192848 2.886531059
数据 X Y

8楼2014-12-17 00:02:03

不倦

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 1663.2
帖子: 89
在线: 214.6小时
虫号: 1626308
注册: 2012-02-19
专业: 天然有机化学

引用回帖:

5楼: Originally posted by feixiaolin at 2014-12-16 18:32:44
把数据挂上来，我晚上帮你看一下。

不好意思，白天一天都没有网络。

9楼2014-12-17 00:02:20

feixiaolin

荣誉版主 (文坛精英)

专家经验: +518
应助: 942 (博后)
贵宾: 1.275
金币: 2930
散金: 58785
红花: 532
沙发: 11
帖子: 24215
在线: 2601.8小时
虫号: 2139575
注册: 2012-11-21
专业: 光学信息获取与处理
管辖: 数学

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
不倦: 金币+5, ★★★很有帮助 2014-12-17 11:20:40

// 1stop  code
Parameters k[0, ], a, b;
Variable x, y;
Function y=k*(x-a)^2+b;
data;
0.003469813 2.356978577
0.003410641 2.47182141
0.003353454 2.613079964
0.003298153 2.69577517
0.003244646 2.790551423
0.003192848 2.886531059
-----------------------------------

迭代数: 23
计算用时(时:分:秒:微秒): 00:00:00:234
优化算法: 麦夸特法(Levenberg-Marquardt) + 通用全局优化法
计算结束原因: 达到收敛判断标准
均方差(RMSE): 0.0119199347122682
残差平方和(SSE): 0.000852509061268412
相关系数(R): 0.997817931981908
相关系数之平方(R^2): 0.995640625384652
决定系数(DC): 0.995640625380019
卡方系数(Chi-Square): 0.000164782518530058
F统计(F-Statistic): 913.562804976051

参数最佳估算
---------- -------------
k       660.053623227762
a       1.44584738126851
b       -1370.8443344474

====== 结果输出 =====
No 实测值y 计算值y
1 2.356978577 2.3662388
2 2.47182141 2.4789101
3 2.613079964 2.5878061
4 2.69577517 2.6931149
5 2.790551423 2.7950112
6 2.886531059 2.8936566

tu.PNG

10楼2014-12-17 05:12:10