有没有大神能把这个程序的框图画出来啊?救救孩子吧!!!!!
function coeff=quasi_syn(period,sample_N)
N=sample_N;
n=period;
coeff=zeros(1,n*N+1);
temp = 0;
value = 0;
for k1=1:1:N+1
g(1,1) = 0.5;
g(1,k1) = 1;
g(1,N+1) = 0.5;
end
for h=2:1:n
for k2 =1:1:(h*N+1)
for k2 =1:1:(h*N+1)
if ( k2 == 1)
g(h,k2) = g(h-1,1) * 0.5;
else if ( k2 < N+1)
for i = 1:1:(k2-1)
temp = temp + g(h-1,i);
end
g(h,k2) = temp + g(h-1,k2) * 0.5;
temp =0;
else if( k2 <= h*(N+1)/2)
for j = (k2-N+1):1:(k2-1)
value = value + g(h-1,j);
end
g(h,k2) = g(h-1,k2-N) * 0.5 + value + g(h-1,k2)*0.5;
value=0;
else
g(h,k2) = g(h,h*N+2-k2);
end
end
end
end
end
end
coeff=g(n,:)/N^n;
end
function coeff=quasi_syn(period,sample_N)
N=sample_N;
n=period;
coeff=zeros(1,n*N+1);
temp = 0;
value = 0;
for k1=1:1:N+1
g(1,1) = 0.5;
g(1,k1) = 1;
g(1,N+1) = 0.5;
end
for h=2:1:n
for k2 =1:1:(h*N+1)
for k2 =1:1:(h*N+1)
if ( k2 == 1)
g(h,k2) = g(h-1,1) * 0.5;
else if ( k2 < N+1)
for i = 1:1:(k2-1)
temp = temp + g(h-1,i);
end
g(h,k2) = temp + g(h-1,k2) * 0.5;
temp =0;
else if( k2 <= h*(N+1)/2)
for j = (k2-N+1):1:(k2-1)
value = value + g(h-1,j);
end
g(h,k2) = g(h-1,k2-N) * 0.5 + value + g(h-1,k2)*0.5;
value=0;
else
g(h,k2) = g(h,h*N+2-k2);
end
end
end
end
end
end
coeff=g(n,:)/N^n;
end