求助：绘制等值线图和梯度场

yxy

绘制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);
[x,y]=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)
大家帮帮我吧，谢谢了！！

多情清秋

1. zz=1*ones(1,n);plot3(x,x,zz),grid off,hold on
   
2.     mesh(x,y,z)
   
3.     xlabel('x'),ylabel('y'),zlabel('z'),box on

复制代码

这部分不要不就行了吗？

yxy

这部分是抄的书上的
那么梯度场怎么绘制呢?

yejet

用jacobian求出梯度，然后绘图

yxy

我们书上用的是gradient函数
问题是我不会用啊
试了很久了都

bainhome

like this?

1. function ContourF
   
2. clc
   
3. clear
   
4. clf
   
5. [x,y]=meshgrid(-1:.2:4);
   
6. z=1./(x.^2-2*x+4)+1./(y.^3-2*y+4);
   
7. [px,py]=gradient(z,.2,.2);
   
8. v=-1:.2:4;
   
9. [C,h]=contour(v,v,z);
   
10. set(gca,'nextplot','add');
    
11. quiver(v,v,px,py);

复制代码

[ 本帖最后由 bainhome 于 2006-10-7 20:24 编辑 ]

fxp

you are best!

lizhenjie17

在气象数据处理中用到的风玫瑰图，怎么用MTLAB做出来呢？有谁知道啊？

xuexixiaozi123

画风玫瑰别用matlab,用Grapher比较好，比较专业，画出来也挺美观。

你去网上找一个，之后在图表——2D图表——玫瑰图里，弄好数据，直接就可以画出来了，很简单。

wusemm

请教6楼的程序ContourF at 9有问题
我的matlab运行出现的问题??? Attempt to execute SCRIPT contour as a function.

[ 本帖最后由 ChaChing 于 2010-2-6 21:04 编辑 ]

ChaChing

请使用which contour -all检查下, 看是否有同名自定script!

程瑞岩

来支持一下，刚来看得不太明白。

程瑞岩

这太有才了。