yiban.zip_弹簧_弹簧 算法_悬臂梁_振动_梁

syms x b c E J k1 k2 L N P; E=3.35*10^11; J=31^3*172*10^(-36)/12; k1=1; k2=1; L=34.87*10^(-6); M=[-sin(x),-cos(x),sinh(x),cosh(x); -cos(x),sin(x),cosh(x),sinh(x); 1,E*J*(x/L)/(k2),1,-E*J*(x/L)/(k2); -E*J*((x/L)^3)/k1,1,E*J*((x/L)^3)/k1,1]; N=det(M) P=simplify(N)%将这里得到的多项式复制粘贴到下面紫色文字位置 ezplot('(72307464630957915344615351492200448*c^4)/261220194723978830517110887124311255848370493187636016781640625 + 2*cos(c)*cosh(c) - (5941917132907099*c*cos(c)*sinh(c))/7242215658284376846129147412480000 + (5941917132907099*c*cosh(c)*sin(c))/7242215658284376846129147412480000 - (72307464630957915344615351492200448*c^4*cos(c)*cosh(c))/261220194723978830517110887124311255848370493187636016781640625 + (1753744394526836458266664733973151744*c^3*cos(c)*sinh(c))/2599052504656833999404912699662505469054716875 + (1753744394526836458266664733973151744*c^3*cosh(c)*sin(c))/2599052504656833999404912699662505469054716875 + 2',[-5,5]) %line([-5,5],[0,0])%绘制横轴 %a1=fzero('(5941917132907099*x*cosh(x)*sin(x))/724221565828437684612914741248 - (18510710945525226328221529982003314688*x^4)/668723698493385806123803871038236814971828462560348202961 - (5941917132907099*x*cos(x)*sinh(x))/724221565828437684612914741248 - 2*cos(x)*cosh(x) + (18510710945525226328221529982003314688*x^4*cos(x)*cosh(x))/668723698493385806123803871038236814971828462560348202961 + (28059910312429383332266635743570427904*x^3*cos(x)*sinh(x))/4158484007450934399047860319460008750487547 + (28059910312429383332266635743570427904*x^3*cosh(x)*sin(x))/4158484007450934399047860319460008750487547 - 2',5) %a10=fzero('(5941917132907099*x*cosh(x)*sin(x))/7242215658284376846129147412480 - 2*cos(x)*cosh(x) - (5941917132907099*x*cos(x)*sinh(x))/7242215658284376846129147412480 - (4627677736381306582055382495500828672*x^4)/16718092462334645153095096775955920374295711564008705074025 + (4627677736381306582055382495500828672*x^4*cos(x)*cosh(x))/16718092462334645153095096775955920374295711564008705074025 + (14029955156214691666133317871785213952*x^3*cos(x)*sinh(x))/20792420037254671995239301597300043752437735 + (14029955156214691666133317871785213952*x^3*cosh(x)*sin(x))/20792420037254671995239301597300043752437735 - 2',5) %a100=fzero('(5941917132907099*x*cosh(x)*sin(x))/72422156582843768461291474124800 - (1156919434095326645513845623875207168*x^4)/417952311558366128827377419398898009357392789100217626850625 - (5941917132907099*x*cos(x)*sinh(x))/72422156582843768461291474124800 - 2*cos(x)*cosh(x) + (1156919434095326645513845623875207168*x^4*cos(x)*cosh(x))/417952311558366128827377419398898009357392789100217626850625 + (7014977578107345833066658935892606976*x^3*cos(x)*sinh(x))/103962100186273359976196507986500218762188675 + (7014977578107345833066658935892606976*x^3*cosh(x)*sin(x))/103962100186273359976196507986500218762188675 - 2',5) %eig(M) vpa(x,15)