磁偶极子源仿真

clear;
clc;
close all;
%format long;
%Para_Meas;
%% definate forword model
% measurement plane
%format long;
MeasPlane = [002020];
MeasPlaneZ = 0;
MeasNum =36;    %6*6 measuring points MeasNum is for changing
NM = sqrt（MeasNum）;
% meshgrid on measurement plane 
[MMX  MMY]= meshgrid（MeasPlane（1）:（MeasPlane（3）-MeasPlane（1））/（NM-1）:MeasPlane（3）...
    MeasPlane（2）:（MeasPlane（4）-MeasPlane（2））/（NM-1）:MeasPlane（4））;
SrcPara1 = [001775];  %scource parametersfrom left to right are : mxmymz，x0y0z0for changing
SrcPara2 = [00110165]; 
SrcPara3 = [00116105]
[nnxnny] = size（MMX）;
%% 源1模型
lx1 = reshape（MMX‘nnx*nny1）-SrcPara1（4）;
ly1= reshape（MMY‘nnx*nny1）-SrcPara1（5）;
lz1 = MeasPlaneZ*ones（nnx*nny1）-SrcPara1（6）;
R1 = [lx1 ly1 lz1];
S1 = repmat（SrcPara1（1:3）nnx*nny1）;
BZ1_1=3*（10e-7）*R1（:1）.*R1（:3）.*S1（:1）;BZ1_2=3*（10e-7）*R1（:2）.*R1（:2）.*S1（:2）;BZ1_3=（10e-7）*（2*R1（:3）.^2-R1（:1）.^2-R1（:2）.^2）.*S1（:3）;
R_Norm1=（R1（:1）.^2+R1（:2）.^2+R1（:3）.^2）.^（5/2）
BZ1=（BZ1_1+BZ1_2+BZ1_3）./R_Norm1;
BZ1=reshape（BZ1NMNM）;
% Bz = 3*（10e-7）*m.*repmat（-lz13）./repmat（（lx.^2+ly.^2+lz.^2）.^（5/2）13）;
% BZ_1=3*（10e-7）*R（:1）.*R（:2）.*S1（:1）;BZ_2=3*（10e-7）*R（:2）.*R（:2）.*S1（:2）;BZ_3=（10e-7）*（2*R（:3）.^2-R（:1）.^2-R（:2）.^2）.*S（:3）;
% R_Norm=（R（:1）.^2+R（:2）.^2+R（:3）.^2）.^（5/2）
% BZ=BZ_1+BZ_2+BZ_3;
%% 源2模型
lx2 = reshape（MMX‘nnx*nny1）-SrcPara2（4）;
ly2= reshape（MMY‘nnx*nny1）-SrcPara2（5）;
lz2 = MeasPlaneZ*ones（nnx*nny1）-SrcPara2（6）;
R2 = [lx2 ly2 lz2];
S2 = repmat（SrcPara2（1:3）nnx*nny1）;
BZ2_1=3*（10e-7）*R2（:1）.*R2（:3）.*S2（:1）;BZ2_2=3*（10e-7）*R2（:2）.*R2（:2）.*S2（:2）;BZ2_3=（10e-7）*（2*R2（:3）.^2-R2（:1）.^2-R2（:2）.^2）.*S2（:3）;
R_Norm2=（R2（:1）.^2+R2（:2）.^2+R2（:3）.^2）.^（5/2）
BZ2=（BZ2_1+BZ2_2+BZ2_3）./R_Norm2;
BZ2=reshape（BZ2NMNM）;
%% 源3模型
lx3 = reshape（MMX‘nnx*nny1）-SrcPara3（4）;
ly3= reshape（MMY‘nnx*nny1）-SrcPara3（5）;
lz3 = MeasPlaneZ*ones（nnx*nny1）-SrcPara3（6）;
R3 = [lx3 ly3 lz3];
S3 = repmat（SrcPara3（1:3）nnx*nny1）;
BZ3_1=3*（10e-7）*R3（:1）.*R3（:3）.*S3（:1）;BZ3_2=3*（10e-7）*R3（:2）.*R3（:2）.*S3（:2）;BZ3_3=（10e-7）*（2*R3（:3）.^2-R3（:1）.^2-R3（:2）.^2）.*S3（:3）;
R_Norm3=（R3（:1）.^2+R3（:2）.^2+R3（:3）.^2）.^（5/2）
BZ3=（BZ3_1+BZ3_2+BZ3_3）./R_Norm3;
BZ3=reshape（BZ3NMNM）;

% delta=0.2/（NM-1）;
% [xy]=meshgrid（0:delta:0.2）;
% [XIYI]=meshgrid（0:delta/10:0.2）;
% ZI1=interp2（xyBZ1XIYI‘spline‘）;
% BZ1=ZI1;
% ZI2=interp2（xyBZ2XIYI‘spline‘）;
% BZ2=ZI2;
% BZ=BZ1+BZ2+BZ3;
% BZ=reshape（BZ136）
% BZ=repmat（BZ60701）;
% BZ=BZ1;
% BZ=BZ1+BZ2;
BZ=BZ1+BZ2+BZ3;
BZ