24小时热门版块排行榜

>论坛更新日志 (3447)
>虫友互识 (396)
>休闲灌水 (238)
>文献求助 (230)
>导师招生 (174)
>考博 (87)
>考研 (75)
>硕博家园 (54)
>博后之家 (53)
>论文投稿 (46)
>基金申请 (39)
>公派出国 (35)
>论文道贺祈福 (32)
>教师之家 (32)
>招聘信息布告栏 (28)
>找工作 (27)

返回列表

lenienter

至尊木虫 (知名作家)

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

回复此楼

» 猜你喜欢

售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急已经有5人回复
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急已经有6人回复
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急已经有3人回复
为什么中国大学教授们水了那么多所谓的顶会顶刊，但还是做不出宇树机器人？已经有3人回复
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急已经有4人回复
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急已经有3人回复
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急已经有4人回复
“人文社科而论，许多学术研究还没有达到民国时期的水平” 已经有5人回复
过年走亲戚时感受到了所开私家车的鄙视链已经有11人回复
什么是人一生最重要的？已经有4人回复

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

查看全部散金贴

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;