一个多元回归的题目,几个地方搞不太清楚,可以帮我看下么?%x(:,1)-冻融循环次数%x(:,2)-围压%x(:,3)-塑性指数x=[010010.7110010.7210010.7310010.7410010.7510010.7610010.7710010.7020010.7120010.7220010.

Question

一个多元回归的题目,几个地方搞不太清楚,可以帮我看下么?
%x(:,1)-冻融循环次数
%x(:,2)-围压
%x(:,3)-塑性指数
x=[0 100 10.7
1 100 10.7
2 100 10.7
3 100 10.7
4 100 10.7
5 100 10.7
6 100 10.7
7 100 10.7
0 200 10.7
1 200 10.7
2 200 10.7
3 200 10.7
4 200 10.7
5 200 10.7
6 200 10.7
7 200 10.7
0 300 10.7
1 300 10.7
2 300 10.7
3 300 10.7
4 300 10.7
5 300 10.7
6 300 10.7
7 300 10.7
0 100 15.98
1 100 15.98
2 100 15.98
3 100 15.98
4 100 15.98
5 100 15.98
6 100 15.98
7 100 15.98
0 200 15.98
1 200 15.98
2 200 15.98
3 200 15.98
4 200 15.98
5 200 15.98
6 200 15.98
7 200 15.98
0 300 15.98
1 300 15.98
2 300 15.98
3 300 15.98
4 300 15.98
5 300 15.98
6 300 15.98
7 300 15.98
]
%y-抗剪强度
y=[110.80301
104.54847
84.84117
82.52244
92.32921
83.87712
83.96642
75.53389
123.73915
113.16948
97.96604
100.56338
102.22341
101.01701
99.90815
104.6684
132.8442
122.48814
116.48945
112.85501
116.04544
125.10151
116.39448
122.91425
135.56231
132.88212
77.77271
107.49568
108.52345
98.69237
111.2504
108.052
179.24035
180.45719
112.78747
123.59022
123.94181
121.77485
127.45383
127.55048
229.56915
233.42255
134.94563
133.17741
138.31906
137.60706
147.51052
153.53207
]
%选择这个函数不知道是不是合理
myfunc=inline('beta(1)*x(:,2).*exp((beta(2)*x(:,3))./(beta(3)*x(:,1)+beta(4)))','beta','x')
beta=nlinfit(x,y,myfunc,[1 0.1 1 2])
a=beta(1),b=beta(2),c=beta(3),d=beta(4)
yy=a*x(:,2).*exp(b*x(:,3))./(c*x(:,1)+d)
%计算复相关系数 ,这个小程序好像哪里有问题
sum = 0.0;
absolute y1 = abs(yy-y);
for i = 1:length(absolute(y1))
sum = sum+(absolute(y1(i))).^2;
end
sum;
mean y = mean(y);
absolute y2 = abs(y-mean y);
sum1 = 0.0;
for i = 1:length(absolute(y2))
sum1 = sum1+(absolute(y2(i))).^2;
end
sum1
r = sqrt(1-sum/sum1)
后面计算相关系数的小程序好像确实有问题,我把前面计算得数导至excel,手算相关系数才0.58左右,确实不好,呵呵,我是不是应该换几个拟合公式试一下呢?“数”海茫茫,实在不好选择啊!

Answer 1

拟合的效果不怎么样!
你后面的
absolute
我也看不懂啊.
你做完前面的,你就
plot(y,yy,'*');hold on;plot(y,y,'r-')
如果点在线两侧附近,应当就是拟合成功了.