24小时热门版块排行榜    

查看: 1637  |  回复: 15

qq295118688

银虫 (小有名气)

[求助] 求助高人:使用Matlab或Origin进行非线性拟合的方法已有1人参与

求助,现在已知公式如下:
y=m*k*c*x/((1-k*x)*(1-k*x+c*k*x))
其中m,k,c为待求参数,
请问怎么使用Matlab或者Origin求的以上三个参数?
感谢!

数据如下
x.      y

0        0
5        2.45198
10        3.41962
15        4.25807
20        4.98196
25        5.69367
30        6.31764
40        7.78736
50        9.24003
60        10.73413
70        12.10637
75        13.16662
80        14.33655
85        16.40343
90        19.0894
95        22.68451
98        25.90914


感恩!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

hollya

至尊木虫 (文坛精英)

葛城?{美

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
qq295118688: 金币+30, ★★★★★最佳答案 2019-07-23 04:08:20
独孤神宇: 金币+5, 鼓励交流 2019-07-23 06:59:42
a=[0 0;
5 2.45198;
10 3.41962;
15 4.25807;
20 4.98196;
25 5.69367;
30 6.31764;
40 7.78736;
50 9.24003;
60 10.73413;
70 12.10637;
75 13.16662;
80 14.33655;
85 16.40343;
90 19.0894;
95 22.68451;
98 25.90914]
f=@(k,x) k(1)*k(2)*k(3)*x./((1-k(2)*x).*(1-k(2)*x+k(3)*k(2)*x))
% k(1)=m,k(2)=k,k(3)=c
k0=[1 1 1]
for ii=1:10
    k=lsqcurvefit(f,k0,a(:,1),a(:,2))
    k0=k
end
plot(a(:,1),a(:,2),'o',a(:,1),f(k,a(:,1)))

» 本帖已获得的红花(最新10朵)

でも私たちはきっと1000回もメールをやりとりして、たぶん心は1センチくらいしか近づけませんでした。
5楼2019-07-23 03:45:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hollya

至尊木虫 (文坛精英)

葛城?{美

引用回帖:
6楼: Originally posted by qq295118688 at 2019-07-23 03:55:37
你好啊,请问怎么看出来R^2的对应数值呢...

后面再加一行:
R_sq=1-sum((a(:,2)-f(k,a(:,1))).^2)/sum((a(:,2)-mean(a(:,1))).^2)

算完R^2=0.9994
再给朵花花呗:)

» 本帖已获得的红花(最新10朵)

でも私たちはきっと1000回もメールをやりとりして、たぶん心は1センチくらいしか近づけませんでした。
7楼2019-07-23 04:03:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖
2楼2019-07-22 22:38:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

southpark007

禁虫 (著名写手)

本帖内容被屏蔽

3楼2019-07-22 22:42:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

465090264

木虫 (正式写手)

MATLAB好久没用了,Python应该可解

发自小木虫Android客户端
4楼2019-07-23 00:41:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qq295118688

银虫 (小有名气)

送红花一朵
引用回帖:
5楼: Originally posted by hollya at 2019-07-23 03:45:31
a=
f=@(k,x) k(1)*k(2)*k(3)*x./((1-k(2)*x).*(1-k(2)*x+k(3)*k(2)*x))
% k(1)=m,k(2)=k,k(3)=c
k0=
for ii=1:10
    k=lsqcurvefit(f,k0,a(:,1),a(:,2))
    k0=k
end
plot(a(:,1),a(:,2),'o',a(:,1),f(k, ...

你好啊,请问怎么看出来R^2的对应数值呢
6楼2019-07-23 03:55:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qq295118688

银虫 (小有名气)

送红花一朵
引用回帖:
7楼: Originally posted by hollya at 2019-07-23 04:03:17
后面再加一行:
R_sq=1-sum((a(:,2)-f(k,a(:,1))).^2)/sum((a(:,2)-mean(a(:,1))).^2)

算完R^2=0.9994
再给朵花花呗:)...

感谢!

发自小木虫IOS客户端
8楼2019-07-23 04:07:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hollya

至尊木虫 (文坛精英)

葛城?{美

引用回帖:
8楼: Originally posted by qq295118688 at 2019-07-23 04:07:01
感谢!
...

多谢啦

» 本帖已获得的红花(最新10朵)

でも私たちはきっと1000回もメールをやりとりして、たぶん心は1センチくらいしか近づけませんでした。
9楼2019-07-23 04:16:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qq295118688

银虫 (小有名气)

引用回帖:
9楼: Originally posted by hollya at 2019-07-23 04:16:51
多谢啦

应该多谢你!以后可能还会再找您咨询!(我发现初始值的设置对结果影响特别特别大,而且不是0.99以上的数据就可靠。)
10楼2019-07-23 15:54:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 qq295118688 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 工材01送了吗? +7 xiaopang8958 2024-05-25 12/600 2024-05-28 19:25 by zhuiming111
[基金申请] 入职高校3年发表10+SCI,尽人事听天命 +21 kaoyan250 2024-05-27 24/1200 2024-05-28 19:23 by ljq158
[基金申请] E10开始送了,希望有好运 +5 sail 2024-05-27 5/250 2024-05-28 18:36 by 芝小芝
[找工作] 找工作如此之难 +5 探123 2024-05-25 5/250 2024-05-28 16:50 by auvauv
[教师之家] 中年 +6 459582015 2024-05-28 7/350 2024-05-28 16:15 by 豫椒
[论文投稿] EI学报,一审返修后,为啥不再送审,直接终审中? 10+3 qweasd12345 2024-05-27 5/250 2024-05-28 14:55 by topedit
[论文投稿] 真急着毕业,CPB主编终审17天了,邮件催稿了两次,就是一点动静没有 5+3 kkkk夏 2024-05-28 5/250 2024-05-28 14:43 by topedit
[有机交流] 机理求助 200+4 15147165026 2024-05-26 10/500 2024-05-28 14:42 by 江东闲人
[硕博家园] 要不要读博 +13 王乔木 2024-05-24 14/700 2024-05-28 13:31 by 905452934
[考博] 招聘催化、储能博士 +6 dfjsf 2024-05-24 10/500 2024-05-28 13:01 by 来自星星jx
[有机交流] D-阿拉伯糖-1,4-内酯的合成 10+4 Leeu55 2024-05-24 10/500 2024-05-28 11:34 by disguiseme
[硕博家园] 又想换工作 +11 brightmj 2024-05-27 14/700 2024-05-28 10:32 by chemhua
[基金申请] 基金上会 +14 mrKiller 2024-05-25 20/1000 2024-05-28 10:11 by bnullh
[电化学] 2200mA锂离子电池设计 +5 小乐987 2024-05-23 5/250 2024-05-28 10:02 by 打工艺术家
[材料综合] 求助,刚入门MOF合成,实验总是失败 20+3 张月半 2024-05-25 8/400 2024-05-28 09:20 by 我想更懂你29
[基金申请] 面上基金会评专家,有回避机制吗? +4 huang1991js 2024-05-27 4/200 2024-05-27 19:08 by 星火12
[硕博家园] 我是很理想化一人 +6 hahamyid 2024-05-26 6/300 2024-05-27 18:13 by 大飞鱼鱼鱼
[基金申请] 感觉自然基金限制通过比例就是有点扯,学学B口,化学学部,不限制比例。 +10 wsjing 2024-05-26 14/700 2024-05-27 11:57 by kanmiaolucky
[硕博家园] 周日 +6 1加油哦棒 2024-05-26 9/450 2024-05-27 10:30 by hahamyid
[硕博家园] 人生 +15 暮色恋伊人 2024-05-22 15/750 2024-05-26 08:23 by elainzai
信息提示
请填处理意见