[求助]求两曲线的交点坐标,帮忙看看啊

wilddrog

t1=input('input t1:');
t2=input('input t2:');
t3=input('input t3:');
t4=input('input t4:');
B='black';
R='red';
G='green';
c=30;
v=34;
d=c*(t2-t1);
K=v*(t4-t3);
f=25;
if t1<t2;
x=-5*pi:0.1:0;
elseif t1==t2
x=0;
else t3>t4
x=0:0.1:5*pi;
end
w=((16*f^2-4*d^2)*x.^2-4*f^2*d^2+d^4)/4*d^2; %f=5;
y1=sqrt(w);
y2=-sqrt(w);
subplot(2,2,1);
plot(x,y1,x,y2,'linewidth',2,'color',R);
hold on;

K=v*(t4-t3);
%y=-5*pi:pi/50:5*pi;
F=10;

%m=((16*F^2-4*K^2)*y.^2-4*F^2*K^2+K^4)/4*K^2;
m=(4*K^2*x.^2-K^4+4*F^2*K^2)/16*F^2-4*K^2
y3=sqrt(m);
y4=-sqrt(m);
subplot(2,2,2);
plot(x,y3,-x,y3,'linewidth',2,'color',B);
hold off;
subplot(2,2,3);
plot(x,y1,x,y2,x,y3,-x,y3,'linewidth',2,'color',G);

x='x\';
y1=sqrt(w);
y3=sqrt(m);
plot(x,y1,x,y3);
hold on;
r0=abs(y3-y1)<=0.02;
yy=r0.*y1;xx=r0.*x;
plot(xx(r0~=0),yy(r0~=0),'r.\')


===============eight==================
更多细节请参阅以下帖子：
[求助]曲线求交点？
====================================

[此贴子已经被cdwxg于2006-6-3 11:47:29编辑过]

[ 本帖最后由 eight 于 2007-1-11 14:06 编辑 ]

cdwxg

http://forum.vibunion.com/thread-15709-1-1.html看下该帖子。