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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 【求助】二元函数拟合
131/2返回列表12下一页
查看: 2363  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +12 Ldrop2023 2026-08-13 13/650 2026-08-17 12:33 by gltch
[基金申请] 今天系统多次维护，明天很可能放榜！ +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 i7NFEVbjQMM5 2026-08-16 4/200 2026-08-17 12:10 by 4GBAYCdVQoK3
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +32 医学老男孩 2026-08-13 74/3700 2026-08-16 21:05 by 医学老男孩
[论文投稿] 岩土工程学报什么时候才能终审完呐 25+3 yeager111 2026-08-10 4/200 2026-08-16 20:42 by tfang
[基金申请] 时间戳又变了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 奶牛小黑
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[精细化工] 招聘 金属平磨液，抛光液研发工程师 +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 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 重要来源：本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] Filecode 又变了，巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 结合人工智能，周易传统文化，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 撸猫猫
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 什么时候出结果，有咨询渠道？？？ +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
信息提示
关闭
请填处理意见
关闭