格林函数的公式matlab仿真

%单位点源对应全空间格林函数（时间域）
clc
clear
close all
u=4*pi*10^（-7）;%取为真空磁导率
o=0.01; %电导率S/m
%r确定t变化
r=100;%观测点距点源距，以电源为坐标原点
v=-7:0.05:2;
t=10.^v;
Gt=（（u*o）^.5）./（8*（pi*t）.^1.5）.*exp（-u*o*r^2./（4.*t））;
set（gcf‘Position‘[160 100 1000 500]）;
set（gca‘Position‘[.13 .17 .8 .74]）;
 subplot（121）
loglog（t*1000Gt‘-*r‘‘MarkerSize‘3）;
title（‘g-t‘‘fontsize‘13‘Fontname‘‘Times New Roman‘）;
xlabel（‘t（ms）‘‘fontsize‘13）;
ylabel（‘g‘‘fontsize‘13）
axis（[10^-4 10^1 10^-4 10^1]）;

%t确定r变化
t=1; %ms
x=0;
y=0;
z=1:3:100;
r=（x.^2+y.^2+z.^2）.^（0.5）;%观测点距点源距，以电源为坐标原点
Gr=（（u*o）^.5）*exp（-u*o*r.^2/（4*t））/（8*（pi*t）^（1.5））*（t>=0）;
subplot（122）
plot（rGr‘-*r‘）;
title（‘g-r‘‘fontsize‘13‘Fontname‘‘Times New Roman‘）
xlabel（‘r（m）‘‘fontsize‘13）
ylabel（‘g‘‘fontsize‘13）
grid;

%三维输出
u=4*pi*10^（-7）;%取为真空磁导率
o=0.01; %电导率S/m
r=1:0.5:91;%观测点距点源距，以电源为坐标原点
v=-7:0.05:2;
t=10.^v;
[rt]=meshgrid（rt）;
G=（（u*o）^.5）./（8*（pi*t）.^1.5）.*exp（-u*o*r^2./（4.*t））;
figure（‘Numbertitle‘‘off‘‘name‘‘3D-g（tr）‘）;