【磁场】基于Matlab实现磁场正演

​✅作者简介：热爱科研的Matlab仿真开发者，修心和技术同步精进，matlab项目合作可私信。

🍎个人主页：Matlab科研工作室

🍊个人信条：格物致知。

更多Matlab仿真内容点击👇

智能优化算法  神经网络预测 雷达通信  无线传感器

信号处理 图像处理 路径规划 元胞自动机 无人机

⛄ 内容介绍

本文主要是利用MATLAB软件对磁性体的磁场进行正演模拟,了解磁性体的一般分布规律以及影响磁性体磁场的主要因素(如物性参数,走向等),掌握磁性体磁场正演的计算方法,为后续更为复杂的研究奠定一定的基础.

⛄ 部分代码

clear;

clf;

clc;

% 测点范围

x=-100:5:100;          % X方向范围

y=-1000:5:1000;        % Y方向范围

[X,Y]=meshgrid(x,y);   % 转化为排列

% 板状参数

i2=90/180*pi;      %板状体倾斜角度      

b=10;             %板状体宽度 m

j=0.5  ;          %磁极化强度 Wb/m^2

h=30;             %顶部埋深m

% Za理论板状体磁异常

Za=2*j*sin(i2).*(atan((X+b)./h)-atan((X-b)./h));

% Hax理论板状体磁异常

Hax=log((X.^2+b^2-X.*2*b+h^2)./(X.^2+b^2+X.*2*b+h^2)).*sin(i2)*j;

%ΔT理论板状体异常

T=Hax*cos(i2)+Za*sin(i2);

%绘平面异常等值线图(二维)

figure(1),clf,

subplot(221),

contourf(X,Y,Hax);xlabel('X(m)'),ylabel('Y(m)'),title('理论板状体Hax异常');

axis equal,axis([-50 50 -50 50]),colorbar;

subplot(222),

contourf(X,Y,Za);xlabel('X(m)'),ylabel('Y(m)'),title('理论板状体Za异常');

axis equal,axis([-50 50 -50 50]),colorbar;

subplot(223),

contourf(X,Y,T);xlabel('X(m)'),ylabel('Y(m)'),title('理论板状体ΔT异常');

axis equal,axis([-50 50 -50 50]),colorbar;

%绘制曲面图（三维）

figure(2),clf,

subplot(221),mesh(X,Y,Hax),shading interp,xlabel('X(m)'),ylabel('Y(m)'),zlabel('理论板状体Hax异常'),colorbar;

subplot(222),mesh(X,Y,Za),shading interp,xlabel('X(m)'),ylabel('Y(m)'),zlabel('理论板状体Za异常'),colorbar;

subplot(223),mesh(X,Y,T),shading interp,xlabel('X(m)'),ylabel('Y(m)'),zlabel('理论板状体ΔT异常'),colorbar;

⛄ 运行结果

⛄ 参考文献

[1]廖一鹏. 基于MATLAB的磁性体磁场正演模拟研究[J]. 当代化工研究, 2018(8):2.

❤️ 关注我领取海量matlab电子书和数学建模资料

❤️部分理论引用网络文献，若有侵权联系博主删除