✅作者簡介:熱愛科研的Matlab仿真開發者,修心和技術同步精進,matlab項目合作可私信。
🍎個人首頁:Matlab科研工作室
🍊個人信條:格物緻知。
⛄ 内容介紹
本文扼要的叙述了大地電磁的基本原理,簡要介紹了大地電磁的正反演方法和流程,重點闡述了基于MATLAB軟體平台的實作,并給出了此種方法可擴充的幾個方面。
⛄ 部分代碼
close all
clc
clear
%%% 受Matlab計算精度所限,觀測時間範圍設定在10^(-7)至10^(-2)s以内。
%%% 受Matlab計算精度所限,各層的測試電阻設定在500歐姆以下。
%%% FORTRAN語言才是電磁波場正反演的經典程式設計語言,MATLAB可用于簡單測試
t=logspace(-7,-2,100);
%兩層模型測試
[V1]=TEMforward4([500],[]);
[V2]=TEMforward4([100],[]);
[V3]=TEMforward4([500,100],[50]);
%三層模型測試
[V4]=TEMforward4([200],[]);
[V5]=TEMforward4([50],[]);
[V6]=TEMforward4([500],[]);
[V7]=TEMforward4([200,50,500],[50,100]);
%四層模型測試
[V20]=TEMforward4([500,100,1000,200],[50,200,500]);
%一層模型繪圖
V11=[V1;V2];
figure(1);
loglog(t,V11,'.-');
legend('500歐姆','100歐姆');
xlabel('時間 / s')
ylabel('電壓 / V')
%兩層模型繪圖
V=[V1;V2;V3];
figure(2);
loglog(t,V,'.-');
legend('第一層','第二層','兩層');
xlabel('時間 / s')
ylabel('電壓 / V')
%三層模型繪圖
V12=[V4;V5;V6;V7];
figure(3);
loglog(t,V12,'.-');
legend('第一層','第二層','第三層','三層');
xlabel('時間 / s')
ylabel('電壓 / V')
%四層模型繪圖
figure(4);
loglog(t,V20,'.-');
legend('{\rho}=[500,100,1000,200] H=[50,200,500]');
xlabel('時間 / s')
ylabel('電壓 / V')