Zero forcing EQ question
CASTalk.com Forum Index CASTalk.com
Discussion of DSP, FPGA, storage and embedded system.
 
 FAQFAQ   MemberlistMemberlist     RegisterRegister 
 ProfileProfile   Log in to check your private messagesLog in to check your private messages   Log inLog in 
 
Google
 
Web castalk.com
Zero forcing EQ question

 
Post new topic   Reply to topic    CASTalk.com Forum Index -> DSP
Author Message
August Reno
Guest





Posted: Fri Dec 09, 2005 1:16 am    Post subject: Zero forcing EQ question Reply with quote

Hello All!
I obtained a matlab code which do the zero forcing equalization on
an awgn channel, and plots the BER vs. SNR
Please look at this small code if it is correct, my question is:
where are the filter coefficients here or how can I obtain this, as I don't
see them ???
Thanks for you effort,
Reno


N: number of symbols we transmit
h: impulse response
snrDB: SNR in dB, actually:
snrdB_min = 1; snrdB_max = 10;
snrdB = snrdB_min:1:snrdB_max;
-----

function z = zf(N,h,snrdB)

lg=length(h);
y=sign(randn(1,N));
sk=filter(h,1,y);
snr=10.^(snrdB/10);

nu=normrnd(0,(1/sqrt(snr)),1,N);
x=sk+nu;
ykegalZF=filter(1,h,x);
aklinZF=sign(ykegalZF(1:N));
PelinZF=sum(abs(y-aklinZF)/2);
errorZF=PelinZF/N;

%semilogy(RSBdB, PelinZF,'r*'); hold on;

-------
Back to top
Fred Marshall
Guest





Posted: Fri Dec 09, 2005 9:16 am    Post subject: Re: Zero forcing EQ question Reply with quote

"August Reno" <a_reno@tiscali.it> wrote in message
news:dna3i9$5d0$1@namru.matavnet.hu...
Quote:
Hello All!
I obtained a matlab code which do the zero forcing equalization on
an awgn channel, and plots the BER vs. SNR
Please look at this small code if it is correct, my question is:
where are the filter coefficients here or how can I obtain this, as I
don't see them ???
Thanks for you effort,
Reno


N: number of symbols we transmit
h: impulse response
snrDB: SNR in dB, actually:
snrdB_min = 1; snrdB_max = 10;
snrdB = snrdB_min:1:snrdB_max;
-----

function z = zf(N,h,snrdB)

lg=length(h);
y=sign(randn(1,N));
sk=filter(h,1,y);
snr=10.^(snrdB/10);

nu=normrnd(0,(1/sqrt(snr)),1,N);
x=sk+nu;
ykegalZF=filter(1,h,x);
aklinZF=sign(ykegalZF(1:N));
PelinZF=sum(abs(y-aklinZF)/2);
errorZF=PelinZF/N;

%semilogy(RSBdB, PelinZF,'r*'); hold on;

Just read the definition for "filter" and you'll see how the coefficients
are used.
The coefficients are "h" and are *inputs* to the function.....

Fred
Back to top
August Reno
Guest





Posted: Fri Dec 09, 2005 5:16 pm    Post subject: Re: Zero forcing EQ question Reply with quote

Thank you I figured it out!
August
"Fred Marshall" <fmarshallx@remove_the_x.acm.org> az alábbiakat írta a
következő hírüzenetben: QrmdnfzrdvW0mgTenZ2dnUVZ_tudnZ2d@centurytel.net...
Quote:

"August Reno" <a_reno@tiscali.it> wrote in message
news:dna3i9$5d0$1@namru.matavnet.hu...
Hello All!
I obtained a matlab code which do the zero forcing equalization on
an awgn channel, and plots the BER vs. SNR
Please look at this small code if it is correct, my question is:
where are the filter coefficients here or how can I obtain this, as I
don't see them ???
Thanks for you effort,
Reno


N: number of symbols we transmit
h: impulse response
snrDB: SNR in dB, actually:
snrdB_min = 1; snrdB_max = 10;
snrdB = snrdB_min:1:snrdB_max;
-----

function z = zf(N,h,snrdB)

lg=length(h);
y=sign(randn(1,N));
sk=filter(h,1,y);
snr=10.^(snrdB/10);

nu=normrnd(0,(1/sqrt(snr)),1,N);
x=sk+nu;
ykegalZF=filter(1,h,x);
aklinZF=sign(ykegalZF(1:N));
PelinZF=sum(abs(y-aklinZF)/2);
errorZF=PelinZF/N;

%semilogy(RSBdB, PelinZF,'r*'); hold on;

Just read the definition for "filter" and you'll see how the coefficients
are used.
The coefficients are "h" and are *inputs* to the function.....

Fred
Back to top
 
Post new topic   Reply to topic    CASTalk.com Forum Index -> DSP All times are GMT
Page 1 of 1

 
You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot vote in polls in this forum




VoIP Electronics Powered by phpBB