Matlab中dblquad作积分时出现的问题 - 未名空间MITBBS历史存档

国际科技财经博客移民网络热点娱乐民生时事公众号

Redian新闻

>未名空间

>Computation - 科学计算

Matlab中dblquad作积分时出现的问题

Matlab中dblquad作积分时出现的问题# Computation - 科学计算

r*a2008-12-05 08:12

1 楼

我在做下面那个积分时老是告诉我“??? Attempted to access y(1); index out of
bounds because numel(y)=0.”，请各位老大帮我看看我什么地方搞错了，式中除积分
变量z,r外都是赋有初值的变量。
p1=dblquad(@(z,r) ((z-(y-q/cos(theta)+r*sin(theta))/thegema)/omegaC)^(1.5)*
thegema*exp(-0.5*(z/thegemaS)^2)/((2*pi)^0.5*thegemaS), (y-q/cos(theta)+r*
sin(theta))/thegema,(y-q/cos(theta)+r*sin(theta))/thegema+omegaC,-0.5/
thegema,0.5/thegema,1e-12);
先谢了！

g*y2008-12-05 08:12

2 楼

try
y = 1;
theta = 0.5;
r= 2;
q =3;
thegemaS = 4;
thegema = 3;
omegaC=1;
dblquad(@(z,r) (((z-(y-q./cos(theta)+r.*sin(theta))./thegema)./omegaC).^(1.5
).*...
thegema.*exp(-0.5.*(z./thegemaS).^2)/((2*pi).^0.5.*thegemaS)), (y-q./cos(
theta)+r.*...
sin(theta))./thegema,(y-q./cos(theta)+r.*sin(theta))./thegema+omegaC,-0.5./.
..
thegema,0.5./thegema,1e-12)
answer is: 0.0773
my matlab version is R2007b

【在 r****a 的大作中提到】

: 我在做下面那个积分时老是告诉我“??? Attempted to access y(1); index out of
: bounds because numel(y)=0.”，请各位老大帮我看看我什么地方搞错了，式中除积分
: 变量z,r外都是赋有初值的变量。
: p1=dblquad(@(z,r) ((z-(y-q/cos(theta)+r*sin(theta))/thegema)/omegaC)^(1.5)*
: thegema*exp(-0.5*(z/thegemaS)^2)/((2*pi)^0.5*thegemaS), (y-q/cos(theta)+r*
: sin(theta))/thegema,(y-q/cos(theta)+r*sin(theta))/thegema+omegaC,-0.5/
: thegema,0.5/thegema,1e-12);
: 先谢了！

r*a2008-12-05 08:12

3 楼

多谢明月清风,未名的高人啊

.5

【在 g****y 的大作中提到】

: try
: y = 1;
: theta = 0.5;
: r= 2;
: q =3;
: thegemaS = 4;
: thegema = 3;
: omegaC=1;
: dblquad(@(z,r) (((z-(y-q./cos(theta)+r.*sin(theta))./thegema)./omegaC).^(1.5
: ).*...

r*a2008-12-05 08:12

4 楼

明月清风,刚发现你给r赋了一个初值，但r是一个积分变量啊？？呵呵

.5

【在 g****y 的大作中提到】

: try
: y = 1;
: theta = 0.5;
: r= 2;
: q =3;
: thegemaS = 4;
: thegema = 3;
: omegaC=1;
: dblquad(@(z,r) (((z-(y-q./cos(theta)+r.*sin(theta))./thegema)./omegaC).^(1.5
: ).*...

g*y2008-12-05 08:12

5 楼

hehe... 搞错了，不好意思，呵呵。你的r既是函数变量又是上限？dblquad是求定积分
的，要给定上下限。

【在 r****a 的大作中提到】

: 明月清风,刚发现你给r赋了一个初值，但r是一个积分变量啊？？呵呵
:
: .5