我想用matlab做个磁法成像的程序来实现地下磁性球体的二维成像并作出平面等值线图！但存在问题！

num2=0;
num3=0;
num4=0;
u0=4pi10^(-7); %u0是绝对磁导率
m=09; %质量m赋初值
z0=5; %地质体埋深
yq=1;
for zq=1:1:10;
for xq=1:1:20;
for i=1:1:20;
xi=i; % 再考虑一下 xi yi zi 的赋值
yi=1;
zi=1;

r=((xq-xi)^2+(yq-yi)^2+(zq-zi)^2)^(1/2); % %%matlab里括号只用圆括号，方括号和大括号有其他意义。 下2行同
T(xi,yi,zi)=u0/(4pi)m/((xi^2+yi^2+(zi-z0)^2)^(5/2))(2(zi-z0)^2-xi^2-yi^2);
%T是测点磁异常
Bq(xi,yi,zi)=(zq-zi)/((xq-xi)^2+(yq-yi)^2+(zq-zi)^(3/2));
%Bq是基函数
num1=num1+T(xi,yi,zi)Bq(xi,yi,zi);
num2=num2+T(xi,yi,zi)^2;
num3=num3+Bq(xi,yi,zi)^2;
num4=(num4num3)^(1/2);
cq=num1/num4; % cq 应该有索引的吧
end
end
end
Plot(xq,yq,cq);

以上就是关于我想用matlab做个磁法成像的程序来实现地下磁性球体的二维成像并作出平面等值线图！但存在问题！全部的内容，包括:我想用matlab做个磁法成像的程序来实现地下磁性球体的二维成像并作出平面等值线图！但存在问题！、、等相关内容解答，如果想了解更多相关内容，可以关注我们，你们的支持是我们更新的动力！