程序错了吧,解对的肯定在图上三:
clear;clc;
y00=883.3048 ;
x00=110.0967 ;
bb=54.1413 ;
aa= 31.4013 ;
tt= 0.4041 ;
k0= 1.9626 ;
syms x y
ezplot([(x-x00)*cos(tt)+(y-y00)*sin(tt)]^2/(aa^2)+[(x-x00)*sin(tt)-(y-y00)*cos(tt)]^2/(bb^2)-1,[50 180],[750 1000]);
hold on
ezplot(y-y00+100 - k0*(x-x00+100),[50 180],[750 1000]);
[jx jy] = solve('[(x-x00)*cos(tt) + (y-y00)*sin(tt)]^2/(aa^2) + [(x-x00)*sin(tt) - (y-y00)*cos(tt)]^2/(bb^2) - 1',' y-y00+100 - k0*(x-x00+100)','x,y');
jx=eval(jx);
jy=eval(jy);
scatter(jx,jy)