如題，運行我的程序，最后的等值線圖出不了，提示：ZData 為非有限時，無法顯示等高線

附上我的代碼，是一個數(shù)值積分：

Sztau=@(x,z)1/4/pi*integral2(@(theta,r)(1./4*r).*(1-1./((1+(L./r).^2).^(3/2)).*(x-r.*cos(theta)).*z./((x.^2-2*x.*r.*cos(theta)+r.^2+z.^2).^(3/2))),eps,2*pi,rmin,inf);

[x,z]=meshgrid(linspace(0,15,150));

f=-reshape(arrayfun(@(k)Sztau(x(k),z(k)),1:length(x(:))),150,150);

%%

%mesh(x,z,f);

[C,h]=contour(x,z,f);

%set(h,'ShowText','on','TextStep',get(h,'LevelStep')*2,'Linewidth',0.5,'LineColor','black')

%set(h,'ShowText','on','LevelStep',0.005,'Linewidth',0.5,'LineColor','black')

set(h,'ShowText','on','LevelList', [-0.5 -0.4 -0.3 -0.1 -0.03 -0.02 -0.01 -0.005 0 0.005 0.01 0.015 0.02 0.03 0.1 0.2 0.3 0.4 0.5],'Linewidth',1.4,'LineColor','black')

set(gca,'DataAspectRatio',[1 1 1],'XDir','reverse','YDir','reverse','FontName','Times New Roman','FontSize',12,'xLim',[0 6],'XTick',0:1:6,'zLim',[0 15],'ZTick',0:1:15)

clabel(C,h,'FontSize',12,'FontName','times new roman','Color','black')

xlabel('x/R','FontSize',16,'FontName','times new roman','FontAngle','italic','Color','black')

ylabel('z/R','FontSize',16,'FontName','times new roman','FontAngle','italic','Color','black')

最后附上我要編的數(shù)值積分式子：

2017-5-15 08:33 上傳

總結(jié)

以上是生活随笔為你收集整理的matlab等高线不显示,Matlab提示ZData 为非有限时，无法显示等高线图的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。