最小二乘法matlab输出系数,最小二乘法求9个系数

Xs = [-1604466.2193, -1487333.5778, -1365029.0194, -1238288.6635, -1107862.7319, -974511.4255, -839000.8017, -702098.6769, -564570.5119, -427175.3011, -290661.5563];

Ys = [6736124.5647, 6654963.7330, 6545886.7267, 6409291.3003, 6245693.7579, 6055726.8467, 5840137.1510, 5599782.0105, 5335625.8311, 5048735.8084, 4740277.2913];

Zs = [1447001.5648, 1879963.1946, 2305278.4178, 2721217.2671, 3126088.2190, 3518245.2085, 3896094.4039, 4258100.8139, 4602794.6913, 4928777.5156, 5234727.4850];

t = [55860,55920,55980,56040,56100,56160,56220,56280,56340,56400,56460];

VSX = [ 1905.09424, 1997.33184, 2077.44098, 2145.15256, 2200.26579, 2242.64880, 2272.23814, 2289.03871, 2293.12414, 2284.63607, 2263.78183]

VSY = [ -1118.33538, -1586.25409, -2048.54081, -2503.20530, -2948.28775, -3381.86793, -3802.07274, -4207.08518, -4595.15508, -4964.60727, -5313.84566]

VSZ = [ 7270.00003, 7157.16003, 7015.20665, 6844.72106, 6646.40309, 6421.06745, 6169.64011, 5893.15596, 5592.75369, 5269.66737, 4925.21945]

CX = Ys.*VSZ - Zs.*VSY;

CY = Zs.*VSX - Xs.*VSZ;

CZ = Xs.*VSY - Ys.*VSX;

C = sqrt(CX.^2 + CY.^2 + CZ.^2);

rho = sqrt(Xs.^2 + Ys.^2 + Zs.^2);

I = acos(CZ./C);

w = asin(Zs./(rho.*sin(I)));

o = atan(-CX./CY);

tbl = table(rho',I',w',o',t','VariableNames',{'rho','I','w','o','t'})

mdl1 = fitlm(tbl,'rho ~ t + t^2');

mdl2 = fitlm(tbl,'I ~ t');

mdl3 = fitlm(tbl,'w ~ t');

mdl4 = fitlm(tbl,'o ~ t');

coef1 = table2cell(mdl1.Coefficients(:,1));

[rho0,rho1,rho2] = deal(coef1{:})

coef2 = table2cell(mdl2.Coefficients(:,1));

[I0,I1] = deal(coef2{:})

coef3 = table2cell(mdl3.Coefficients(:,1));

[w0,w1] = deal(coef3{:})

coef4 = table2cell(mdl4.Coefficients(:,1));

[o0,o1] = deal(coef4{:})