Prev: Call for Papers: World Congress on Engineering and Computer Science WCECS 2010
Next: Recalculating difference equation coefficients when the time base changes
From: tmnhuda09 on 4 Jun 2010 09:50 The real and imaginary have the same PDF(Gaussain PDF)after IFFT of complex signal. But in my matlab code it looks different. why? can anybody explain it please? here is my code N=10; ml=2; for p=0:1000-1 x=randint(N*ml,1)*2-1; ich = x((1:N),1); qch = x((N+1:ml*N),1); x1=ich+qch*i; y=ifft(x1,10); yr=real(y); yi=imag(y); env=sqrt(yr.^2+yi.^2); for i=1:10 real_p(p*10+i)=yr(i); ima_p(p*10+i)=yi(i); st(p*10+i)=env(i); end end Best regards T.M. Nazmul Huda
From: Greg Heath on 4 Jun 2010 16:51
On Jun 4, 9:50 am, "tmnhuda09" <tmnhuda09(a)n_o_s_p_a_m.gmail.com> wrote: > The real and imaginary have the same PDF(Gaussain PDF)after IFFT of complex > signal. But in my matlab code it looks different. why? can anybody explain > it please? here is my code > > N=10; > ml=2; > for p=0:1000-1 > x=randint(N*ml,1)*2-1; > ich = x((1:N),1); > qch = x((N+1:ml*N),1); > x1=ich+qch*i; > y=ifft(x1,10); > yr=real(y); > yi=imag(y); > env=sqrt(yr.^2+yi.^2); > for i=1:10 > real_p(p*10+i)=yr(i); > ima_p(p*10+i)=yi(i); > st(p*10+i)=env(i); > end > end > > Best regards > T.M. Nazmul Huda I don't have randint. However, I interpret the documentation to mean that randint(m,n) is the same as round(rand(m,n)). Both cases are distributed uniformly, not normally. Hope this helps. Greg |