matlab 符号运算的化简问题

shouce

首先来个简单
0 N- X4 Y$ a7 K
) }  Q, w$ r2 A! B* S/ q>> syms t r a b q
) L. Y7 G5 F3 [* h/ b9 m; X>> x=r*cos(t)+a;
; q' {; C3 C5 p>> y=r*sin(t)+b;
  N: N0 d; j+ g" d9 j+ I. c>> X=-x*cos(q)-y*sin(q);
>> simplify(X)

ans =

! a" O$ C5 d( m, M, R) T- a*cos(q) - b*sin(q) - r*cos(q - t)   
. t; l& Y, g+ C0 a4 T" Y* B% T
* {" M' b. |( \/ p* C) ?$ r复杂点   就不化简了
>> syms t r a b k q A
( T- V2 I: M$ B5 D; _>>  x=r*cos(t)+a;
>> y=r*sin(t)+b;
, y9 _# r7 o& @9 R& ~' m>>  N=[-cos(k*q) -sin(k*q) cos(q);-sin(k*q) cos(k*q) sin(q)];
0 B' R; U' G2 E1 e9 C; ]>> M=[x;y;A];
  n" L- T: L% P1 ^) |4 g" ?. M>> NM=N*M
/ J: G7 l. a; h& _3 g: ], |
, g. `) }. B* L$ b6 t) F& t/ JNM =

0 c$ s- D8 P8 q A*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t))
cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q)

>> x1=A*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t));
% f; `# z) `5 J& o; t>> simplify(x1)

ans =
2 h$ }5 U$ r  Z- V& V! ?
7 I$ K* ]8 [" b1 EA*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t))       化简无效果    按正常应该是  

7 a$ r& K6 T) U; WA*cos(q) - a*cos(k*q) - b*sin(k*q)-r*cos(k*q-t)

' p) e% N* _9 Z' ]) s>> y1=cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q);
" ]' \) U& c& z' g$ n>> simplify(y1)
$ ?; k. C! O( d& s9 G1 m6 t/ o/ e, L8 m
ans =

6 k; V+ F6 l. o# A  Vcos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q)
  |) N0 n  f' i6 }8 ?) [9 S
* P1 t# [3 u( |& Z

明月山河

simplifty函数要有控制参数才会按照你想像的运算。

shouce

明月山河 发表于 2016-2-20 18:27
+ p/ r% x2 Y3 Fsimplifty函数要有控制参数才会按照你想像的运算。

7 v6 U: W: \1 x7 v大侠   能举个例子吗   我这个不懂

明月山河

shouce 发表于 2016-2-20 18:33
3 [& c4 G& Y7 t8 Y: V+ o7 d大侠   能举个例子吗   我这个不懂

http://blog.sina.com.cn/s/blog_636a8b120100jnyt.html
" c" c/ \- F2 S" F1 Q* z
3.带有附加条件的化简
, V3 i# h4 S* @# D, h>simplify(expr,{x*y=1});
>siderel:=x^2+y^2=1;
# s" B( L7 K1 B9 G>simplify(expr,{siderel},[y,x]);
#maple在表达式中作代换y^2=1-x^2，然后再试图对x^2做代换，由于找不到，就停止了
4.展开与组合
% ~: \: c" O2 @1 Z# Cexpand &combine
combine是与expand相反的命令，它可以按照数学规则将表达式中的某些项组合在一起



7 H6 L0 i* [/ ^) J7 Smatlab是用的maple的符号运算内核，使用规则一样。
" e  V& M; R" o, J9 M/ w& c5 s

shouce

明月山河 发表于 2016-2-20 18:38
  q4 t& [. p: r# _& shttp://blog.sina.com.cn/s/blog_636a8b120100jnyt.html

7 ?5 ~: L% p( P3.带有附加条件的化简

多谢了   等下去试试

shouce

明月山河 发表于 2016-2-20 18:38
http://blog.sina.com.cn/s/blog_636a8b120100jnyt.html
1 Z$ A$ S- `7 c
3.带有附加条件的化简

  J3 ]' G. q2 b# A2 d4 \0 @/ ]试了下还是不行   可能matlab和maple 使用不一样
+ d3 X! A9 F7 _# P8 [8 z2 o>> syms x
>> combine(sin(x)^2+cos(x)^2)   

4 m# c6 n- G7 x* b6 kans =

8 G* p( K& `6 A( P5 zcos(x)^2 + sin(x)^2          这个不化简单
) Q2 A$ x9 Q* k- N
>> simplify(sin(x)^2+cos(x)^2)
: o, D& G; H& G/ V& @, h) ?
ans =
4 G6 D( W# X4 S
) H! A! l/ M! k8 D0 s; W1
& \& q+ ]! \( L4 G1 b
- o$ a: u- @$ ?7 o  W6 b2 ~, _
>> syms x
, H. S/ h+ v7 L& b>> A=sym('ln(3*x)+sin(x)^2+cos(x)^2')

A =

log(3*x) + cos(x)^2 + sin(x)^2
' W# s1 ~$ U3 i( R5 j, C4 ^! C
>> simplify(A)

- N! F# Y- }: A2 K+ H* Wans =
& M2 u) U, w# P; B
, ?, C3 D4 |6 i6 t8 {: alog(3*x) + 1
# ^  |; J1 R8 s/ f% j5 ^8 a) [
, D9 `8 d5 ~8 d9 U0 x# K' L
' j( x  H6 l" r
>> simplify(A,x)
错误使用 sym/simplify (line 99)
2 m. J* n! B0 Q6 u8 c( }Too many input arguments.

393810298

matlab大学里学完了，我就在也没有用过了，已经看不懂了

shouce

>>  expand(x1)
8 ~! }1 F0 j, u2 Q1 ~' D* uans =
1 v* L: x5 |( R4 y5 E0 V' dA*cos(q) - b*sin(k*q) - a*cos(k*q) - r*cos(k*q)*cos(t) - r*sin(k*q)*sin(t)
# M# Y5 n4 w5 N>> KMM=sym('A*cos(q) - b*sin(k*q) - a*cos(k*q) - r*cos(k*q)*cos(t) - r*sin(k*q)*sin(t)')
0 W! X1 J1 k  P5 aKMM =
- m6 c# J1 L' q$ L- {% C7 I: _  l# @A*cos(q) - b*sin(k*q) - a*cos(k*q) - r*cos(k*q)*cos(t) - r*sin(k*q)*sin(t)
>> simplify(KMM)
ans =
A*cos(q) - a*cos(k*q) - b*sin(k*q) - r*cos(t - k*q)

% s2 P; p: w7 L2 L要先展开然后合并