JFF昨天问题的MATLAB程序

送交者: Amsel 于 2008-12-16, 14:03:07:

mesh = 500;
prob_A = [ones(mesh,1)/mesh; 0];
prob_S = 1;
N = 10;

for count = 1:N
　　prob_S =conv(prob_S,prob_A);
　　disp([count,trapz(prob_s(count*mesh/2+mesh*0.5:end-1))])
　　plot(prob_S), hold on
end

（trapz求积分有些误差）

积分那行我写的是prob_S怎么copy过来就变成prob_s了 (无内容) - Amsel (0 bytes) 2008-12-16, 14:54:37 (284552)
THX! (无内容) - JFF (0 bytes) 2008-12-16, 14:43:43 (284550)
- 不用谢，练练手。画图那句改成这样更直观 - Amsel (44 bytes) 2008-12-16, 14:48:23 (284551)