function V=mdp_event_spare(p,s1,s2,L1,L2,E1,E2,N)
%P为1*4矩阵，s1,s2,L1,L2,E1,E2,N均为实数，N为迭代次数%
V=zeros(s1+1,s2+1,N+1);
for n=0:N-1
  for i=1:s1
    for j=1:s2
     if i==1&j==1
         V(i,j,N-n)=p*[s1*V(i+1,j,N-n+1);s2*V(i,j+1,N-n+1);V(i,j,N-n+1)+E1;V(i,j,N-n+1)+E2];
     elseif i>=2&i<=&j==1
        V(i,j,N-n)=p*[(s1-i+1)*V(i+1,j,20-n+1)+(i-1)*V(i,j,N-n+1);s2*V(i,j,N-n+1);min([V(i-1,j,N-n+1) V(i,j,N-n+1)+E1]);min([V(i-1,j,N-n+1)+2 V(i,j,N-n+1)+E2])];
     elseif j>=2&j<=9&i==1
       V(i,j,N-n)=p*[s1*V(i,j,N-n+1);(s2-j+1)*V(i,j+1,N-n+1)+(j-1)*V(i,j,N-n+1);min([V(i,j,N-n+1)+L1 V(i,j,N-n+1)+E1]);min([V(i,j-1,N-n+1) V(i,j,N-n+1)+E2])];
     elseif i>=2&i<=7&j==9+1
       V(i,j,N-n)=p*[(s1-i+1)*V(i+1,j,20-n+1)+(i-1)*V(i,j,N-n+1);s2*V(i,j,N-n+1);min([V(i-1,j,N-n+1) V(i,j,N-n+1)+L1 V(i,j,N-n+1)+E1]);min([V(i,j-1,N-n+1) V(i-1,j,N-n+1)+L2 V(i,j,N-n+1)+E2])];
    elseif i==7+1&j>=2&j<=9
      V(i,j,N-n)=p*[s1*V(i,j,N-n+1);(s2-j+1)*V(i,j+1,N-n+1)+(j-1)*V(i,j,N-n+1);min([V(i-1,j,N-n+1) V(i,j,N-n+1)+L1 V(i,j,N-n+1)+E1]);min([V(i,j-1,N-n+1) V(i-1,j,N-n+1)+L2 V(i,j,N-n+1)+E2])];
    elseif i==7+1&j==9+1
     V(i,j,N-n)=p*[7*V(i,j,N-n+1);9*V(i,j,N-n+1);min([V(i-1,j,N-n+1) V(i,j,N-n+1)+L1 V(i,j,N-n+1)+E1]);min([V(i,j-1,N-n+1) V(i-1,j,N-n+1)+L2 V(i,j,N-n+1)+E2])];
     end
  end
end
end
比如比如输入P=[0.1 0.1 0.2 0.2];s1=4;s2=4;L1=10;E1=25;L2=2;E2=15;N=10;
mdp_event_spare(p,s1,s2,L1,L2,E1,E2,N)
显示??? Undefined function or method 'mdp_event_spare_p' for input arguments of type 'double'.