这是我上次问的问题
这是答案
我现在想尝试把几个图放一块类似于下图但我好菜啊
p=200e3;
z=0:0.1:6;
figure(1)
for r=[0,0.25,0.5,1]
R=sqrt(r*r+z.*z);
sigma=-3*p*z.^3./2./pi./R.^5;
semilogx(abs(sigma),z), hold on %x轴用log表示
end
hold off
xlabel("sigma")
ylabel("z")
axis ij; %翻转y轴坐标
上面这个稍微容易理解一点吧....
下面这个纯用matrix计算了。
p=200e3;
r=[0,0.25,0.5,1];
z=transpose(0:0.1:6);
R=sqrt(r.*r+z.*z);
sigma=-3*p*z.^3./2./pi./R.^5;
figure(1)
semilogx(abs(sigma),z)
xlabel("sigma")
ylabel("z")
axis ij; %翻转y轴坐标
重点不在于画图,重点在于你要知道你要的数据在哪里,要怎么制造你的数据。画图的话搜一搜就行了。