My Matlab script is as follows (which works perfectly but i have to change it to GUI for user readiness)
R_pp=[]; Vf_pp=[];
for i=1:Nmeas_pp Vf_pp=vertcat(Vf_pp,x_pp(i),y_pp(i)); k=B_pp(i)/(2*rigidity); % m^(-1)
M11_pp=cos(k*l)^2*(cos(k*l)^2 - D*k*cos(k*l)*sin(k*l)) + k*sin(k*l)^2*(sin(k*l)^2/k + D*cos(k*l)*sin(k*l)) - cos(k*l)*sin(k*l)*(cos(k*l)*sin(k*l) - D*k*sin(k*l)^2) - k*cos(k*l)*sin(k*l)*(D*cos(k*l)^2 + (cos(k*l)*sin(k*l))/k); M12_pp=cos(k*l)^2*(D*cos(k*l)^2 + (cos(k*l)*sin(k*l))/k) - cos(k*l)*sin(k*l)*(sin(k*l)^2/k + D*cos(k*l)*sin(k*l)) - (sin(k*l)^2*(cos(k*l)*sin(k*l) - D*k*sin(k*l)^2))/k + (cos(k*l)*sin(k*l)*(cos(k*l)^2 - D*k*cos(k*l)*sin(k*l)))/k; M13_pp=cos(k*l)^2*(cos(k*l)*sin(k*l) - D*k*sin(k*l)^2) - k*sin(k*l)^2*(D*cos(k*l)^2 + (cos(k*l)*sin(k*l))/k) + cos(k*l)*sin(k*l)*(cos(k*l)^2 - D*k*cos(k*l)*sin(k*l)) - k*cos(k*l)*sin(k*l)*(sin(k*l)^2/k + D*cos(k*l)*sin(k*l)); M14_pp=cos(k*l)^2*(sin(k*l)^2/k + D*cos(k*l)*sin(k*l)) + cos(k*l)*sin(k*l)*(D*cos(k*l)^2 + (cos(k*l)*sin(k*l))/k) + (sin(k*l)^2*(cos(k*l)^2 - D*k*cos(k*l)*sin(k*l)))/k + (cos(k*l)*sin(k*l)*(cos(k*l)*sin(k*l) - D*k*sin(k*l)^2))/k; M31_pp=k*sin(k*l)^2*(D*cos(k*l)^2 + (cos(k*l)*sin(k*l))/k) - cos(k*l)^2*(cos(k*l)*sin(k*l) - D*k*sin(k*l)^2) - cos(k*l)*sin(k*l)*(cos(k*l)^2 - D*k*cos(k*l)*sin(k*l)) + k*cos(k*l)*sin(k*l)*(sin(k*l)^2/k + D*cos(k*l)*sin(k*l)); M32_pp= -cos(k*l)^2*(sin(k*l)^2/k + D*cos(k*l)*sin(k*l)) - cos(k*l)*sin(k*l)*(D*cos(k*l)^2 + (cos(k*l)*sin(k*l))/k) - (sin(k*l)^2*(cos(k*l)^2 - D*k*cos(k*l)*sin(k*l)))/k - (cos(k*l)*sin(k*l)*(cos(k*l)*sin(k*l) - D*k*sin(k*l)^2))/k; M33_pp=cos(k*l)^2*(cos(k*l)^2 - D*k*cos(k*l)*sin(k*l)) + k*sin(k*l)^2*(sin(k*l)^2/k + D*cos(k*l)*sin(k*l)) - cos(k*l)*sin(k*l)*(cos(k*l)*sin(k*l) - D*k*sin(k*l)^2) - k*cos(k*l)*sin(k*l)*(D*cos(k*l)^2 + (cos(k*l)*sin(k*l))/k); M34_pp=cos(k*l)^2*(D*cos(k*l)^2 + (cos(k*l)*sin(k*l))/k) - cos(k*l)*sin(k*l)*(sin(k*l)^2/k + D*cos(k*l)*sin(k*l)) - (sin(k*l)^2*(cos(k*l)*sin(k*l) - D*k*sin(k*l)^2))/k + (cos(k*l)*sin(k*l)*(cos(k*l)^2 - D*k*cos(k*l)*sin(k*l)))/k; Rlinex_pp= [M11_pp, M12_pp, M13_pp, M14_pp]; Rliney_pp= [M31_pp, M32_pp, M33_pp, M34_pp]; R_pp=vertcat(R_pp,Rlinex_pp,Rliney_pp); endR_mp=[];Vf_mp=[];for i=1:Nmeas_mp Vf_mp=vertcat(Vf_mp,x_mp(i),y_mp(i)); k=B_mp(i)/(2*rigidity); % m^(-1) M11_mp=cos(k*l)^2*(cos(k*l)^2 - D*k*cos(k*l)*sin(k*l)) - k*sin(k*l)^2*(sin(k*l)^2/k + D*cos(k*l)*sin(k*l)) + cos(k*l)*sin(k*l)*(cos(k*l)*sin(k*l) - D*k*sin(k*l)^2) - k*cos(k*l)*sin(k*l)*(D*cos(k*l)^2 + (cos(k*l)*sin(k*l))/k); M12_mp=cos(k*l)^2*(D*cos(k*l)^2 + (cos(k*l)*sin(k*l))/k) + cos(k*l)*sin(k*l)*(sin(k*l)^2/k + D*cos(k*l)*sin(k*l)) + (sin(k*l)^2*(cos(k*l)*sin(k*l) - D*k*sin(k*l)^2))/k + (cos(k*l)*sin(k*l)*(cos(k*l)^2 - D*k*cos(k*l)*sin(k*l)))/k; M13_mp=cos(k*l)^2*(cos(k*l)*sin(k*l) - D*k*sin(k*l)^2) + k*sin(k*l)^2*(D*cos(k*l)^2 + (cos(k*l)*sin(k*l))/k) - cos(k*l)*sin(k*l)*(cos(k*l)^2 - D*k*cos(k*l)*sin(k*l)) - k*cos(k*l)*sin(k*l)*(sin(k*l)^2/k + D*cos(k*l)*sin(k*l)); M14_mp=cos(k*l)^2*(sin(k*l)^2/k + D*cos(k*l)*sin(k*l)) - cos(k*l)*sin(k*l)*(D*cos(k*l)^2 + (cos(k*l)*sin(k*l))/k) - (sin(k*l)^2*(cos(k*l)^2 - D*k*cos(k*l)*sin(k*l)))/k + (cos(k*l)*sin(k*l)*(cos(k*l)*sin(k*l) - D*k*sin(k*l)^2))/k; M31_mp=cos(k*l)*sin(k*l)*(cos(k*l)^2 - D*k*cos(k*l)*sin(k*l)) - k*sin(k*l)^2*(D*cos(k*l)^2 + (cos(k*l)*sin(k*l))/k) - cos(k*l)^2*(cos(k*l)*sin(k*l) - D*k*sin(k*l)^2) + k*cos(k*l)*sin(k*l)*(sin(k*l)^2/k + D*cos(k*l)*sin(k*l)); M32_mp=cos(k*l)*sin(k*l)*(D*cos(k*l)^2 + (cos(k*l)*sin(k*l))/k) - cos(k*l)^2*(sin(k*l)^2/k + D*cos(k*l)*sin(k*l)) + (sin(k*l)^2*(cos(k*l)^2 - D*k*cos(k*l)*sin(k*l)))/k - (cos(k*l)*sin(k*l)*(cos(k*l)*sin(k*l) - D*k*sin(k*l)^2))/k; M33_mp=cos(k*l)^2*(cos(k*l)^2 - D*k*cos(k*l)*sin(k*l)) - k*sin(k*l)^2*(sin(k*l)^2/k + D*cos(k*l)*sin(k*l)) + cos(k*l)*sin(k*l)*(cos(k*l)*sin(k*l) - D*k*sin(k*l)^2) - k*cos(k*l)*sin(k*l)*(D*cos(k*l)^2 + (cos(k*l)*sin(k*l))/k); M34_mp=cos(k*l)^2*(D*cos(k*l)^2 + (cos(k*l)*sin(k*l))/k) + cos(k*l)*sin(k*l)*(sin(k*l)^2/k + D*cos(k*l)*sin(k*l)) + (sin(k*l)^2*(cos(k*l)*sin(k*l) - D*k*sin(k*l)^2))/k + (cos(k*l)*sin(k*l)*(cos(k*l)^2 - D*k*cos(k*l)*sin(k*l)))/k; Rlinex_mp= [M11_mp, M12_mp, M13_mp, M14_mp]; Rliney_mp= [M31_mp, M32_mp, M33_mp, M34_mp]; R_mp=vertcat(R_mp,Rlinex_mp,Rliney_mp); endlsqr(R_pp, Vf_pp)lsqr(R_mp, Vf_mp)
Best Answer