 2.*(flm - (1+sqrt(-1))./2);
0065 
0066 disp(<span class="string">'Constraint on flms to generate real signal'</span>)
0067 <span class="keyword">for</span> el = 0:L-1
0068    ind = el*el + el + 1;
0069    flm(ind,1) = real(flm(ind,1));
0070    <span class="keyword">for</span> m = 1:el
0071       ind_pm = el*el + el + m + 1;
0072       ind_nm = el*el + el - m + 1;
0073       flm(ind_nm,1) = (-1)^m * conj(flm(ind_pm,1));
0074    <span class="keyword">end</span>  
0075 <span class="keyword">end</span>
0076 
0