搜搜考试网请问这段程序出了什么问题,本人菜鸟,
来源：学生作业帮 编辑：搜搜考试网作业帮 分类：综合作业 时间：2024/07/08 11:46:20
请问这段程序出了什么问题,本人菜鸟,
clear;
clc;
X=0:1:180;
y=0:1:100;
r=2; %半径
rd=2; %密度
rm=(4/3)*pi*r^3*rd; %质量
x1=75;
y1=30;
h1=5;
x2=95;
y2=60;
h2=5;
x3=85;
y3=40;
h3=8;
x4=105;
y4=70;
h4=8;
G=66.7;
for m=1:1:101
for n=1:1:181
gg1(m,n)=G*rm*h1/((x(n)-x1)^2+(y(m)-y1)^2+h1^2)^(3/2);
gg2(m,n)=G*rm*h2/((x(n)-x2)^2+(y(m)-y2)^2+h2^2)^(3/2);
gg3(m,n)=G*rm*h3/((x(n)-x3)^2+(y(m)-y3)^2+h3^2)^(3/2);
gg4(m,n)=G*rm*h4/((x(n)-x4)^2+(y(m)-y4)^2+h4^2)^(3/2);
gg(m,n)=gg1(m,n)+gg2(m,n)+gg3(m,n)+gg4(m,n);
end
end
contour(gg);
clear;
clc;
X=0:1:180;
y=0:1:100;
r=2; %半径
rd=2; %密度
rm=(4/3)*pi*r^3*rd; %质量
x1=75;
y1=30;
h1=5;
x2=95;
y2=60;
h2=5;
x3=85;
y3=40;
h3=8;
x4=105;
y4=70;
h4=8;
G=66.7;
for m=1:1:101
for n=1:1:181
gg1(m,n)=G*rm*h1/((x(n)-x1)^2+(y(m)-y1)^2+h1^2)^(3/2);
gg2(m,n)=G*rm*h2/((x(n)-x2)^2+(y(m)-y2)^2+h2^2)^(3/2);
gg3(m,n)=G*rm*h3/((x(n)-x3)^2+(y(m)-y3)^2+h3^2)^(3/2);
gg4(m,n)=G*rm*h4/((x(n)-x4)^2+(y(m)-y4)^2+h4^2)^(3/2);
gg(m,n)=gg1(m,n)+gg2(m,n)+gg3(m,n)+gg4(m,n);
end
end
contour(gg);
亲,matlab大小写要分清的,第三行~~x=0:1:180;就好了~~