怎么样才能得出想要的结果啊
l0=177.8;R=175;l=482.6;%xitamax=40;%xc=R*l/xitamaxxc0=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;
=solve('(xa-xd)^2+(ya-yd)^2-ad^2=0','(xa-xc)^2+(xa-yc)^2-ac^2=0');
%
syms xb yb;
=solve('(xb-xa)^2+(yb-ya)^2-ab^2=0','(xa-xc)^2+(ya-yc)^2-ac^2=0');
syms 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 编辑 ]
回复 楼主 的帖子
说清楚你的问题想要得到什么样的结果
以及有什么样的错误提示
如何求管中物质径向密度分布函数
新人发帖前请先看置顶的会员守则
页:
[1]