24小时热门版块排行榜    

查看: 894  |  回复: 4

zy020845

银虫 (小有名气)

[求助] Matlab三元非线性回归,急求 已有3人参与

自变量X1,X2,X3,系数a,b,c,方程y=a*x1+b*x2+c*x3,x1,x2,x3和y数据如下,怎么编写Matlab非线性回归程序,计算出系数a,b,c
-0.00359        -4.12343E-18        -1.77469E-23        2.01725E-23
0.00208        -3.60585E-18        -1.10372E-23        1.04121E-23
-2.0234E-4        -3.16446E-18        -8.94233E-25        5.34096E-24
-0.00107        -2.75194E-18        -1.85101E-24        1.52903E-24
6.72658E-4        -2.23653E-18        -1.08442E-24        -9.10219E-24
0.00165        -1.81359E-18        2.74972E-24        -1.65717E-23
0.00107        -1.36763E-18        7.73421E-24        -2.00725E-23
0.00133        -9.56332E-19        6.28137E-24        -2.24341E-23
-0.00103        -4.95781E-19        1.07942E-23        -2.91833E-23
-6.53972E-4        -6.06294E-20        4.0151E-24        -3.24459E-23
-5.90757E-4        1.81729E-19        8.4007E-24        -3.24577E-23
-8.7934E-4        4.21581E-19        1.0951E-23        -3.21475E-23
6.68916E-4        5.43246E-19        1.42788E-23        -3.53133E-23
4.69911E-4        6.3032E-19        6.34717E-24        -3.50421E-23
-7.82425E-4        7.41742E-19        1.40619E-23        -3.29354E-23
2.97379E-4        8.60789E-19        1.23828E-23        -2.91212E-23
6.94096E-4        9.47386E-19        1.46529E-23        -2.62905E-23
0.00157        9.67946E-19        7.15769E-24        -2.12477E-23
0.00135        1.03696E-18        1.63978E-23        -1.47819E-23
1.33784E-5        1.12469E-18        -3.88451E-24        -1.19203E-23
-1.06037E-4        1.20089E-18        1.16233E-23        -1.10687E-23
4.55055E-5        1.27435E-18        9.21188E-24        -8.30528E-24
-6.07391E-4        1.37333E-18        -1.77745E-25        -4.22855E-24
-0.00104        1.37709E-18        5.00915E-25        -2.66514E-24
-0.00101        1.44639E-18        -3.09529E-25        -2.35639E-24
-0.00141        1.48302E-18        4.9184E-24        -1.51231E-24
-0.00126        1.55521E-18        -4.71789E-24        1.13818E-24
-0.00111        1.60539E-18        -1.6555E-24        3.76909E-24
-0.00119        1.64962E-18        -8.81191E-24        5.35461E-24
-0.00112        1.66058E-18        -1.57947E-23        6.27912E-24
-0.00106        1.69684E-18        -1.64222E-23        6.67755E-24
-0.00123        1.71129E-18        -1.6715E-23        9.15294E-24
-5.36041E-4        1.73407E-18        -1.07468E-23        1.15802E-23
-4.78845E-4        1.72189E-18        -1.2255E-23        1.25628E-23
-5.98937E-5        1.76237E-18        -7.41704E-24        1.32947E-23
2.08466E-4        1.78578E-18        -2.02457E-23        1.41029E-23
0.00104        1.81716E-18        -1.24537E-23        1.51631E-23
0.00145        1.80832E-18        -1.45689E-23        1.66509E-23
0.00119        1.79953E-18        -6.80905E-24        1.76409E-23
9.4855E-4        1.81684E-18        -1.11621E-23        1.95819E-23
9.97514E-4        1.76832E-18        -1.52392E-23        2.40777E-23
6.92467E-4        1.73818E-18        -1.16444E-23        2.70048E-23
8.66976E-4        1.69993E-18        -1.40758E-23        3.37399E-23
4.11891E-4        1.64462E-18        -9.901E-24        3.75995E-23
3.44032E-4        1.58625E-18        -4.10636E-24        4.15453E-23
3.22188E-4        1.50451E-18        -1.20865E-23        4.70289E-23
5.11915E-4        1.41122E-18        1.8058E-24        5.41276E-23
5.42298E-4        1.29946E-18        -1.07575E-23        6.4919E-23
2.26713E-4        1.21333E-18        -1.57445E-24        7.24342E-23
-9.31989E-5        1.11966E-18        -6.44683E-24        8.17031E-23
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

CelestialCYJ

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
zy020845: 金币+30, ★★★很有帮助 2015-07-16 13:41:03
clc
clear
D=[-0.00359        -4.12343E-18        -1.77469E-23        2.01725E-23
0.00208        -3.60585E-18        -1.10372E-23        1.04121E-23
-2.0234E-4        -3.16446E-18        -8.94233E-25        5.34096E-24
-0.00107        -2.75194E-18        -1.85101E-24        1.52903E-24
6.72658E-4        -2.23653E-18        -1.08442E-24        -9.10219E-24
0.00165        -1.81359E-18        2.74972E-24        -1.65717E-23
0.00107        -1.36763E-18        7.73421E-24        -2.00725E-23
0.00133        -9.56332E-19        6.28137E-24        -2.24341E-23
-0.00103        -4.95781E-19        1.07942E-23        -2.91833E-23
-6.53972E-4        -6.06294E-20        4.0151E-24        -3.24459E-23
-5.90757E-4        1.81729E-19        8.4007E-24        -3.24577E-23
-8.7934E-4        4.21581E-19        1.0951E-23        -3.21475E-23
6.68916E-4        5.43246E-19        1.42788E-23        -3.53133E-23
4.69911E-4        6.3032E-19        6.34717E-24        -3.50421E-23
-7.82425E-4        7.41742E-19        1.40619E-23        -3.29354E-23
2.97379E-4        8.60789E-19        1.23828E-23        -2.91212E-23
6.94096E-4        9.47386E-19        1.46529E-23        -2.62905E-23
0.00157        9.67946E-19        7.15769E-24        -2.12477E-23
0.00135        1.03696E-18        1.63978E-23        -1.47819E-23
1.33784E-5        1.12469E-18        -3.88451E-24        -1.19203E-23
-1.06037E-4        1.20089E-18        1.16233E-23        -1.10687E-23
4.55055E-5        1.27435E-18        9.21188E-24        -8.30528E-24
-6.07391E-4        1.37333E-18        -1.77745E-25        -4.22855E-24
-0.00104        1.37709E-18        5.00915E-25        -2.66514E-24
-0.00101        1.44639E-18        -3.09529E-25        -2.35639E-24
-0.00141        1.48302E-18        4.9184E-24        -1.51231E-24
-0.00126        1.55521E-18        -4.71789E-24        1.13818E-24
-0.00111        1.60539E-18        -1.6555E-24        3.76909E-24
-0.00119        1.64962E-18        -8.81191E-24        5.35461E-24
-0.00112        1.66058E-18        -1.57947E-23        6.27912E-24
-0.00106        1.69684E-18        -1.64222E-23        6.67755E-24
-0.00123        1.71129E-18        -1.6715E-23        9.15294E-24
-5.36041E-4        1.73407E-18        -1.07468E-23        1.15802E-23
-4.78845E-4        1.72189E-18        -1.2255E-23        1.25628E-23
-5.98937E-5        1.76237E-18        -7.41704E-24        1.32947E-23
2.08466E-4        1.78578E-18        -2.02457E-23        1.41029E-23
0.00104        1.81716E-18        -1.24537E-23        1.51631E-23
0.00145        1.80832E-18        -1.45689E-23        1.66509E-23
0.00119        1.79953E-18        -6.80905E-24        1.76409E-23
9.4855E-4        1.81684E-18        -1.11621E-23        1.95819E-23
9.97514E-4        1.76832E-18        -1.52392E-23        2.40777E-23
6.92467E-4        1.73818E-18        -1.16444E-23        2.70048E-23
8.66976E-4        1.69993E-18        -1.40758E-23        3.37399E-23
4.11891E-4        1.64462E-18        -9.901E-24        3.75995E-23
3.44032E-4        1.58625E-18        -4.10636E-24        4.15453E-23
3.22188E-4        1.50451E-18        -1.20865E-23        4.70289E-23
5.11915E-4        1.41122E-18        1.8058E-24        5.41276E-23
5.42298E-4        1.29946E-18        -1.07575E-23        6.4919E-23
2.26713E-4        1.21333E-18        -1.57445E-24        7.24342E-23
-9.31989E-5        1.11966E-18        -6.44683E-24        8.17031E-23];
X=D(:,1:3);
Y=D(:,4);
V=X\Y;
2楼2015-07-15 18:03:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangzhiguang

至尊木虫 (文坛精英)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
zy020845: 金币+45, ★★★★★最佳答案 2015-07-16 13:41:17
建立一个函数,另外再建立一个脚本,如下:
函数为:
function y=fxnhx(beta,x)
A=beta(1);
B=beta(2);
C=beta(3);
x1=x(:,1);
x2=x(:,2);
x3=x(:,3);
y=A.*x1+B.*x2+C.*x3;%此处为函数形式
脚本为:
clear;
clc;
data=xlsread('3.xlsx');%所有数据
x=data(:,1:3);%自变量
y=data(:,4);%因变量
beta0=[1,1,1]';%三个系数的初值
[b,r]=nlinfit(x,y,'fxnhx',beta0);
其中,b就是三个系数。不过,matlab的非线性拟合依赖于初值的选择,所以最好知道大体的范围,在此,取值为1,1,1.
好好学习,天天向上!
3楼2015-07-15 18:57:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangzhiguang

至尊木虫 (文坛精英)

我把你的数据合到一个文件3.xlsx中了,前三列为自变量,最后一列为因变量。
好好学习,天天向上!
4楼2015-07-15 18:59:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
zy020845: 金币+5, 有帮助 2015-07-16 13:41:25
MATLAB里面有个拟合工具箱,简单,方便,适合比较简单的函数和一些特殊的函数拟合

输入: cftool  回车
数值计算
5楼2015-07-15 19:43:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zy020845 的主题更新
信息提示
请填处理意见