convert.m

航天工程工具箱

function y=convert(x,u1,u2)
%CONVERT  Convert between units from different systems.
%   Y = CONVERT(X,U1,U2) converts number X from unit U1
%   to unit U2 if possible.
%
%   For example,
%      CONVERT(15,'gallon_US_liq','barrel_UK')
%   converts the weight 15 US gallons to UK barrels.
%
%   See also CONSTANT.

% Copyright (c) 2001-12-07, B. Rasmus Anthin.

%Length [m]
m=1;
A=1e-10;
XU=1.00208e-13;
fermi=1e-15;
au_l=.529177*A;
AU=1.49597870e11;
ly=9.46055e15;
pc=3.0857e16;
mile_n=1852;
mile=1609.344;
yd=.9144;
ft=.3048;
in=2.54e-2;
fot=.29690;

%Area [m^2]
m2=1;
barn=1e-28;
tunnland=4936;
acre=4046.86;
ha=1e4;
sqin=.64516e-3;
in2=sqin;
sqft=92.90304e-3;
ft2=sqft;
sqyd=.83612736;
yd2=sqyd;
sqmile=2.5899881e6;
mile2=sqmile;

%Volume [m^3]
m3=1;
l=1/1000;
barrel_UK=163.66*l;
barrel_US_pet=158.98*l;
barrel_US_liq=119.24*l;
gallon_UK=4.54609*l;
gallon_US_liq=3.785411784*l;
gallon_US=gallon_US_liq;
pint_UK=.56825*l;
pint_US_dry=.55060*l;
pint_US_liq=.47316*l;
cuin=.016387064*l;
in3=cuin;
cuft=28.316847*l;
ft3=cuft;
cuyd=.76455486;
yd3=cuyd;

%Angle [rad]
rad_a=1;
deg=pi/180;
o=deg;
min_a=deg/60;
sec=min_a/60;
gon=pi/200;
grad=gon;

%Time [s]
s=1;
min_t=60;
h=3600;
day=24*h;
d=day;
yr_trop=365.24219878*day;
yr=yr_trop;
yr_sid=365.25637*day;
yr_cal=365*day;
yr_leap=366*day;
au_t=2.418884e-17;

%Frequency [Hz]
Hz=1;
s01=1;
rpm=1/60;

%Speed [m/s]
m_s=1;
mph=.44704;
knot=.51444;
knop=knot;
km_h=1/3.6;
fps=.3048;

%Mass [kg]
kg=1;
u=1.66054e-27;
lb=.45359237;
lbm=6.48e-5;
slug=14.593903;
ton=1000;
tn=1016.0469088;
shtn=907.18474;
cwt=50.80234544;
shcwt=45.359237;
grain=.0000648;
oz=.02835;

%Density [mol/m^3]
mol_m3=1;
amagat=.04096/1000;

%Density [kg/m^3]
kg_m3=1;
lb_ft3=16.0185;
lb_in3=27.6799e3;

%Temperature [K]

%Energy [J]
J=1;
eV=1.60217733e-19;
Ry=2.1799e-18;
au_W=2*Ry;
erg=1e-7;
kWh=3.6e6;
cal=4.1868;
kcal=1000*cal;
Btu=1055.06;
quad=10^15*Btu;
Q=1000*quad;

%Power [W]
W=1;
hk=735.5;
hp=745.7;

%Angular momentum [Js]
Js=1;
au_am=1.0545727e-34;

%Force [N]
N=1;
dyn=1e-5;
kp=9.80665;
lbf=4.4482;

%Pressure [Pa]
Pa=1;
torr=1.33322e2;
atm=760*torr;
bar=1e5;
at=9.80665e4;
psi=6.8948e3;

%Current [A]
%A=1;
%abamp=10;
%statamp=3.3356e-10;

%Ampl/Att [Np]
Np=1;
B=log(10)/2;
dB=B/10;

%Electric dipole moment [C m]
Cm=1;
au_de=8.47828e-30;
D=3.3356e-30;
esu=1e18*D;

%Magnetic flux density [T]
T=1;
Vs_m2=1;
Vsm02=1;
G=1e-4;
gamma=1e-9;
gauss=1e-4;

%Magnetic flux [Wb]=[V s]
Wb=1;
Vs=1;
Mx=1e-8;

%Magnetizing field [A/m]
A_m=1;
oersted=1e3/4/pi;

%Magnetic dipole moment [Am^2]=[J/T]
Am2=1;
J_T=1;
au_dm=1.8548e-23;

%Activity [Bq]
Bq=1;
Ci=3.7e10;
Rd=1e6;

%Exposure [C/kg]
C_kg=1;
R=2.58e-4;

%Absorbed dose [J/kg]
J_kg=1;
Gy=1;
rad_d=1e-2;
%Dose equivalent [J/kg]
Sv=1;
rem=1e-2;

%Luminous intensity [cd]
cd=1;
hefner=.9;

%Luminance [cd/m^2]
cd_m2=1;
sb=1e4;
asb=1/pi;
cd_ft2=10.76391;
lambert=1e4/pi;

%Gas exposure [Pa s]
Pas=1;
torrs=1.33e2;
L=1.33e-4;

%Effect flux
W_m2Hz01=1;
Wm02Hz01=1;
W_m2_Hz=1;
W_m2s=1;
Wm02s=1;
Jy=1e-26;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

u1=strrep(u1,' ','');
u2=strrep(u2,' ','');
uu={u1 u2};
uu=strrep(uu,'-','0');
uu=strrep(uu,'^','');
uu=strrep(uu,'*','');
uu=strrep(uu,'/','_');
uu=strrep(uu,'.','_');
uu=strrep(uu,'