24小时热门版块排行榜

返回列表

lenienter

至尊木虫 (知名作家)

应助: 1 (幼儿园)
金币: 35150.5
帖子: 6631
在线: 604.7小时
虫号: 125191

[交流] 【求助】帮忙看个小程序。

我写了小程序，用来做数值积分的。老是出错。看不出是怎么回事。请各位高手帮忙看看。
p1=0.04;
p0=0.96;
D=64.03;
n=@(x) 0.00578*512.36*x.^0.5;
N=quadl(n,0.0,128.06);
m=@(x) x.^2.5*0.00578*512.36*(p1*p0-0.5*p1*p0*(exp(-x/44.58)+exp(-x/1.98)));
M=quadl(m,0.0,128.06);
n2=p1-M/N;

错误是：
??? Error using ==> mtimes
Inner matrix dimensions must agree.

Error in ==>
@(x)x.^2.5*0.00578*512.36*(p1*p0-0.5*p1*p0*(exp(-x/44.58)+exp(-x/1.98)))

Error in ==> quadl at 70
y = feval(f,x,varargin{:}); y = y(

.';

Error in ==> simple at 7
M=quadl(m,0.0,128.06);

回复此楼

» 猜你喜欢

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2011-03-02 22:23:53

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

xiegangmai

版主 (职业作家)

lenienter(金币+10): 多谢，我现在手头上没有MATLAB这个软件，不知道你运行下来，n2的值是多少？ 2011-03-03 00:31:20

引用回帖:

Originally posted by lenienter at 2011-03-02 22:23:53:
我写了小程序，用来做数值积分的。老是出错。看不出是怎么回事。请各位高手帮忙看看。
p1=0.04;
p0=0.96;
D=64.03;
n=@(x) 0.00578*512.36*x.^0.5;
N=quadl(n,0.0,128.06);
m=@(x) x.^2.5*0.00578*512.36*( ...

少了三个点号，与自变量x相关的地方，都需要加上。

CODE:

p1 = 0.04;

p0 = 0.96;

D = 64.03;

n = @( x ) 0.00578 * 512.36 * x .^ 0.5;

N = quadl( n, 0.0, 128.06 );

m = @( x ) x .^ 2.5 * 0.00578 * 512.36 .* ( p1 * p0 - 0.5 * p1 *p0 * ( exp( -x ./ 44.58 ) + exp( -x ./ 1.98 ) ) );

M = quadl( m, 0.0, 128.06 );

n2 = p1 - M / N;