由买买提看人间百态

boards

本页内容为未名空间相应帖子的节选和存档,一周内的贴子最多显示50字,超过一周显示500字 访问原贴
Computation版 - Matlab中dblquad作积分时出现的问题
相关主题
[合集] 请问大家在matlab里面画的图存为.bmp 还是.jpg 后再放到wor[转载] 求一个方程的解析解
[合集] 问问做学生时搞计算的 (转载)求解两个联立的二阶常微分方程 (我的方程)
fortran给数组赋初值一问[转载] 我的论文标题的问题
matlab可以画三维球坐标图吗?请教求解二阶ODE方程
为啥ode45没办法解这个PDEhelp! Fit envelope for a signal spectrum
请教怎样把一个Matlab画的图压小?a question about Theta(nlgn)
matlab中的robust fit可以画confidence bound for function吗?真诚请教解方程组的问题!!
怎么求这个转角?[转载] 请maple 或者 mathematica 的高手帮个忙!
相关话题的讨论汇总
话题: thegema话题: theta话题: omegac话题: thegemas话题: cos
进入Computation版参与讨论
1 (共1页)
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

1 (共1页)
进入Computation版参与讨论
相关主题
[转载] 请maple 或者 mathematica 的高手帮个忙!为啥ode45没办法解这个PDE
[转载] 这叫什么方向?请教怎样把一个Matlab画的图压小?
如何在maple下定义一个i×j矩阵,其中元素是ai,j?matlab中的robust fit可以画confidence bound for function吗?
求助:downhill simplex minimization怎么求这个转角?
[合集] 请问大家在matlab里面画的图存为.bmp 还是.jpg 后再放到wor[转载] 求一个方程的解析解
[合集] 问问做学生时搞计算的 (转载)求解两个联立的二阶常微分方程 (我的方程)
fortran给数组赋初值一问[转载] 我的论文标题的问题
matlab可以画三维球坐标图吗?请教求解二阶ODE方程
相关话题的讨论汇总
话题: thegema话题: theta话题: omegac话题: thegemas话题: cos