| 查看: 3324 | 回复: 7 | |||
kkx1986金虫 (正式写手)
|
[求助]
应用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）这种显示表示的。 |
回复此楼» 收录本帖的淘帖专辑推荐
 考研 |
» 猜你喜欢
今天务委会开完了，明天出结果吗
已经有21人回复
-
2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？
已经有14人回复
-
在坚冰还盖着北海的时候，我看到了怒放的梅花。
已经有9人回复
-
明天应该可查了！？
已经有6人回复
-
有没有大神帮我看看基金代码
已经有10人回复
-
放榜前的不淡定
已经有19人回复
-
某些机构，以效率低为荣，以效率低作为存在感
已经有10人回复
-
2026国自然函评费到账
已经有23人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有22人回复
-
什么时候开奖？
已经有13人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于Matlab画三维曲面图的问题，求助高手，感激感激
已经有3人回复
- 求助：matlab3D曲面方程拟合（急）
已经有2人回复
- 【求助】matlab曲面拟合表达式
已经有10人回复
- 【求助】在MATLAB中如何将离散点拟合成曲面图
已经有18人回复
715211229
木虫 (正式写手)
- 应助: 52 (初中生)
- 金币: 3652.8
- 散金: 20
- 红花: 6
- 帖子: 339
- 在线: 249.2小时
- 虫号: 1357002
- 注册: 2011-07-27
- 性别: GG
- 专业: 传动机械学
2楼2013-03-22 07:52:31
kkx1986
金虫 (正式写手)
- 应助: 5 (幼儿园)
- 金币: 1484.3
- 散金: 253
- 红花: 22
- 帖子: 937
- 在线: 272.1小时
- 虫号: 896506
- 注册: 2009-11-07
- 性别: GG
- 专业: 机械动力学
3楼2013-03-22 10:31:33
4楼2013-03-22 13:16:16
kkx1986
金虫 (正式写手)
- 应助: 5 (幼儿园)
- 金币: 1484.3
- 散金: 253
- 红花: 22
- 帖子: 937
- 在线: 272.1小时
- 虫号: 896506
- 注册: 2009-11-07
- 性别: GG
- 专业: 机械动力学
5楼2013-03-22 13:24:32
6楼2013-03-22 14:15:24
wgdxidname
木虫 (著名写手)
- 应助: 32 (小学生)
- 金币: 3350.1
- 散金: 878
- 红花: 9
- 帖子: 1375
- 在线: 284.2小时
- 虫号: 476192
- 注册: 2007-12-09
- 性别: GG
- 专业: 计算数学与科学工程计算
7楼2013-03-22 14:28:57
kkx1986
金虫 (正式写手)
- 应助: 5 (幼儿园)
- 金币: 1484.3
- 散金: 253
- 红花: 22
- 帖子: 937
- 在线: 272.1小时
- 虫号: 896506
- 注册: 2009-11-07
- 性别: GG
- 专业: 机械动力学
|
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)'); |
8楼2013-03-23 11:11:12
28