怎么样才能得出想要的结果啊

skyriver2008

l0=177.8;R=175;l=482.6;%xitamax=40;%xc=R*l/xitamax
xc0=1000;yc0=200;xg0=900;yg0=500;M=1100;g=9.8;
xd=800;yd=150;ad=900;ac=850;ab=400;ye=-150;xe=900;
for xita=40*pi/180%0:5*pi/180:35*pi/180
    xc=xc0*cos(xita)-yc0*sin(xita);
    yc=xc0*sin(xita)+yc0*cos(xita);
    xg=xg0*cos(xita)-yg0*sin(xita);
    yg=xg0*sin(xita)+yg0*cos(xita);
    syms xa ya;
    [xa,ya]=solve('(xa-xd)^2+(ya-yd)^2-ad^2=0','(xa-xc)^2+(xa-yc)^2-ac^2=0');
     %[eval('(xa-xd).^2+(ya-yd).^2-ad.^2'),eval('(xa-xc).^2+(xa-yc).^2-ac.^2')]
    syms xb yb;
    [xb yb]=solve('(xb-xa)^2+(yb-ya)^2-ab^2=0','(xa-xc)^2+(ya-yc)^2-ac^2=0');
    syms xf yf;
    [xf,yf]=solve('(yd-ya)*xf+(xa-xd)*yf+ya*(xd-xa)-(yd-ya)*xa=0',...
    '(ye-yb)*xf+yf*(xb-xe)+yb*(xe-xb)-xb*(ye-yb)=0 ');
    Dofc=abs(yf*(xc-xf)-xf*(yc-yf))/sqrt((yf-yc)^2+(xf-xc)^2);
    Wxg=M*g*yg;
    Ffc=Wxg/Dofc
    Dabe=abs(xa.*(ye-yb)+ya.*(xb-xe)+yb.*(xe-xb)-xb.*(ye-yb))./sqrt((yb-ye).^2+(xb-xe).^2);
    Dafc=abs(xa.*(yc-yf)+ya.*(xf-xc)+ye.*(xc-xf)-xf.*(yc-yf))./sqrt((yf-ya).^2+(xf-xc).^2);
    Fbe=Ffc*Dafc/Dabe;
end
以上是我的程序,我想算出Fbe和Ffc这两个力来,但是就没得我想要的结果,这是怎么回事啊?
高手来帮帮忙啊!

[ 本帖最后由 eight 于 2008-4-14 20:45 编辑 ]

sigma665

说清楚你的问题
想要得到什么样的结果
以及有什么样的错误提示

eight

新人发帖前请先看置顶的会员守则