早教吧作业答案频道 -->其他-->
matlab int函数没有显式积分 vpa两个例子:1.y=x*log(x^4)*asin(1/x^2);int(y,x,1,3);提示没有解析解vpa(ans)得到数值结果,和quad函数的结果比较也是一样的.2.syms x yvpa(int( int( sqrt(1+x^2*sin(y)),x,0,1),y,0,2*pi)
题目详情
matlab int函数没有显式积分 vpa
两个例子:
1.y=x*log(x^4)*asin(1/x^2);
int(y,x,1,3);
提示没有解析解
vpa(ans)
得到数值结果,和quad函数的结果比较也是一样的.
2.
syms x y
vpa(int( int( sqrt(1+x^2*sin(y)),x,0,1),y,0,2*pi))
也提示没有显示解,得到数值结果
ans=3.4441;
但如果用dblquad函数求数值:
ans=6.1879;
这是什么原因呢?
我想知道,在没有解析解时,vpa是通过什么算法得到数值解的?那么,这个结果什么时候是正确的?
两个例子:
1.y=x*log(x^4)*asin(1/x^2);
int(y,x,1,3);
提示没有解析解
vpa(ans)
得到数值结果,和quad函数的结果比较也是一样的.
2.
syms x y
vpa(int( int( sqrt(1+x^2*sin(y)),x,0,1),y,0,2*pi))
也提示没有显示解,得到数值结果
ans=3.4441;
但如果用dblquad函数求数值:
ans=6.1879;
这是什么原因呢?
我想知道,在没有解析解时,vpa是通过什么算法得到数值解的?那么,这个结果什么时候是正确的?
▼优质解答
答案和解析
楼主是用什么版本求出第2题的解是3.4441?
我用6.5、R2007b两个版本,等了好久都没有得到结果(R2007b算了好半天之后甚至干脆崩溃了),用R2012b计算,得到下面的结果:
intlib::intOverSet(piecewise([not 0 < sin(y) and (y/pi in Z_ or arg(sin(y)) == pi), int((sin(y)*x^2 + 1)^(1/2), x == 0..1)]), y, {0} union Dom::Interval([pi], [2*pi])) + 3.4441419676288028273509895776497
显然这个结果并不是3.4441,而是另外加上前面一大串的积分表达式.
6.5、R2007b符号计算内核是用Maple,R2012b符号计算内核是用MuPad,但不论用哪个内核,其内部的计算细节基本上是不可能知道的.
到目前为止,我使用vpa,在能得到具体结果时都是正确的,但有可能就像上面这样,得不到具体的数值,而是一个表达式.
用dblquad函数计算的结果是可信的,其实你用
ezmesh('sqrt(1+x^2*sin(y))',[0,1],[0,2*pi])
画出曲面来,可以看到,曲面的大部分高度值都接近1,其它突出的部分比起凹陷的部分应该稍微小一点,所以算出来的体积(该重积分的几何意义)应该比2*pi稍小,但不至于差太多.
我用6.5、R2007b两个版本,等了好久都没有得到结果(R2007b算了好半天之后甚至干脆崩溃了),用R2012b计算,得到下面的结果:
intlib::intOverSet(piecewise([not 0 < sin(y) and (y/pi in Z_ or arg(sin(y)) == pi), int((sin(y)*x^2 + 1)^(1/2), x == 0..1)]), y, {0} union Dom::Interval([pi], [2*pi])) + 3.4441419676288028273509895776497
显然这个结果并不是3.4441,而是另外加上前面一大串的积分表达式.
6.5、R2007b符号计算内核是用Maple,R2012b符号计算内核是用MuPad,但不论用哪个内核,其内部的计算细节基本上是不可能知道的.
到目前为止,我使用vpa,在能得到具体结果时都是正确的,但有可能就像上面这样,得不到具体的数值,而是一个表达式.
用dblquad函数计算的结果是可信的,其实你用
ezmesh('sqrt(1+x^2*sin(y))',[0,1],[0,2*pi])
画出曲面来,可以看到,曲面的大部分高度值都接近1,其它突出的部分比起凹陷的部分应该稍微小一点,所以算出来的体积(该重积分的几何意义)应该比2*pi稍小,但不至于差太多.
看了 matlab int函数没有...的网友还看了以下:
如图所示,设A重10N,B重20N,A、B间的动摩擦因数为0.1,B与地面的摩擦因数为0.2.问: 2020-04-09 …
写出符合条件的一元一次不等式(1)-2,-1,0,1都是不等式的解(2)不等式的整数解只有-1,0 2020-04-27 …
在标准状况下,0.0672LH2S气体通入含有1.0×10-3molX2O72-离子的酸性溶液中, 2020-05-13 …
质量相同的两个铜导线其长度之比为3:5则他们的电阻之比?一个带电金属球当它带的电荷量增加后(稳定) 2020-06-04 …
第一题:12/(-3-4分之1+1又3分之1)第二题:(-2又5分之4*1又6分之5-0.25*5 2020-06-10 …
大一逻辑代数中变量和常量的关系比如A+1=1,以前高中只有1,0这懂的,而大学这里加入了变量A,还 2020-06-11 …
已知4阶实对称矩阵A只有两个不同的特征值λ1,λ2,且A的属于λ1的特征向量仅有(1,0,0,1) 2020-06-16 …
某电解池中,在2s内各有1.0×1019个二价正离子和2.0×1019个一价负离子通过该截面,那么 2020-06-23 …
求初中虚词的意思:之,其,而,然,则,乃,以,者.为,夫,于每种要有1个例句(初中文言文,人教版) 2020-06-25 …
一根铜导线,横截面积为1.0mm2,载有1.0安电流,已知铜导线内自由电子的密度n=8.5×102 2020-06-26 …