求助求助，这个如何实现？用什么代码？

一个曲面，被一个垂直于xOy面的平面所截，如何绘制截面处的曲线？平面和xOy面的投影上的两个点就是曲线的定义区间的左右端点。我语言表述不清楚，大家看看图吧，顺便一提，曲面没方程，是通过数值方法得到的偏微分方程的解。私以为如果AB投影是垂直或平行于坐标轴的直线还好，毕竟这时平面的表达式就是x=a或y=a，问题在于如果平面的形式是ax+by=c这种情况。

我的背景其实就类似于这样，我说的“曲面”其实就是海底地形的高程，“截曲面的平面”其实就是这张图里的M1、M2、M3。

不感兴趣

开通SVIP免广告

大家当得到我说的曲面上的点的方式是命令“u=A\b”，就当A和b都是已知的矩阵和向量，然后u中的各个元素用v4等插值法拟合成曲面，如何实现这个的命令我是知道的。

你可以查一下meshgrid，这个命令，它可以生成两个坐标的矩阵，分别对应着x,y的数值，然后再从这个矩阵中找出满足ax+by=c这样的所有坐标，然后记得取整，之后就可以通过这些去找到对应的数值了，就可以画出你后面要的线了。如果你前面的二维图片精度不够或者大小不够大，建议你用插值试试。

griddata就可以，要插值到的点坐标(xi,yi)就是平面和xoy交线

给你个例子：
clc; clear;
n = 201;
x = linspace(-2, 2, n);
y = x;
[X, Y] = meshgrid(x, y);
z = sin(X.^2+Y.^3);
x1 = -2; y1 = -2; x2 = 2; y2 = 2;
n = 201;
xx = linspace(x1, x2, n);
k = (y2-y1) / (x2-x1);
b = y1 - k * x1;
yy = k * xx + b;
v = zeros(1, n);
d0 = 0.05;
for i = 1: n
xi = xx(i); yi = yy(i);
d = sqrt((xi-X).^2+(yi-Y).^2);
vv = z(d<=d0);
v(i) = mean(vv);
end
figure('position', [100, 500,1000, 350]);
subplot(121);
colormap(jet);
[~, c] = contourf(x, y, z, 51);
hold on;
plot(xx, yy, 'linewidth', 2, 'color', 'w');
set(c, 'color', 'none');
colorbar();
xlabel('x'); ylabel('y');
subplot(122);
plot(xx, v, 'b');
xlabel('x');
ylabel('v');

扫二维码下载贴吧客户端

下载贴吧APP
看高清直播、视频！

6回复贴，共1页

<<返回matlab吧

分享到:

日	一	二	三	四	五	六