[m,n]=size(X11); X22=[]; for i=2:n for j=2:n if i<=j
X22=([X22,X11(:,i).*X11(:,j)]); else continue end end end
X=[X11,X22];
[B,BINT,R,RINT,STATS]=regress(Y,X) plot(Y,X*B,'*'),
hold on,line([min(y),max(y)],[min(y),max(y)]) axis([min(y) max(y) min(y) max(y)]) legend('曲线拟合')
xlabel('实际值');ylabel('计算值') X*B
B = 0 -3.64267562479661 1.72707370186922 4.43441169883136 0.0020463400088433 -0.00012494765204425 -0.000212405933828659 -0.00102542042789814 -0.00297221560425269 4.32021484482528e-005
BINT = 0 28
0
-6.3474049893683 -0.937946260224914 -4.51967556506294 7.97382296880138 -7.05846578314345 15.9272891808062 0.000903069409478543 0.00318961060820805 -0.00305282303978543 0.00280292773569693 -0.000570448243772087 0.000145636376114768 -0.00642060480746999 -0.0109578976295441 -7.17448030178289e-005 R = -16.1681949759024 73.3055843187963 72.0505282654747 -89.5115373178149 -96.9855516920856 -160.117234856269 253.496816551423 -16.1681949759024 347.799442047917 615.295272307765 308.546117078223 248.43155808637 -168.88939426284 -154.354122622331 -16.1681949759024 969.100745240424 740.616730694779 -789.552534864793 -1653.51824656714 1652.97678305466 -319.685940236184
0.00436976395167372 0.00501346642103875 0.000158149099914335 29
2.288818359375e-005
RINT = -1833.98851270957 1801.65212275777 -1727.10219947926 1873.71336811685 -1509.28146003113 1653.38251656208 -1546.8825666164 1367.85949198077 -1547.56718902067 1353.5960856365 -1532.406532561 -724.623539173622 -1833.98851271943 -1441.30861331399 -1028.14577349689 -1226.67640227228 -1089.6972575665 -1364.35297784757 -847.731301140414 -1833.98851271943 -594.158284313735 -484.69142159573 -1621.37243897002 -2156.98713298104 1321.82244939608 -467.382651321028 -Inf STATS = Columns 1 through 3
0.992303471936571 Column 4
654400.225184767 ans = 17.1681949759024 -70.0055843187963 -61.3905282654747
1212.17206284846 1231.61717227647 1801.65212276763 2136.90749740982 2258.73631811242 1843.76863642873 1586.56037373924 1026.57418932189 539.023055895752 1801.65212276763 2532.35977479458 1965.92488298529 42.2673692404347 -1150.04936015324 1984.13111671323 -171.98922915134 Inf 211.692266702313 1.68772296199683e-012 30
101.241537317815 119.395551692086 195.257234856269 -202.346816551423 17.1681949759024 -341.079442047917 -596.515272307765 -287.946117078223 152.97844191363 799.71939426284 786.354122622331 17.1681949759024 -952.100745240424 -694.056730694779 849.582534864793 9080.75824656714 14502.4632169453 18590.1359402362 25999.9999771118 可得图像如下：
31