m=0.75;L=0.35;k2=300;x2=150;a2=1e-03;N=24;g=9.8;Lj=[0.35,1,0.35];f=0.95;
ze1=zeros(1,8);ze2=zeros(1,4);
syms x;L1=0.35;L2=1;L3=0.35;
y1=det([0,1,0,1,ze1;sin(x.*L1),cos(x.*L1),sinh(x.*L1),cosh(x.*L1),ze1;0,-1,0,1,ze1;cos(x.*L1),-sin(x.*L1),cosh(x.*L1),sinh(x.*L1),-1,0,-1,0,ze2;
-sin(x.*L1),-cos(x.*L1),sinh(x.*L1),cosh(x.*L1),0,1,0,-1,ze2;ze2,0,1,0,1,ze2;ze2,sin(x.*L2),cos(x.*L2),sinh(x.*L2),cosh(x.*L2),ze2;
ze2,cos(x.*L2),-sin(x.*L2),cosh(x.*L2),sinh(x.*L2),-1,0,-1,0;ze2, -sin(x.*L2),-cos(x.*L2),sinh(x.*L2),cosh(x.*L2),0,1,0,-1;
ze1,0,1,0,1;ze1,sin(x.*L3),cos(x.*L3),sinh(x.*L3),cosh(x.*L3);ze1,-sin(x.*L3),-cos(x.*L3),sinh(x.*L3),cosh(x.*L3)]);
fun=@(x) 8*cos((7*x)/10)+8*cosh((7*x)/10)+8*cos(x)*cosh(x)-8*cos((17*x)/10)*cosh(x)-8*cosh((17*x)/10)*cos(x)-8*s
matlab薛之谦,用integral函数求积分出现警告:已达到正在使用的区间最大数目的限制。误差的近似范围为 5.8e-07。积分可能不存在,或者可能很难在数值上接近要......
最新推荐文章于 2024-06-06 17:39:58 发布