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

onesimus

6 a/ B7 h( x/ o: p6 e7 `) T牛顿法优化程序代码如下：
: `1 _0 ?4 X; _: B& V7 }( l9 N+ @0 }function [x,minf] = minNewton(f,x0,eps)
format long;
if nargin == 2
& x! X) e! e: s2 \- j- n    eps = 1.0e-6;
: Q: _* A4 Y' M. ^4 Yend
3 ]! |. K' o) u  G
) y: r8 l9 S) }df = diff(f);
- A7 J  w' h& q5 L* }d2f = diff(df);
1 j* ~3 t; ~% @9 X, }k = 0;
tol = 1;
. f1 l! E$ g; _. ^* m9 o8 V% g
while tol>eps
    dfx = subs(df,findsym(df),x0);
    if diff(d2f) == 0
6 o' {- |8 p, w1 u1 l# }* N+ Y        d2fx = double(d2f);
    else
0 H: [9 W* R: g        d2fx = subs(d2f,findsym(d2f),x0);
    end
4 D+ z* [; J; a% z9 B    x1 = x0 - dfx/d2fx;
3 Y+ K6 j5 p9 w2 a    k = k + 1;
% q8 |0 \- A9 o7 T    tol = abs(dfx);
    x0 = x1;
! J% J$ D  f( e% N1 F; s& @7 iend

" W# A; `9 A0 H1 @4 Hx = x1;
* s0 {6 t$ ]3 @' M9 S& Bminf =  subs(f,findsym(f),x);
format short;
0 A  ^9 r' w) \. ^/ Z! j复制代码

' d, c' Y; j  y  v主程序代码如下：
+ L# u" p2 T+ |syms t
f=t^2-log(t)-5
! y3 Z4 |! b# X' s! ~' C5 F) jx=minNewton(f,2)
( r8 e& ^4 U; D8 @4 ^复制代码
. V  h2 ?! i: A$ C/ v
* f3 w/ s5 @# B' K; b/ j& a; ~( k运行主程序后，得到的结果如下：
* R0 f8 N- [& s* E4 j>> minNewton_ex
! t$ H3 m+ |* ?, X
f =

t^2 - log(t) - 5

) _" D0 I% [6 O* h4 O) r2 e
8 z1 x# b9 D5 a% i" O) [" u+ ex =

7273770929600966997213123375563852928/10286665498236842695784281251177659137
, W& l- U: i* c/ `/ D3 w8 t. U* i
# c) W; S: K1 [- J! ]1 U
1 M& t+ g8 o1 F/ kx的值怎么会是这样，不应该是一个小数吗?

shelby

默认求得是符号解

shelby

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

IRESSEMT

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

3 o7 r/ K! W) `$ p9 u# ?7 M嗯嗯，楼主可以试试这个方法。

tuoni

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

) }2 F/ e( Q: [5 `2 o: \9 W
8 g" K: Q+ r0 @3 h" J

gaoxings

学习一下。