matlab的regress如何使用,[转载]Matlab regress用法

clear

clc

data = [5 5 5 5 5 5 5 10 10 10 10 10 10 10 15 15 15 15 15 15 15 20 20 20 20 20 20 20 25 25 25 25 25 25 25 30 30 30 30 30 30 30

144 176 208 240 272 320 512 144 176 208 240 272 320 512 144 176 208 240 272 320 512 144 176 208 240 272 320 512 144 176 208 240 272 320 512 144 176 208 240 272 320 512

135.346 138.193 147.369 137.871 140.782 143.022 146.985 179.075 230.787 236.99 240.18 232.734 248.473 249.812 218.834 242.56 276.038 282.54 293.582 320.17 279.204 215.68 286.505 313.645 334.782 325.526 332.145 337.282 231.54 261.884 324.214 343.436 367.167 372.446 340.181 259.584 302.132 313.351 347.63 340.478 414.03 443.249

]';

x1 = data(:,1);

x2 = data(:,2);

y = data(:,3);

X = [ones(size(x1)) x1.*x1 x1 x2.*x2 x2 x1.*x2];

[b,bint] = regress(y,X)

scatter3(x1,x2,y,'filled')

hold on

x1fit = min(x1):0.5:max(x1);

x2fit = min(x2):0.5:max(x2);

[X1FIT,X2FIT] = meshgrid(x1fit,x2fit);

YFIT = b(1)+ b(2)*X1FIT.*X1FIT+b(3)*X1FIT + b(4)*X2FIT.*X2FIT + b(5)*X2FIT + b(6)*X1FIT.*X2FIT;

mesh(X1FIT,X2FIT,YFIT)

xlabel('x1')

ylabel('x2')

zlabel('Y')

view(140,30)

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值