一、问题:最近在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条回答 默认 最新
悬赏问题
- ¥15 求差集那个函数有问题,有无佬可以解决
- ¥15 【提问】基于Invest的水源涵养
- ¥20 微信网友居然可以通过vx号找到我绑的手机号
- ¥15 寻一个支付宝扫码远程授权登录的软件助手app
- ¥15 解riccati方程组
- ¥15 display:none;样式在嵌套结构中的已设置了display样式的元素上不起作用?
- ¥15 使用rabbitMQ 消息队列作为url源进行多线程爬取时,总有几个url没有处理的问题。
- ¥15 Ubuntu在安装序列比对软件STAR时出现报错如何解决
- ¥50 树莓派安卓APK系统签名
- ¥65 汇编语言除法溢出问题