24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

kkx1986

金虫 (正式写手)

应助: 5 (幼儿园)
金币: 1484.3
散金: 253
红花: 22
帖子: 937
在线: 272.1小时
虫号: 896506
注册: 2009-11-07
性别: GG
专业: 机械动力学

[求助] 应用matlab求两个曲面的交线

%第一曲面三个顶点的坐标
%第一点轨迹

x11=A*sin(w*t);
y11=vw*t+(ds/2)*sin(ws*t);
z11=0.5*ds*(1-cos(ws*t));
%plot3(x11,y11,z11)
%第二点轨迹
x12=x11+H*tan(theta);
y12=y11;
z12=z11+H;
%第三点轨迹
x13=x11-H*tan(theta);
y13=y11;
z13=z11+H;

%第二个曲面三个顶点的坐标
%第一点轨迹
x21=A*sin(w*(t+5*L/vs));
y21=vw*(t+5*L/vs)+(ds/2)*sin(ws*t);
z21=0.5*ds*(1-cos(ws*t));
%hold on
%plot3(x21,y21,z21,'r')
%第二点轨迹
x22=x21+H*tan(theta);
y22=y21;
z22=z21+H;
%第三点轨迹
x23=x21-H*tan(theta);
y23=y21;
z23=z21+H;

X1=[x12; x11; x13];
Y1=[y12; y11; y13];
Z1=[z12; z11; z13];

X2=[x22; x21; x23];
Y2=[y22; y21; y23];
Z2=[z22; z21; z23];
mesh(X1,Y1,Z1);
hold on
mesh(X2,Y2,Z2);

怎么求这两个曲面的交线，我在网上搜了一下，都是要求方程是z=f（x，y）这种显示表示的。

回复此楼

» 收录本帖的淘帖专辑推荐

考研

» 猜你喜欢

上海工程技术大学【激光智能制造】课题组招收硕士已经有6人回复
带资进组求博导收留已经有11人回复
自荐读博已经有5人回复
求个博导看看已经有16人回复
上海工程技术大学张培磊教授团队招收博士生已经有4人回复
求助院士们，这个如何合成呀已经有4人回复
临港实验室与上科大联培博士招生1名已经有9人回复
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？已经有6人回复
最近几年招的学生写论文不引自己组发的文章已经有11人回复
中科院杭州医学所招收博士生一名（生物分析化学、药物递送）已经有3人回复

» 本主题相关价值贴推荐，对您同样有帮助:

千头万绪净心第一不净心全是废话，二六时中精进为先不精进一事无成。

1楼 2013-03-21 09:58:18

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

kkx1986

金虫 (正式写手)

应助: 5 (幼儿园)
金币: 1484.3
散金: 253
红花: 22
帖子: 937
在线: 272.1小时
虫号: 896506
注册: 2009-11-07
性别: GG
专业: 机械动力学

引用回帖:

6楼: Originally posted by leedobb at 2013-03-22 14:15:24
只有一个参数坐标，只能是曲线，而不可能是曲面。...

clear all; clc; close all
ds=0.05;
vw=-0.6702;
vs=18.3;
a=0.001;
f=20000;
h=0.001;
A=0.000016;
w=2*pi*f;
ws=2*vs/ds;
angle=acos(1-2*h/ds);
L=0.001;
H=0.0001;
theta=pi/9;
T=angle/ws;
t=-T:0.000001:T;
%第一个曲面三个顶点轨迹
%第一点轨迹

x11=A*sin(w*t);
y11=vw*t+(ds/2)*sin(ws*t);
z11=0.5*ds*(1-cos(ws*t));
%plot3(x11,y11,z11)
%第二点轨迹
x12=x11+H*tan(theta);
y12=y11;
z12=z11+H;
%第三点轨迹
x13=x11-H*tan(theta);
y13=y11;
z13=z11+H;

%第二个三个顶点轨迹
%第一点轨迹
x21=A*sin(w*(t+5*L/vs));
y21=vw*(t+5*L/vs)+(ds/2)*sin(ws*t);
z21=0.5*ds*(1-cos(ws*t));
%hold on
%plot3(x21,y21,z21,'r')
%第二点轨迹
x22=x21+H*tan(theta);
y22=y21;
z22=z21+H;
%第三点轨迹
x23=x21-H*tan(theta);
y23=y21;
z23=z21+H;
%点连成线
%X=[x11; x12; x13; x11];
%Y=[y11; y12; y13; y11];
%Z=[z11; z12; z13; z11];
%plot3(X,Y,Z);
%点连成面

X1=[x11; x12; x13; x11];
Y1=[y11; y12; y13; y11];
Z1=[z11; z12; z13; z11];

X2=[x21; x22; x23; x21];
Y2=[y21; y22; y23; y21];
Z2=[z21; z22; z23; z21];
C1=2*ones(size(Z1));
mesh(X1,Y1,Z1,C1)

hold on
C2=4*ones(size(Z2));
patch(X2,Y2,Z2,C2)

grid
axis equal

xlabel('x(t)');
ylabel('y(m)');
zlabel('z(m)');