zoeppritz(zoeppritz)
zoeppritz()所属R语言包：zoeppritz
Zoeppritz Equations
佐普里兹方程
译者：生物统计家园网 机器人LoveR
描述----------Description----------
Calculate the P and S-wave scattering amplitudes for a plane wave at an interface.
计算出的P波和S波散射振幅平面波在接口。
用法----------Usage----------
zoeppritz(icoef, vp1, vp2, vs1, vs2, rho1, rho2, incw)
参数----------Arguments----------
参数：icoef
type of out put Amplitude=1, Potential=2, Energy=3
类型的输出振幅= 1，潜在= 2，能量= 3
参数：vp1
P-wave Velocity of Upper Layer, km/s
上层的P波速度，km / s的
参数：vp2
P-wave Velocity of Lower Layer, km/s
下层的P波速度，km / s的
参数：vs1
S-wave Velocity of Upper Layer, km/s
上层的S波速度，km / s的
参数：vs2
S-wave Velocity of Lower Layer, km/s
km / s的S波速度的较低层，
参数：rho1
Density of Upper Layer, kg/m3
上层的密度，kg/m3的
参数：rho2
Density of Lower Layer, kg/m3
下层的密度，kg/m3的
参数：incw
integer,Incident Wave: P=1, S=2
整数，入射波：P = 1，S = 2
Details
详细信息----------Details----------
Coeficiants are calculated at angles from 0-90 degrees. Zero is returned where coefficients are imaginary.
Coeficiants计算的角度从0-90度。返回零系数是虚构的。
值----------Value----------
List: <table summary="R valueblock"> <tr valign="top"><td>angle</td> <td> Incident angles (degrees)</td></tr> <tr valign="top"><td>rmat</td> <td> Matrix of 4 by n reflection coefficients for each angle</td></tr> <tr valign="top"><td>rra</td> <td> Matrix of 4 by n real part of scattering matrix</td></tr> <tr valign="top"><td>rra</td> <td> Matrix of 4 by n imaginary part of scattering matrix</td></tr> <tr valign="top"><td>ang</td> <td> Matrix of 4 by n phase angle</td></tr> <tr valign="top"><td>incw</td> <td> integer, from input parameter </td></tr> <tr valign="top"><td>icoef</td> <td> integer, from input parameter </td></tr> </table>
列表：<table summary="R valueblock"> <tr valign="top"> <TD> angle </ TD> <TD>入射角（度）</ TD> </ TR> <TR VALIGN =“”> <TD> rmat </ TD> <TD> 4矩阵每个角度的反射系数由n </ TD> </ TR> <tr valign="top"> <TD> rra </ TD> <TD> 4矩阵由n个散射矩阵的实部</ TD> </ TR> <tr valign="top"> <TD> rra</ TD> < TD> 4矩阵由n </ TD> </ TR> <tr valign="top"> <TD>ang </ TD> <TD>由n矩阵的4相位角散射矩阵的虚部</ TD> </ TR> <tr valign="top"> <TD>incw </ TD> <TD>整数，从输入参数</ TD> </ TR> <TR VALIGN =“顶部“> <TD> icoef </ TD> <TD>整数，从输入参数</ TD> </ TR> </ TABLE>
注意----------Note----------
Based on the fortran algorithm in Young and Braile. Uses a linear approximation by Aki and Richards.
基于FORTRAN算法在年轻和Braile的。使用线性逼近Aki和理查兹。
（作者）----------Author(s)----------
Jonathan M. Lees<jonathan.lees@unc.edu>
参考文献----------References----------
Young, G.B., Braile, L. W. 1976. A computer program for the application of Zoeppritz's amplitude equations and Knott's energy equations, Bulletin of the Seismological Society of America, vol.66, no.6,1881-1885.
K. Aki and P.G. Richards.Quantitative seismology. University Science Books, Sausalito, Calif., 2nd edition, 2002.
参见----------See Also----------
pzoeppritz, plotzoeppritz
pzoeppritz，plotzoeppritz
实例----------Examples----------
######### set up 2-layer model[＃＃＃＃＃＃＃＃成立2层模型]
alpha1 = 4.98
beta1 = 2.9
rho1 = 2.667
alpha2 = 8.0
beta2 = 4.6
rho2 = 3.38
################### P-wave incident = 1[＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃P波事件= 1]
incw=1;
icoef=1
A = zoeppritz(icoef, alpha1, alpha2, beta1, beta2, rho1,rho2, incw)
plot(A$angle, A$rmat[,1], xlab="Incident Angle", ylab="Ratio of Amplitudes",
main="P-wave incident/P-wave Reflected" )
plot(A$angle, A$rmat[,2], xlab="Incident Angle", ylab="Ratio of Amplitudes",
main="P-wave incident/S-wave Reflected" )
plot(A$angle, A$rmat[,3], xlab="Incident Angle", ylab="Ratio of Amplitudes",
main="P-wave incident/P-wave Refracted" )
plot(A$angle, A$rmat[,4], xlab="Incident Angle", ylab="Ratio of Amplitudes",
main="P-wave incident/S-wave Refracted" )
################### S-wave incident = 2[＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃S-波入射= 2]
incw=2
icoef=1
A = zoeppritz(icoef, alpha1, alpha2, beta1, beta2, rho1,rho2, incw)
plot(A$angle, A$rmat[,1], xlab="Incident Angle", ylab="Ratio of Amplitudes",
main="S-wave incident/P-wave Reflected" )
plot(A$angle, A$rmat[,2], xlab="Incident Angle", ylab="Ratio of Amplitudes",
main="S-wave incident/S-wave Reflected" )
plot(A$angle, A$rmat[,3], xlab="Incident Angle", ylab="Ratio of Amplitudes",
main="S-wave incident/P-wave Refracted" )
plot(A$angle, A$rmat[,4], xlab="Incident Angle", ylab="Ratio of Amplitudes",
main="S-wave incident/S-wave Refracted" )
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注：
注1：为了方便大家学习，本文档为生物统计家园网机器人LoveR翻译而成，仅供个人R语言学习参考使用，生物统计家园保留版权。
注2：由于是机器人自动翻译，难免有不准确之处，使用时仔细对照中、英文内容进行反复理解，可以帮助R语言的学习。
注3：如遇到不准确之处，请在本贴的后面进行回帖，我们会逐渐进行修订。
|
发表于 2012-10-2 07:51:47
