请教，特定函数的非先线性拟合

fxp

对于离散数据y =[0.5010    0.4021     0.3033     0.2048    0.1063     0.0066    -0.0955    -0.2020    -0.3147  ];
x=[1.2000     1.4000     1.6000     1.8000     2.0000     2.2000     2.4000     2.6000     2.8000];  
   想用l函数

                               
登录/注册后可看大图

在x变换范围内进行拟合，其中a等于0。不知我的算法错在了那里？请高手指点一个。

function f=myfun2(t,a)
f=sqrt(2./(t.*pi))*cos(t-pi/4)-a;

for i=1:length(y)
    cn=y(i)
    [zz,err]=lsqnonlin(@(t) myfun2(t,cn),x(i));
    R(i)=real(zz);
end

按理说结果所求的R(i)应该和y(i)的值相当，但实际计算结果，却差异很大。不知错误在哪呢人，请高手指点！