24小时热门版块排行榜    

查看: 1769  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 21/1050 2026-08-24 14:21 by refreshing11
[基金申请] 能否退出参与的面上项目解除限项 +15 koalala 2026-08-24 17/850 2026-08-24 14:16 by hhs666
[基金申请] 我面上完蛋了 +8 且听虎啸 2026-08-20 9/450 2026-08-24 13:52 by zzuzxg
[基金申请] 2026国自然函评费到账 +17 羊腰板 2026-08-21 18/900 2026-08-24 12:03 by iaeyuan
[基金申请] 建议基金发布提前给出明确的时间点 +12 kulium 2026-08-21 15/750 2026-08-24 12:02 by iaeyuan
[基金申请] 范进中举一文的中心思想 +6 炎黄贵胄 2026-08-22 7/350 2026-08-24 11:58 by 6543yes
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 17/850 2026-08-24 10:20 by echo8914667
[基金申请] 什么时候开奖? +10 CrisMessi 2026-08-18 11/550 2026-08-24 06:50 by 开心的小狮子
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 今天基金会出结果吗?20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见