examp.m

来自「这是在网上下的一个东东」· M 代码 · 共 24 行

%
% This script simply plots the posterior distributions for Example 11.1
%
m=(5.0:0.1:15.0)';
for i=1:length(m)
  p1(i)=(1/sqrt(2*pi))*exp(-(10.3-m(i))^2/2);
  p2(i)=(1/sqrt(pi))*exp(-(10.2-m(i))^2);
end
figure(1);
bookfonts;
plot(m,p1,'k');
xlabel('m');
ylabel('q(m|d=10.3 ug)');
axis([5 15 0 0.7]);
%print -deps post1.eps
figure(2);
bookfonts;
plot(m,p2,'k');
xlabel('m');
ylabel('q(m|d_{1}=10.3 ug, d_{2}=10.1 ug)');
axis([5 15 0 0.7]);
%print -deps post2.eps