﻿功能:力学求惯性矩(Mechanical moment of inertia)

格式:mechmoi(x,y,dx,dy);mechmoi(x,y)

说明:x,y是同维数的矩阵变量,其里面的数据依次存放围成凸多边形【对于是凹多边形,可以先把图形变成凸多边形,利用本函数求变化后的凸多边形的惯性矩,再减去凸多边形与凹多边形求差集的图形区域的惯性矩即可,对于有曲线边界的封闭图形,本函数无能为力】的点的x与y坐标.因为要组成平面图形,因此x,y里面的数据不能少于3个.dx、dy是一条直线上的2个点,dx里是存储2个x坐标,dy存储2个点的坐标.例如dx=[x1,x2],dy=[y1,y2]那么求得的惯性矩是以(x1,y1)为原点方向(x1-x2,y1-y2)为我们认为惯性矩Ix的方向.dx与dy默认组成的直线为(0,0)与(0,10),即默认是当前的x轴,就是只接收2个参数的格式时的默认值.本函数执行完毕返回4个参数,其依次代表Ix,Iy,Ip

例子:

x =
[  -4.0000000000000  4.00000000000000  4.00000000000000  -4.0000000000000  ]
y =
[  -5.0000000000000  -5.0000000000000  5.00000000000000  5.00000000000000  ]
dx =
[  0.00000000000000  1.00000000000000  ]
dy =
[  0.00000000000000  0.00000000000000  ]

c=mechmoi(x,y,dx,dy)//注意这里的惯性矩是默认以dx与dy组成的直线为新的x轴,对应的轴求得的是Ix,而dx,dy里的第一点变换后为原点,程序进而确定变换后的y轴,从而求出Iy,Ixy,而Ip=Ix+Iy
c =
[ 666.666666666667//=Ix
  426.666666666667//=Iy
  1093.33333333333  ]//=Ip

//By 2011/10/29