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

thinkfunny

- S/ B9 h3 S3 y

matlab中有些计算结果比较长，直接查看有些困难，下面介绍pretty和symdisp函数优化输出结果，使结果更为直观。
5 \4 }+ o2 j% U0 I3 m6 `
, E" s+ ?, ~8 K4 U, {演示示例1
) C  H9 u( }$ P0 J有一个计算结果如下：

4 Z& u# h3 i, ^: j' c$ L

• >> 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
• 
  

   

8 O' w5 S- F0 e- n5 I  S4 o
% R6 t! V3 L( C1. 使用pretty函数美化输出
. s, A+ X  [  J# ]% N/ V/ x! a% N

• >> 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
• 
    c4 u9 {5 O2 r9 m

3 N% `& o. V# c9 Y( a' p/ y该函数可使输出更接近数学格式。


2. 使用symdisp函数美化输出
. X' z: \9 p5 i" L9 g/ L
$ C& j7 M5 K  t4 Y* \$ \9 U$ c

• symdisp(f1);
  

1 ^; C9 _1 n. H



演示示例2
, R5 }8 X& S$ r7 l
有一个计算结果如下：
5 A' n+ y; k0 w0 x9 q6 D% e

• >> 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
  

  


1. 使用pretty函数美化输出
. D' }9 l2 O$ r8 m0 E% ^; J5 m; l

• >> 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
  

# a/ I# Y. b; H( k, ^1 o3 ]
3 R1 ^: O7 A1 |3 g/ v) d* l+ m* ]该函数可使输出更接近数学格式。

) Y- J* a/ k  N
7 H. M( x. I9 j/ g0 T0 ^2. 使用symdisp函数美化输出
- \4 a0 r: h) {0 K6 |0 _& m6 ~6 T
+ i1 O( \$ j4 d" `; P% g0 o

• symdisp(F(3));
  . E: ^4 _9 W2 ~( D- j: g! z

& h& u) V- L7 o4 x( p- _
7 |3 A. q+ ?% q: L/ w/ [
8 A8 Y' `) ]/ d  @5 n- P3 y

; Y) _! }  N, c6 h演示示例3
. k7 o* X9 a( N0 X
6 \! B4 i- |8 ]5 p) T8 B* v6 y: L有一个计算结果如下：
7 ?- B3 X$ t% a' X8 J3 [, ]( l
/ W. h; b, Y- `- V( N) R! F4 @

• >> 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)]
  

  
! s6 M+ y4 v5 R2 R
1. 使用pretty函数美化输出
' G) w! e: z6 ~) ^. M" J3 n

• >> 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
  ( W8 l& f. d; l$ n3 H' S

3 h* a$ L( M3 H$ `+ s5 P- p
该函数可使输出更接近数学格式。
8 m! v& q! i9 q  |

- U" ~: H7 M8 M8 C0 r3 S; O' c2. 使用symdisp函数美化输出

4 n0 H( w' C& S, w* w; @

• symdisp(n);
  $ q! e) l5 ?% ~' u6 B& R/ J

6 t6 n) R4 v' _
9 W/ t5 m# l) n2 y
- j  q! `0 {$ f6 t; y
总结
- T2 O$ v$ w* B' V$ D
经过以上实验，发现symdisp函数可将输出结果转化为更易读的格式，且效果较好
! N- E( _2 m  _/ i. b; N

附录：symdisp函数源码

• 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');
  4 q4 {& _1 s$ o7 `0 w- X

% I$ @5 _5 B% ]  B/ |; \. e6 U

senlenced3

厉害了