function?[rho_aphase]=MT1D_FWD（rhoh）
%解析法，層狀介質的大地電磁響應
%rho：各層電阻率
%h：層厚度
%rho_a：視電阻率
%phase：相位
mu=（4e-7）*pi;
T=logspace（-3440）;
k=zeros（size（rho2）size（T2））;
for?N=1:size（rho2）
????k（N:）=sqrt（-i*2*pi*mu./（T.*rho（N）））;
end
m=size（rho2）;
z=-（i*mu*2*pi）./（T.*k（m:））;
for?n=m-1:-1:1
????A=-（i*mu*2*pi）./（T.*k（n:））;
????B=exp（-2*k（n:）*h（n））;
????z=A.*（A.*（1-B）+z.*（1+B））./（A.*（1+B）+z.*（1-B））;
end
rho_a=（T./（mu*2*pi））.*（abs（z）.^2）;
phase=-atan（imag（z）./real（z））.*180/pi;
%?subplot（211）;
%?semilogx（Trho_a‘-*‘）
%?xlabel（‘T（s）‘）
%?ylabel（‘\rho_a（\Omega\cdotm）‘）
%?subplot（212）;
%?semilogx（Tphase‘-*‘）
%?xlabel（‘T（s）‘）
%?ylabel（‘\circ）‘）

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件?????????484??2013-10-13?11:31??MT_forward.m
?????文件????????1452??2013-10-13?11:11??nonlin_inv.m
?????文件?????????709??2013-10-12?19:43??MT1D_FWD.m