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