clear all;
close all;
clc;
h0=0.01;
beta=0.000391*10^(-8);
loss=0.00871;
density=1.150;
v=0.989;
delta2=0.002;
gam=3.25;
c=1500;
f=20*10^3;
lamda=c/f;
k=2*pi/lamda;
theta=(10:1:90)*pi/180;
[delta_whole,sb_kr,sb_cr,sb_mr,sb_v,sb]=Sb(h0,beta,loss,density,v,delta2,gam,k,theta);
plot(theta*180/pi,sb);
hold on;

h0=0.01;
beta=0.002166*10^(-8);
loss=0.01327;
density=1.49;
v=0.989;
delta2=0.005;
gam=3.18;
c=1500;
f=20*10^3;
lamda=c/f;
k=2*pi/lamda;
theta=(10:1:90)*pi/180;
[delta_whole,sb_kr,sb_cr,sb_mr,sb_v,sb]=Sb(h0,beta,loss,density,v,delta2,gam,k,theta);
plot(theta*180/pi,sb,'--');
grid on;