spbs3d_ybbh.m

在matlab里实现三维仿真布里渊散射信号随温度变化的趋势

%三维图SPBS应变变化
z1=0:30:3000;
z2=3000:5:3500;
z3=3500:5:4000;
f=11100:1:11300;
fB1=11200;
fB2=11250;
fB3=11200;
fBi=35;
W=150e-9;
Pp0=4e-3;
a=0.046e-3;
c=3e+8;  

[F,Z1]=meshgrid(f,z1);
n1=1.45; 
v1=c/n1;
S1=14.297510106e-4; 
a_b1=1.07714800197098e-6;  
P1=a_b1*S1*Pp0*exp(-2*a*Z1)*v1*W/2;
Q1=1./(1+((F-fB1)/(fBi/2)).^2);
I1=P1.*Q1;

[F2,Z2]=meshgrid(f,z2);
n2=1.45;
v2=c/n2;
S2=14.2912606e-4;
a_b2=1.0753884726965e-6;                                                      
P2=a_b2*S2*Pp0*exp(-2*a*Z2)*v2*W/2;
Q2=1./(1+((F-fB2)/(fBi/2)).^2);
I2=P2.*Q2;

[F3,Z3]=meshgrid(f,z3);
n3=1.45;
v3=c/n3;
S3=14.297510106e-4; 
a_b3=1.07714800197098e-6; 
P3=a_b3*S3*Pp0*exp(-2*a*Z3)*v3*W/2;
Q3=1./(1+((F-fB3)/(fBi/2)).^2);
I3=P3.*Q3;

mesh(Z1,F,I1);
hidden off
hold on
mesh(Z2,F,I2);
hidden off
hold on
mesh(Z3,F,I3);
hidden off
hold on