运行不出来啊 郁闷 。 这是等误差查补程序 高手过目啊
x1=input('输入起点横坐标 X\n x1 = ');y1=input('输入起点纵坐标 Y\n y1 = ');
x2=input('输入终点横坐标 X\n x2 = ');
y2=input('输入起点纵坐标 Y\n y2 = ');
xa=x1;
ya=y1;
x=x1;
y=y1;
xi=x1;
yi=y1;
while(xa<=x2)
yp =ya-0.04;
while(yp-ya<=0.04)
m=(yp- ya)*(yp- ya)
xp=xa-sqrt(0.04*0.04-m);
k=sqrt(0.04*0.04/m -1);
xt=(k-1)/2;
yt =xt*xt+xt+1;
k1=(yt-yp)/(xt-xp);
if abs(k-k1)<=0.01
n=(1-k)*(1-k);
u=4*(1+k*xa-ya);
t=sqrt(n-u);
xi=k/2-0.5+t/2;
yi=-k*xa+k*xi+ya;
xa=xi;
ya=yi;
xi=;
yi=;
yp = yp +0.01;
else
yp = yp +0.01;
end
end
end
plot(xi,yi);
grid;
不懂，高手快来。。。 我表示看不懂你的程序！能不能把你的原题拿出来看下！
页:
[1]