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

回复此楼

» 猜你喜欢

26/27申博自荐已经有10人回复
东北林业大学材料科学与工程学院“一流”A+学科国家级人才团队课题组招收2026级博士生已经有3人回复
医学类期刊求推荐已经有5人回复
生活琐事由它去已经有4人回复
提交了我也来说说感想已经有12人回复
青B发送上会通知了吗已经有9人回复
西安交大新媒学院副院长用撤稿论文结题已经有6人回复
论文撤稿了已经有8人回复
化学专业申博已经有4人回复
某211大学教师把个人教师官方主页改成:我跑了我跑了我跑了!官宣跑路！已经有5人回复

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

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

pdl9527

专家顾问 (小有名气)

专家经验: +8
应助: 100 (初中生)
金币: 2110.3
红花: 23
帖子: 282
在线: 112小时
虫号: 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);