x的值怎么会是这样，不应该是一个小数吗?

onesimus · 发表于 2020-5-7 11:26

EDA365欢迎您登录！

您需要登录才可以下载或查看，没有帐号？注册

x

牛顿法优化程序代码如下：
function [x,minf] = minNewton(f,x0,eps)
format long;
if nargin == 2
eps = 1.0e-6;
end

df = diff(f);
d2f = diff(df);
k = 0;
tol = 1;

while tol>eps
dfx = subs(df,findsym(df),x0);
if diff(d2f) == 0
      d2fx = double(d2f);
else
      d2fx = subs(d2f,findsym(d2f),x0);
end
x1 = x0 - dfx/d2fx;
k = k + 1;
tol = abs(dfx);
x0 = x1;
end

x = x1;
minf =  subs(f,findsym(f),x);
format short;
复制代码

主程序代码如下：
syms t
f=t^2-log(t)-5
x=minNewton(f,2)
复制代码

运行主程序后，得到的结果如下：
>> minNewton_ex

f =

t^2 - log(t) - 5

x =

7273770929600966997213123375563852928/10286665498236842695784281251177659137

x的值怎么会是这样，不应该是一个小数吗?