Gruppeoppgaver for uke 4
Tilbake til hovedsidenTilbake til oversikten over l?sningsforslag
F?rst, en liten modifikasjon av matlabkoden i Allens appendix for bruk i oppgave 6:
clear all
set(0,'DefaultAxesFontSize', 18);
time=500;
b = 0.025;
d = b;
N=20;
r=0.004;
en=5;
T=zeros(N+1,N+1);
p=zeros(time+1,N+1);
p(1,6)=1;
v=linspace(0,N,N+1);
T(1,1)=1;
T(N+1,N+1)=1-d*N;
T(N,N+1)=d*N;
for i=2:N
??? T(i,i) = 1-i*(b+d);
??? T(i-1,i) = i*d;
??? T(i+1,i) = i*b;
end
for t=1:time
??? y=T*p(t,:)';
??? p(t+1,:)=y';
end
pm(1,:) = p(1,:);
for t=1:time/en;
??? pm(t+1,:)=p(en*t,:);
end
mesh([0:1:N],[0:en:time], pm);
xlabel('State');
ylabel('Time, n');
zlabel('Probability');
view(140,30)
Nils legger ut en R-l?sning p? ekstraoppgave B. For ekstraoppgave A finnes matlabkode i boken som sv?rt lett kan modifiseres.
