QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 4937|回复: 1
打印 上一主题 下一主题

[代码资源] fprintf函数问题

[复制链接]
字体大小: 正常 放大

6

主题

11

听众

268

积分

升级  84%

  • TA的每日心情
    奋斗
    2016-5-1 21:35
  • 签到天数: 82 天

    [LV.6]常住居民II

    邮箱绑定达人 社区QQ达人

    群组2015司守奎matlab培训

    群组C 语言讨论组

    群组Matlab讨论组

    群组LINGO

    跳转到指定楼层
    1#
    发表于 2015-7-23 09:56 |只看该作者 |倒序浏览
    |招呼Ta 关注Ta |邮箱已经成功绑定
    2体力
    syms x y! \: T8 d/ H6 }0 n
    f=x^3-y^3+3*x^2+3*y^2-9*x;
    " D; y% V8 T; R" ?3 Ldf=jacobian(f);  %求一阶偏导数6 l# x! S9 g; |+ g
    d2f=jacobian(df); %求Hessian阵0 h: \& M* D9 H2 X# T1 P
    [xx,yy]=solve(df)  %求驻点
    # j# k, X* B" s5 g: A9 ^* {- `! nxx=double(xx);yy=double(yy); %转化成双精度浮点型数据,下面判断特征值的正负,必须是数值型数据
      C& A( D9 U2 @& F4 G: ofor i=1:length(xx)
      S! U" e, W) m, U( }    a=subs(d2f,{x,y},{xx(i),yy(i)});  
    ) E: B* i! W! Y6 ]2 C3 @" U! K: B    b=eig(a);  %求矩阵的特征值
    % R6 O- S, K( X- O- u1 M    f=subs(f,{x,y},{xx(i),yy(i)});* z) E* `7 H  r( r2 A
        if all(b>0)
    + z/ }: Z. e# [' v        fprintf('(%f,%f)是极小值点,对应的极小值为%f\n',xx(i),yy(i),f);
    5 D9 f! \5 k3 Q* r    elseif all(b<0)
    + X' O8 r; O  v4 ~* ~/ L        fprintf('(%f,%f)是极大值点,对应的极大值为%f\n',xx(i),yy(i),f);3 l6 V* {8 @' }+ p- T: l' L
        elseif any(b>0) & any(b<0)
    " r' X( U, x: X. J- s        fprintf('(%f,%f)不是极值点\n',xx(i),yy(i));! a" V: B8 p" l8 A
        else
    : G, H, s2 n4 t3 a) U. `5 m        fprintf('无法判断(%f,%f)是否是极值点\n',xx(i),yy(i));  * E8 ]" R8 K1 D! A9 q
        end
    1 V9 o. h( B& Q* ^' p0 Eend
    3 i% ^9 c8 c: M6 G9 W, Z0 Q8 s$ E5 \3 i  N$ Y' F8 Y1 S

    * A% P* ?1 w- }0 N5 u5 r" Bxx =2 A3 X; R# E7 x8 I) J. a; I
    4 F* d! S, Y) }8 N' N( Q# e: l  W' U
      1  ]; W4 Z- ^" I; x0 Q" f6 o6 Y
    -3
      d9 [& y; ]5 H) E: p  1
    ! f8 s! [6 C% x, m0 @! u5 M -3+ D: l' ~; Y2 z5 R) y3 y

    8 U" l( m1 `" F7 k, ~6 A4 f# h
    & W. O8 x2 T& T' R6 `; Zyy =2 X# }7 \4 q$ N. Q

    5 x" a; @* K1 ^ 0" B4 s& A) h) A  P9 y+ Q. x
    0
    + X. |( z9 W' h5 v$ ` 2/ Y7 D8 o* j- ]7 o
    2, l! g) c, l  G4 r+ u

    * S$ _6 t0 u: D7 x- M错误使用 fprintf
    % S5 n4 x1 Y: C2 J+ K6 l: [+ _没有为 'sym' 输入定义函数。
    & C+ x& x- n* r! y% Z5 A
    $ c5 h- X% @/ }. ?, o能告诉我哪里出问题了吗?
    # B! ?: n3 q+ b# O我用的是2014aMATLAB
    * v( \: f0 e  N3 R6 O, B

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

    142

    主题

    50

    听众

    2万

    积分

    升级  0%

  • TA的每日心情
    开心
    2022-8-28 18:40
  • 签到天数: 1187 天

    [LV.10]以坛为家III

    网络挑战赛参赛者

    国际赛参赛者

    邮箱绑定达人 社区QQ达人 元老勋章 发帖功臣 新人进步奖 原创写作奖 最具活力勋章 风雨历程奖

    群组Matlab讨论组

    群组MATLAB技术交流

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

    群组学术交流A

    群组数学建模

    1. syms x y
      5 P\" ^4 Z, I1 X5 e
    2. f=x^3-y^3+3*x^2+3*y^2-9*x;
      ! A& l* Z) M' }3 J( z( N& Y6 ^4 h
    3. df=jacobian(f);  %求一阶偏导数2 K5 Z- x1 Q! U( I8 K  K9 Z- z9 `
      . B8 p' `# Y\" n6 F& u0 k
    4. d2f=jacobian(df); %求Hessian阵
      1 A8 A' h2 P\" U/ A
    5. [xx,yy]=solve(df)  %求驻点/ M0 M: y) }5 N3 e1 M( J
    6. xx=double(xx);yy=double(yy); %转化成双精度浮点型数据,下面判断特征值的正负,必须是数值型数据  d+ C1 y" F/ b+ b# ^! g7 o( B; p- C+ O
    7. for i=1:length(xx)
      * A6 o* B1 _8 s
    8.     a=subs(d2f,{x,y},{xx(i),yy(i)});
      8 G$ v* ]2 b' P/ `3 P! Y2 M+ }$ A
    9.     b=eig(a);  %求矩阵的特征值0 + ]7 C. K3 h& v\" M  i: [, B
    10.     f=subs(f,{x,y},{xx(i),yy(i)});2 |& S/ j$ U9 E0 r, L
    11.     if all(b>0)2 q8 T' i( o+ T
    12.         fprintf('(%f,%f)是极小值点,对应的极小值为%f\n',xx(i),yy(i),double(f));
      2 ~/ K8 N# g' `
    13.     elseif all(b<0)
      6 [8 z* ^' ?% s4 }& d( t
    14.         fprintf('(%f,%f)是极大值点,对应的极大值为%f\n',xx(i),yy(i),double(f));
      % M, s' O) C  H
    15.     elseif any(b>0) & any(b<0)
      \" w+ {4 i7 C% C\" A
    16.         fprintf('(%f,%f)不是极值点\n',xx(i),yy(i));' K& X3 M- [# [* {+ w* M' y
    17.     else
      ! w% V; @/ Z- [/ V  @; V4 c2 Y
    18.         fprintf('无法判断(%f,%f)是否是极值点\n',xx(i),yy(i));  ( Z$ H/ O, _6 J( O6 ^. M
    19.     end
      ( I3 p$ o! E+ O& j4 i+ f
    20. end
    复制代码

    * @! J! m  |4 P3 [+ j2 F
    新浪微博:http://weibo.com/yzhsky
    Q        Q:
    736466568
    欢迎来数学中国交流,您的问题我们尽量回答,谢谢。
    回复

    使用道具 举报

    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-4-10 19:15 , Processed in 0.304795 second(s), 58 queries .

    回顶部