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

thinkfunny

matlab中有些计算结果比较长，直接查看有些困难，下面介绍pretty和symdisp函数优化输出结果，使结果更为直观。
2 Y5 l* W/ K! K  b5 n
0 B+ z: N2 Y9 |0 H/ v演示示例1
% c# R' }. x% x* l. C$ Y% m有一个计算结果如下：

0 G& @1 \& R% }/ M% L7 w; R

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

   
0 K+ I6 D1 Y! `
1 Q9 A5 h4 c9 b. _( m
1. 使用pretty函数美化输出
  F6 t6 c7 y) K- w' ^4 Z1 f1 l

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

/ d5 |4 [8 z# v9 d8 z# i9 T该函数可使输出更接近数学格式。

9 T- X9 U5 k' z( y7 |0 @; q
2. 使用symdisp函数美化输出

$ \1 x" N$ x7 f$ {+ j! @

• symdisp(f1);
  $ S& z! t, ~4 Q3 N( h- X: Y7 X: f) I8 p

* V  I3 }/ E1 y; ?6 S( ]



* h; O  d7 i  T5 i2 Q& h  n& T演示示例2

有一个计算结果如下：
" I! T! F5 d+ w0 ]
. i! y8 y# X$ A9 j" T

• >> 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
  2 s9 }/ a4 i. B* ~5 [1 c; l: _

  
% i" o- {( d  \9 s5 y! e/ p
2 B; m1 _8 T* K! t0 c9 [
1. 使用pretty函数美化输出

3 \  y& n1 j1 y1 g2 c

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

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


2. 使用symdisp函数美化输出
! a0 O( y. p0 H& B2 ~6 v# ]
7 c$ Z3 B9 h3 @4 ?* m, }

• symdisp(F(3));
  4 w# o; Z! h9 r/ L

6 d- \; s1 U2 P
9 l/ ]- R4 G) ^; r( q/ ~: B8 t

+ ]! O- k7 q) [演示示例3
# R' {1 n5 _: |5 {1 Z. V$ y
4 o( P1 D2 q  v  A5 W) B有一个计算结果如下：

• >> 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)]
  + C6 G) v/ r  @5 t/ ]

  
' w. o0 W$ C- B% i+ Q
$ `4 e. Z2 a9 m0 A" N. z' ?1. 使用pretty函数美化输出
  q# q# s7 c( Y% L; }; l

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

( D* B$ C- q) r% s. W
0 F8 s* y2 E1 v1 s: K该函数可使输出更接近数学格式。

0 e) o& m$ b" N
2. 使用symdisp函数美化输出

* B0 U' G9 U! d

• symdisp(n);
  % o$ d6 C4 D3 x- W, _

7 w3 o( Y2 c( D9 `
1 ^& E9 N* D* v3 `5 q

7 S5 s1 V/ \% v$ [" l
总结
" k1 R) N0 Y7 @& W/ X
经过以上实验，发现symdisp函数可将输出结果转化为更易读的格式，且效果较好
5 M$ M5 d. ?% P/ r
6 a2 T) Z2 r4 K6 y+ u0 p# f
附录：symdisp函数源码

! {6 i7 w& K9 Z% J7 i$ e

• 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');
  

senlenced3

厉害了