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

onesimus

牛顿法优化程序代码如下：
function [x,minf] = minNewton(f,x0,eps)
2 B+ N# F. D2 k5 H  Tformat long;
/ `, D4 y2 j8 Rif nargin == 2
    eps = 1.0e-6;
end
* V( X/ u8 ?. E3 y$ z% U' D# f
" j& g, G& ]9 n$ H3 ?# _& ldf = diff(f);
, W$ z) P7 B% k) \) @, i9 J& sd2f = diff(df);
6 B# @4 g: n# S3 A$ nk = 0;
tol = 1;

while tol>eps
% V0 F, c: o" X    dfx = subs(df,findsym(df),x0);
    if diff(d2f) == 0
        d2fx = double(d2f);
! H6 ?; c, `$ K) h8 a3 H    else
        d2fx = subs(d2f,findsym(d2f),x0);
2 h, L, H3 p! c& ^  d' `    end
) t: b6 v- g& J& [" J8 ~6 i    x1 = x0 - dfx/d2fx;
    k = k + 1;
9 }( d2 a1 K+ P) e4 g0 v1 f8 H1 N    tol = abs(dfx);
) R& ]3 X% X5 m$ p) D) x7 M. c- o    x0 = x1;
end

, P$ C2 }0 H1 Fx = x1;
' P  i4 f# f- r8 S$ P) x" V9 ~minf =  subs(f,findsym(f),x);
$ d6 C% J* C" W+ m) n# vformat short;
复制代码
3 S8 W4 p8 I' q; _' G8 U
主程序代码如下：
# G, _! D1 f3 ]4 J9 G5 {7 hsyms t
: n3 y6 J7 Q; C- A+ L3 tf=t^2-log(t)-5
* p4 [; |1 x# i+ M5 tx=minNewton(f,2)
; A- n6 ?0 P5 u, F复制代码
/ {9 X' M: Z6 h+ W1 V4 e
- U7 d  M/ |5 S7 S( {6 F  i* U运行主程序后，得到的结果如下：
4 ?) ~! q, E' n5 a7 L1 x>> minNewton_ex

f =

t^2 - log(t) - 5


7 [! w& U, l( i, b: f5 c3 `. \) lx =

7273770929600966997213123375563852928/10286665498236842695784281251177659137

9 M. Q1 D2 m- u! v7 W& u0 R
1 c. U7 q; j9 [& l: I% z5 H  _x的值怎么会是这样，不应该是一个小数吗?

gaoxings

学习一下。

tuoni

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

; C3 Z' r3 _; }3 F/ H

IRESSEMT

shelby 发表于 2020-5-7 13:49
- i) w' S- B6 R2 `; N7 U你可以用vpa(x)或者double(x)转成小数。

; s% O, m7 @' K( ]) X5 ]- q5 M嗯嗯，楼主可以试试这个方法。
$ Y4 a3 F6 D' N' I/ X" e* \! e! U/ J

shelby

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

shelby

默认求得是符号解