clear?all?
??
a=150e3;????????????????????%?晶格常數(shù)?
rhoA=2730;???????????????????%?A鋁的密度?
rhoB=1180;???????????????????%?B環(huán)氧樹脂的密度?
EA=7.76e10;??????????????????%?A鋁的楊氏模量?
EB=4.35e9;???????????????????%?B環(huán)氧樹脂的楊氏模量????????????????????????????????????????????????????????????
cA=（EA/rhoA）^（1/2）;??????????%?A鋁中的聲速?
cB=（EB/rhoB）^（1/2）;??????????%?B環(huán)氧樹脂中的聲速????????????????????
?
NP=5.5;??????????????????????%?
fmax=3e4;????????????????????%?最大頻率??
?
??
dx=5e-3;?????????????????????%?空間步長?
dt=5e-7;?????????????????????%?時間步長?
?
nx=a/dx;?????????????????????%?一個晶格中的網(wǎng)格數(shù)?
Nx=NP*nx+1;??????????????????%?劃分的網(wǎng)格總數(shù)?
?
E=EB*ones（1Nx）;?????????????%?1行Nx列的EB??????定義所有網(wǎng)格的楊氏模量為樹脂的?
rho=rhoB*ones（1Nx）;?????????%?1行Nx列的rhoB????定義所有網(wǎng)格的密度為樹脂的?
?????????????????????????
for?ii=0:（NP-1）??????????????%?定義聲子晶體?
??for?jj=（nx/2+2）:nx?
??????E（jj+ii*nx）=EA;?????????????????????????????????????????????????????????????
??????rho（jj+ii*nx）=rhoA;?
??end??
end?
?
for?ii=1:（NP*2-1）????????????%?定義兩種材料交界處密度與楊氏模量為平均值???????????????????????????????????????????????????????
????E（ii*nx/2+1）=（EB+EA）/2;??????
????rho（ii*nx/2+1）=（rhoB+rhoA）/2;??
end?