作业帮 > 综合 > 作业

MATLAB进行对数回归分析求助(高分)

来源:学生作业帮 编辑:搜搜做题作业网作业帮 分类:综合作业 时间:2024/07/20 23:59:29
MATLAB进行对数回归分析求助(高分)
y=[720 725 760 695 745 755 720 750 770 715 715 705 695 715 675 725 715 700 695 745 750 745 720 740 690 710 765 720 745 700 720 720 750 750 730 725 725 720 780 800 760 735 785 770 800
]';
x1=[13.43 11.72 12.23 11.65 12.30 13.08 12.72 13.17 13.29 14.00 12.82 12.04 11.84 12.48 13.24 14.00 13.93 12.58 11.60 12.95 13.25 12.29 12.15 12.54 12.58 12.25 13.38 12.09 11.87 13.17 12.68 13.49 13.70 12.20 11.80 12.73 13.43 12.81 13.65 13.69 13.36 13.79 13.03 12.92 13.54]';
x2=[3.95 3.75 3.69 3.26 2.96 4.00 3.55 3.87 4.00 3.74 3.05 3.39 3.11 3.57 3.37 3.57 3.63 2.88 3.16 3.01 3.49 3.20 3.71 3.24 3.50 3.22 3.95 3.08 2.82 3.22 3.77 3.36 3.45 3.21 3.13 3.27 3.27 3.24 3.37 3.10 3.21 3.07 3.52 3.36 2.72]';
x3=[1.81 1.87 2.50 1.85 2.22 1.80 1.92 1.80 1.82 1.83 1.96 1.96 1.98 1.95 1.85 2.00 1.87 2.14 1.96 2.17 1.93 1.80 2.31 1.98 2.16 2.14 2.39 2.31 2.24 1.88 2.06 1.94 2.06 2.28 1.95 2.06 1.93 2.04 1.91 1.84 1.95 2.49 2.26 2.07 2.01]';
要求进行对数回归分析(非线性回归分析)
求MATLAB代码,
对数线性回归
LOG Y=b0+b1*log x1+b2*log x2+b3*log x3
MATLAB进行对数回归分析求助(高分)
对数回归与线性回归并没有本质区别,把原数据先取log再做线性回归则可.
代码:
y=[720 725 760 695 745 755 720 750 770 715 715 705 695 715 675 725 715 700 695 745 750 745 720 740 690 710 765 720 745 700 720 720 750 750 730 725 725 720 780 800 760 735 785 770 800
]';
x1=[13.43 11.72 12.23 11.65 12.30 13.08 12.72 13.17 13.29 14.00 12.82 12.04 11.84 12.48 13.24 14.00 13.93 12.58 11.60 12.95 13.25 12.29 12.15 12.54 12.58 12.25 13.38 12.09 11.87 13.17 12.68 13.49 13.70 12.20 11.80 12.73 13.43 12.81 13.65 13.69 13.36 13.79 13.03 12.92 13.54]';
x2=[3.95 3.75 3.69 3.26 2.96 4.00 3.55 3.87 4.00 3.74 3.05 3.39 3.11 3.57 3.37 3.57 3.63 2.88 3.16 3.01 3.49 3.20 3.71 3.24 3.50 3.22 3.95 3.08 2.82 3.22 3.77 3.36 3.45 3.21 3.13 3.27 3.27 3.24 3.37 3.10 3.21 3.07 3.52 3.36 2.72]';
x3=[1.81 1.87 2.50 1.85 2.22 1.80 1.92 1.80 1.82 1.83 1.96 1.96 1.98 1.95 1.85 2.00 1.87 2.14 1.96 2.17 1.93 1.80 2.31 1.98 2.16 2.14 2.39 2.31 2.24 1.88 2.06 1.94 2.06 2.28 1.95 2.06 1.93 2.04 1.91 1.84 1.95 2.49 2.26 2.07 2.01]';
P=[ones(size(x1)),log(x1),log(x2),log(x3)];
Q=log(y);
B=P\Q;
b0=B(1)
b1=B(2)
b2=B(3)
b3=B(4)
logyy=b0+b1*log(x1)+b2*log(x2)+b3*log(x3);
yy=exp(logyy)