Matlab中，如何获得这个曲面与x=1，v3，e=300平面的交线并显示出来？

Matlab中，如何获得这个曲面与x=1，v=3，e=300平面的交线并显示出来？~

1、题主所给 x 的范围比较奇怪，0.173-0.987 这个区间并不包括x=1，怎么和这个平面相交？另，x和e取的过于密集，这样会导致资源消耗很大（内存、CPU、图形各方面），但对最终效果没什么明显帮助。

2、基本思路：以x=1为例，可以让v在给定取值范围内变化，然后求解方程 f(x,v,e)=0 得到e。解方程可以用 fzero，初值可以取 e 的均值。

参考代码：
f = @(x,v,e)50.111.*v.^2.*(3.55-x).^-2+2.159.*v.^2-e;X = 0.1:.01:1.5;V = 0:.4:8;E = 0:10:700;[x,v,e] = meshgrid(X,V,E);f1 = f(x,v,e);clfh = patch(isosurface(x,v,e,f1,0)); isonormals(x,v,e,f1,h) set(h,'FaceColor','r','EdgeColor','none');xlabel('x');ylabel('v');zlabel('e'); alpha(0.7) grid on; view([1,1,1]); axis; camlight; lighting gouraud% 与 x=1 平面交线X1 = V*0 + 1;E1 = V*0;for ii = 1:length(V) E1(ii) = fzero(@(e)f(1,V(ii),e), mean(E));endhold allh(2) = plot3(X1, V, E1);% 与 v=3 平面交线V1 = X*0 + 3;E1 = X*0;for ii = 1:length(X) E1(ii) = fzero(@(e)f(X(ii),3,e), mean(E));endh(3) = plot3(X, V1, E1);% 与 e=300 平面交线V1 = X*0;E1 = X*0 + 300;for ii = 1:length(X) V1(ii) = fzero(@(v)f(X(ii),v,300), mean(V));endh(4) = plot3(X, V1, E1);rotate3dview(55,30)set(h(2:4), 'linewidth',2)legend(h(2:4),'x=1','v=3','e=300')

[X,Y,Z]=sphere(50);
Z(Z<0)=-(Z(Z<0));
mesh(X,Y,Z)%画上半球面
hold on;
ezmesh('1.5-x-y',[-1 1])%画平面
[x,y,z]=meshgrid(linspace(-1,1));
contourslice(x,y,z,x+y+z-1.5,X,Y,Z,[0 0])%画交线
axis equal
view(135,20)

r0=find(abs(z1-z2)<=0.1); % 找到对应z1与z2差值小于0.1的序号，即近似z1=z2的位置，也就是找到交点
px=x(r0); % 取出交点的横坐标
py=y(r0); % 取出交点的纵坐标
pz=z1(r0); % 取出交点的竖坐标
plot3(px,py,pz,'*') % 根据取出的交点坐标画图，样式为*表示的点
colormap(gray) % 采用灰度色标，也就是你看到的曲面图是以灰色深浅表示z坐标大小的

---

洪山区17223623401： matlab中,已知大量点的三维坐标x,y,z ,如何绘制出这个曲面,要网格的曲面,希望写出具体代码 谢谢 - ？
苗志小儿： 看这些点是不是很规矩的网格点的数据. 如果是网格点数据,适当做点重组,就可以用surf函数得到曲面图. 如果不是网格点数据,那么请用griddata函数插值到网格点数据,然后再surf 如: data=[ 1 1 5.14 1 2 9.56 1 3 8.66 2 1 12.55 2 2 10.51 ...

洪山区17223623401： matlab如何将离散点拟合成曲面 - ？
苗志小儿： % 方法1 T = [1:100];D = [1:100]; K = rand(1,100);% 插值 [X,Y,Z]=griddata(T,D,K,linspace(min(T),max(T))',linspace(min(D),max(D)),'v4'); figure,surf(X,Y,Z);%方法2 x = [1:100];y = [1:100]; z = rand(1,100); xlin = linspace(min(x),max(x),33); ylin = ...

洪山区17223623401： Matlab已知三维坐标,怎么得到相应的曲面图 - ？
苗志小儿： 可以通过matlab提供的view函数来获取三维坐标图的视角信息.其调用格式为:[az,el]=view;%返回当前的方位角和仰角az是azimuth(方位角)的缩写,EL是elevation(仰角)的缩写.它们均以度为单位.系统缺省的视点定义为方位角-37.5°,...

洪山区17223623401： Matlab 中如何得到用isosurface 形成的曲面的点集？
苗志小儿： 举个例子 [x y z v] = flow; [fc vt] = isosurface(x, y, z, v, -3);%vt就是曲面点集 p=patch('Faces',fc,'Vertices',vt); isonormals(x,y,z,v, p) set(p, 'FaceColor', 'red', 'EdgeColor', 'none'); daspect([1 1 1]) view(3) camlight; lighting phong figure plot3(vt(:,1),vt(:,2),vt(:,3),'.')

洪山区17223623401： 已知x,y,z的值,怎么在matlab中绘出曲面图 - ？
苗志小儿： 如果 x,y,z 是规格相同的矩阵,可以用 surf(x,y,z) 但如果 x,y,z 是行向量或列向量时,用 scatter3(x,y,z)

洪山区17223623401： 已知x,y,z的值,怎么在matlab中绘出曲面图?我有x,y,z的值,想在matlab中绘出三维曲面图?例如x=[1 2 3 4 5 6];y=[1 2 3 4 5 6];z=[23 43 51 14 16 35];[X,Y]=... - ？
苗志小儿：[答案] 这样试试. z=[23 43 51 14 16 15]; [X,Y]=meshgrid([1 2 3 4 5 6]); Z=meshgrid(z); surf(X,Y,Z)

洪山区17223623401： MATLAB中如果已经有X、Y、Z的一大堆数值,想要输入数值,得到曲线,并得到表达式,应该怎么做. - ？
苗志小儿： plot(x,y,z)%就能画图了(x,y,z三个数组,元素个数相同) 光有数据是不能凭空得到表达式滴~ 你只能指定一个固定的表达式形式,比如z=ax+by,然后根据你的数据用matlab找到一组最适合的a,b 这样你的表达式就有了,这个过程就是所谓的“拟合”. 具体怎么做你百度下曲线拟合,看两个简单的例子吧

洪山区17223623401： 已知一组三维散乱点,怎么用matlab拟合出曲面,并得到每个点的法线 - ？
苗志小儿： 嘿嘿,这是个数学问题,每点的法线的求取没那么麻烦,不求切平面,使用surfnorm可求取每一点法向在三维空间xyz三方向的分量,用这三个分量可以使用quiver3成图,也可直接使用surfnorm成图.在matlab的帮助文档中有很详细的例程,你查查吧.

洪山区17223623401： matlab xyz数据如何进行曲面拟合获得曲面方程 - ？
苗志小儿： 用regress函数吧. REGRESS(y,X),其中,y就是你的Z,X就是你的X和Y. REGRESS Multiple linear regression using least squares. b = returns the vector of regression coefficients, b, in the linear model y = Xb, (X is an nxp matrix, y is the nx1 vector of observations).

洪山区17223623401： 已知三维数据点,用matlab求拟合的曲面及其方程表达 - ？
苗志小儿： 可以根据你熟知的曲面方程,去进行拟合,如拟合后的相关系数R^2接近于0.9以上,就可以基本确定该曲面方程了.