No BSD License  

Highlights from
MAP Algorithm

4.66667

4.7 | 6 ratings Rate this file 13 Downloads (last 30 days) File Size: 3.09 KB File ID: #7883
image thumbnail

MAP Algorithm

by Bahadir Gedikbey

 

21 Jun 2005 (Updated 23 Jun 2005)

Implementation of MAP Algorithm.

| Watch this File

File Information
Description

This program includes:
[5 7] convolutional code (encoder) + BPSK + AWGN + MAP (decoder). It evaluates Bit Error Rate and plots it versus SNR(signal to Noise Ratio).

Required Products Communications System Toolbox
MATLAB release MATLAB 6.5 (R13)
Tags for This File  
Everyone's Tags
Tags I've Applied
Add New Tags Please login to tag files.
Comments and Ratings (8)
20 Aug 2005 Vu Martin

How is about in Simulink ?

26 Oct 2005 Ugo Seux

Limited, doesn't seem to accept รก priori information. Does it work with other than R=1/2 codes?

12 May 2006 h xd  
24 Dec 2006 li lini  
30 Jul 2007 ali sabery

hello That was very useful for me.Thanks BR

21 Jun 2008 Hammad Dilpazir

Dear Bahadir Gedikbey it is a good method but i think it has some mistake, The formula used to compute "gama" in your code is
gama11(k)=exp((-(r(i)-x111)^2-(r(i+1)-x112)^2)/2/sigma);
            gama12(k)=exp((-(r(i)-x121)^2-(r(i+1)-x122)^2)/2/sigma);
            gama23(k)=exp((-(r(i)-x231)^2-(r(i+1)-x232)^2)/2/sigma);
            gama24(k)=exp((-(r(i)-x241)^2-(r(i+1)-x242)^2)/2/sigma);
            gama31(k)=exp((-(r(i)-x311)^2-(r(i+1)-x312)^2)/2/sigma);
            gama32(k)=exp((-(r(i)-x321)^2-(r(i+1)-x322)^2)/2/sigma);
            gama43(k)=exp((-(r(i)-x431)^2-(r(i+1)-x432)^2)/2/sigma);
            gama44(k)=exp((-(r(i)-x441)^2-(r(i+1)-x442)^2)/2/sigma);

but i think it should be like

gama11(k)=exp((-(r(i)-x111)^2-(r(i+1)-x112)^2)/(2*sigma)); %see moons book page 595 equation or eq 14.21
            gama12(k)=exp((-(r(i)-x121)^2-(r(i+1)-x122)^2)/(2*sigma));
            gama23(k)=exp((-(r(i)-x231)^2-(r(i+1)-x232)^2)/(2*sigma));
            gama24(k)=exp((-(r(i)-x241)^2-(r(i+1)-x242)^2)/(2*sigma));
            gama31(k)=exp((-(r(i)-x311)^2-(r(i+1)-x312)^2)/(2*sigma));
            gama32(k)=exp((-(r(i)-x321)^2-(r(i+1)-x322)^2)/(2*sigma));
            gama43(k)=exp((-(r(i)-x431)^2-(r(i+1)-x432)^2)/(2*sigma));
            gama44(k)=exp((-(r(i)-x441)^2-(r(i+1)-x442)^2)/(2*sigma));

the sigma should be divided instead of multiplied
Thanks and to check it
hammaddilpazir@yahoo.com

01 Apr 2009 Shahriar Rahman

Dear Bahadir,

I want to implement Systematic encoder like kafes = poly2trellis(3,[7 5], 7). For that your code has to changed. And I changed it accordingly. But the performance is poor. Can you please help me?

shahriar52@gmail.com

11 Jun 2009 dragos paslaru

Hello!I want to introduce a non-recursive and nonsistematic code first, then i want a systematic and nonrecursive code.what should i modify in this : kafes=poly2trellis(3,[5 7]);

Please login to add a comment or rating.
Updates
23 Jun 2005

the Communications Blockset is not needed to run this demo. Actually, only the Communications Toolbox is necessary

Tag Activity for this File
Tag Applied By Date/Time
coding theory Bahadir Gedikbey 22 Oct 2008 07:50:53
information theory Bahadir Gedikbey 22 Oct 2008 07:50:53
map Bahadir Gedikbey 22 Oct 2008 07:50:53
coding Bahadir Gedikbey 22 Oct 2008 07:50:54
convolutional Bahadir Gedikbey 22 Oct 2008 07:50:54
bit error rate Bahadir Gedikbey 22 Oct 2008 07:50:54
noise radio Bahadir Gedikbey 22 Oct 2008 07:50:54
snr Bahadir Gedikbey 22 Oct 2008 07:50:54
snr AMER Arif 06 Nov 2010 17:48:34
bit error rate AMER Arif 06 Nov 2010 17:48:44
coding theory AMER Arif 06 Nov 2010 17:48:55
noise radio AMER Arif 06 Nov 2010 17:49:02

Contact us at files@mathworks.com