查看: 4939|回复: 1

[代码资源] fprintf函数问题

字体大小: 正常放大

6 主题	11 听众	268 积分

升级 84%

TA的每日心情

	奋斗 2016-5-1 21:35

签到天数: 82 天

[LV.6]常住居民II

群组: 2015司守奎matlab培训

群组: C 语言讨论组

群组: Matlab讨论组

群组: LINGO

电梯直达

^#

发表于 2015-7-23 09:56 |只看该作者 |正序浏览

|招呼Ta 关注Ta |邮箱已经成功绑定

2体力

syms x y
f=x^3-y^3+3*x^2+3*y^2-9*x;
df=jacobian(f);  %求一阶偏导数
d2f=jacobian(df); %求Hessian阵
[xx,yy]=solve(df)  %求驻点
xx=double(xx);yy=double(yy); %转化成双精度浮点型数据，下面判断特征值的正负，必须是数值型数据
for i=1:length(xx)
a=subs(d2f,{x,y},{xx(i),yy(i)});
b=eig(a);  %求矩阵的特征值
f=subs(f,{x,y},{xx(i),yy(i)});
if all(b>0)
      fprintf('(%f,%f)是极小值点,对应的极小值为%f\n',xx(i),yy(i),f);
elseif all(b<0)
      fprintf('(%f,%f)是极大值点,对应的极大值为%f\n',xx(i),yy(i),f);
elseif any(b>0) & any(b<0)
      fprintf('(%f,%f)不是极值点\n',xx(i),yy(i));
else
      fprintf('无法判断(%f,%f)是否是极值点\n',xx(i),yy(i));
end
end

xx =

  1
-3
  1
-3

yy =

0
0
2
2

错误使用 fprintf
没有为 'sym' 输入定义函数。

能告诉我哪里出问题了吗？
我用的是2014aMATLAB

zan

fprintf问题

转播0 淘帖0 分享0 收藏0 支持0 反对0 微信

使用道具举报

士心之约

142 主题	50 听众	2万积分

升级 0%

TA的每日心情

	开心 2022-8-28 18:40

签到天数: 1187 天

[LV.10]以坛为家III

网络挑战赛参赛者

国际赛参赛者

群组: Matlab讨论组

群组: MATLAB技术交流

群组: 西南交通大学数学建模

群组: 学术交流A

群组: 数学建模

1^#

发表于 2015-8-23 09:20 |只看该作者 |招呼Ta 关注Ta |邮箱已经成功绑定

syms x y
% T0 i) K) A/ U5 U% `
f=x^3-y^3+3*x^2+3*y^2-9*x;. S\" s) Y0 u2 T- B% @
df=jacobian(f); %求一阶偏导数2 K5 Z- x1 Q! U( I8 K K9 Z- z9 `+ ^$ w( @. F8 [2 M8 ~
d2f=jacobian(df); %求Hessian阵
- H a. }7 }; a& }8 `) ]) u
[xx,yy]=solve(df) %求驻点/ b& c9 X# V- k Z a8 ]. O I9 }, n
xx=double(xx);yy=double(yy); %转化成双精度浮点型数据，下面判断特征值的正负，必须是数值型数据 d+ C1 y" F/ b+ b: i: W1 H$ j4 p. ^3 x# B* x( X* A5 J
for i=1:length(xx)
. [; M9 `$ v9 n3 X6 p6 g) d8 H
a=subs(d2f,{x,y},{xx(i),yy(i)});
, a) V6 {3 \- _9 k% v/ R
b=eig(a); %求矩阵的特征值0 $ \0 B5 c8 ?; d: k1 X
f=subs(f,{x,y},{xx(i),yy(i)});0 \' {2 r6 r9 D, {6 }9 n5 e8 f
if all(b>0)
) g$ \+ A6 I- y$ b( U0 a5 H1 p
fprintf('(%f,%f)是极小值点,对应的极小值为%f\n',xx(i),yy(i),double(f));
0 B! p% B* m/ f' A6 x: V8 a) N, H
elseif all(b<0)/ H) l& ^6 f' ]- V\" D. W9 M
fprintf('(%f,%f)是极大值点,对应的极大值为%f\n',xx(i),yy(i),double(f));1 V7 N3 C, z1 x+ d
elseif any(b>0) & any(b<0)
}* p7 s; I. ?( Z9 y, ?2 M7 a7 _
fprintf('(%f,%f)不是极值点\n',xx(i),yy(i));4 N. G\" M: B% x
else J- B) ?& ]* w& ]( _8 T* h8 }
fprintf('无法判断(%f,%f)是否是极值点\n',xx(i),yy(i)); 8 F& M0 k; p( [7 _+ v! M2 ?
end
4 {- D) j4 a; _* {5 I/ `+ o1 z- w
end