声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 829|回复: 2

[综合讨论] 怎么样才能得出想要的结果啊

[复制链接]
发表于 2008-4-11 21:17 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
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 编辑 ]
回复
分享到:

使用道具 举报

发表于 2008-4-12 09:35 | 显示全部楼层

回复 楼主 的帖子

说清楚你的问题
想要得到什么样的结果
以及有什么样的错误提示
发表于 2008-4-14 20:46 | 显示全部楼层

如何求管中物质径向密度分布函数

新人发帖前请先看置顶的会员守则
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-26 10:52 , Processed in 0.061524 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表