24小时热门版块排行榜    

查看: 1624  |  回复: 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的回帖

hollya

至尊木虫 (文坛精英)

葛城?{美

引用回帖:
12楼: Originally posted by qq295118688 at 2019-07-24 18:04:50
请问,如果把我的公式改成这个,m*K*x*(1-k*(1-w)*x)./((1+K*x).*(1-k*x)),代码要怎么写呢?
原始数据如下:
x.      y

0        0
5        2.45198
10        3.41962
15        4.25807
20        4. ...

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)*x.*((1-k(4)*(1-k(3))*x))./((1+k(2)*x).*(1-k(4)*x))
% k(1)=m,k(2)=K,k(3)=w,k(4)=k
k0=[1 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)))
R_sq=1-sum((a(:,2)-f(k,a(:,1))).^2)/sum((a(:,2)-mean(a(:,1))).^2)
でも私たちはきっと1000回もメールをやりとりして、たぶん心は1センチくらいしか近づけませんでした。
14楼2019-07-25 18:32:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 16 个回答
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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 求推荐博导 +6 想读博的我 2024-05-05 8/400 2024-05-07 09:24 by 丁jiayou
[论文投稿] Optics Express一直卡在Final Decision +3 MissPro 2024-05-05 8/400 2024-05-07 09:23 by MissPro
[高分子] 求推荐研究U V热熔胶的高手 +4 千年一虫 2024-05-04 4/200 2024-05-07 09:16 by Sunrisepay
[考研] 关于要不要三战 +8 煎饼果子007 2024-05-06 12/600 2024-05-07 08:32 by 煎饼果子007
[论文投稿] 求各位大神推荐复合材料类中文核心 50+3 ayj080121 2024-05-01 5/250 2024-05-07 07:56 by 修为01
[基金申请] 2024国家自然科学基金资助额度 +13 i6277 2024-05-05 14/700 2024-05-07 00:24 by dxcharlary
[基金申请] 祈福今年优青或者QB能有个好消息 66+3 5011fudan 2024-05-04 12/600 2024-05-06 23:48 by dxcharlary
[有机交流] 酰胺合成 +3 Spiralup 2024-05-06 6/300 2024-05-06 19:39 by Spiralup
[教师之家] 七级岗,到手不到8000,带本科论文12教学分/生,额定教学分280 +25 河西夜郎 2024-05-02 28/1400 2024-05-06 19:30 by jhhcooi
[教师之家] 博士毕业找工作选择! +11 gsq嘿嘿 2024-04-30 15/750 2024-05-06 16:46 by 半生梦君
[论文投稿] 林学有投过总环的朋友吗? 50+3 New…… 2024-05-04 4/200 2024-05-06 16:15 by wuyuanzhao
[考博] 考博求推荐24-25都行 +3 godweiliu 2024-05-01 5/250 2024-05-06 15:43 by 愿你快乐啊
[考博] 2024还有没有补录的博士 +7 crazy peng 2024-04-30 8/400 2024-05-06 15:37 by 愿你快乐啊
[基金申请] 基金函评开始了吗 +5 fuzy 2024-05-05 5/250 2024-05-06 11:57 by hdzw9071
[硕博家园] 等盲审 +10 暮色恋伊人 2024-05-03 18/900 2024-05-06 11:44 by 暮色恋伊人
[基金申请] 多大算大龄? +5 wyjecho666 2024-05-01 6/300 2024-05-06 11:06 by 梦想成真2014
[考博] 求博导收留 +4 黑森林nono派 2024-05-06 4/200 2024-05-06 08:02 by 安塔瓦拉多
[考博] 西安交通大学转化医学研究院蔡文课题组招收2025年入学的申请-考核制博士生 +3 caiwen2016 2024-05-03 3/150 2024-05-05 09:39 by 宋mlk
[基金申请] 关于公布2024年度国家自然科学基金项目申请初审结果的通告 (金币+2) +12 hdzw9071 2024-04-30 18/900 2024-05-04 17:49 by felicity6056
[基金申请] 今年面上预计又打酱油了 +15 speedxyl 2024-04-30 17/850 2024-05-04 02:07 by jbb0523
信息提示
请填处理意见