sbs3d_wdbh.m

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

%三维图SBS温度变化
z1=0:30:3000;
z2=3000:5:3500;
z3=3500:5:4000;
f=11100:1:11300;
fB1=11200;
fB2=11260;
fB3=11200;
fBi=35;
W=150e-9;
Pp0=4e-3;
a=0.046e-3;
c=3e+8;  
S=14.1737321199e-4;                                                          
a_b=1.032121625e-6; 
a=0.046e-3;
g1=3.973662998e-12;
g2=3.950567766e-12; 
g3=3.973662998e-12;
A=6.4e-11;
L=4000;
c=3e+8;  
n=1.45;  
v=c/n;   
W=150e-9;  
Pp0=4e-3;
Pcw0=4e-6;
%Psp_z=a_b*S*Pp0*exp(-2*a*z)*v*W/2;
Psp_L1=0.7620e-10;  %Z=L自发散射功率强度
Psp_L2=0.0878e-9;
Psp_L3=0.7620e-10;
[F,Z1]=meshgrid(f,z1);
P1=Psp_L1*exp(a*Z1).*exp(+(g1/A)*Pcw0*(exp(-a*Z1)-exp(-a*L))/a);
Q1=1./(1+((F-fB1)/(fBi/2)).^2);
I1=P1.*Q1;

[F,Z2]=meshgrid(f,z2);
P2=Psp_L2*exp(a*Z2).*exp(+(g2/A)*Pcw0*(exp(-a*Z2)-exp(-a*L))/a);
Q2=1./(1+((F-fB2)/(fBi/2)).^2);
I2=P2.*Q2;

[F,Z3]=meshgrid(f,z3);
P3=Psp_L3*exp(a*Z3).*exp(+(g3/A)*Pcw0*(exp(-a*Z3)-exp(-a*L))/a);
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
mesh(Z3,F,I3);
hidden off