quantities��荅������帥��
���:http://packages.python.org/quantities/�� 荐潟��������


1.莎桁��

hercules:~% python
Python 2.7.3 (default, Jan  2 2013, 16:53:07)
[GCC 4.7.2] on linux2
Type "help", "copyright", "credits" or "license" for more information.


>>> import quantities as pq
pq����physical quantities�� or ��python quantities�� = ����������憗���������ゃ�純�渇���(�����純�潟�����㏍�違�����潟�域�茯�)����

m(=�＜�若����)��ft(=���ｃ�若��)�����ｃ�����������篏���茯㏍�粋昭����?


2.膩雁�

>>> q = 1 * pq.m����������:紊���q��� 臂���q=1�(pq��m)
>>> print q���������� ����:紊���q���阪��
1.0 m������������������������:紊���q��1.0 m(�＜�若����)


3.紊���

>>> q.units = pq.ft��������:紊���q����篏���ft(���ｃ�若��)�����
>>> print q
3.28083989501 ft������������:紊���q��3.28083989501 ft(���ｃ�若��)

50kpc��
>>> q = 50000 * pq.pc
>>> q.units = pq.m
>>> q
array(1.542840125e+21) * m

1.542840125�10^21 m���������


J(=joule)��W(=watt)����������
>>> q = 10 * pq.joule
>>> q.units = pq.watts
ValueError: Unable to convert between units of "J" and "W"
>>> print q
10.0 J


4.荐�膊�

>>> q = (10 * pq.meter)**3 ��������:鐚�鐚������
>>> q2 = q/(5*pq.meter)**2
>>> print q2
40 m

��篏������ｃ��������
>>> q = (10 * pq.meter)**3
>>> q2 = q/(5*pq.ft)**2
>>> print q2������������������������:�����障�障�с������
40.0 m**3/ft**2
>>> q3 = q2.simplified��������������:��篏���輝筝�(��膣���)
>>> print q3
430.556416668 m


5.��篏�

>>> print pq.J.simplified�� ��������: J����膣�������������������鐚�
1.0 kg*m**2/s**2

>>> pq.set_default_units('cgs') ����:荐㊤�紊��眼���с����
>>> print pq.J.simplified
10000000.0 g*cm**2/s**2

��篏�絨���
>>> uK = pq.UnitQuantity('microkelvin', pq.degK/1e6, symbol='uK')
>>> print uK
1 uK (microkelvin)
>>> q = 1000*uK
>>> print q.simplified
0.001 K



篏��ゅ���������ｃ������綺������荐�膊�

hercules:~> ipython ������������������������������������������������������:ipython�с��莎桁���с����
In [1]: import quantities as pq
In [2]: f = 1e-12 * pq.erg / pq.sec / pq.cm**2 ��������������������������:flux���ュ����f=1*10^-12 erg/s/cm^2
In [3]: r = 1 * pq.kilo * pq.parsec ��������������������������������������:莊��≪���ュ����r=1 kpc
In [4]: f, r ����������������������������������������������������������������:腆肴�
Out[4]: (array(1e-12) * erg/(cm**2*s), array(1000.0) * pc)
In [5]: luminosity = 4 * pq.pi * r**2 * f ��������������������������������: ��綺����臂���L=4�r^2f
In [6]: luminosity ��������������������������������������������������������:腆肴�
Out[6]: array(1.2566370614359172e-05) * pc**2*erg/(cm**2*s)
In [7]: luminosity.units = pq.erg / pq.sec ����������������������������:��篏���������
In [8]: luminosity ��������������������������������������������������������:腆肴�
Out[8]: array(1.1964972523338386e+32) * erg/s��������������������:腟���