24小时热门版块排行榜    

查看: 1758  |  回复: 11

apple880518

新虫 (小有名气)

[求助] 求拟合一组数据

请大家帮帮忙拟合一下,我要得到Ca-t曲线各点处切线的斜率,到底用origin还是matlab拟合比较好?谢谢!
t        Ca
10        0.34
30        0.29
60        0.25
120        0.17
180        0.09
240        0.03
回复此楼

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

apple880518

新虫 (小有名气)

用matlab的话是不是需要拟合公式的形式?公式为Ca=k1[1-exp(-k2t)]
2楼2013-05-31 16:17:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lauysche

禁虫 (著名写手)

★ ★
感谢参与,应助指数 +1
jjdg: 编辑内容 2013-06-02 23:08
apple880518: 金币+2 2013-08-27 09:28:51
本帖内容被屏蔽

3楼2013-05-31 16:27:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

apple880518

新虫 (小有名气)

引用回帖:
3楼: Originally posted by lauysche at 2013-05-31 16:27:53
看你对哪个更加熟悉了,这组数据不错,很有规律性,猜测为一元回归:
Y=a*X+b
采用最小二号才呢过回归,代码:
xx =

   10.0000    0.3400
   30.0000    0.2900
   60.0000    0.2500
  120.0000    0.17 ...

这是线性的吧?我想要曲线的,因为其他的数据就是曲线的
这组数据
t               ca
0        0.0871
10        0.0402
20        0.0335
30        0.0315
40        0.0293
4楼2013-05-31 17:24:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangzhiguang

至尊木虫 (文坛精英)

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
xzhdty: 金币+1, 谢谢参与 2013-06-02 19:58:45
apple880518: 金币+2 2013-08-27 09:28:55
使用拟合工具箱获得的结果:
General model Exp1:
     f(x) = a*exp(b*x)
Coefficients (with 95% confidence bounds):
       a =      0.3739  (0.3221, 0.4257)
       b =   -0.007664  (-0.009944, -0.005385)

Goodness of fit:
  SSE: 0.001605
  R-square: 0.9777
  Adjusted R-square: 0.9721
  RMSE: 0.02003
好好学习,天天向上!
5楼2013-05-31 21:16:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangzhiguang

至尊木虫 (文坛精英)

【答案】应助回帖

自己通过编程,使用你给的那个方程始终出错,当然楼主也可以自己编一下,使用曲线拟合函数lsqcurvefit或者nlinfit
好好学习,天天向上!
6楼2013-05-31 21:19:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

apple880518

新虫 (小有名气)

引用回帖:
6楼: Originally posted by zhangzhiguang at 2013-05-31 21:19:29
自己通过编程,使用你给的那个方程始终出错,当然楼主也可以自己编一下,使用曲线拟合函数lsqcurvefit或者nlinfit

恩,谢谢,我用过lsqcurvefit,总运行出错
7楼2013-06-01 09:12:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

忆如沉

银虫 (小有名气)

【答案】应助回帖


感谢参与,应助指数 +1
apple880518(xzhdty代发): 金币+1, 谢谢参与 2013-06-02 19:58:16
给你一个链接,很好的一个拟合
已知一个曲线的方程,如何才能精确拟合其系数?程序代码如下:
%生成待回归信号-------------------------------
%生成总信号
clear;clc;
t = 1:10000;
factor_aci(1)=2.35;factor_aci(2) = 10;factor_aci(3) = 780e-6; factor_aci(4)=33;factor_aci(5) = 0;
data = myfun(factor_aci,t);
%取一段进行回归
tt=2000:4000;
ydata = data(tt);
xdata = 1:length(ydata);
%------------------------------------------
%回归系数项
factor0 = [1,1,0.0001,1,500];
Lb =[0,0,0,0,500];
Ub =[5,30,0.01,50,4000];
[a,residual] = lsqcurvefit(@myfun,factor0,xdata,ydata,Lb ,Ub);
%采用回归系数计算回归曲线
regression = myfun(a,xdata);
figure(1);
plot(regression,'b');
hold on;
plot(ydata,'r--');
legend('拟合','实测');
你若不离不弃,我便生死相依。
8楼2013-06-02 18:55:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

apple880518

新虫 (小有名气)

引用回帖:
8楼: Originally posted by 忆如沉 at 2013-06-02 18:55:33
给你一个链接,很好的一个拟合
已知一个曲线的方程,如何才能精确拟合其系数?程序代码如下:
%生成待回归信号-------------------------------
%生成总信号
clear;clc;
t = 1:10000;
factor_aci(1)=2.35;fac ...

谢谢!
9楼2013-06-02 19:41:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangzhiguang

至尊木虫 (文坛精英)

这个拟合使用matlab拟合最麻烦之处在于系数初值的计算
好好学习,天天向上!
10楼2013-06-03 08:21:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 apple880518 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +36 医学老男孩 2026-08-13 85/4250 2026-08-17 21:23 by lcmdlut
[基金申请] 感觉是下周放榜了 +5 angus9576 2026-08-17 9/450 2026-08-17 20:57 by ssxclkj
[基金申请] filecode=后面第一个是大写字母 +8 wangze12014 2026-08-14 10/500 2026-08-17 17:05 by xter9665
[基金申请] 时间戳变了,能看出什么问题? +4 基诺咪客 2026-08-17 4/200 2026-08-17 16:10 by Vivilian
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +12 Ldrop2023 2026-08-13 15/750 2026-08-17 15:02 by 小豌豆_发芽
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[基金申请] 时间戳又变了8-15 +13 archvillain 2026-08-15 25/1250 2026-08-16 20:13 by zhaosm1982
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 6/300 2026-08-16 20:07 by 奶牛小黑
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +3 Tide man 2026-08-14 3/150 2026-08-16 17:47 by jurkat.1640
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +10 majunge000 2026-08-11 12/600 2026-08-16 08:18 by xli1984
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
信息提示
请填处理意见