f.m

应用模拟退火方法求解非线形方程

function [s]=f(db)
%global x y z yd a b c
x=1;y=10;z=0;yd=sqrt(x^2+y^2+z^2);a=pi/2;b=pi/2;c=0;
s=(((y/yd)-cos(a)*cos(b+db)*(z/yd)/(sin(b+db)))^2)-((sin(a))^2)*(1-((z/yd)^2)/(sin((b+db))^2));