1. 球体磁异常正演模拟
% 球体参数
R = 10; % 球体半径 (m)
chi = 0.01; % 球体磁化率
inc = 30; % 磁倾角 (度)
dec = 0; % 磁偏角 (度)% 地球磁场参数
B0 = 50000; % 地球磁场强度 (nT)
inc0 = 60; % 地球磁场倾角 (度)
dec0 = 0; % 地球磁场偏角 (度)% 计算磁场分布
[x, y, z] = meshgrid(-100:5:100, -100:5:100, -100:5:100);
[Bx, By, Bz] = sphere_mag(R, chi, inc, dec, x, y, z, B0, inc0, dec0);% 计算磁异常
[Bx0, By0, Bz0] = sphere_mag(R, 0, inc, dec, x, y, z, B0, inc0, dec0);
dBx = Bx - Bx0;
dBy = By - By0;
dBz = Bz - Bz0;
dB = sqrt(dBx.^2 + dBy.^2 + dBz.^2);% 绘制磁异常图像
figure;
slice(x, y, z, dB, [], [], [-50, 0, 50]);
colorbar;
xlabel('X');
ylabel('Y');
zlabel('Z');
title('Sphere Magnetic Anomaly');
2. 圆柱体磁异常正演模拟
% 圆柱体参数
R = 10; % 圆柱体半径 (m)
L = 50; % 圆柱体长度 (m)
chi = 0.01; % 圆柱体磁化率
inc = 30; % 磁倾角 (度)
dec = 0; % 磁偏角 (度)% 地球磁场参数
B0 = 50000; % 地球磁场强度 (nT)
inc0 = 60; % 地球磁场倾角 (度)
dec0 = 0; % 地球磁场偏角 (度)% 计算磁场分布
[x, y, z] = meshgrid(-100:5:100, -100:5:100, -100:5:100);
[Bx, By, Bz] = cylinder_mag(R, L, chi, inc, dec, x, y, z, B0, inc0, dec0);% 计算磁异常
[Bx0, By0, Bz0] = cylinder_mag(R, L, 0, inc, dec, x, y, z, B0, inc0, dec0);
dBx = Bx - Bx0;
dBy = By - By0;
dBz = Bz - Bz0;
dB = sqrt(dBx.^2 + dBy.^2 + dBz.^2);% 绘制磁异常图像
figure;
slice(x, y, z, dB, [], [], [-50, 0, 50]);
colorbar;
xlabel('X');
ylabel('Y');
zlabel('Z');
title('Cylinder Magnetic Anomaly');
3. 板状体磁异常正演模拟
% 板状体参数
Lx = 50; % 板状体长度 (m)
Ly = 30; % 板状体宽度 (m)
Lz = 10; % 板状体厚度 (m)
chi = 0.01; % 板状体磁化率
inc = 30; % 磁倾角 (度)
dec = 0; % 磁偏角 (度)% 地球磁场参数
B0 = 50000; % 地球磁场强度 (nT)
inc0 = 60; % 地球磁场倾角 (度)
dec0 = 0; % 地球磁场偏角 (度)% 计算磁场分布
[x, y, z] = meshgrid(-100:5:100, -100:5:100, -100:5:100);
[Bx, By, Bz] = plate_mag(Lx, Ly, Lz, chi, inc, dec, x, y, z, B0, inc0, dec0);% 计算磁异常
[Bx0, By0, Bz0] = plate_mag(Lx, Ly, Lz, 0, inc, dec, x, y, z, B0, inc0, dec0);
dBx = Bx - Bx0;
dBy = By - By0;
dBz = Bz - Bz0;
dB = sqrt(dBx.^2 + dBy.^2 + dBz.^2);% 绘制磁异常图像
figure;
slice(x, y, z, dB, [], [], [-50, 0, 50]);
colorbar;
xlabel('X');
ylabel('Y');
zlabel('Z');
title('Plate Magnetic Anomaly');
参考代码 基于matlab编写的板状体、球体和圆柱体的磁异常正演模拟 youwenfan.com/contentcno/80730.html
附加函数
以下是用于计算球体、圆柱体和板状体磁场的附加函数示例：
% 球体磁场计算函数
function [Bx, By, Bz] = sphere_mag(R, chi, inc, dec, x, y, z, B0, inc0, dec0)% 球体磁场计算公式% ...
end% 圆柱体磁场计算函数
function [Bx, By, Bz] = cylinder_mag(R, L, chi, inc, dec, x, y, z, B0, inc0, dec0)% 圆柱体磁场计算公式% ...
end% 板状体磁场计算函数
function [Bx, By, Bz] = plate_mag(Lx, Ly, Lz, chi, inc, dec, x, y, z, B0, inc0, dec0)% 板状体磁场计算公式% ...
end