Question

我是matlab小白，急求大佬解答啊，怎么求解这个方程啊

N=250000;k=1416.59;k0=1*(10^7);EI=3.518*(10^7);D=0.32;M1=247.79;M3=164.95;W1=804.7;W2=1616.51;p=184.32;
N1=N+0.5*W1*L*sin(25*pi/180);E=(N1/(2*EI)+((N1*0.5/EI)^2+k/EI))^0.5;F=(((N1*0.5/EI)^2+k/EI)^0.5-N*0.5/EI)^0.5;
c=N/EI;d=(W2*sin(25*pi/180)-p);e=(5.717*M3-k0*D)/EI;q=((k0*D-5.717*M3)/(4*EI))^0.25;
a9=d*p-e-4*p^4;a10=d*p-2*c*p^2;a11=a9^2-p*a10;a12=a9*a10+(a9^2-a10^2)/2;a13=a10^2+a9*a10;a14=a10^2-a9*a10;a15=(a9^2-a10^2)/2-a9*a10;a16=a9^2+a9*a10;a17=(2*c^2+d^2)/(4*c*(c^2+d^2));a18=d*0.5/(c^2+d^2);a19=d^2/(4*c*(c^2+d^2));a26=a11*a17+a12*a18+a13*a19;a27=a14*a17+a15*a18+a16*a19;
a28=4*E*F*q^4*a26;a29=2*E*F^3*q^2*a27;a30=2*E^3*F*q^2*a26;a31=E^3*F^3*a27;a32=2*E*F^3*q^2*a26;a33=2*E^3*F*q^2*a27;a34=2*E^3*F*q^2*(a26+a27);a35=2*E^3*F^3*a27;a36=8*E*F*q^4*a26;a37=2*E*F^3*q^2*a26;a38=2*E*F^2*q^3*a26;a39=E*F^4*q*(a27-a26);a40=E^3*F^2*q*(a27-a26);a41=2*E^3*F^3*(a26-a27);a42=2*E*F^2*q^3*a27;a43=2*F^3*q^3*(a27-a26);a44=2*E^2*F*q^3*(a27-a26);a45=E^2*F^3*q*(a27-a26);a46=E^4*F*q*(a27-a26);a47=4*E^2*q^4*a26;a48=E^2*F^4*a27;a49=4*F^2*q^4*a26;a50=E^4*F^2*a27;a51=2*F^4*q^2*(a26+a27);a52=2*E^4*q^2*(a26+a27);
a53=a28+a29-a30-a31;a54=a28-a31+a32-a33;a55=a34+a35-a36-a37-a29;a56=a38+a39+a40+a41-a42;a57=a43+a44+a45+a46;a58=a47+a48-a49-a50-a51-a52;
a53*((cosh(E*L))^2-(sinh(E*L))^2)+a54+a55*cosh(E*L)*sinh(F*L)+a56*cosh(E*L)*sin(F*L)+a57*sinh(E*L)*cos(F*L)+a58*sinh(E*L)*sin(F*L)=0

Answer 1

题主的问题不明确。44个变量，而方程有45个。无法求解。请完善问题。

追问

非常感谢您的回答，前面的方程都是可通过已知的数据表达出来的，相当于迭代出最后一个方程的系数，最后面那个才是需要求解的方程，只有L是未知的，像这样的方程能求解吗，大佬求指点啊，非常感谢！

追答

最后的方程还有a53，a54， a55，a56，a57，a58，L等未知变量，还是不能求解。

追问

a53，a54， a55，a56，a57，a58应该可以算出一个具体的值啊，为什么求解的时候没代进去，崩溃了，懂不起