杭州锐达数字技术有限公司
查看: 16952|回复: 28

[综合讨论] [求助]曲线求交点?

  [复制链接]
发表于 2006-4-5 19:18 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有帐号?我要加入

x
怎么求两条时间序列曲线的交点???急啊

[ 本帖最后由 eight 于 2007-1-11 14:05 编辑 ]
回复
分享到:

使用道具 举报

发表于 2006-4-6 09:50 | 显示全部楼层
转自《工程计算可视化与MATLAB实现》尚涛等编著 武汉大学出版社

两直线相交

function [X,Y]=pll(X1,Y1,X2,Y2)

% 直线相交求交点

A1=Y1(1)-Y1(2);
B1=X1(2)-X1(1);
C1=Y1(2)*X1(1)-Y1(1)*X1(2);
A2=Y2(1)-Y2(2);
B2=X2(2)-X2(1);
C2=Y2(2)*X2(1)-Y2(1)*X2(2);
D=det([A1,B1;A2,B2]);
X=det([-C1 B1;-C2 B2])/D;
Y=det([A1 -C1;A2,-C2])/D;

调用格式:

x1=[1 5];y1=[1 5];x2=[1 5];y2=[5,1];
[x,y]=pll(x1,y1,x2,y2);
plot(x1,y1,\'r\');
hold on
plot(x2,y2,\'b\');
plot(x,y,\'ko\');



直线与多条直线相交

xi=[1 2 3 4 5];yi=[2 6 3 6 1];
plot(xi,yi);hold on
x1=[1 5];y1=[4 5];line(x1,y1);
x=zeros(size(xi));
y=x;
for i=1:5-1
x2=xi([i i+1]);y2=yi([i i+1]);
[x,y]=pll(x1,y1,x2,y2);
plot(x,y,\'ro\')
end



直线与曲线相交

x=-8:0.1:8;y=x;[X,Y]=meshgrid(x,y);
R=sqrt(X.^2+Y.^2)+eps;Z=sin(R)./R;
contour(Z,3);hold on
c=contour(Z,3);
x=[0 360];y=[0 400];
y=(y(2)-y(1))/(x(2)-x(1))*(x-x(1))+y(1);z=[0 0];
line(x,y,z);c=c\';
X=c(:,1);Y=c(:,2);
r0=abs(Y-(y(2)-y(1))/(x(2)-x(1))*(X-x(1))+y(1))<=.93;
zz=0;yy=r0.*Y;xx=r0.*X;
plot(xx(r0~=0),yy(r0~=0),\'r\')



曲线与曲线相交

x=0:pi/400:2*pi;
x=x\';
y1=sin(pi*x);y2=cos(pi*x);plot(x,y1,x,y2);hold on
r0=abs(y2-sin(pi*x))<=0.02;
yy=r0.*y1;xx=r0.*x;plot(xx(r0~=0),yy(r0~=0),\'r.\')



直线与曲面相交

x=-8:0.3:8;y=x;[X,Y]=meshgrid(x,y);
Z=X.^2+Y.^2;
mesh(X,Y,Z);hold on
x=[-10 10];y=[-10 3];z=[30 35];line(x,y,z);
r0=(abs(Y-y(1)-(y(2)-y(1))/(x(2)-x(1))*(X-x(1)))<=0.45)&...
(abs(Z-z(1)-(z(2)-z(1))/(x(2)-x(1))*(X-x(1)))<0.45)&...
(abs(Y-y(1)-(y(2)-y(1))/(z(2)-z(1))*(Z-z(1)))<=0.45);
zz=r0.*Z;yy=r0.*Y;xx=r0.*X;
plot3(xx(r0~=0),yy(r0~=0),zz(r0~=0),\'r*\');




平面与曲面相交

x=-8:0.1:8;
y=x;
[X,Y]=meshgrid(x,y);
Z1=2*ones(size(X));
Z2=X.^2-Y.^2;
mesh(X,Y,Z1);
hold on
mesh(X,Y,Z2);
r0=(abs(Z1-Z2)<=.65);
zz=r0.*Z1;yy=r0.*Y;xx=r0.*X;
plot3(xx(r0~=0),yy(r0~=0),zz(r0~=0),\'k*\')
clc
disp(\'观察曲面后,按任意键画交线\');
pause
clf
plot3(xx(r0~=0),yy(r0~=0),zz(r0~=0),\'k*\');

%曲面与多个截平面相交
y=-10:0.5:10;
z=y;
[Z,Y]=meshgrid(z,y);
X=Z;
X1=0*ones(size(Z));
X2=3*ones(size(Z));
X3=-3*ones(size(Z));
Z4=(X.^2-Y.^2)/10;
mesh(X1,Y,Z);hold on
mesh(X2,Y,Z)
mesh(X3,Y,Z);
mesh(X,Y,Z4);
r1=(abs(X1-X)<0.05);
r2=(abs(X2-X)<0.05);
r3=(abs(X3-X)<0.05);
zz1=r1.*Z4;yy1=r1.*Y;xx1=r1.*X;
zz2=r2.*Z4;yy2=r1.*Y;xx2=r1.*X;
zz3=r3.*Z4;yy3=r1.*Y;xx3=r1.*X;
plot3(xx1(r1~=0),yy1(r1~=0),zz1(r1~=0),\'k*\');
plot3(xx2(r2~=0),yy2(r2~=0),zz2(r2~=0),\'k*\');
plot3(xx3(r3~=0),yy3(r3~=0),zz3(r3~=0),\'k*\');
colormap(hsv)
clc;
disp(\'观察曲面后,按任意键画交线\');
hold off



平面与曲面相交

y=-8:0.4:8;
z=y;
[Z,Y]=meshgrid(z,y);
X=Z;
X1=zeros(size(Z));
Z2=zeros(size(Z));
Z3=(X.^2-Y.^2)/10;
mesh(X1,Y,Z);
hold on
mesh(X,Y,Z2);
mesh(X,Y,Z3);
r1=(abs(X1-X)<0.05);
r2=(abs(Z3-Z2)<0.05);
r3=(abs(X1-X)<0.05)&(abs(Z-Z2)<=0.05);
zz1=r1.*Z3;yy1=r1.*Y;xx1=r1.*X;
zz2=r2.*Z3;yy2=r2.*Y;xx2=r2.*X;
zz3=r3.*Z;yy1=r3.*Y;xx1=r3.*X1;

plot3(xx1(r1~=0),yy1(r1~=0),zz1(r1~=0),\'k*\');
plot3(xx2(r2~=0),yy2(r2~=0),zz2(r2~=0),\'k*\');
plot3(xx3(r3~=0),yy3(r3~=0),zz3(r3~=0),\'k*\');

colormap(hsv);



曲面与曲面求交

[x,y]=meshgrid(-2:.1:2);
z1=x.^2-2*y.^2;
z2=x.^2+y.^2-5;
mesh(x,y,z1);
hold on
mesh(x,y,z2);
r0=(abs(z1-z2)<=.1);
zz=r0.*z1;yy=r0.*y;xx=r0.*x;
plot3(xx(r0~=0),yy(r0~=0),zz(r0~=0),\'k*\')
colormap(cool)
hold off

[ 本帖最后由 happy 于 2006-10-24 20:00 编辑 ]

点评

赞成: 5.0
赞成: 5
  发表于 2014-8-13 10:02

评分

1

查看全部评分

 楼主| 发表于 2006-4-6 12:42 | 显示全部楼层
非常感谢!仔细研究一下
发表于 2006-7-26 10:40 | 显示全部楼层
fzero
fsolve
发表于 2006-9-14 17:22 | 显示全部楼层
这是怎么回事啊~~
X1=[5 1];Y1=[1 5];
>> plot(X1,Y1,\'r\');
??? plot(X1,Y1,\'r\');
               |
Error: Missing variable or function
发表于 2006-9-14 17:24 | 显示全部楼层
plot(X1,Y1,\'r\');改成plot(X1,Y1,'r');

评分

1

查看全部评分

发表于 2007-9-13 21:00 | 显示全部楼层

回复 #2 happy 的帖子

直线与多条直线相交

xi=[1 2 3 4 5];yi=[2 6 3 6 1];
plot(xi,yi);hold on
x1=[1 5];y1=[4 5];line(x1,y1);
x=zeros(size(xi));
y=x;
for i=1:5-1
x2=xi([i i+1]);y2=yi([i i+1]);
[x,y]=pll(x1,y1,x2,y2);
plot(x,y,\'ro\')
end

请问Pll是什么,MATLAB 无法识别?
发表于 2007-9-13 21:28 | 显示全部楼层
本帖最后由 ChaChing 于 2010-8-25 01:02 编辑

请认真看看二楼的程序
发表于 2007-9-14 12:30 | 显示全部楼层

回复 #13 eight 的帖子

谢谢!我太大意了,竟然没有看到pll是一个自编的函数:@L
发表于 2007-11-11 14:45 | 显示全部楼层
contour可以返回数据点,已知函数也可以solve
如果是 plot 得到的不规则曲线又该怎么做呢?
发表于 2010-8-24 22:10 | 显示全部楼层
真是好东西,辛苦楼主了!
发表于 2010-9-6 11:26 | 显示全部楼层
本帖最后由 xiaoxiaobeifeng 于 2010-9-6 11:29 编辑

两条曲线求交点吧 我觉得可以根据定义来嘛
abs(x(i)-y(i))<=abs(x(i-1)-y(i-1)) $$ abs(x(i)-y(i))<=abs(x(i+1)-y(i+1)) 应该可以吧
发表于 2010-9-7 11:24 | 显示全部楼层
回复 happy 的帖子
不错,谢谢


   
发表于 2010-10-17 21:33 | 显示全部楼层
找到答案了!
发表于 2010-10-18 05:37 | 显示全部楼层
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

快速回复 返回顶部 返回列表