matlab 符号运算的化简问题

shouce

首先来个简单
9 E* s& t1 V; {
( W/ z2 c4 m$ k" L9 m- `>> syms t r a b q
: [8 z& S$ }0 C/ W1 g8 _>> x=r*cos(t)+a;
>> y=r*sin(t)+b;
7 _0 p! I! k  k8 b" a+ a6 c8 n>> X=-x*cos(q)-y*sin(q);
+ q+ i: H& u+ N5 C>> simplify(X)

2 h, F* }+ r0 I9 S0 _ans =

" N9 {$ M; |* ?. c- a*cos(q) - b*sin(q) - r*cos(q - t)   

复杂点   就不化简了
>> syms t r a b k q A
3 h3 H1 B0 x& e* q) }>>  x=r*cos(t)+a;
% R/ `. R$ _1 |3 P! N4 D( j>> y=r*sin(t)+b;
' a8 ^; W1 t; d( }2 R  K>>  N=[-cos(k*q) -sin(k*q) cos(q);-sin(k*q) cos(k*q) sin(q)];
>> M=[x;y;A];
>> NM=N*M
4 S7 P" {6 ]! U7 r9 V3 F0 H
NM =

2 G$ w- P5 n# h3 y# ~/ b! R A*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t))
: A( ?$ C# u; {7 ] cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q)
9 {# ?( d( z6 w* c+ x5 C
>> x1=A*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t));
>> simplify(x1)
" ]8 t. b8 L( P/ C! o1 U
ans =
6 A, E# ]' i; P  f% r+ q: h) V: I
# g& p1 F4 v' L) y" ~" @A*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t))       化简无效果    按正常应该是  

5 w. V" w: G9 h6 h% u5 d6 s- ~A*cos(q) - a*cos(k*q) - b*sin(k*q)-r*cos(k*q-t)

>> y1=cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q);
>> simplify(y1)
0 g6 w! s8 C& @
ans =
# N# _4 _( C' U
2 Q' S9 J; a: j) P$ C- D* I. Kcos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q)
4 z2 B7 c8 G9 @) A

明月山河

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

shouce

明月山河 发表于 2016-2-20 18:27
simplifty函数要有控制参数才会按照你想像的运算。

% S. z1 `# a4 |  c* B! a4 d3 ?' `大侠   能举个例子吗   我这个不懂

明月山河

shouce 发表于 2016-2-20 18:33
0 U5 `# O  H: w( F5 K0 f大侠   能举个例子吗   我这个不懂

http://blog.sina.com.cn/s/blog_636a8b120100jnyt.html
5 i9 K" a! |! e0 |5 }6 r0 N7 W
3.带有附加条件的化简
( P' k8 f0 b- J& L, ]>simplify(expr,{x*y=1});
>siderel:=x^2+y^2=1;
& l0 o  o0 A" N7 C0 t/ o" e  k>simplify(expr,{siderel},[y,x]);
#maple在表达式中作代换y^2=1-x^2，然后再试图对x^2做代换，由于找不到，就停止了
+ F* U' \: Y2 D' o% `* C+ }+ n9 ^4 @4.展开与组合
! J! a& t$ U4 o" ^5 d! k" yexpand &combine
6 X: D6 g1 ^3 I* acombine是与expand相反的命令，它可以按照数学规则将表达式中的某些项组合在一起


2 z! D, n+ g; X( s# s
matlab是用的maple的符号运算内核，使用规则一样。
5 r6 g7 k% Y& m* J  h! b9 g
+ H+ N  p1 q( k  g

shouce

明月山河 发表于 2016-2-20 18:38
' T- Q( [; }* Jhttp://blog.sina.com.cn/s/blog_636a8b120100jnyt.html

; ~$ R9 B& @0 z# f8 M$ o3.带有附加条件的化简

多谢了   等下去试试
0 H. m" J2 t! K3 m

shouce

明月山河 发表于 2016-2-20 18:38
http://blog.sina.com.cn/s/blog_636a8b120100jnyt.html

3.带有附加条件的化简

# [* @! f; F6 L2 T( {试了下还是不行   可能matlab和maple 使用不一样
! A# d" E$ `. U7 ]- W! _8 {>> syms x
>> combine(sin(x)^2+cos(x)^2)   
3 y, g" g* X! D. t2 {
- S% x1 d0 R' G; @. n7 }ans =
+ x! Q1 |& ]* M6 ?9 `5 @# d1 {# `
cos(x)^2 + sin(x)^2          这个不化简单
1 G1 s) ]% {- r# g) [8 t1 R
+ _! [0 G) l, G% M/ O0 K' G>> simplify(sin(x)^2+cos(x)^2)
6 B. Q2 D  R, x! ?( \& X3 V! Q
ans =
' |5 H( j! O5 W  R2 P: [7 N# [+ G
1

, f& ?& j- A" B% |% C' g
4 V4 e; R( L$ K$ G4 F% t7 S>> syms x
>> A=sym('ln(3*x)+sin(x)^2+cos(x)^2')
& ?' z4 j& J0 l5 n2 u% w2 g
# g$ F6 y* F& [4 d. ~) j& x- s& L2 UA =
4 j* \0 ]# l" }
2 G; A8 ?0 N8 g; g; J, Rlog(3*x) + cos(x)^2 + sin(x)^2
0 q! d. e5 B0 L6 c( P7 V+ {
$ q$ [/ U% U" T& w& p, I" S; \# U% A/ A>> simplify(A)

4 M$ b- c0 z6 L" K- o5 vans =
6 V: Q: _$ r3 }. s* r
log(3*x) + 1
7 p+ _8 A; }1 K5 a
6 e. p" f% g& R. p, Y
1 P3 _: B5 `) c1 C% r7 W
>> simplify(A,x)
/ I" Y5 W. |9 y/ r9 U错误使用 sym/simplify (line 99)
Too many input arguments.

393810298

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

shouce

>>  expand(x1)
ans =
% N$ y* }. {  U( B6 UA*cos(q) - b*sin(k*q) - a*cos(k*q) - r*cos(k*q)*cos(t) - r*sin(k*q)*sin(t)
>> 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)')
  }- E- D; M  e6 z% AKMM =
' q8 P% |8 w9 \8 s' [A*cos(q) - b*sin(k*q) - a*cos(k*q) - r*cos(k*q)*cos(t) - r*sin(k*q)*sin(t)
' m8 y) {2 O+ v! K* x>> simplify(KMM)
ans =
A*cos(q) - a*cos(k*q) - b*sin(k*q) - r*cos(t - k*q)

要先展开然后合并