對(duì)于求平方根,變成方程模式為f(x)=x^2-a,即求此方程的實(shí)根；
下面編寫(xiě)了兩個(gè)function函數(shù),可以直接調(diào)用.
二分法：
function x=sqrt_bisect(a)
f=@(x)x^2-a;
if a0
xb=x;
elseif f(xa)*f(x)>0
xa=x;
else
break
end
end
end
x;
牛頓迭代法：
function x=sqrt_newton(a)
f=@(x)x^2-a;
df=diff(sym('x^2-a'));
if a1e-6
x0=x1;
x1=x0-f(x0)/subs(df,x0);
end
end
x=x1;
調(diào)用格式為：
sqrt_bisect(3)
ans =
1.7321
或者
sqrt_newton(2)
ans =
1.4142