求前辈指点迷津，感觉椭流线表达式有问题，又不知道哪里不对 K~
VUD(
'+iqbcUd,
程序如下 iYz!:TxP
4fau
9bW
clearclc R |K#nh
h=1000;N=100; A{ Ejk|
R=100;for i=2:101 i]IZ0 .?Y
theata(i)=asin(sqrt((i-1)/N));end +$D~?sk
theata(1)=0;for i=1:101 #6za
% r(i)=R*sqrt((i-1)/N);r(i)=sqrt(R^2/i) oT[8Iu
end% r(1)=0; T0lbMp
a=zeros(101,1);c=zeros(101,1); BZRC0^-C@
x=zeros(101,1);y=zeros(101,1); `AxhA.&V
a(1)=550;x(1)=-50; (Q6}N' T
y(1)=0; 6Ej@;]^^-
Y?V>%eBu
for i=2:100 k=2*a(i-1); av}pT)]\
F=[h,r(i-1)/2-r(i)/2]; XY=((k^2-norm(F)^2)/(2*k-2*norm(F)*cos(pi-theata(i))))*[cos(pi-theata(i)),sin(theata(i))]; gf U!sYZ
x(i)=XY(1); y(i)=XY(2); P*6&0\af|
FX=F-XY; a(i)=(norm(FX)+norm(XY))/2; &g5+ |g (
i x(i) ~
H $q
y(i)end Pt8 U0)i)
plot(x，y）
