电磁分析matlab,瞬变电磁一维解析解MATLAB代码

【实例简介】研究了大回线圈激励水平分层媒质时,非零偏移距处磁场响应的数值计算以及实现拟地震的偏移成像技术的关键问题-提高波场变换方程数值解的分辨能力.该文首先讨论了半径较大的细圆环天线在自由空间中的瞬态辐射场脉冲宽度分布.在线圈直径远小于激励电流脉冲宽度与光速乘积的条件下,将大回线圈的辐射等效为由磁偶极子所构成的口径面的辐射.

【实例截图】

【核心代码】

% 中心回线中心点感应电动势解析解

% H=0是地面瞬变电磁法 均匀大地解析解

clear;clc;

a=70/(pi^0.5);

I0=1;

u0=4*pi*10^-7;

p1=50;

ConE=1./p1;

HzA=zeros(5000,1);

t=logspace(-8,-1,5000);

for m=1:5000

% u1=a*(u0*ConE/2/t(m))^0.5;

u=a*(u0*ConE/t(m))^0.5/2.;

Hzt(m)=(I0*p1/a.^3)*(3*erf(u)-2*u*(3 2*u^2)*exp(-u^2)/sqrt(pi)); %均匀大地磁场对时间的导数

Hz(m)=I0*(3/sqrt(pi)/u*exp(-u^2) (1-3/2/u^2)*erf(u))/2/a; %均匀大地磁场

end

figure(1);

loglog(t,Hz,'-b','LineWidth',2) % 解析解

hold on;

% tt=load('ttt.dat'); % 有限差分数值解

% hz=load('hhh.dat');

% loglog(tt,hz,'-r','LineWidth',2)