| 查看: 259 | 回复: 0 | |||
duandian12木虫 (小有名气)
|
[求助]
求拟合指点
|
|
烦大神们帮帮忙，我有数据x=2.822E-4，0.00495，0.00995，0.0197，0.03985，0.06047，0.07989，0.10009，0.14979，0.20018，0.24984，0.29978，0.39984，0.5003，0.60034，0.70045，0.79981，0.89997，0.99995，1.09996，1.19987，1.29989，1.40027，1.50004，1.59983，1.69939，1.79978，1.89995 y=7.5E-6，6.71875E-4，0.00391，0.0106，0.023，0.0358，0.04649，0.05719，0.08382，0.10749，0.13196，0.15581，0.20194，0.24542，0.2914，0.33464，0.37335，0.42031，0.46305，0.50221，0.54497，0.58705，0.62944，0.66904，0.71459，0.75343，0.79638，0.83817 想用方程y=y0*x*((1-y/y0)^a)*b*exp(-c/2311.3)进行拟合，求出参数y0,a,b,c。跪求各位支下招，小弟不胜感激.......... clc; clear A=[... 2.822E-4 7.5E-6 0.00495 6.71875E-4 0.00995 0.00391 0.0197 0.0106 0.03985 0.023 0.06047 0.0358 0.07989 0.04649 0.10009 0.05719 0.14979 0.08382 0.20018 0.10749 0.24984 0.13196 0.29978 0.15581 0.39984 0.20194 0.5003 0.24542 0.60034 0.2914 0.70045 0.33464 0.79981 0.37335 0.89997 0.42031 0.99995 0.46305 1.09996 0.50221 1.19987 0.54497 1.29989 0.58705 1.40027 0.62944 1.50004 0.66904 1.59983 0.71459 1.69939 0.75343 1.79978 0.79638 1.89995 0.83817]; format short g x=A(:,1);y=A(:,2); X=[x,y]; %y=y0*k0*x*((1-y/y0)^a)*exp(-H/2311.3) %y0=a(1),k0=a(2),a=a(3),H=a(4) fun=inline('a(1)*a(2)* x(:,1)*((1- x(:,2)/a(1)).^a(3))*exp(-a(2)/2311.3)-x(:,2)','a','x') a=nlinfit(x,y,fun,[1 1 1 1]) %{ figure(1) plot(x,y) figure(2) semilogy(x,y) figure(3) loglog(x,y) %} 不知怎么回事就是不对，麻烦指点一下 |
回复此楼» 猜你喜欢
不应该看fileCode
已经有8人回复
-
FileCode能看出啥？
已经有26人回复
-
关于Filecode分析方法
已经有11人回复
-
重要来源：本周末出结果
已经有9人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有14人回复
-
结合人工智能，周易传统文化，filecode打分制来了，3分以上希望很大。
已经有4人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有11人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有36人回复
-
好奇怪的filecode
已经有6人回复
-
filecode
已经有15人回复
