菜单

MATLAB求解二重积分案例。MATLAB求解二重积分案例。

2018年9月19日 - 中甲报道

凯鲁嘎吉 – 博客园

凯鲁嘎吉 – 博客园

 

 

http://www.cnblogs.com/kailugaji/

http://www.cnblogs.com/kailugaji/

 

 

 
定积分解决的凡同样维连续量求和的题目,而解决多维连续量之求与题材将用到又积分了。重积分是树立以定积分的根底及之,它的中坚思想也是将重积分化为定积分来计算,其中要是积分限的规定,这也是重新积分的难关所在。正是因为重积分从计算上来说仍是采用的定积分的法子,MATLAB系统并无提供特别的命函数来拍卖又积分,因此当咱们规定了积分限后循是动int()命令来拍卖又积分问题。有些积分区间形状比较复杂,为了便利表达积分的上下限常常将比较复杂的间距分割成几个相对简便易行的距离,然后针对两样之间隔分别积分,最后把各个积分结果相互加起。

 
定积分解决的凡同一维连续量求和之题目,而解决多维连续量之请与题材即将用到重积分了。重积分是白手起家以定积分的根底及之,它的主导思想为是用重积分化也定积分来算,其中重大是积分限的确定,这也是重复积分的难关所在。正是为重积分从计算上来说仍是运的定积分的艺术,MATLAB系统并不曾提供特别的吩咐函数来处理还积分,因此于我们规定了积分限后以是采取int()命令来处理还积分问题。有些积分区间形状比较复杂,为了方便表达积分的上下限常常把比较复杂的区间分割成多单相对简单的间隔,然后对两样的间距分别积分,最后将各个积分结果相互加起来。

计算图片 1

计算图片 2

其中,

其中,

图片 3

图片 4

图片 5

图片 6

图片 7

图片 8

wrr.m:

function tt=wrr(r,v,c,b,y,w)
syms u t
p=-(r/c)^v;
g=exp(p);
h=(b/y)*(t/y)^(b-1);
hu=(b/y)*(u/y)^(b-1);
f=int(g*hu,u,0,t);
tt=int(t*g*h*exp(-f),t,0,w);
tt=vpa(tt,10);

//结果为:
>> clear
>> tt=wrr(0.5,0.787,0.033,1.097,0.241,2)

tt =

0.002183207532
wrr.m:

function tt=wrr(r,v,c,b,y,w)
syms u t
p=-(r/c)^v;
g=exp(p);
h=(b/y)*(t/y)^(b-1);
hu=(b/y)*(u/y)^(b-1);
f=int(g*hu,u,0,t);
tt=int(t*g*h*exp(-f),t,0,w);
tt=vpa(tt,10);

//结果为:
>> clear
>> tt=wrr(0.5,0.787,0.033,1.097,0.241,2)

tt =

0.002183207532

于第二头函数的号积分,可以先转化成为逐次积分形式,利用int函数进行求解。  

对此第二冠函数的标志积分,可以先行转化成逐次积分形式,利用int函数进行求解。  

相关文章

发表评论

电子邮件地址不会被公开。 必填项已用*标注

网站地图xml地图