求助:绘制等值线图和梯度场
绘制z=1/(x^2-2x+4)+1/(y^3-2y+4)的等值线和梯度场我的程序:
clear,clf,
n=20;
x=linspace(-8,8,20);y=linspace(-8,8,20);
=meshgrid(x,y);
for i=1:n
for j=1:n
z(i,j)=1/(x(i,j).^2-2.*x(i,j)+4)+1/(y(i,j).^3-2.*y(i,j)+4);
end
end
zz=1*ones(1,n);plot3(x,x,zz),grid off,hold on
mesh(x,y,z)
xlabel('x'),ylabel('y'),zlabel('z'),box on
contourf(x,y,z,n)
大家帮帮我吧,谢谢了!! zz=1*ones(1,n);plot3(x,x,zz),grid off,hold on
mesh(x,y,z)
xlabel('x'),ylabel('y'),zlabel('z'),box on
这部分不要不就行了吗? 这部分是抄的书上的
那么梯度场怎么绘制呢? 用jacobian求出梯度,然后绘图 我们书上用的是gradient函数
问题是我不会用啊
试了很久了都 like this?
function ContourF
clc
clear
clf
=meshgrid(-1:.2:4);
z=1./(x.^2-2*x+4)+1./(y.^3-2*y+4);
=gradient(z,.2,.2);
v=-1:.2:4;
=contour(v,v,z);
set(gca,'nextplot','add');
quiver(v,v,px,py);
[ 本帖最后由 bainhome 于 2006-10-7 20:24 编辑 ]
Thank you!
you are best!有哪位能做气象上面的风玫瑰图啊?
在气象数据处理中用到的风玫瑰图,怎么用MTLAB做出来呢?有谁知道啊?回复 8楼 lizhenjie17 的帖子
画风玫瑰别用matlab,用Grapher比较好,比较专业,画出来也挺美观。你去网上找一个,之后在图表——2D图表——玫瑰图里,弄好数据,直接就可以画出来了,很简单。
6楼的程序有问题
请教6楼的程序ContourF at 9有问题我的matlab运行出现的问题??? Attempt to execute SCRIPT contour as a function.
[ 本帖最后由 ChaChing 于 2010-2-6 21:04 编辑 ]
回复 10楼 wusemm 的帖子
请使用which contour -all检查下, 看是否有同名自定script! 来支持一下,刚来看得不太明白。 这太有才了。
页:
[1]