求助 关于隐函数的画图问题

zxchinasun

最近在编一个程序
就是h是一个三乘三的矩阵，其中包括an和a 两个未知量。
然后我想用solve解det（h）=0得出an和a的关系，然后在画图，可程序总是出错。求各位帮帮忙，不胜感激。
程序如下：
clear;
c11=8.674;c33=10.72;c44=5.794;c12=0.699;c13=1.191;c14=-1.791;md=26.51;c66=(c11-c12)/2;
syms a  an  
lx=0; ly=sin(an); lz=cos(an);     
b=c11*(lx^2)+c66*(ly^2)+c44*(lz^2)+2*c14*(ly*lz);
c=c66*(lx^2)+c11*(ly^2)+c44*(lz^2)-2*c14*ly*lz;
d=c44*(lx^2+ly^2)+c33*(lz^2);
e=(c12+c66)*lx*ly+2*c14*lx*lz;
f=(c13+c44)*lx*lz+2*c14*lx*ly;
g=(c13+c44)*ly*lz+c14*(lx^2-ly^2);
h=[(a^2)*b-md (a^2)*e (a^2)*f;
    (a^2)*e (a^2)*c-md (a^2)*g;
    (a^2)*f (a^2)*g (a^2)*d-md];
l=det(h);
solve('l','a');
an=(0:360)*pi/180;
ezplot(det(h));
grid on

[ 本帖最后由 ChaChing 于 2010-6-15 15:30 编辑 ]

xiezhh

solve('l','a');改为solve(l,a);即可

zxchinasun

还是有一点问题啊
程序运行得到的图像为

而应该要得到的是

会不会是坐标选取不对啊。
求各位的帮助啊。多谢。

[ 本帖最后由 zxchinasun 于 2009-5-21 21:18 编辑 ]

xiezhh

至少运行上是没有问题，把你的问题说明白，不要说一半留一半。

zxchinasun

已经说了 在2楼 希望能再看看  多谢多谢。
在这里我是画了一个慢度图。

[ 本帖最后由 zxchinasun 于 2009-5-21 22:20 编辑 ]

xiezhh

不好意思，刚才没有看到，我又看了看，不了解问题的背景，也不知道怎样产生你所要的图形。

zxchinasun

没关系 我再试试吧  还是谢谢了