基本信息
源码名称:瞬变电磁一维解析解MATLAB代码
源码大小:0.64KB
文件格式:.m
开发语言:MATLAB
更新时间:2019-06-04
友情提示:(无需注册或充值,赞助后即可获取资源下载链接)
嘿,亲!知识可是无价之宝呢,但咱这精心整理的资料也耗费了不少心血呀。小小地破费一下,绝对物超所值哦!如有下载和支付问题,请联系我们QQ(微信同号):813200300
本次赞助数额为: 2 元×
微信扫码支付:2 元
×
请留下您的邮箱,我们将在2小时内将文件发到您的邮箱
源码介绍
研究了大回线圈激励水平分层媒质时,非零偏移距处磁场响应的数值计算以及实现拟地震的偏移成像技术的关键问题-提高波场变换方程数值解的分辨能力.该文首先讨论了半径较大的细圆环天线在自由空间中的瞬态辐射场脉冲宽度分布.在线圈直径远小于激励电流脉冲宽度与光速乘积的条件下,将大回线圈的辐射等效为由磁偶极子所构成的口径面的辐射.
研究了大回线圈激励水平分层媒质时,非零偏移距处磁场响应的数值计算以及实现拟地震的偏移成像技术的关键问题-提高波场变换方程数值解的分辨能力.该文首先讨论了半径较大的细圆环天线在自由空间中的瞬态辐射场脉冲宽度分布.在线圈直径远小于激励电流脉冲宽度与光速乘积的条件下,将大回线圈的辐射等效为由磁偶极子所构成的口径面的辐射.
% 中心回线中心点感应电动势解析解 % 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)