fj2是一个5212*180的已知矩阵
我想用最小二乘法求超定方程的最优解、但是好像最后一步跑不出，想问一下代码哪里错了？？

syms t th d k;
f=@(t,th,d,k)sqrt((-55*sin(2*th*abs(t-k*d)))*(-55*sin(2*th*abs(t-k*d)))...
    -(4*55*sin(th)*sin(th)+9)*(9*abs(t-k*d).^2+55*abs(t-k*d).^2*cos(th).^2-64*9*25))...
    /(4*55*sin(th)*sin(th)+9);
n=8;
y=find(fj2(:,1));
y=y(1:n,1);
R=fj2(y,1);%获取所有R
c=cell(n,1);
for k=1:n
    c{k,1}=(f(t,th,d,k)-R(k,1)).^2;
end
F=c{1,1}+c{2,1}+c{3,1}+c{4,1}+c{5,1}+c{6,1}+c{7,1}+c{8,1};
[a,fval]=fminsearch(F,[0 0 0]');