24小时热门版块排行榜    

查看: 526  |  回复: 2

竢实扬华臻

铁虫 (小有名气)

[求助] Matlab求解连续梁方程组

function b3monment(q,l)
syms F1 F2 z l1 l2
l12=l1+l2;l3=l-l1-l2;
a=q/l*(l1*int('1',0,l1)-int('x',0,l1));            %反弯矩M1
b=q/l*(int('x^2',l12,l)-(l1+l2)*int('x',l12,l));%反弯矩M2
c=q/l*int('x',0,l);                                   %合力
d=q/l*int('x',0,z);                                  %跨中极大值左侧剪力
e=q/l*int('x',z,l);                                      %跨中极大值右侧剪力——确定极值点
f=q/l*int('x^2',0,z);                                 %正弯矩Mmax
g=q/l*int('x^2',0,l12);                             %反正前面是定义函数表达式。
%求解单跨连续梁条件:两个反弯矩相等,以及等于正弯矩,分布是三角分布,q,l需要指定值



[F1,F2,l1,l2,x]=solve('F1+F2=c','a+b=0','d-F1-e+F2=0','a+f=0','g-F1*l2-b=0',...
    'F1','F2','l1','l2','z');

可是最后求解时出先无解析解,无解析解释什么意思呢,是无解还是方程有误?
求高手
运行:
> > bmonment(1,1)
Warning: Explicit solution could not be found.
> In solve at 140
  In bmonment at 11
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

一种无序的文化繁荣和混乱的思想活跃,带给人们的不是知识和智慧,而是精神的疲惫和理性的迷茫。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

梁学武lxw

木虫 (小有名气)

【答案】应助回帖


感谢参与,应助指数 +1
xiegangmai: 金币+1, 谢谢参与! 2013-05-14 21:30:03
求解的时候应该在命令窗口输入的是:b3monment(1,1)吧?
2楼2013-05-14 09:19:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

竢实扬华臻

铁虫 (小有名气)

引用回帖:
2楼: Originally posted by 梁学武lxw at 2013-05-14 09:19:52
求解的时候应该在命令窗口输入的是:b3monment(1,1)吧?

对的,就是,但是总是说无解析解
一种无序的文化繁荣和混乱的思想活跃,带给人们的不是知识和智慧,而是精神的疲惫和理性的迷茫。
3楼2013-05-14 10:06:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 竢实扬华臻 的主题更新
信息提示
请填处理意见