求前辈指点迷津，感觉椭流线表达式有问题，又不知道哪里不对 J W"`i
' <'5BeU
程序如下 ?xuWha@:
h-x~:$Z,
clearclc , eSpt #M
h=1000;N=100; -j1]H"-
R=100;for i=2:101 UzW]kY[A<
theata(i)=asin(sqrt((i-1)/N));end qRz /$|.
theata(1)=0;for i=1:101 A\v53AT
% r(i)=R*sqrt((i-1)/N);r(i)=sqrt(R^2/i) olKM0K
end% r(1)=0; /m i&7C(6
a=zeros(101,1);c=zeros(101,1); PEaZ3{-
x=zeros(101,1);y=zeros(101,1); OzR<jCOS
a(1)=550;x(1)=-50; 2*]
[M,L0c
y(1)=0; /FjdcH=
\jZ)r>US"
for i=2:100 k=2*a(i-1); KOoV'YSC[(
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))]; o_5[}d
x(i)=XY(1); y(i)=XY(2); = J]M#6N0
FX=F-XY; a(i)=(norm(FX)+norm(XY))/2; Z$%!H7w
i x(i) /%)(Uz
y(i)end 1H -~+lf
plot(x，y）
