张继民
2021-05-23 10:17
采纳率: 0%
浏览 164

matlab中int函数无法得出解

一、问题:最近在matlab中用int函数求积分时一直得不出解,软件左下角一直显示正忙。 二、已尝试过如下解决办法: 1、更换积分函数:可以得出解,怀疑原函数本身过于复杂; 2、更换积分上下限为常数:同样显示正忙。 综上,认为原函数的存在问题 三、求助各位大神: 这种情况如何处理?跪谢! 四、代码: syms x h u=1*10^-6; a_33=12.56*10^(-9); F=h; L=50*h; b=h; G_E=1.41*(10^10)*h^4; k_w=1.12*10^4; k_p=2.8*(10^7)*h^2; s_1=k_p/G_E; s_2=k_w/G_E; O_1=((s_1+(s_1^2-4*s_2)^0.5)/2)^0.5; O_2=((s_1-(s_1^2-4*s_2))^0.5/2)^0.5; m=O_1^3*cosh(O_1*L)-O_1*(O_2^2)*cosh(O_2*L); n=O_1^3*sinh(O_1*L)-(O_2^3)*sinh(O_2*L); t=O_1^3*sinh(O_1*L)-(O_1^2)*O_2*sinh(O_2*L); Z=2*G_E*(m^2-n*t); u_1=cosh(O_1*L)- cosh(O_2*L); u_2=O_2*sinh(O_2*L)-O_1*sinh(O_1*L); A31=O_1*O_1*sinh(O_1*x)-O_1*O_2*sinh(O_2*x); B31=O_2*O_2*cosh(O_2*x)-O_1*O_1*cosh(O_1*x); O3=(h*F*u*(m*O_1*u_1+t*u_2))/(2*(u^2)*b*O_1*h*(O_1*n*u_1+m*u_2)-Z*a_33*L); C31=(2*u*b*O_1*O3*n-F*m)/Z; D31=(2*u*b*O_1*O3*m-F*t)/Z; W3(2)=(C31*A31+D31*B31)^2; W03(2)=((-F*m/Z)*A31+(-F*t/Z)*B31)^2; Y31=int(W3(2),x,-h/2,h/2); Y32=int(W03(2),x,-h/2,h/2); Y33=Y32/Y31;

  • 写回答
  • 好问题 提建议
  • 追加酬金
  • 关注问题
  • 邀请回答

2条回答 默认 最新

相关推荐 更多相似问题