数学建模社区-数学中国

标题: 运行不出来啊郁闷。这是等误差查补程序高手过目啊 [打印本页]

作者: 梦想成真 时间: 2011-2-27 10:24
标题: 运行不出来啊郁闷。这是等误差查补程序高手过目啊
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=[x1,xi];
yi=[y1,yi];
yp = yp +0.01;
else
yp = yp +0.01;
end
end
end
plot(xi,yi);
grid;

作者: gaoshanliu水 时间: 2011-2-27 10:28
不懂，高手快来。。。

作者: 数模鸿雁 时间: 2011-2-27 14:05
我表示看不懂你的程序！能不能把你的原题拿出来看下！

欢迎光临数学建模社区-数学中国 (http://www.madio.net/)