1和D2是从1到11的数,他们满足下面的条件:
D1(0)+D4(0)=12;
D2(k)=sk/45向正无穷取整;
D3(k)=sk1/45向正无穷取整;
D1(k)=D2(k)+D1(k+1)+D3(k-45);
D4(K)=D1(K-45)+D3(K)+D4(K+1);
D3(k)<=D4(K);
D2(k)<=D1(K)
我写了下面的程序,结果出来不能理解,错在哪里?
k=1:1046;
n=length(k);
sk1=zeros(n,1);
sk2=zeros(n,1);
for i=1:11
c=45:n;
m=length(c);
a=ones(m,1);
for k=1:45
D2(k)=0;
D3(k)=0;
D1(k)=i;
D4(k)=12-i;
end
for k=45:n
D2(k)=ceil(sk1(k)/45);
D3(k)=ceil(sk2(k)/45);
D1(k+1)=D1(k)-D2(k)+D3(k+1-45);
D4(k+1)=D4(k)-D3(k)+D2(k+1-45);
end
for k=45:n
if (D1(k)>=D2(k)) &&(D4(k)>=D3(k))
fff(k)=1;
else fff(k)=0;
end
end