clear all;
format short;
format compact;
syms n w ds df d kd vic b M11 M12 M22 L N A R Q P
syms a1 a2 a3
syms kp1 kp2 kp1re kp1im kp2re kp2im ks ksre ksim vp1 vp2 vs
syms k kR0 kR0re kR0im cR0 kRre kRim c cR
f=100;
w=2*pi*f; % 频率
n=0.1; % 孔隙率
ds=2.65*10^3; % 固体颗粒密度
df=1.0*10^3; % 流体密度
d=(1-n)*ds+n*df; % 土体单元密度
kd=10^(-10); % 动力渗透系数
vic=8*10^(-3); % 流体粘滞系数(viscosity)
b=n^2*vic/kd; % 耗散系数
M11=(1-n)*ds-1i*b/w;
M12=1i*b/w;