Matlab int函数积分

代码如下，错在哪里？？？symsxhA=83.8805+53.4*asin((h-x*0.07168-6)/6)+1.483*(h-x*0.07168-6)*(12*(h... 代码如下，错在哪里？？？
syms x h
A=83.8805+53.4*asin((h-x*0.07168-6)/6)+1.483*(h-x*0.07168-6)*(12*(h-x*0.07168)-(h-x*0.07168)^2)^0.5
V=int(A,x,-4,20.5) 展开

 我来答

1个回答

#热议# 生活中有哪些实用的心理学知识？

dbb627
推荐于2016-11-15 · TA获得超过1.2万个赞

知道大有可为答主

回答量：2127

采纳率：88%

帮助的人：1476万

我也去答题访问个人页

关注

展开全部

没错，只是表达式比较复杂，matlab求不出解析式，可以带入h求数值积分
syms x h
A=83.8805+53.4*asin((h-x*0.07168-6)/6)+1.483*(h-x*0.07168-6)*(12*(h-x*0.07168)-(h-x*0.07168)^2)^0.5
V=int(A,x,-4,20.5)

A =

(267*asin(h/6 - (112*x)/9375 - 1))/5 - (12*h - (2688*x)/3125 - (h - (224*x)/3125)^2)^(1/2)*((41524*x)/390625 - (1483*h)/1000 + 4449/500) + 5902565445994545/70368744177664

Warning: Explicit integral could not be found.

V =

int((267*asin(h/6 - (112*x)/9375 - 1))/5 - (12*h - (2688*x)/3125 - (h - (224*x)/3125)^2)^(1/2)*((41524*x)/390625 - (1483*h)/1000 + 4449/500) + 5902565445994545/70368744177664, x = -4..41/2)

追问

如何带入h的值求数值积分？ 而且数值积分是用quad函数求吧。与int有何区别

追答

h如果不是变量，很好求的
syms x h 
A=83.8805+53.4*asin((h-x*0.07168-6)/6)+1.483*(h-x*0.07168-6)*(12*(h-x*0.07168)-(h-x*0.07168)^2)^0.5;
h=5;
A=subs(A);
V=int(A,x,-4,20.5)
vpa(V,10)
 
V =
 
(333675*atan(346594059^(1/2)/2229))/896 - (333675*atan(292010411^(1/2)/7717))/896 + (578421*asin(743/6250))/640 - (1415733*asin(7717/18750))/640 + (801*38510451^(1/2))/1120 - (1131401685487*292010411^(1/2))/6562500000000 - (24476142357*346594059^(1/2))/312500000000 + 289225706853732705/140737488355328

ans =
 
1371.957937