永发信息网

matlab中怎么把f=@(x,y,z)x.^2+y.^2+4*z.^4+(4*y-4).*z.^2这个方程化成立体的形式,类似下图

答案:2  悬赏:0  手机版
解决时间 2021-01-11 08:44
matlab中怎么把f=@(x,y,z)x.^2+y.^2+4*z.^4+(4*y-4).*z.^2这个方程化成立体的形式,类似下图
最佳答案
你这个是个三元函数,要是画图就是四维的了
按你的意思,你说要画满足f(x,y,z)=0方程的曲面吧
那么比较麻烦,先要解出方程z=fz(x,y),再根据fz画图
由于是四次方程,所以有四个解,还要考虑在实数范围
根据以上种种,写出了程序


f=@(x,y,z)x.^2+y.^2+4*z.^4+(4*y-4).*z.^2;
%从f(x,y,z)=0 解出来的z(x,y),有四组解分别是z1(x,y) z2(x,y) -z1(x,y) -z2(x,y)
z1=@(x,y) sqrt((1-sqrt(- x.^2 - 2*y + 1) - y)/2);
z2=@(x,y) sqrt((1+sqrt(- x.^2 - 2*y + 1) - y)/2);


[xx,yy]=meshgrid(-3:0.1:3,-3:0.1:0.6);%x范围-1~1 y范围-1.5~0.5(实数范围内y不会超过0.5)
zz1=z1(xx,yy);
zz1(abs(imag(zz1))>1e-6)=nan;%考虑在实数范围内,去除虚数部分
zz1=real(zz1);
zz2=z2(xx,yy);
zz2(abs(imag(zz2))>1e-6)=nan;
zz2=real(zz2);
zz3=-zz1;
zz4=-zz2;
mesh(xx,yy,zz1);hold on;
mesh(xx,yy,zz2);
mesh(xx,yy,zz3);
mesh(xx,yy,zz4);hold off;
xlabel('x');
ylabel('y');
zlabel('z');




得到的图像如下



像是两个筒其中有一点连起来
由于在边界处,计算的点出现虚数被删除,曲面边缘出现有点不连续
追问你好厉害,但是我不是想画这个图,这个图我能画出来,我只是想在它上面立体感强,如我发的图上有等势线,然后是立体的感觉,好的话我再追加追答不是太明白你的意思
全部回答
不会这个呢
我要举报
如以上问答信息为低俗、色情、不良、暴力、侵权、涉及违法等信息,可以点下面链接进行举报!
大家都在看
阴历93.11.30 阳历94.1.11 19:40 到底是属什
胭脂的主要成份是什么
西服190/112c是什么意思
手被特别热的铁盆烫伤了,很疼,有什么好办
6T管理是什么?
育雷小学怎么样
改型地带(石桥铺店)地址在什么地方,我要处理
跪求陈道明版《女巡按》!
求教如何定制一双AJ鞋
如何出塑料零件图纸
以阅读为对象,就莎士比亚这一段仿写,10悬赏
四柱分别,丙申,辛卯,己亥,甲子,五行欠缺什么
1060钢含碳量,它和60号碳钢!那个硬度更大?
买股票要下载什么软件吗
这是欧派的吗?我怕自己买到假的。5000买这三
推荐资讯
农业银行etc可以在超市消费吗
男孩名字大全2015马文什么
洛奇英雄传杀不杀盖联
刚考完驾照的新司机应该如何买车
光明大陆ios和安卓能一起玩吗 互通吗
windows10系统windows defender怎么打开?
在安博教育培训,主要是培训什么课程
西游记86版开始是丢丢丢 后面是什么
一个正方体的棱长总和是36cm,它的表面积是__
在文海投资被骗了,有人根我一样吗
在男人心中挣钱和陪老婆哪个重要
我爸出轨之后,现在更离谱了,我该怎么办。
正方形一边上任一点到这个正方形两条对角线的
阴历怎么看 ?