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

onesimus

$ _" a9 S# F8 D6 L4 v9 r牛顿法优化程序代码如下：
function [x,minf] = minNewton(f,x0,eps)
5 j& z9 A* }$ ~% _format long;
* E* L0 _7 g6 {- Z, zif nargin == 2
    eps = 1.0e-6;
end

; g$ V% u1 v$ k2 J$ Udf = diff(f);
2 J, ?& @9 X! f; Sd2f = diff(df);
k = 0;
/ {4 D" {2 |% Gtol = 1;

6 k# v$ ?* F/ R. i# w2 Uwhile tol>eps
, L! c) g, |# s/ H- }    dfx = subs(df,findsym(df),x0);
) c* a' H4 {  A  X6 x    if diff(d2f) == 0
        d2fx = double(d2f);
    else
        d2fx = subs(d2f,findsym(d2f),x0);
- F7 b9 s/ {- J8 t2 I1 `$ Q# O    end
    x1 = x0 - dfx/d2fx;
    k = k + 1;
    tol = abs(dfx);
    x0 = x1;
end
$ Z2 a, A) D* j2 I0 z7 e% U
x = x1;
minf =  subs(f,findsym(f),x);
1 o( D& v, S5 k: Z5 ^format short;
8 t# y9 t% U6 _  r+ b复制代码

) i$ T, F; x  Q  i! X  m主程序代码如下：
syms t
f=t^2-log(t)-5
x=minNewton(f,2)
复制代码

运行主程序后，得到的结果如下：
>> minNewton_ex
; ?- F$ Z. P1 @2 L5 p; Z: G2 f
f =

t^2 - log(t) - 5

& I1 w$ ]9 o; @8 l* ?8 [2 X
x =

7273770929600966997213123375563852928/10286665498236842695784281251177659137
" y8 r; @1 o& X: m0 ?4 v5 j7 l

+ J+ h9 E7 ?8 p3 ^4 F3 `! X8 cx的值怎么会是这样，不应该是一个小数吗?

shelby

默认求得是符号解

shelby

你可以用vpa(x)或者double(x)转成小数。

IRESSEMT

shelby 发表于 2020-5-7 13:49
你可以用vpa(x)或者double(x)转成小数。

嗯嗯，楼主可以试试这个方法。
% c6 R; `4 J. P0 W9 E

tuoni

shelby 发表于 2020-5-7 13:49
你可以用vpa(x)或者double(x)转成小数。

9 J3 F! F  ~9 B

gaoxings

学习一下。