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

onesimus

- [; c# v. x. B  D8 t& u* R/ T. g( F牛顿法优化程序代码如下：
function [x,minf] = minNewton(f,x0,eps)
) @0 ~5 \9 l5 ]) q9 hformat long;
if nargin == 2
6 R- w" U9 s4 c7 `2 Q    eps = 1.0e-6;
" z) A$ {4 C5 O3 b0 Q  {end

df = diff(f);
d2f = diff(df);
k = 0;
- m' a" X; p$ L5 M% ]( Ltol = 1;
9 l* ]7 I5 v9 H& ?+ M4 {( b" i4 M4 {
0 p( Y3 o, o; B; q  V1 m. f! m3 bwhile tol>eps
    dfx = subs(df,findsym(df),x0);
    if diff(d2f) == 0
  U0 L1 I6 P; j* X        d2fx = double(d2f);
1 `" o8 I8 X5 _4 u0 f. E& w    else
0 x6 B2 t; G! k9 n3 ]' w        d2fx = subs(d2f,findsym(d2f),x0);
    end
7 Y+ ]7 i: u0 \. ^0 G    x1 = x0 - dfx/d2fx;
2 A  L0 G! J+ w# L# u1 ~1 y    k = k + 1;
0 ], \4 Q$ L8 L4 f0 g9 C/ D9 m2 U0 F/ m    tol = abs(dfx);
' k; o, k/ L( \/ j    x0 = x1;
end
' T6 y) \4 Z* M7 ?; }( P8 S7 G
x = x1;
minf =  subs(f,findsym(f),x);
format short;
复制代码

* p4 a& f) G( U2 r7 j主程序代码如下：
/ N( O  n' S4 ?; g; `4 Qsyms t
f=t^2-log(t)-5
x=minNewton(f,2)
复制代码
) `) _+ d" M, P- X5 ^" V
运行主程序后，得到的结果如下：
, E8 d* z9 s0 D/ R# o>> minNewton_ex

, p" D: i3 S4 C0 W# \, p4 Qf =
* ^  T7 R; Q3 T( W6 Y/ L  Z- O; [
" A9 H% x$ D# D& r2 ?0 Wt^2 - log(t) - 5


x =
! ?+ S# `9 [8 U1 ?
7273770929600966997213123375563852928/10286665498236842695784281251177659137
* {' H% O! V( l( N( {+ ~
% @1 n+ d( R: R: `
4 M- v4 m, Y* G# g- n& s2 {x的值怎么会是这样，不应该是一个小数吗?

shelby

默认求得是符号解

shelby

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

IRESSEMT

shelby 发表于 2020-5-7 13:49
7 w: [6 f7 k; @你可以用vpa(x)或者double(x)转成小数。

) _' X7 p4 r. O, l8 X" G5 t嗯嗯，楼主可以试试这个方法。
* J+ T7 J' f3 n' |; T' M

tuoni

shelby 发表于 2020-5-7 13:49
/ F( ~& _$ p2 O6 P% n! A" P  S你可以用vpa(x)或者double(x)转成小数。

$ s* n. v4 l- V$ `3 Z& G

gaoxings

学习一下。