求助高人：使用Matlab或Origin进行非线性拟合的方法 - 计算模拟 - 程序代码

でも私たちはきっと1000回もメールをやりとりして、たぶん心は1センチくらいしか近づけませんでした。

11楼2019-07-23 16:36:19

qq295118688

银虫 (小有名气)

应助: 2 (幼儿园)
金币: 476.2
散金: 75
红花: 3
帖子: 279
在线: 23.9小时
虫号: 3442403
注册: 2014-09-26
性别: GG
专业: 高分子材料结构与性能

送红花一朵

引用回帖:

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

请问，如果把我的公式改成这个，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.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
求四个参数及其R2

12楼2019-07-24 18:04:50

dingd

铁杆木虫 (职业作家)

计算强帖: 4
应助: 1641 (讲师)
金币: 15356.8
散金: 101
红花: 234
帖子: 3409
在线: 1222.7小时
虫号: 291104
注册: 2006-10-28

试下1stOpt，不用猜初值：

CODE:

Function y=m*K*x*(1-k1*(1-w)*x)/((1+K*x)*(1-k1*x));

Data;

x=0,5,10,15,20,25,30,40,50,60,70,75,80,85,90,95,98;

y=0,2.45198,3.41962,4.25807,4.98196,5.69367,6.31764,7.78736,9.24003,10.73413,12.10637,13.16662,14.33655,16.40343,19.0894,22.68451,25.90914;

Root of Mean Square Error (RMSE): 0.2146670130891
Sum of Squared Residual: 0.783392750646129
Correlation Coef. (R): 0.999552613989238
R-Square: 0.999105428132719
Adjusted R-Square: 0.998977632151679
Determination Coef. (DC): 0.999092348635769
Chi-Square: 0.0921118108747331
F-Statistic: 4798.54852809769

Parameter Best Estimate
---------- -------------
m 9.74847626007366
k 0.0485707895587777
k1 0.00867325606751355
w 0.391166540974436

1stOpt不区分大小写，原公式中“k”由现在的“k1”替代。

» 本帖已获得的红花（最新10朵）

qq295118688

13楼2019-07-25 10:01:32

hollya

至尊木虫 (文坛精英)

葛城?{美

应助: 744 (博后)
金币: 43440.6
散金: 1143
红花: 127
沙发: 779
帖子: 24904
在线: 834.8小时
虫号: 458790
注册: 2007-11-13
性别: MM
专业: 计算机软件

引用回帖:

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

qq295118688

银虫 (小有名气)

应助: 2 (幼儿园)
金币: 476.2
散金: 75
红花: 3
帖子: 279
在线: 23.9小时
虫号: 3442403
注册: 2014-09-26
性别: GG
专业: 高分子材料结构与性能

引用回帖:

14楼: Originally posted by hollya at 2019-07-25 18:32:39
a=
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=
for ii=1:10
k=lsqcurvefit(f,k0,a(:,1),a(:,2))
k0=k
end
plot(a(:,1),a(:,2),' ...

感谢！！！

15楼2019-07-25 23:00:03

qq295118688

银虫 (小有名气)

应助: 2 (幼儿园)
金币: 476.2
散金: 75
红花: 3
帖子: 279
在线: 23.9小时
虫号: 3442403
注册: 2014-09-26
性别: GG
专业: 高分子材料结构与性能

送红花一朵

引用回帖:

13楼: Originally posted by dingd at 2019-07-25 10:01:32
试下1stOpt，不用猜初值：

Function y=m*K*x*(1-k1*(1-w)*x)/((1+K*x)*(1-k1*x));
Data;
x=0,5,10,15,20,25,30,40,50,60,70,75,80,85,90,95,98;
y=0,2.45198,3.41962,4.25807,4.98196,5.69367,6.31764,7.7873 ...

非常感谢！我周末回去了尝试一下这个软件，之前没有接触过呢！

16楼2019-07-25 23:00:36