題目

f(x)=-2x^4+3x^3+2x^2-x+5=0

程式碼

clear all
z=-2:0.01:1
A=unidrnd(5)-3
x=A
p=[1,3,2,-1,5]
y=polyval(p,z)
q=polyder(p)
r = roots(p)
for i=1:10000 
y1=polyval(q,x)
x1=x-y1/1000 %%1000這個數當作步伐
x=x1
end
plot(z,y,x,polyval(p,x),'x')

結果：

數據1：

當ｘ的初始值為２時，起始點在右半邊，所以求到的局部最佳點為右邊的凹槽x=0.1754

數據2：

當ｘ的初始值為-２時，起始點在左半邊，所以求到的局部最佳點為左邊的凹槽x=-1.4254