matlab符号编程real类型数据 real imag expand使用实例

alljoyland

1. syms m1 m2 c k1 k2 Y w  den num  real
   
2. %这个real很重要的  如果在maple里面应该是 assume指令相当的
   
3. % sym(m1,'real')
   
4. % sym(m2,'real')
   
5. % sym(c,'real')
   
6. % sym(k1,'real')
   
7. % sym(k2,'real')
   
8. % sym(Y,'real')
   
9. % sym(w,'real')
   
10. % sym(den,'real')
    
11. % sym(num,'real')
    
12. % sym(,'real')
    
13. % sym(,'real')
    
14. M =[ m1 0; 0 m2]
    
15. C =[ c -c ; -c c]
    
16. K = [ k1,-k1;-k1,k1+k2]
    
17. U0 = -w^2*M + i*w*C +K
    
18. F=[ 0;k2*Y]
    
19. U = inv(U0)*F
    
20. U = inv(U0)*F
    
21. U1=U(1,:)
    
22. [num,den]=numden(det(U1))
    
23. % collect(den,i)
    
24. % U2=U(2,:)
    
25. Imden=imag(-den)
    
26. Imnum=imag(-num)
    
27. imden = expand(Imden)
    
28. imnum = expand(Imnum)
    
29. realden = expand(real(-den))

复制代码

本题是 两自由度系统的一个例子的求解

sigma665

lz发这段程序是什么意思

ch_j1985

原帖由 sigma665 于 2008-8-7 22:14 发表

                               
登录/注册后可看大图

lz发这段程序是什么意思

LZ发这段程序应该是要告诉大家real在求解两自由度系统时的重要性

sigma665

syms m1 m2 c k1 k2 Y w  den num  real
这个，把real当成变量名了，是怎么回事

ch_j1985

原帖由 sigma665 于 2008-8-8 09:46 发表

                               
登录/注册后可看大图

syms m1 m2 c k1 k2 Y w  den num  real
这个，把real当成变量名了，是怎么回事

real在这里不是变量名，而是用来声明m1 m2 c k1 k2 Y w  den num这些变量为实变量（实数），而不是虚变量（虚数）