%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % Creative Commons % Attribution-Noncommercial 2.5 India % You are free: % to Share — to copy, distribute and transmit the work % to Remix — to adapt the work % Under the following conditions: % Attribution. You must attribute the work in the manner % specified by the author or licensor (but not in any way % that suggests that they endorse you or your use of the work). % Noncommercial. You may not use this work for commercial purposes. % For any reuse or distribution, you must make clear to others the % license terms of this work. The best way to do this is with a % link to this web page. % Any of the above conditions can be waived if you get permission % from the copyright holder. % Nothing in this license impairs or restricts the author's moral rights. % http://creativecommons.org/licenses/by-nc/2.5/in/ % Checked for proper operation with Octave Version 3.0.0 % Author : Krishna % Email : [email protected] % Version : 1.1 % Date : 25 May 2008 % Details : Corrected the % Version : 1.0 % Date : 30 April 2008 % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % Script for plotting the eye diagram where transmit filtering % is performed by raised cosine filtering with alpha=0.5, alpha=1. clear N = 10^3; % number of symbols am = 2*(rand(1,N)>0.5)-1 + j*(2*(rand(1,N)>0.5)-1); % generating random binary sequence fs = 10; % sampling frequency in Hz % defining the sinc filter sincNum = sin(pi*[-fs:1/fs:fs]); % numerator of the sinc function sincDen = (pi*[-fs:1/fs:fs]); % denominator of the sinc function sincDenZero = find(abs(sincDen) < 10^-10); sincOp = sincNum./sincDen; sincOp(sincDenZero) = 1; % sin(pix/(pix) =1 for x =0 % raised cosine filter alpha = 0.5; cosNum = cos(alpha*pi*[-fs:1/fs:fs]); cosDen = (1-(2*alpha*[-fs:1/fs:fs]).^2); cosDenZero = find(abs(cosDen)<10^-10); cosOp = cosNum./cosDen;