使用matlab求多元非线性模型

zzz.dan

下面是我的数据，是x1,x2,x3,y1,y2类型，回归的模型对于y1目测应该是三元二次的形式，y2的形式不清楚。y=p1+p2*a+p3*b+p4*c+p5*a*b+ p6*b*c+p7*c*a+p8*a^2+p9*b^2+p10*c^2，自己在matlab中操作不太会，请教大佬们如何进行回归在matlab中，最好能附上代码，感谢给位大佬了。

x1	x2	x3	y1	y2
10	0	50	0.283980	3.825798
10	0	100	0.283020	3.837966
10	0	150	0.282200	3.850134
10	0	200	0.281500	3.862302
10	0.05	50	0.284250	3.82519
10	0.05	100	0.283530	3.83675
10	0.05	150	0.282930	3.84831
10	0.05	200	0.282440	3.85987
10	0.1	50	0.284560	3.824582
10	0.1	100	0.284130	3.835534
10	0.1	150	0.283800	3.846485
10	0.1	200	0.283550	3.857437
10	0.15	50	0.284910	3.823974
10	0.15	100	0.284810	3.834317
10	0.15	150	0.284790	3.844661
10	0.15	200	0.284840	3.855005
10	0.2	50	0.285300	3.823366
10	0.2	100	0.285570	3.833101
10	0.2	150	0.285890	3.842837
10	0.2	200	0.286270	3.852573
10	0.25	50	0.285730	3.822757
10	0.25	100	0.286420	3.831885
10	0.25	150	0.287140	3.841013
10	0.25	200	0.287900	3.85014
10	0.3	50	0.286200	3.822149
10	0.3	100	0.287350	3.830669
10	0.3	150	0.288520	3.839188
10	0.3	200	0.289710	3.847708
8	0	50	0.342120	2.985156
8	0	100	0.339560	3.034751
8	0	150	0.337230	3.085751
8	0	200	0.335130	3.13764
8	0.05	50	0.342610	2.984759
8	0.05	100	0.340490	3.033937
8	0.05	150	0.338580	3.084498
8	0.05	200	0.336840	3.135927
8	0.1	50	0.343190	2.984362
8	0.1	100	0.341580	3.033122
8	0.1	150	0.340140	3.083245
8	0.1	200	0.338840	3.134214
8	0.15	50	0.343820	2.983965
8	0.15	100	0.342810	3.032308
8	0.15	150	0.341910	3.081992
8	0.15	200	0.341120	3.132501
8	0.2	50	0.344510	2.983568
8	0.2	100	0.344140	3.031494
8	0.2	150	0.343850	3.080739
8	0.2	200	0.343630	3.130788
8	0.25	50	0.345270	2.983171
8	0.25	100	0.345630	3.030679
8	0.25	150	0.346020	3.079486
8	0.25	200	0.346440	3.129074
8	0.3	50	0.346090	2.982774
8	0.3	100	0.347220	3.029865
8	0.3	150	0.348370	3.078233
8	0.3	200	0.349530	3.127361
6	0	50	0.438460	2.228357
6	0	100	0.433220	2.269241
6	0	150	0.428450	2.310923
6	0	200	0.424080	2.353301
6	0.05	50	0.439390	2.228132
6	0.05	100	0.434970	2.268775
6	0.05	150	0.430940	2.310199
6	0.05	200	0.427240	2.352303
6	0.1	50	0.440460	2.227906
6	0.1	100	0.437010	2.268308
6	0.1	150	0.433850	2.309476
6	0.1	200	0.430950	2.351305
6	0.15	50	0.441650	2.22768
6	0.15	100	0.439280	2.267841
6	0.15	150	0.437140	2.308752
6	0.15	200	0.435160	2.350308
6	0.2	50	0.442920	2.227455
6	0.2	100	0.441760	2.267375
6	0.2	150	0.440710	2.308028
6	0.2	200	0.439770	2.34931
6	0.25	50	0.444340	2.227229
6	0.25	100	0.444490	2.266908
6	0.25	150	0.444690	2.307304
6	0.25	200	0.444940	2.348313
6	0.3	50	0.445830	2.227004
6	0.3	100	0.447410	2.266441
6	0.3	150	0.448980	2.30658
6	0.3	200	0.450540	2.347315
4	0	50	0.624890	1.587742
4	0	100	0.614770	1.620349
4	0	150	0.604650	1.653529
4	0	200	0.595640	1.687265
4	0.05	50	0.627990	1.58764
4	0.05	100	0.618480	1.620134
4	0.05	150	0.609960	1.653191
4	0.05	200	0.602240	1.686791
4	0.1	50	0.630300	1.587538
4	0.1	100	0.622880	1.619919
4	0.1	150	0.616160	1.652852
4	0.1	200	0.610050	1.686317
4	0.15	50	0.632900	1.587436
4	0.15	100	0.627780	1.619704
4	0.15	150	0.623170	1.652513
4	0.15	200	0.618960	1.685844
4	0.2	50	0.635700	1.587333
4	0.2	100	0.633190	1.619489
4	0.2	150	0.630870	1.652175
4	0.2	200	0.628760	1.68537
4	0.25	50	0.638700	1.587231
4	0.25	100	0.639090	1.619274
4	0.25	150	0.639480	1.651836
4	0.25	200	0.639870	1.684896
4	0.3	50	0.642000	1.587129
4	0.3	100	0.645490	1.619059
4	0.3	150	0.648780	1.651498
4	0.3	200	0.651980	1.684423

) C# n2 b$ W/ Y. Y- c$ j% w/ |( k: [
! K, E- f: S8 H$ `! [) c3 J1 g

小小鲁班

( b- m1 S3 B, L5 b* m  oy1和y2都用这个公式试试：
y = 1/(p1+p2*x1+p3*x2+p4*x3)+p5;
( ?" q% _% n% y# w( N8 ?
" l4 J2 E9 f7 cy1：
8 V, k/ m7 v( ]" i9 [均方差(RMSE): 0.00303582928390663
残差平方和(SSR): 0.0010322210573948
) }* C" _8 A" y' K1 K5 }8 x: j相关系数(R): 0.999727213832403
9 \+ ]& y! p$ n! w' G/ |' u相关系数之平方(R^2): 0.9994545020771
修正R平方(Adj. R^2): 0.999434109631384
$ S. w6 i0 M) Y3 K& z4 p: o6 E确定系数(DC): 0.9994545020771
卡方系数(Chi-Square): 0.000950753215952914
F统计(F-Statistic): 49011.0168734231
- Z1 ~/ e& w( R" _3 z
" C( `- t( g) f, U! P1 a7 _) A# X参数        最佳估算
$ K  k6 g* r7 Y' W% G+ O1 d----------        -------------
/ q1 o' ]# N8 w: U( ep1        0.088856700939011
0 |, _" @: G& R& }5 y/ I2 M: Yp2        0.417023913591413
p3        -0.367124775736504
* a! M! K8 M, j+ c7 R  Mp4        0.000224185510016508
4 K' j7 T+ q% n/ bp5        0.0486268290232504
$ h0 _1 K. }: ~
9 |) I! j( @( Z2 p1 _# dy2：
均方差(RMSE): 0.0221863589448616
4 |: X3 L; u+ P残差平方和(SSR): 0.0551302666017869
相关系数(R): 0.999639775152374
相关系数之平方(R^2): 0.999279680066688
6 U! c$ S1 o" r" P. T修正R平方(Adj. R^2): 0.999252752218714
确定系数(DC): 0.999279680066687
卡方系数(Chi-Square): 0.00923570718418041
% |1 t" n0 z5 A0 r9 R% QF统计(F-Statistic): 37109.5262254675
8 l: x2 a1 [! h  Y
% q3 C: i) P; T5 J/ o5 M0 @( f参数        最佳估算
----------        -------------
p1        0.0641980629886512
& z& \# [1 t3 tp2        -0.00114927294188735
6 F* T/ {3 F1 C; ip3        5.5884605145893E-5
p4        -2.01579028865093E-6
p5        -15.2184718761159

大小的小

楼上大神，来学习一下

zzz.dan

小小鲁班 发表于 2021-1-21 13:10
# _" c# q9 w, o2 r( s, Ky1和y2都用这个公式试试：
y = 1/(p1+p2*x1+p3*x2+p4*x3)+p5;

谢谢大佬！！！