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

thinkfunny

& s* Q3 X4 K# `' E: B( n2 g
+ ?; ~4 l% A% d. K/ i5 i3 B4 Hmatlab中有些计算结果比较长，直接查看有些困难，下面介绍pretty和symdisp函数优化输出结果，使结果更为直观。

演示示例1
. ?7 m4 b4 p- P4 ?有一个计算结果如下：

* |# l8 \7 i7 D7 w1 ]

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

   


1. 使用pretty函数美化输出
3 h- h1 n6 Y/ w, c9 c

• >> 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
• 
  & L' m3 ]: O2 q

3 y) E8 i0 L& W3 a% b2 C$ X该函数可使输出更接近数学格式。
* B9 a3 f+ U3 Z' }  C  @: g
" r: i( M* k7 M1 C
9 V+ z4 k, e6 }8 V! H, D2. 使用symdisp函数美化输出

0 i% l# \5 h% q  x! E

• symdisp(f1);
  , z9 M7 T/ V+ Z) x

& L; y' l0 m) k$ F( X" ?  e* ]


1 L+ q* @, y3 c; |6 ~
8 {/ ~6 I, E6 n' A& i# H演示示例2
9 ?3 k' d1 K* x. [/ N
有一个计算结果如下：
& i, X: m' j: m, e6 V
. b5 i9 \  P% p+ A. c0 |

• >> 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
  & n5 ?0 K# V5 h: G' U- a

  

3 h% ^2 y" d3 B/ A/ W
1. 使用pretty函数美化输出

• >> 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
  4 Z% I3 g2 w5 D( o

2 n7 Z  o1 r" ^2 y& E
0 h2 Z7 h0 G! ~( O* T! q- Y  R& q该函数可使输出更接近数学格式。


2. 使用symdisp函数美化输出

( s, [! j7 \/ z4 P+ i) X% g! y5 F

• symdisp(F(3));
  $ M* b6 H% [  {& f( Q

5 m6 y- s; M/ y( g3 ~% ?" o

8 t0 W# I0 r9 Q' D$ R: M
演示示例3

有一个计算结果如下：

- }+ m; p1 e2 Z# F' {  M6 _0 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)]
  * G4 k/ Z$ |# v

  
/ u! F* i  J0 l, V' B6 z
. i( A" x; E1 U, x) X7 t$ c1. 使用pretty函数美化输出
2 L/ y3 p, s6 X

• >> 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
  2 E7 a% g5 A" w3 O5 @8 B0 s0 T' t

该函数可使输出更接近数学格式。
+ d' s+ U# L9 r

* g- Y7 y, X3 n+ U4 ~4 `; ?2. 使用symdisp函数美化输出

• symdisp(n);
  

: B: x! Z  D0 T6 t& F, _
3 T1 x. Y7 W- i' B$ Q
: C! w) ?$ `- _" @5 _
& _. t4 x2 O% K" y4 O; u5 g
总结

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

6 ?6 o. p* B+ W. D' G/ j1 H" }9 M8 M
附录：symdisp函数源码
- K6 @. s8 `; h% t7 t/ s( y3 r

• 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');
  0 J2 x6 M8 G/ N# D+ [# n; h7 r

2 d9 L$ Q, [- I8 G

senlenced3

厉害了