題目

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:

數據1

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

結果圖1

數據2:

數據2

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

結果圖2

Last modification:December 4, 2018
If you think my article is useful to you, please feel free to appreciate