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

onesimus

0 B" L  B8 m( M% R; s4 p牛顿法优化程序代码如下：
9 O/ G6 t+ T& `8 t( wfunction [x,minf] = minNewton(f,x0,eps)
( k  v% {  r2 }# [5 Iformat long;
; H# f- T. O% mif nargin == 2
    eps = 1.0e-6;
7 n; o, g4 _9 ~" oend
7 u4 P; G$ Y* _* K0 ?( E
: k$ |; z4 X: i. pdf = diff(f);
% F6 U' _3 p  T1 P- }) K9 vd2f = diff(df);
  D3 h/ |5 M) i) G) n  I( C. uk = 0;
tol = 1;
1 a' _' _' a/ Z3 X- }) d8 j/ k* c
6 v) K6 |; g- }0 vwhile tol>eps
    dfx = subs(df,findsym(df),x0);
9 N; x$ }5 L+ `- @- E) k    if diff(d2f) == 0
        d2fx = double(d2f);
5 n4 j+ ^. J; C' u" m2 D    else
1 k$ _5 K' t' h6 _        d2fx = subs(d2f,findsym(d2f),x0);
1 t4 x; m  ^2 O% S2 B    end
6 c, m% `  i3 e' V0 O' q    x1 = x0 - dfx/d2fx;
* S9 n& r% S- T    k = k + 1;
    tol = abs(dfx);
    x0 = x1;
end
3 t* i* \8 \' N' H  l
x = x1;
- o3 k5 f; r. n# J2 L/ Iminf =  subs(f,findsym(f),x);
* \7 u+ Z4 G' A# m2 P2 jformat short;
7 l( ]5 z( K, h: ?复制代码

主程序代码如下：
/ x8 B5 w7 i5 o. {# m' W" R3 }9 Asyms t
f=t^2-log(t)-5
x=minNewton(f,2)
- X9 ~1 w+ C" F+ n复制代码

运行主程序后，得到的结果如下：
5 H6 A# n$ O7 }4 z- ^1 G>> minNewton_ex
0 S$ c9 E6 r( ~, {
9 ?8 c* V8 |0 R7 }; Y! Hf =
* f2 Y. X" Q( x  R
; U8 h# M7 o- c0 T2 K: D  ct^2 - log(t) - 5


x =
* F1 V: d, @+ x/ p; i
! Z7 G- {% T- V  g2 K( a- T9 f7273770929600966997213123375563852928/10286665498236842695784281251177659137
7 A) W/ G: Y5 b: P6 |
4 m3 c: f& L4 y2 j, b  i
0 A: X/ ~3 {8 P- v& wx的值怎么会是这样，不应该是一个小数吗?

shelby

默认求得是符号解

shelby

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

IRESSEMT

shelby 发表于 2020-5-7 13:49
) G' J. |, v$ I( c# R& ]你可以用vpa(x)或者double(x)转成小数。

嗯嗯，楼主可以试试这个方法。
3 b, p0 V. \5 z& v: K# d2 h- h

tuoni

shelby 发表于 2020-5-7 13:49
% ]$ Y. e+ D* E$ R4 a' {) U你可以用vpa(x)或者double(x)转成小数。

gaoxings

学习一下。