24小时热门版块排行榜

返回列表

lenienter

至尊木虫 (知名作家)

应助: 1 (幼儿园)
金币: 33817.7
帖子: 6435
在线: 604.4小时
虫号: 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);

回复此楼

» 猜你喜欢

拟解决的关键科学问题还要不要写已经有7人回复
请教限项目规定已经有3人回复
存款400万可以在学校里躺平吗已经有15人回复
Materials Today Chemistry审稿周期已经有6人回复
基金委咋了？2026年的指南还没有出来？已经有10人回复
基金申报已经有6人回复
推荐一本书已经有13人回复
国自然申请面上模板最新2026版出了吗？已经有17人回复
纳米粒子粒径的测量已经有8人回复
疑惑？已经有5人回复

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

查看全部散金贴

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;