- 在线时间
- 48 小时
- 最后登录
- 2017-1-23
- 注册时间
- 2015-6-7
- 听众数
- 11
- 收听数
- 2
- 能力
- 0 分
- 体力
- 558 点
- 威望
- 0 点
- 阅读权限
- 30
- 积分
- 268
- 相册
- 0
- 日志
- 0
- 记录
- 0
- 帖子
- 195
- 主题
- 6
- 精华
- 0
- 分享
- 0
- 好友
- 29
升级   84% TA的每日心情 | 奋斗 2016-5-1 21:35 |
|---|
签到天数: 82 天 [LV.6]常住居民II
 群组: 2015司守奎matlab培训 群组: C 语言讨论组 群组: Matlab讨论组 群组: LINGO |
2体力
syms x y2 v, H& f* M' m0 y; u1 W
f=x^3-y^3+3*x^2+3*y^2-9*x;
2 \! O4 Y7 }. H: Pdf=jacobian(f); %求一阶偏导数3 B5 @) ]/ k4 R) f; l. H
d2f=jacobian(df); %求Hessian阵+ p/ h& v8 J1 L0 X: ] R1 t/ o
[xx,yy]=solve(df) %求驻点; L6 w% z6 D# M( d; x- r
xx=double(xx);yy=double(yy); %转化成双精度浮点型数据,下面判断特征值的正负,必须是数值型数据
1 g2 {% {" G0 X6 Ufor i=1:length(xx)
0 y8 o- E3 ?0 \$ @" G) { a=subs(d2f,{x,y},{xx(i),yy(i)}); + M; A& e; x d6 ~# m5 ]
b=eig(a); %求矩阵的特征值- w" ^$ W6 F: |( h8 x2 W2 Z
f=subs(f,{x,y},{xx(i),yy(i)});8 d8 F7 Q$ e" L- b; b' i" _
if all(b>0)3 K' }( w; C' j7 @
fprintf('(%f,%f)是极小值点,对应的极小值为%f\n',xx(i),yy(i),f);2 q* M% E d" ?8 k% Z8 }4 `8 B
elseif all(b<0)
& Z; Q, ~! f" W& H5 _4 l fprintf('(%f,%f)是极大值点,对应的极大值为%f\n',xx(i),yy(i),f);" S- |1 W8 w" d) M/ q# \: x
elseif any(b>0) & any(b<0)9 Q1 T" I. t: a8 |+ x0 K& o- P* w
fprintf('(%f,%f)不是极值点\n',xx(i),yy(i));
# [4 q; [$ _- w9 N else
: V$ u* A# F) p, e& r0 n fprintf('无法判断(%f,%f)是否是极值点\n',xx(i),yy(i));
/ g+ ?( w. R# S2 e/ `0 A end
( u7 b5 s$ k8 Lend
* L( Q3 Q" y. K8 s9 B1 l8 F/ }$ q
7 {: v, \$ ]! S" Y' ]3 a6 ] , I5 }$ ?/ Z5 X6 o0 [: p* k
xx =3 ?6 K# Y) g8 |% m# T; J$ D! l6 b5 [2 T
+ c( |" C9 `8 w. t# r6 T O
1
, E4 S' {4 w4 n# h/ E' R -30 Y3 G% P9 D+ j3 ~3 t, Q) Y
1
2 T% y+ n& m+ _ -3
( M2 `0 z( S$ A4 E' r3 x% {
0 B2 k* Q _, u- C
* C h5 v4 z( O0 U9 L2 Ryy =
! N3 d' j6 z) v2 H . S1 }/ v6 L, U( K
0
; |& I4 F' `* s4 @1 k 0% A- X& [5 u* \8 c# a1 ~
2; @5 ^0 l8 y% Q V5 Y/ W* ]
2" u4 H/ j( u: J F: M0 \
; b: @4 k: {$ L5 q2 y0 A1 o错误使用 fprintf
) ^7 @1 s `& e没有为 'sym' 输入定义函数。8 [. ^9 j8 p4 Z0 I2 N" s6 r) q
5 g% G+ C4 t: E- N/ l
能告诉我哪里出问题了吗?: D8 Z1 I- e# }) K5 F3 b7 n/ p
我用的是2014aMATLAB
3 ~6 j4 x2 t$ k# Q |
zan
|