Question

matlab怎么解四元二次方程组

x1=x3*cos(t)+y3*sin(t);
y1=y3*cos(t)+x3*sin(t);
x2=x4*cos(t)+y4*sin(t);
y2=y4*cos(t)+x4*sin(t);
tan(t)=(y4-y3)/(x4-x3);
其中，x3,x4,y3,y4是未知数，x1,x2,y1,y2已知；
x1=19/64;y1=0;x2=1;y2=57/64;
怎么求出x3,x4,y3,y4的值
斜椭圆C：((x-p)*cos(t)+(y-q)*sin(t)).^2/a^2+((y-q)*cos(t)-(x-p)*sin(t)).^2/d^2=1；
p=(0+x4)/2,
q=(y3+y4)/2;
t=atan((y4-y3)/(x4-0))；
a=0.5*sqrt((0-x4)^2+(y3-y4)^2)；
x1=19/64;y1=0;x2=1;y2=57/64;在椭圆上
怎么求y3,x4,y4的值，x3=0；(x3,y3),(x4,y4)是斜椭圆两个端点
我是想求这个，帮一下忙吧

Answer 1

用solve（）函数求解四元一次方程组。代码如下
syms x3 x4 y3 y4 t
x1=19/64;y1=0;x2=1;y2=57/64;
eq1=x3*cos(t)+y3*sin(t)-x1;
eq2=y3*cos(t)+x3*sin(t)-y1;
eq3=x4*cos(t)+y4*sin(t)-x2;
eq4=y4*cos(t)+x4*sin(t)-y2;
eq5=(y4-y3)/(x4-x3)-tan(t);
s=solve(eq1,eq2,eq3,eq4,eq5,'x3','x4','y3','y4','t')
double(s.t)
double(s.x3)
double(s.x4)
double(s.y3)
double(s.y4)

追问

斜椭圆C：((x-p)*cos(t)+(y-q)*sin(t)).^2/a^2+((y-q)*cos(t)-(x-p)*sin(t)).^2/d^2=1；p=(0+x4)/2,q=(y3+y4)/2;
t=atan((y4-y3)/(x4-0))；a=0.5*sqrt((0-x4)^2+(y3-y4)^2)；
x1=19/64;y1=0;x2=1;y2=57/64;在椭圆上；