24小时热门版块排行榜    

查看: 1733  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:27 by oEVWOejN9taj
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 5/250 2026-08-08 17:02 by oEVWOejN9taj
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 qTvzQBRAHjCi 2026-08-07 5/250 2026-08-08 16:47 by oEVWOejN9taj
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 16:47 by oEVWOejN9taj
[考博] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:27 by oEVWOejN9taj
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:19 by oEVWOejN9taj
[博后之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:02 by oEVWOejN9taj
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 5/250 2026-08-08 15:59 by oEVWOejN9taj
[论文投稿] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 6/300 2026-08-08 15:07 by oEVWOejN9taj
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 WQyTGMbfH7Cx 2026-08-07 4/200 2026-08-08 15:02 by oEVWOejN9taj
[基金申请] 这样的filecode谁见过 +5 布布和一二 2026-08-08 8/400 2026-08-08 14:54 by archvillain
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:42 by oEVWOejN9taj
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 HEQlVqMTIA7d 2026-08-07 5/250 2026-08-08 14:22 by oEVWOejN9taj
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 KXLV3nuBVBY7 2026-08-07 6/300 2026-08-08 14:01 by oEVWOejN9taj
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] filecode与中标关系的预测 +3 布布和一二 2026-08-07 3/150 2026-08-07 15:09 by gltch
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
信息提示
请填处理意见