瞬变电磁TEM层状介质正演（TEM forward）

瞬变电磁TEM层状介质正演.基于Matlab平台。地球物理学

function [V1] = TEMforward4（RHOH）
% RHO为各层电阻，H为各层厚度
% V为观测电压
%%%%层状大地磁偶源瞬变电磁场正演   
%%% Liuweiqiang CSU
a=100;              %发送回线半径
r=1;                %接收回线半径
q=r.^2*pi*100;      %接收线圈的有效面积
miu0=4*pi*10^（-7）;
N=length（RHO）;
w0=0.000001:200:2000000;%原则上最大值越大，精度越大越好。受计算速度限制，如此设置
N1=length（w0）;
for jjj=1:N1
    w=w0（jjj）;
for nn=1:47
a0=-3.05078187595e+00;
s0=1.10599010095e-01;   
lamda=（1/a）*10^（a0+（nn-1）*s0）;  
ramta=lamda;
I0=10;%电流
hh=0;%线圈高度
if N==1;
u（N）=（ramta*ramta-j*w/RHO（N）*miu0）.^（1/2）;
z（N）=-j*w*miu0/u（N）;
Z（N）=z（N）;
elseif N~=1
for ii=1:N-1    
u（N）=（ramta*ramta-j*w/RHO（N）*miu0）.^（1/2）;
z（N）=-j*w*miu0/u（N）;
Z（N）=z（N）;
u（N-ii）=（ramta*ramta-j*w/RHO（N-ii）*miu0）.^（1/2）;
z（N-ii）=-j*w*miu0/u（N-ii）;
Z（N-ii）=z（N-ii）*（Z（N-ii+1）+z（N-ii）*tanh（u（N-ii）*H（N-ii）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        2500  2015-01-29 22:43  TEM层状磁偶源正演\TEMforward4.m
     文件        1121  2014-12-04 13:54  TEM层状磁偶源正演\zhengyanceshi.m
     文件       26355  2014-12-10 11:31  TEM层状磁偶源正演\代码说明.docx
     目录           0  2015-01-29 22:43  TEM层状磁偶源正演\