matlab莱斯特模型

建模书code实现
f = [ 0 , 0.2 , 1.8, 0.8, 0.2];
c = [0.5, 0.8 , 0.8, 0.1];
x = zeros(30,5);
x = double(x);
for i = 1 : 1
for j = 1 : 5
x(i,j) = 100;
end
end
%%
for i = 2 : 30
for j = 1 : 5
x(i,1) = x(i,1) + x(i-1,j)*f(j);
t = j + 1;
if t ~= 6
x(i,t) = x(i - 1,j) * c(j);
end
if t == 6
break;
end
end
end
s  = sum(x(2,:));
%%比值
y = zeros(30,5);
y = double(y);
for  i = 1 : 30
     for j = 1 : 5
         y(i,j) = x(i,j) / sum(x(i,:));
     end
end