ch6_foreman.m

these are the Mtlab functions to plot the figures in the book Communication system (by Haykin)

%generating psfiles for the foreman images
show('foreman69.Y',352,288);
print -deps2 foreman69.eps

show('foreman72.Y',352,288);
print -deps2 foreman72.eps

show('foreman_mvshalf',352,288);
print -deps2 foreman_EBMA_MV.eps

show('foreman_EBMA_pimg_72.Y',352,288);
%show('foreman_reconshalf',352,288);
print -deps2 foreman_EBMA.eps


show('foreman_HBmvs9.raw',352,288);
print -deps2 foreman_HBMA_MV.eps

show('foreman_HBrecons9.raw',352,288);
print -deps2 foreman_HBMA.eps

show('foreman_mesh_69.mesh',352,288);
print -deps2 foreman_HMMA_MV.eps

show('foreman_mesh_pimg_72.Y',352,288);
print -deps2 foreman_HMMA.eps


orig=fread(fopen('foreman72.Y'),[352,288]);
%ebma=fread(fopen('foreman_reconshalf'),[352,288]);
ebma=fread(fopen('foreman_EBMA_pimg_72.Y'),[352,288]);
ebmamse=sum(sum((orig-ebma).*(orig-ebma)))/352/288
psnrebma=10*log10(255*255/ebmamse)

hbma=fread(fopen('foreman_HBrecons16.raw'),[352,288]);
hbmamse=sum(sum((orig-hbma).*(orig-hbma)))/352/288
psnrhbma=10*log10(255*255/hbmamse)


hmma=fread(fopen('foreman_mesh_pimg_72.Y'),[352,288]);
hmmamse=sum(sum((orig-hmma).*(orig-hmma)))/352/288
psnrhmma=10*log10(255*255/hmmamse)