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

thinkfunny

5 ?, N  \1 u  G, o: q% Y

  Z9 h; }* n5 z' r2 u" Lmatlab中有些计算结果比较长，直接查看有些困难，下面介绍pretty和symdisp函数优化输出结果，使结果更为直观。
& V, j( _* S- m/ e! N* d
演示示例1
6 D0 q9 Z8 K0 R/ I有一个计算结果如下：

( \% L7 Y3 n8 `' K. G

• >> 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 ~5 c' s; w! W# e: Q  [

   
9 O! X- }* o- L" w) F4 l" D. R

# w3 f7 ?5 O) }  j# t" t1. 使用pretty函数美化输出
% i  }/ _( d9 N8 K# v7 J/ t

• >> 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
• 
  & O% M$ G7 K9 s

+ Y$ S% \' R4 Z' _6 I
该函数可使输出更接近数学格式。
2 h8 ^$ k' z4 k1 l
7 W% h# |1 j  ?7 `8 a  b7 o
2. 使用symdisp函数美化输出
8 I* T7 I/ ~0 K( ]0 |
) e$ w2 B3 b7 H0 a' |

• symdisp(f1);
  

0 E; D( t8 x) N) R  q& q# K
% l5 a' W% ?7 n/ W5 x, J* V; \
+ ~9 R: q1 I6 i+ [
$ U( m9 }% L2 G/ I- V3 S
演示示例2

: [% F" W$ Y8 w+ _' g7 V7 w有一个计算结果如下：
/ l! p( n( b1 h  M, ]4 ^

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

  


% W7 d7 O( Y1 _& n$ @7 ]) v  ?# d1. 使用pretty函数美化输出
6 v: p. Y8 z# ~+ ]% V$ h6 l. h& y
! D+ a3 a" A! t+ X, e) d) O

• >> 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函数美化输出

• symdisp(F(3));
  9 o" j) u& B( Q9 D4 G

" O# b- n5 l% O
- c9 P: k8 @! y( ~% T
- x# r/ R: g% D$ q4 A- }. f

  \* ~+ j% @% z( T. w* X演示示例3

有一个计算结果如下：

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

  

1. 使用pretty函数美化输出
, F2 v7 B: b8 R" `
6 F& K9 A# g0 s: n% D. P

• >> 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
  1 C; z8 ]0 g  @: x

$ E4 e3 F- k7 b9 b
) o6 ?, ~' b: C5 b3 F该函数可使输出更接近数学格式。

* F7 @* v) |+ o- ?( R- B8 G  `) L
3 h9 e# O- y" |  v2 I% e# @2. 使用symdisp函数美化输出

• symdisp(n);
  

- e( X5 s1 Q" n) U* I; O. J/ ?' p
$ b* D! o3 J% E
) z' z4 P) f% _" C2 Q0 u  e总结

经过以上实验，发现symdisp函数可将输出结果转化为更易读的格式，且效果较好
- _: L  a4 X$ [' u' z( R1 Z, y; |
- Q2 f0 {! z$ x& h3 ~5 o
% W9 `" d: W7 z  C" v2 A% r附录：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');
  

6 l2 v9 b! M; P2 S& c; u' C

senlenced3

厉害了