Matlab 版 (精华区)

发信人: zjliu (秋天的萝卜), 信区: Matlab
标  题: [合集] 谁有对符号矩阵求逆的实例以及如何建立符号矩阵
发信站: 哈工大紫丁香 (Mon Nov 22 11:39:32 2004), 站内

────────────────────────────────────────
  damutou (大木头)                 于  (Wed Sep 22 10:01:01 2004)  说道:

这个矩阵:                                                                     
  
[k1 k2 k3                                                                     
  
 l  m  n                                                                      
  
    p  q]                                                                     
  
求它的逆矩阵，结果如何？

如题,请各位大侠不吝赐教,本菜鸟在此谢谢大家



────────────────────────────────────────
  zjliu (秋天的萝卜)               于  (Wed Sep 22 10:05:43 2004)  说道:

例子
syms x y z
A=[x,y,z;z,z,z;y,z,x];
inv(A)
【 在 damutou (大木头) 的大作中提到: 】
: 如题,请各位大侠不吝赐教,本菜鸟在此谢谢大家




────────────────────────────────────────
  damutou (大木头)                 于  (Wed Sep 22 10:07:19 2004)  说道:

这个矩阵:
[k1 k2 k3
 l  m  n
    p  q]
求它的逆矩阵，结果如何？
【 在 damutou (大木头) 的大作中提到: 】
: 如题,请各位大侠不吝赐教,本菜鸟在此谢谢大家




────────────────────────────────────────
  zjliu (秋天的萝卜)               于  (Wed Sep 22 10:11:25 2004)  说道:

赫赫，你是不是漏了一个数？
【 在 damutou (大木头) 的大作中提到: 】
: 这个矩阵:
: [k1 k2 k3
:  l  m  n
: ...................



────────────────────────────────────────
  damutou (大木头)                 于  (Wed Sep 22 10:12:22 2004)  说道:

没有漏掉，a31的位置是零，我怕写成o，所以留了空格
【 在 zjliu (秋天的萝卜) 的大作中提到: 】
: 赫赫，你是不是漏了一个数？




────────────────────────────────────────
  damutou (大木头)                 于  (Wed Sep 22 10:13:19 2004)  说道:

我用手算了半天刚算好一行，太慢了
【 在 zjliu (秋天的萝卜) 的大作中提到: 】
: 赫赫，你是不是漏了一个数？




────────────────────────────────────────
  billlang (天也...也天)           于  (Wed Sep 22 12:04:36 2004)  说道:

matlab下运算如下：

>> syms k1 k2 k3 l m n p q
>> A=[k1 k2 k3;l m n;0 p q];
>> inv(A)
 
ans =
 
[   -(-m*q+n*p)/(k1*m*q-k1*n*p-l*k2*q+l*k3*p),  (-k2*q+k3*p)/(k1*m*q-k1*n*p-l*
k2*q+l*k3*p), -(-k2*n+k3*m)/(k1*m*q-k1*n*p-l*k2*q+l*k3*p)]
[          -l*q/(k1*m*q-k1*n*p-l*k2*q+l*k3*p),          k1*q/(k1*m*q-k1*n*p-l*
k2*q+l*k3*p),  (-k1*n+k3*l)/(k1*m*q-k1*n*p-l*k2*q+l*k3*p)]
[           l*p/(k1*m*q-k1*n*p-l*k2*q+l*k3*p),         -k1*p/(k1*m*q-k1*n*p-l*
k2*q+l*k3*p),   (k1*m-k2*l)/(k1*m*q-k1*n*p-l*k2*q+l*k3*p)]
 
>> 


【 在 damutou (大木头) 的大作中提到: 】
: 我用手算了半天刚算好一行，太慢了


结果好复杂，手算很容易出错的 :)


────────────────────────────────────────
  damutou (大木头)                 于  (Wed Sep 22 23:08:18 2004)  说道:

看来结果还很有规律,分母都相同,谢谢你了,大侠
【 在 billlang (天也...也天) 的大作中提到: 】
: matlab下运算如下：
: ans =
: [   -(-m*q+n*p)/(k1*m*q-k1*n*p-l*k2*q+l*k3*p),  (-k2*q+k3*p)/(k1*m*q-k1*n*p-l*
: ...................



────────────────────────────────────────
  rtm (南神之火)                   于  (Wed Sep 22 23:31:49 2004)  说道:

用^(-1)也行
【 在 damutou (大木头) 的大作中提到: 】
: 看来结果还很有规律,分母都相同,谢谢你了,大侠




────────────────────────────────────────
  wolfooo (wolfooo)                于  (Fri Sep 24 07:55:45 2004)  说道:

如何验证这个结果是否正确
请各位大侠各抒高见
【 在 billlang (天也...也天) 的大作中提到: 】
: matlab下运算如下：
: ans =
: [   -(-m*q+n*p)/(k1*m*q-k1*n*p-l*k2*q+l*k3*p),  (-k2*q+k3*p)/(k1*m*q-k1*n*p-l*
: ...................