有人帮忙么。我贴出程序
k1=2;k2=3;k3=4;k4=5;k5=6;k6=7;
a=0.45;
x1=-6.5;x2=-5.5;x3=-4.5;x4=-3.5;x5=0;
S11=k1*cos(a)*x1;S12=k1*cos(a)*x2;S13=k1*cos(a)*x3;S14=k1*cos(a)*x4;
S21=k2*cos(a)*x1;S22=k2*cos(a)*x2;S23=k2*cos(a)*x3;S24=k2*cos(a)*x4;
S31=k3*cos(a)*x1;S32=k3*cos(a)*x2;S33=k3*cos(a)*x3;S34=k3*cos(a)*x4;
S41=k4*cos(a)*x1;S42=k4*cos(a)*x2;S43=k4*cos(a)*x3;S44=k4*cos(a)*x4;
S51=k5*cos(a)*x1;S52=k5*cos(a)*x2;S53=k5*cos(a)*x3;S54=k5*cos(a)*x4;
S61=k6*cos(a)*x1;S62=k6*cos(a)*x2;S63=k6*cos(a)*x3;S64=k6*cos(a)*x4;
y11=2;y12=3;y13=4;y14=5;
y21=2.1;y22=3.1;y23=4.1;y24=5.1;
y31=2.2;y32=3.2;y33=4.2;y34=5.2;
y41=2.3;y42=3.3;y43=4.3;y44=5.3;
y51=2.4;y52=3.4;y53=4.4;y54=5.4;
y61=2.5;y62=3.5;y63=4.5;y64=5.5;
B1=[y11;y12;y13;y14];B2=[y21;y22;y23;y24];B3=[y31;y32;y33;y34];B4=[y41;y42;y43;y44];B5=[y51;y52;y53;y54];B6=[y61;y62;y63;y64];
Z1=[exp(i*S11) exp(i*S11); exp(i*S12) exp(i*S12); exp(i*S13) exp(i*S13); exp(i*S14) exp(i*S14)];
Z2=[exp(i*S21) exp(i*S21); exp(i*S22) exp(i*S22); exp(i*S23) exp(i*S23); exp(i*S24) exp(i*S24)];
Z3=[exp(i*S31) exp(i*S31); exp(i*S32) exp(i*S32); exp(i*S33) exp(i*S33); exp(i*S34) exp(i*S34)];
Z4=[exp(i*S41) exp(i*S41); exp(i*S42) exp(i*S42); exp(i*S43) exp(i*S43); exp(i*S44) exp(i*S44)];
Z5=[exp(i*S51) exp(i*S51); exp(i*S52) exp(i*S52); exp(i*S53) exp(i*S53); exp(i*S54) exp(i*S54)];
Z6=[exp(i*S61) exp(i*S61); exp(i*S62) exp(i*S62); exp(i*S63) exp(i*S63); exp(i*S64) exp(i*S64)];
A1=pinv(Z1.'*Z1)*(Z1.')*B1
k1=2;k2=3;k3=4;k4=5;k5=6;k6=7;
a=0.45;
x1=-6.5;x2=-5.5;x3=-4.5;x4=-3.5;x5=0;
S11=k1*cos(a)*x1;S12=k1*cos(a)*x2;S13=k1*cos(a)*x3;S14=k1*cos(a)*x4;
S21=k2*cos(a)*x1;S22=k2*cos(a)*x2;S23=k2*cos(a)*x3;S24=k2*cos(a)*x4;
S31=k3*cos(a)*x1;S32=k3*cos(a)*x2;S33=k3*cos(a)*x3;S34=k3*cos(a)*x4;
S41=k4*cos(a)*x1;S42=k4*cos(a)*x2;S43=k4*cos(a)*x3;S44=k4*cos(a)*x4;
S51=k5*cos(a)*x1;S52=k5*cos(a)*x2;S53=k5*cos(a)*x3;S54=k5*cos(a)*x4;
S61=k6*cos(a)*x1;S62=k6*cos(a)*x2;S63=k6*cos(a)*x3;S64=k6*cos(a)*x4;
y11=2;y12=3;y13=4;y14=5;
y21=2.1;y22=3.1;y23=4.1;y24=5.1;
y31=2.2;y32=3.2;y33=4.2;y34=5.2;
y41=2.3;y42=3.3;y43=4.3;y44=5.3;
y51=2.4;y52=3.4;y53=4.4;y54=5.4;
y61=2.5;y62=3.5;y63=4.5;y64=5.5;
B1=[y11;y12;y13;y14];B2=[y21;y22;y23;y24];B3=[y31;y32;y33;y34];B4=[y41;y42;y43;y44];B5=[y51;y52;y53;y54];B6=[y61;y62;y63;y64];
Z1=[exp(i*S11) exp(i*S11); exp(i*S12) exp(i*S12); exp(i*S13) exp(i*S13); exp(i*S14) exp(i*S14)];
Z2=[exp(i*S21) exp(i*S21); exp(i*S22) exp(i*S22); exp(i*S23) exp(i*S23); exp(i*S24) exp(i*S24)];
Z3=[exp(i*S31) exp(i*S31); exp(i*S32) exp(i*S32); exp(i*S33) exp(i*S33); exp(i*S34) exp(i*S34)];
Z4=[exp(i*S41) exp(i*S41); exp(i*S42) exp(i*S42); exp(i*S43) exp(i*S43); exp(i*S44) exp(i*S44)];
Z5=[exp(i*S51) exp(i*S51); exp(i*S52) exp(i*S52); exp(i*S53) exp(i*S53); exp(i*S54) exp(i*S54)];
Z6=[exp(i*S61) exp(i*S61); exp(i*S62) exp(i*S62); exp(i*S63) exp(i*S63); exp(i*S64) exp(i*S64)];
A1=pinv(Z1.'*Z1)*(Z1.')*B1