matlab 符号运算的化简问题

shouce

首先来个简单
: L/ G; T/ z, ^# ^% N3 N' A4 x8 z
>> syms t r a b q
) w. T5 ]' A7 ]) Z" ~( ]>> x=r*cos(t)+a;
' I& p' q+ V4 I3 W7 q% C. j>> y=r*sin(t)+b;
>> X=-x*cos(q)-y*sin(q);
>> simplify(X)
+ ^& d) K! V) ~5 g
ans =
  _1 G1 M8 l. Z9 a1 E! `$ s
8 s! W4 K8 u! x; E& ~; _# C4 P- a*cos(q) - b*sin(q) - r*cos(q - t)   
/ h" }7 s- i7 o9 O- F
复杂点   就不化简了
" @' c# k5 ]' J) M# N >> syms t r a b k q A
+ ]. m* x# K4 h7 g>>  x=r*cos(t)+a;
+ v" E$ i: {8 Z>> y=r*sin(t)+b;
  K7 Y0 D5 [& X! V6 ?1 Y>>  N=[-cos(k*q) -sin(k*q) cos(q);-sin(k*q) cos(k*q) sin(q)];
>> M=[x;y;A];
- ^: D" Q3 U1 A+ z; O>> NM=N*M

NM =

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)
9 h4 \9 T; y( ?( }/ ?
& X- e0 }+ M( ~% R* a>> x1=A*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t));
, t) U9 [! i4 `+ Z>> simplify(x1)
8 Z. t/ x/ [5 L, ?
5 Z/ y1 I8 u( T8 Q7 A  Eans =
! z8 }( h* ]+ q; C
7 [; c) f$ J2 }+ }' }A*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t))       化简无效果    按正常应该是  
( \- n% A3 t/ w9 y: o$ d
A*cos(q) - a*cos(k*q) - b*sin(k*q)-r*cos(k*q-t)

9 I: V/ c* o* Z! I>> y1=cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q);
6 |* F- W) c% ^$ k. C! x>> simplify(y1)

, y+ o3 t( X/ S0 M8 L6 P( L) qans =

3 X* u) R+ ^9 B2 A- Dcos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q)

明月山河

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

shouce

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

大侠   能举个例子吗   我这个不懂

明月山河

shouce 发表于 2016-2-20 18:33
大侠   能举个例子吗   我这个不懂

& }% Z# \: J5 U9 g" j7 vhttp://blog.sina.com.cn/s/blog_636a8b120100jnyt.html

  p5 r* z7 V8 X# t$ k8 u3.带有附加条件的化简
>simplify(expr,{x*y=1});
>siderel:=x^2+y^2=1;
0 o3 ~1 {; t  w& `>simplify(expr,{siderel},[y,x]);
7 Y8 I( h: [. `4 }% b#maple在表达式中作代换y^2=1-x^2，然后再试图对x^2做代换，由于找不到，就停止了
2 \7 E8 c* A  b# c( T6 R* V4.展开与组合
3 h; N- D( \1 n$ f( F; \" Eexpand &combine
combine是与expand相反的命令，它可以按照数学规则将表达式中的某些项组合在一起

! |5 W; Q) ]3 L' o7 @! l4 p$ G- K# M2 A
" a1 a  i/ ?( }; w( L
( R0 a& c$ [$ c3 B" D# Jmatlab是用的maple的符号运算内核，使用规则一样。
9 J# Z9 h+ c! d" \& X  K

shouce

明月山河 发表于 2016-2-20 18:38
http://blog.sina.com.cn/s/blog_636a8b120100jnyt.html
" d; R! c4 k& C, r7 T
- }3 M) O2 L4 v- x8 b( B3 b$ D3.带有附加条件的化简

0 h8 k0 U  g7 e多谢了   等下去试试
3 z, Z3 L( J, F0 P* L5 ^$ i

shouce

明月山河 发表于 2016-2-20 18:38
( |+ z# d. v% h+ Q: \  W/ Shttp://blog.sina.com.cn/s/blog_636a8b120100jnyt.html
* N, x: E+ Y- _9 h5 G8 l4 k4 k
5 K+ Y4 y# r5 J' o& [' n$ V3.带有附加条件的化简

试了下还是不行   可能matlab和maple 使用不一样
>> syms x
>> combine(sin(x)^2+cos(x)^2)   

ans =

( a7 Y% n  v2 v# `' s8 Gcos(x)^2 + sin(x)^2          这个不化简单

  a5 l: |" \( z7 w! E' v7 g. \' f0 A>> simplify(sin(x)^2+cos(x)^2)
! P2 K! y7 X9 W/ j' |6 B, K
6 h- ]  ~5 ]- Y, q3 R) r+ x; eans =
! a4 ~4 Y+ ]; G, Q/ z0 r
1
' J2 j, H7 a  F2 Q. E; K1 G5 n1 v
; |! s8 V  b4 T1 z0 M! k
>> syms x
>> A=sym('ln(3*x)+sin(x)^2+cos(x)^2')
5 C" d7 W! T0 J1 h
A =
4 W# T. P  h7 c* h5 ~$ `! ~2 ^
log(3*x) + cos(x)^2 + sin(x)^2

>> simplify(A)
- v$ g! g9 o: B5 h
ans =

log(3*x) + 1

+ x% X5 J1 b8 G1 O2 k% w8 {; V/ d

>> simplify(A,x)
# y6 Q( a5 l$ S/ s- Q错误使用 sym/simplify (line 99)
  v: Z6 G& b% U# E* G3 |Too many input arguments.
; q* e  G& v% H. r

393810298

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

shouce

>>  expand(x1)
# Y/ a$ V: T# E* G0 K: ^ans =
A*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)')
' P. |/ }% c7 V9 Y+ U  ^KMM =
A*cos(q) - b*sin(k*q) - a*cos(k*q) - r*cos(k*q)*cos(t) - r*sin(k*q)*sin(t)
8 a2 m  e# I- M) Q>> simplify(KMM)
  Q' N) S. }/ {1 L. q% Z6 Qans =
A*cos(q) - a*cos(k*q) - b*sin(k*q) - r*cos(t - k*q)

9 H; s1 j5 U7 p) l+ f要先展开然后合并