r****a
发帖数: 3249 | 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****y
发帖数: 199 | 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****a
发帖数: 3249 | 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****a
发帖数: 3249 | 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****y
发帖数: 199 | 5 hehe... 搞错了，不好意思，呵呵。你的r既是函数变量又是上限？dblquad是求定积分
的，要给定上下限。
【在 r****a 的大作中提到】
: 明月清风,刚发现你给r赋了一个初值，但r是一个积分变量啊？？呵呵
:
: .5
|
|
