声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1387|回复: 6

[综合讨论] 求助 关于隐函数的画图问题

[复制链接]
发表于 2009-5-21 19:18 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
最近在编一个程序
就是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 编辑 ]
回复
分享到:

使用道具 举报

发表于 2009-5-21 20:47 | 显示全部楼层
solve('l','a');改为solve(l,a);即可

评分

1

查看全部评分

 楼主| 发表于 2009-5-21 21:06 | 显示全部楼层
还是有一点问题啊
程序运行得到的图像为
2.jpg
而应该要得到的是
1.jpg
会不会是坐标选取不对啊。
求各位的帮助啊。多谢。

[ 本帖最后由 zxchinasun 于 2009-5-21 21:18 编辑 ]
发表于 2009-5-21 21:11 | 显示全部楼层
至少运行上是没有问题,把你的问题说明白,不要说一半留一半。
 楼主| 发表于 2009-5-21 21:19 | 显示全部楼层

回复 地板 xiezhh 的帖子

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

[ 本帖最后由 zxchinasun 于 2009-5-21 22:20 编辑 ]
发表于 2009-5-21 21:24 | 显示全部楼层
不好意思,刚才没有看到,我又看了看,不了解问题的背景,也不知道怎样产生你所要的图形。
 楼主| 发表于 2009-5-21 21:31 | 显示全部楼层

回复 6楼 xiezhh 的帖子

没关系 我再试试吧  还是谢谢了
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-14 17:46 , Processed in 0.073827 second(s), 25 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表