24小时热门版块排行榜

返回列表

837261051

木虫 (正式写手)

应助: 1 (幼儿园)
金币: 3318.4
散金: 76
红花: 3
帖子: 576
在线: 158.9小时
虫号: 2878270
注册: 2013-12-17
专业: 力学中的基本问题和方法

[求助] 画出一个含有x，y，z三个变量的隐函数图像？不知为什么会报错？已有1人参与

该程序中x，y，z是用o1，o2，o3来代替的，运行到f那个函数表达式后就出现“Error: Unexpected MATLAB operator.”的字样，请大神指点一二。
   由于f的表达式太长，为了便于观看，我分成三行写的

代码如下：
clc
clear
f = @(o1,o2,o3)  (1/(((2*o2.-2*o3.)/(o1.-o3.)-1)^2/3 + 1)^(1/2)-((2*o2.-2*o3.)/(o1.-o3.)-1)/(2*(((2*o2.-2*o3.)/(o1.-o3.)-1)^2/3+1)^(1/2))) ...
*((o1.+o2.+o3.)^2/2-(3*o1.*o2.)/2-(3*o1.*o3.)/2-(3*o2.*o3.)/2)^(1/2)-5*3^(1/2)+(3^(1/2)*((2*o2.-2*o3.)/(o1.-o3.)-1) ...
*(o1./3+o2./3+o3./3))/(3*(((2*o2.-2*o3.)/(o1.-o3.)-1)^2/3+1)^(1/2));       %表达式
[o1,o2,o3] = meshgrid(-10:.2:10,-10:.2:10,-10:.2:10);          % 画图范围
v = f(o1,o2,o3);
h = patch(isosurface(o1,o2,o3,v,0));
isonormals(o1,o2,o3,v,h)
set(h,'FaceColor','r','EdgeColor','none');
xlabel('x');ylabel('y');zlabel('z');
alpha(1)
grid on; view([1,1,1]); axis equal; camlight; lighting gouraud

回复此楼

» 猜你喜欢

投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿已经有22人回复
申请26博士已经有5人回复
职称评审没过，求安慰已经有22人回复
垃圾破二本职称评审标准已经有15人回复
EST投稿状态问题已经有7人回复
毕业后当辅导员了，天天各种学生超烦已经有4人回复
聘U V热熔胶研究人员已经有10人回复
求助文献已经有3人回复
投稿返修后收到这样的回复，还有希望吗已经有8人回复
三无产品还有机会吗已经有6人回复

1楼2016-12-14 15:26:50

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

pdl9527

专家顾问 (小有名气)

专家经验: +8
应助: 100 (初中生)
金币: 2107.8
红花: 23
帖子: 282
在线: 111.9小时
虫号: 1227333
注册: 2011-03-09
性别: GG
专业: 分离过程
管辖: 计算模拟

【答案】应助回帖

感谢参与，应助指数 +1

把程序中对应的地方改成下面的方程，另外检查你的方程表达式和o的取值范围，因为v中存在很多复数和NAN的情况，所以作图有问题。

CODE:

f = @(o1,o2,o3)  (1/(((2*o2-2*o3)/(o1-o3)-1)^2/3 + 1)^(1/2)-((2*o2-2*o3)/(o1-o3)-1)/(2*(((2*o2-2*o3)/(o1-o3)-1)^2/3+1)^(1/2))) ...

    *((o1+o2+o3)^2/2-(3*o1*o2)/2-(3*o1*o3)/2-(3*o2*o3)/2)^(1/2)-5*3^(1/2)+(3^(1/2)*((2*o2-2*o3)/(o1-o3)-1) ...

    *(o1/3+o2/3+o3/3))/(3*(((2*o2-2*o3)/(o1-o3)-1)^2/3+1)^(1/2));        %表达式

[o1,o2,o3] = meshgrid(-10:.2:10,-10:.2:10,-10:.2:10);           % 画图范围

v = arrayfun(f,o1,o2,o3);