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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 【求助】二元函数拟合
131/2返回列表12下一页
查看: 2404  |  回复: 12
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

xiexuewang

金虫 (小有名气)

[交流] 【求助】二元函数拟合 已有6人参与

t,c0,c每组35个数据
t=[0 24 48 60 72 84 96 0 24 48 60 72 84 96 0 24 48 60 72 84 96 0 24 48 60 72 84 96 0 24 48 60 72 84 96];
c0=[25 25 25 25 25 25 25 50 50 50 50 50 50 50 75 75 75 75 75 75 75 100 100 100 100 100 100 100 150 150 150 150 150 150 150];
c=[25 24.6 21.9 14.3 9 6 4 50 46 40 26 22 17 14 75 73 69.7 63.6 50.9 31 26.3 100 97.9 92.4 83.3 58.3 41.6 30 150 148 142 120.7 109.5 86.8 54];

拟合函数形式：c=(A(1)*c0+A(2))*(A(3)*exp(A(4)*t)+A(5)*exp(A(6)*t)
求解A1-A6的值

[ Last edited by xiexuewang on 2010-5-13 at 14:05 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
没有不劳而获，任何成就都是汗水浇灌的结果
1楼 2010-05-12 16:09:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaojun376

金虫 (小有名气)
★
zzuwangshilei(金币+1):帮助他人解决问题 2010-05-12 19:25:27
你这好像不是拟合，是参数估计，你是想求A的值是吧，那就用lsqnonlin寻找最优值
赞一下(1人) 回复此楼
在交流中共同进步！
2楼2010-05-12 18:05:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiexuewang

金虫 (小有名气)
是啊，我是想求A的值，就是不会用lsqnonlin函数啊，
赞一下 回复此楼
没有不劳而获，任何成就都是汗水浇灌的结果
3楼2010-05-12 21:04:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiexuewang

金虫 (小有名气)
zzuwangshilei:具体的步骤，一般要自己做的，别人给出的是建议或思路 2010-05-13 15:30:45
用matlab求出A1-A6的值，给出具体的步骤，先谢了
赞一下(1人) 回复此楼
没有不劳而获，任何成就都是汗水浇灌的结果
4楼2010-05-12 21:08:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rewind

至尊木虫 (正式写手)
★
robert2020(金币+1):谢谢提供建议，欢迎常来仿真模拟版交流。 2010-05-13 07:51:44
引用回帖:
Originally posted by xiexuewang at 2010-05-12 16:09:46:
t,c0,c每组35个数据
t=[0        24        48        60        72        84        96 0        24        48        60        72        84        96 0        24        48        60        72        84        96 0        24        48        60        72        84        96 0        24        48        60        72        84        96];
c0=[25        25        25        25        25 25 25 50        50        50        50        50 50 50 75        75        75        75        75 75 75 100         ...

您的C值是什么？要是二元函数拟合的话 结果矩阵应该是35*35的吧
赞一下(1人) 回复此楼
在外种地...
5楼2010-05-12 21:14:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiexuewang

金虫 (小有名气)
zzuwangshilei:鼓励继续讨论问题，欢迎常来仿真模拟版交流 2010-05-13 09:35:17
c=f(c0,t), c是和c0，t相关的数啊，应该是c=function==（A(1)*c0+A(2))*(A(3)*exp(A(4)*t)+A(5)*exp(A(6)*t)
赞一下(1人) 回复此楼
没有不劳而获，任何成就都是汗水浇灌的结果
6楼2010-05-12 22:09:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rewind

至尊木虫 (正式写手)
★ ★
robert2020(金币+2):鼓励交流经验~~ 2010-05-13 12:17:31
zzuwangshilei:积极参与讨论，欢迎常来 2010-05-13 15:31:32
引用回帖:
Originally posted by xiexuewang at 2010-05-12 22:09:23:
c=f(c0,t), c是和c0，t相关的数啊，应该是c=function==（A(1)*c0+A(2))*(A(3)*exp(A(4)*t)+A(5)*exp(A(6)*t)

要是那样的话 就不是拟合问题 可以使用非线性回归求解

matlab里面有个自带例子

调用命令为：

load reaction

beta = nlinfit(reactants,rate,@hougen,beta)

照做改下就可以了 祝一切顺利
赞一下(2人) 回复此楼
在外种地...
7楼2010-05-13 10:20:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

change0618

铁杆木虫 (著名写手)

方丈大师
★
robert2020(金币+1):谢谢提供意见，欢迎常来仿真模拟版交流。 2010-05-13 22:36:40
c=(A(1)*c0+A(2))*(A(3)*exp(A(4)*t)+A(5)*exp(A(6)*t)
这个式子不正确!!
如果是c=(A(1)*c0+A(2))*(A(3)*exp(A(4)*t)+A(5)*exp(A(6)*t)  )
也会有无数个解,因为可以提取因式,直接影响到A的六个数值的大小.
赞一下(1人) 回复此楼
8楼2010-05-13 19:42:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiexuewang

金虫 (小有名气)
zzuwangshilei:积极参与讨论，欢迎常来 2010-05-14 08:07:04
只要满足数据的条件就可以了，最重要的是给出满足条件的函数就行，
赞一下(1人) 回复此楼
没有不劳而获，任何成就都是汗水浇灌的结果
9楼2010-05-13 21:01:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hopingzou

金虫 (著名写手)
★
zzuwangshilei(金币+1):提供问题的解决结果，欢迎参与讨论 2010-05-14 16:04:20
抽空帮你做了一下，仅供参考。觉得这个拟合模型未必是最好的！
赞一下(1人) 回复此楼
10楼2010-05-14 15:01:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiexuewang 的主题更新
131/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 j5mmowWXNWxv 2026-09-05 6/300 2026-09-06 18:41 by eObJPa4gWmp6
[硕博家园] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 18:21 by eObJPa4gWmp6
[博后之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 18:09 by eObJPa4gWmp6
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 jPDp0sc2B7zQ 2026-09-05 4/200 2026-09-06 18:09 by eObJPa4gWmp6
[论文投稿] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 18:01 by eObJPa4gWmp6
[找工作] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 8/400 2026-09-06 17:49 by eObJPa4gWmp6
[文学芳草园] 两块石头 +7 阿美_Lml888 2026-09-03 9/450 2026-09-06 17:48 by myrtle
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 17:21 by eObJPa4gWmp6
[教师之家] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 9/450 2026-09-06 17:09 by eObJPa4gWmp6
[博后之家] 广西大学-广州大学招聘博士后 欢迎广大优秀人才！！！ +5 SCSIOyxguo 2026-09-03 8/400 2026-09-06 15:06 by 可爱小花朵
[考博] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 13:09 by eObJPa4gWmp6
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 12:49 by eObJPa4gWmp6
[考博] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 12:29 by eObJPa4gWmp6
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 7/350 2026-09-06 11:20 by eObJPa4gWmp6
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 8/400 2026-09-06 11:09 by eObJPa4gWmp6
[论文投稿] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 jPDp0sc2B7zQ 2026-09-04 8/400 2026-09-06 11:00 by eObJPa4gWmp6
[考研] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 4/200 2026-09-06 07:20 by E1iiBLMU2XnD
[文学芳草园] 初秋的晨风 +4 阿美_Lml888 2026-09-04 6/300 2026-09-05 13:51 by 阿美_Lml888
[硕博家园] Ei源刊怎么投 +4 Fengshun9711 2026-09-04 4/200 2026-09-04 19:30 by 研途知予
[基金申请] 要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够 +9 瞬息宇宙 2026-08-31 17/850 2026-09-02 14:04 by anjeeshine
信息提示
关闭
请填处理意见
关闭