matlab如何使输出结果更美观(symdisp函数——pretty函数升级版)

thinkfunny

1 ?" ^' I9 Q6 h. f: j
5 U# G/ c$ R, C( r$ ~: S& L
3 ^8 _1 s: q, ~5 O' Umatlab中有些计算结果比较长，直接查看有些困难，下面介绍pretty和symdisp函数优化输出结果，使结果更为直观。

3 P1 _& z8 L3 J. C6 g演示示例1
- b' ^: S/ W$ v" t* L5 k6 W8 K4 I) f有一个计算结果如下：
8 s( t* B! u" {: X: f8 F

• >> f1
• f1 =
• y^5 + (- w - y0)*y^4 + 1800*y^3 + (1498200*w - 1800*y0)*y^2 + (3600*w*y0 + 810000)*y - 1350810000*w - 810000*y0
• 
  9 ]# O# C/ y, i, R& l$ L( ^

   

: u( x- K* l6 I( Z  {) h! @- F
1. 使用pretty函数美化输出

• >> pretty(f1)
• 5               4         3                          2
• y  + (- w - y0) y  + 1800 y  + (1498200 w - 1800 y0) y  + (3600 w y0 + 810000) y - 1350810000 w - 810000 y0
• 
  6 O3 r/ X$ P1 l" f5 u6 _

该函数可使输出更接近数学格式。
6 ^1 A9 Y- ^1 E: Q

+ e, @8 i1 L" N. w2. 使用symdisp函数美化输出

+ w5 A; [' i' G- ?# w; L* {* m

• symdisp(f1);
  ! o8 }, Q( k# H8 S/ `

4 v8 g- n; p! |, Z3 W

. q% W! J1 x2 j0 y% `9 m7 w" K

演示示例2
8 y% W. c& J2 V) l6 {7 @. T3 j
* S1 D9 r6 w" t/ f: Z有一个计算结果如下：

• >> F(3)
• ans =
• (y^2*((w - (2*w*y^2)/(y^2 + 900))/(y - y0) + 1))/60 - (25015*(w - (2*w*y^2)/(y^2 + 900)))/(y - y0) - (60*w*y)/(y^2 + 900) + 15
  5 q+ {7 Y9 n: k9 L

  

7 z6 l$ Q7 S8 m( P# G
1. 使用pretty函数美化输出
4 G! |& R( p+ F# r' y3 w# g9 h* E

• >> pretty(F(3))
•    /           2      \
•    |      2 w y       |
•    | w - --------     |   /           2  \
•    |      2           |   |      2 w y   |
• 2 |     y  + 900     |   | w - -------- | 25015
• y  | ------------ + 1 |   |      2       |
•    \    y - y0        /   \     y  + 900 /          60 w y
• ----------------------- - ---------------------- - -------- + 15
•            60                     y - y0            2
•                                                    y  + 900
  / ^# m# D6 h. @+ s0 O# e

该函数可使输出更接近数学格式。


  t5 \" H! ^- }2. 使用symdisp函数美化输出

• symdisp(F(3));
  

7 h9 E) n, @1 Z! N3 V) i


演示示例3
! f" u( E' l  `8 P$ Y
5 \) t* @, Y6 \8 P' [+ o  q( P  z有一个计算结果如下：

5 I8 A# s+ G  e+ ]' a0 f

• >> n
• n =
• [ -(2*x*(2*r^2 - 4*x^2 + 4*y*(w - y) - 1800))/(4*x^2 + 4*y^2 + 3600), -(2*y*(2*r^2 - 4*x^2 + 4*y*(w - y) - 1800))/(4*x^2 + 4*y^2 + 3600), (60*(2*r^2 - 4*x^2 + 4*y*(w - y) - 1800))/(4*x^2 + 4*y^2 + 3600)]
  

  
- z# @; S6 u% C' Y6 w
. k: O8 j6 z' ?! G% k1. 使用pretty函数美化输出

8 t9 I: H* G% }  y

• >> pretty(n)
• /         x #1 2                y #1 2               60 #1       \
• | - ------------------, - ------------------, ------------------ |
• |      2      2              2      2            2      2        |
• \   4 x  + 4 y  + 3600    4 x  + 4 y  + 3600  4 x  + 4 y  + 3600 /
• where
•             2      2
•    #1 == 2 r  - 4 x  + 4 y (w - y) - 1800
  

( O. a; V( H! r
! g& s& }/ s8 n" E# P# w. {) \# O该函数可使输出更接近数学格式。

/ L0 B0 E3 f% F$ g
2. 使用symdisp函数美化输出
0 \) q  `$ P, r: c1 w
: _" f! u- a6 J. H

• symdisp(n);
  8 h! |6 T. j9 u% ^- v+ p$ D

) O6 Q, @9 S: o1 F5 U+ w+ p

6 [% T, h6 j6 ^( B
' ]& G; h( x9 u3 e9 o5 R
4 z3 I+ q' I% Y' A5 @1 w9 A7 a总结

经过以上实验，发现symdisp函数可将输出结果转化为更易读的格式，且效果较好

0 l1 N, U# p# L
( V& V/ v  L* ^附录：symdisp函数源码

3 d& o  _- a8 R7 ]* ^* \2 S  i

• function h=symdisp(s)
• %//SYMDISP Display a symbolic expression in human readable form.
• %// symdisplay(S) displays the symbolic expression S in a small figure window,
• %// using standard mathematical notation.
• %//
• %// Examples:
• %//   syms x t positive
• %//   f=taylor(cos(x));
• %//   symdisp(f)
• %//   f=int(exp(-t)*t^(x-1),t,0,inf);
• %//   symdisp(f)
• %//
• %// Required toolbox: Symbolic Math
• %//
• %// See also SYMBOLIC PRETTY.
• if ~isa(s,'sym')
•     s=sym(s);
•     %error('输入参数必须是sym类型,请使用 sym() 将你的结果转化为sym类型.')
• end
• S=['$',latex(s),'$'];
• S=strrep(S,'&','& \quad');
• S=strrep(S,'{\it','\mathrm{');
• h=msgbox(S,'字符的数学展示形式');
• h1=get(h,'children');
• h2=h1(1);
• h3=get(h2,'children');
• if isempty(h3)
•     h2=h1(2); h3=get(h2,'children');
• end
• set(h3,'visible','off')
• set(h3,'interpreter','latex')
• set(h3,'string',S)
• set(h3,'fontsize',20)
• w=get(h3,'extent');
• W=get(h,'position');
• W(3)=max(w(3)+10,125);
• W(4)=w(4)+40;
• set(h,'position',W)
• h4=h1(2);
• if ~strcmp(get(h4,'tag'),'OKButton'), h4=h1(1); end
• o=get(h4,'position');
• o(1)=(W(3)-o(3))/2;
• set(h4,'position',o)
• set(h3,'visible','on')
• set(h,'color','w');
  " C+ }+ K. r5 i9 y6 Z# y9 u1 M5 w

senlenced3

厉害了