声子晶体的传输特性

clear all 
  
a=150e3;                    % 晶格常数 
rhoA=2730;                   % A铝的密度 
rhoB=1180;                   % B环氧树脂的密度 
EA=7.76e10;                  % A铝的杨氏模量 
EB=4.35e9;                   % B环氧树脂的杨氏模量                                                            
cA=（EA/rhoA）^（1/2）;          % A铝中的声速 
cB=（EB/rhoB）^（1/2）;          % B环氧树脂中的声速                    
 
NP=5.5;                      % 
fmax=3e4;                    % 最大频率  
 
  
dx=5e-3;                     % 空间步长 
dt=5e-7;                     % 时间步长 
 
nx=a/dx;                     % 一个晶格中的网格数 
Nx=NP*nx+1;                  % 划分的网格总数 
 
E=EB*ones（1Nx）;             % 1行Nx列的EB      定义所有网格的杨氏模量为树脂的 
rho=rhoB*ones（1Nx）;         % 1行Nx列的rhoB    定义所有网格的密度为树脂的 
                         
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