matlab多元非線性回歸及顯著性分析(實(shí)例)

matlab多元非線性回歸及顯著性分析給各位高手：小弟有一些數(shù)據(jù)需要回歸分析（非線性）及顯著性檢驗(yàn)（回歸模型，次要項(xiàng)，誤差及失擬項(xiàng)純誤差，F(xiàn)值和P值），求大俠幫助，給出程序，不勝感激。

模型：

DATA=...

%DATA前三列是影響因子，第四列為響應(yīng)值

[21307548.61;

21107556.43;

21304561.32;

21104565.28;

11104555.80;

11307545.65;

11107550.91;

11304567.94;

1.51206074.15;

1.51206071.28;

1.51206077.95;

1.51206074.16;

1.51206075.20;

1.51208535.65;

1.51406048.66;

1.51203074.10;

1.51006062.30;

0.51206066.00;

2.51206075.10];回歸分析過(guò)程：（1）MATLAB編程步驟1：首先為非線性回歸函數(shù)編程，程序存盤(pán)為user_function.m

functiony=user_function(beta,x)

b0=beta(1);

b1=beta(2);

b2=beta(3);

b3=beta(4);

x0=x(:,1);

x1=x(:,2);

x2=x(:,3);

x3=x(:,4);

y=b0*x0+b1*x1.^2+b2*x2.^2+b3*x3.^2;

（2）MATLAB編程步驟2：編寫(xiě)非線性回歸主程序，程序運(yùn)行時(shí)調(diào)用函數(shù)user_function

x=[121307548.61;

121107556.43;

121304561.32;

121104565.28;

111104555.80;

111307545.65;

111107550.91;

111304567.94;

11.51206074.15;

11.51206071.28;

11.51206077.95;

11.51206074.16;

11.51206075.20;

11.51208535.65;

11.51406048.66;

11.51203074.10;

11.51006062.30;

10.51206066.00;

12.51206075.10];%%第1列全是1，第6列是指標(biāo)變量，其余列是自變量

xx=x(:,1:5);yy=x(:,5);

%%指定響應(yīng)變量yy和自變量xx

beta0=[0.50.40.70.5];

%%設(shè)置初始回歸系數(shù)（如何確定初值？）

[beta_fit,residual]=nlinfit(xx,yy,@user_function,beta0)

%%非線性回歸

結(jié)果

beta_fit=

91.3757

1.2712

-0.0009

-0.0049

residual=

-4.2935

-1.0248

-9.2044

-9.7957

-15.4620

-3.4398

-2.7311

1.2293

11.1889

8.3189

14.9889

11.1989

12.2389

-9.5678

-9.3704

-2.0767

-4.8331

5.5814