✅作者简介：热爱科研的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.
