MatLab

-Transmisin y Modulacin-
- Vct or Ju l i n N egrn P l a t a - 44581
- A i t or V zq u ez Or t ega - 51512
Mayo-2010
NDICE
1-Funciones Desarrolladas (Cdigo y Funcionamiento)
1.1-Transmisor
1.2-Medio
1.3-Receptor
1.4-PlotDEPvp
1.5-PruebaCom
1.6-Pot
2-Pruebas del Correcto Funcionamiento
(Portadoras)
2.1-Seno
2.2-Tringulo
1-Funciones Desarrolladas (Cdigo y Funcionamiento)
1.1-Transmisor
f unct i on [ y, t ] =TRANSMI SOR( TR, x, t , At e, d, N0, Ac, f c, f s, W)
%se i nvoca : >>TRANSMI SOR( TR, x, t , At e, d, N0, Ac, f c, f s, W)
%TR=0 TRANSMI SI ON EN BANDA BASE
%TR=1 TRANSMI SI ON EN DBL
%TR=. . . . .
%x=seal modul ador a
%t =vect or de t i empos de l a seal modul ador a
%At e=At enuaci n del medi o en Db/ Km
%d=di st anci a ent r e Tr ansmi sor y Recept or en km
%N0=Pot enci a de Rui do en W/ Hz
%W=Ancho de Banda de l a t r ansmi si n. [ Hz]
%devuel ve [ y, t ]
%y=seal en Recepci n( despues de l a cadena de t r ansmi si n)
%t =vect or de t i empos de ' y'
i f TR==1
%W=Ancho de banda de t r ansmi si n( depende de l a seal de ent r ada x)

%TR=1 - > Tr ansmi si n de l a seal en DBL
pl ot ( t , x)
t i t l e( ' Seal en Banda Base' )
gr i d on
pause
pl ot DEPvp( x, ' DEP de l a Seal en Banda Base' ) ;
%gr i d on
pause
xmax=max( x) ;
xn=x/ xmax;
xt =Ac*xn. *cos( 2*pi *f c*t ) ;
pl ot ( t , xt ) ;
t i t l e( ' Seal Modul ada en DBL' ) %Aqui se pi nt a l a seal modul ada en DBL
gr i d on
pause
[ Xt , Om] =pl ot DEPvp( xt , ' Espect r o de l a Seal Modul ada en DBL' ) ;
%Aqu se pi nt a el espect r o de l a seal modul ada en dBL
%t i t l e( ' Espect r o de l a Seal Modul ada en DBL' )
pause
[ xm, N0m] =MEDI O( xt , t , At e, d, N0) ;
pot N=N0m*W;
SNRsF=( ( pot ( xm) ) / pot N)
st em( [ 0 1 2 3] , [ 0 pot ( xm) pot N 0] ) ;
t i t l e( ' 0 - Pot . de Seal - Pot . de Rui do - 0 ' ) ;
gr i d on
pause
[ y t ] =RECEPTOR( TR, f s, f c, xm, t , W)
el sei f TR==0 %TR=0 - >
Tr ansmi si n de l a seal en Banda Base
pl ot ( t , x)
t i t l e( ' Seal en banda base' )
gr i d on
pause
pl ot DEPvp( x, ' DEP de l a Seal en Banda Base' ) ;
pause
[ xm, N0m] =MEDI O( x, t , At e, d, N0) ;
pot N=N0m*W;
SNRe=( ( pot ( xm) ) / pot N)
st em( [ 0 1 2 3] , [ 0 pot ( xm) pot N 0] ) ;
t i t l e( ' 0 - Pot . de Seal - Pot . de Rui do - 0 ' ) ;
gr i d on
pause
[ y t ] =RECEPTOR( TR, f s, f c, xm, t , W)
el se
y=( ' Bander a de Sel ecci on de Modul aci n I ncor r ect a' )

end
end
1.2-Medio
f unct i on [ xm, N0m] =MEDI O( x, t , At e, d, N0)
xp1=x/ ( 10^( At e*d/ 10) ) ;
pl ot ( t , xp1)
t i t l e( ' Seal ( por el medi o) , at enuada y si n r ui do' )
gr i d on
pause
pl ot DEPvp( xp1, ' Espect r o de l a Seal ( por el medi o) , at enuada y si n r ui do' ) ;
pause
noi se=sqr t ( N0/ 2) *r andn( 1, l engt h( xp1) ) ;
xm=xp1+noi se;
pl ot ( t , noi se) ;
t i t l e( ' r ui do( t ) ' )
gr i d on
pause
pl ot DEPvp( noi se, ' DEP de Rui do' ) ;
pause
N0m=N0/ 2
pl ot ( t , xm) ; %Aqui se pi nt a l a
seal modul ada una vez pasa el medi o
t i t l e( ' Seal ( ent r a al Recept or ) , at enuada y CON r ui do' )
gr i d on;
pause
pl ot DEPvp( xm, ' Espect r o de l a Seal ( ent r a al Recept or ) , at enuada y CON r ui do' )
%Aqui se pi nt a el espect r o de l a seal modul ada una vez ha
pasado el medi 0
end
1.3-Receptor
f unct i on[ yt , t ] =RECEPTOR( TR, f s, f c, x, t , W)
i f TR==1
%f i l t r o de pr edet ecci n;
N=10;
%Or den 3 del f i l t r o de moment o
OmC=2*f c/ f s;
OmW=2*W/ f s; %Wen Her ci os
Omcor t e=[ OmC- OmWOmC+OmW]
b=f i r 1( N, Omcor t e)
[ xf ] =f i l t er ( b, [ 1] , x)
pl ot ( t , xf ) ; %Aqui se pi nt a l a seal
modul ada pasada por medi o y f i l t r ada paso banda
gr i d on
t i t l e( ' Seal Modul ada pasada por el medi o y f i l t r ada paso banda' )
pause
pl ot DEPvp( xf , ' DEP de Seal Modul ada, pasada por el medi o y f i l t r ada paso banda' ) ;
%Aqui se pi nt a el espect r o de l a seal modul ada pasada por medi o y f i l t r ada paso
banda
gr i d on
pause
D=cos( 2*pi *f c*t )
y=xf . *D;
%Fi l t r o paso baj o
N=10;
b=f i r 1( N, OmW)
[ yt ] =f i l t er ( b, [ 1] , y) ;
pl ot ( t , yt ) ;
t i t l e( ' Seal Reci bi da' )
gr i d on
pause
pl ot DEPvp( yt , ' Espect r o de l a Seal Reci bi da' ) ;
pause
el sei f TR==0
N=10;
b=f i r 1( N, 2*W/ f s, GAUSSWI N( N+1) ) ;
[ yt ] =f i l t er ( b, [ 1] , x) ;
pl ot ( t , yt ) ;
t i t l e( ' Seal Reci bi da' )
gr i d on
pause
pl ot DEPvp( yt , ' Espect r o de l a Seal Reci bi da' ) ;
el se
yt =' ERROR POR DEMODULACI N I NCORRECTA'
end
end

1.4-PlotDEPvp
f unct i on [ Dep, w] = pl ot Depvp ( x, t xt )
%Di buj a el mdul o de l a DFT ( x) ent r e - pi y pi . Con una separ aci n ent r e
%r ayas espect r al es de 2pi / ( 5*l engt h ( x) )
%Devuel ve el mdul o de l a DFT ( X) y el vect or de f r ecuenci as ( w) en que se ha
%medi do
l enX = 5*l engt h ( x) ;
X = f f t shi f t ( ( abs( f f t ( x, l enX) ) ) ) ;
Dep = ( X. ^2) / l enX;
del t a = 2*pi / l enX; %Separ aci n de l as r ayas espect r al es
wMi nus = [ - pi : del t a: - del t a] ;
i f r em( l enX, 2) == 0,
%Sol o par a l ongi t udes ( l enX) par es
wPl us = [ 0: del t a: pi - del t a] ;
el se
%Sol o par a l ongi t udes ( l enX) i mpar es
wPl us = [ 0: del t a: pi ] ;
end
w = [ wMi nus wPl us] ;
pl ot ( w, Dep) , xl abel ( ' r ads' ) ; t i t l e ( t xt ) ;
gr i d on;


1.5-PruebaCom
%Scr i pt de t est eo de si mul ador ' BB' y ' DBL'
f m1=150;
f m2=300;
f s=16000;
f c1=2e3;
Bt dbl 1=2*f m1;
Bt bb1=f m1;
%i mpr esi on de Omegas
Om1=f m1*2*pi / f s
Omc1=f c1*2*pi / f s
Om2=2*pi *f m2/ f s
suma=Om1+Omc1
t 1=[ 0: ( 1/ f s) : 5*( 1/ f m1) - ( 1/ f s) ] ;
%1. seno
x1=si n( 2*pi *t 1*f m1) ;
%2. yuxt aposi ci on de t onos
f c2=10000;
f m3=21;
f m4=146;
f m5=f m2;
%i mpr esi on de Omegas
Om3=2*pi *f m3/ f s
Om4=2*pi *f m4/ f s
Om5=2*pi *f m5/ f s
t 2=[ 0: ( 1/ f s) : 5*( 1/ ( f m3) ) - ( 1/ f s) ] ;
a=si n( 2*pi *t 2*f m3) ;
b=si n( 2*pi *t 2*f m4) ;
c=si n( 2*pi *t 2*f m5) ;
x2=a+b+c;
Bt dbl 2=2*f m5;
Bt bb2=f m5;
%3. t r i angul o
T=- 1: 1/ 10: 1;
f s3=10; %Como se apr eci a en l a l i nea ant er i or , se muest r ea a 10 Hz.
xi 3=t r i pul s( T) ;
x3=[ xi 3 xi 3 xi 3]
t 3=[ 1: 1: l engt h( x3) ] ;
pl ot ( t 3, x3) ;
f m7=1/ l engt h( xi 3) ;
f c3=25*f m7;
Om7=2*pi *f m7/ f s3
Omc3=2*pi *f c3/ f s3
Bt dbl 3=2*f m7;
Bt bb3=f m7;
%pr ueba seno DBL
TRANSMI SOR( 1, x1, t 1, 1, 4, 10e- 5, 1, f c1, f s, Bt dbl 1)
%pr ueba seno Banda Base
TRANSMI SOR( 0, x1, t 1, 1, 4, 10e- 5, 1, f c1, f s, Bt bb1)
%Suma de t onos t r ansmi t i dos en DBL
%TRANSMI SOR( 1, x2, t 2, 1, 4, 10e- 7, 1, f c2, f s, Bt dbl 2)
%suma de t onos t r ansmi t i dos en Banda Base
%TRANSMI SOR( 0, x2, t 2, 1, 4, 10e- 7, 1, f c2, f s, Bt bb2)
%seal t r i angul ar t r ansmi t i da en DBL
TRANSMI SOR( 1, x3, t 3, 1, 4, 10e- 5, 1, f c3, f s3, 1. 05)
%seal t r i angul ar t r ansmi t i da en Banda Base
TRANSMI SOR( 0, x3, t 3, 1, 4, 10e- 5, 1, f c3, f s3, Bt bb3)
%se i nvoca : >>TRANSMI SOR( TR, x, t , At e, d, N0, Ac, f c, f s, W)
1.6-Pot
f unct i on [ p] =pot ( x) ;
p=sum( x. ^2) / ( l engt h( x) ) ;
end
2-Pruebas del Correcto Funcionamiento
2.1-Seno
Banda Base
0 0.005 0.01 0.015 0.02 0.025 0.03 0.035
-1
-0.8
-0.6
-0.4
-0.2
0
0.2
0.4
0.6
0.8
1
Seal en banda base
-0.2 -0.1 0 0.1 0.2 0.3
5
10
15
20
25
rads
DEP de la Seal en Banda Base
0 0.005 0.01 0.015 0.02 0.025 0.03 0.035
-0.4
-0.3
-0.2
-0.1
0
0.1
0.2
0.3
0.4
Seal (por el medio),atenuada y sin ruido
-0.1 -0.05 0 0.05 0.1 0.15
0
0.5
1
1.5
2
2.5
3
3.5
4
rads
Espectro de la Seal(por el medio),atenuada y sin ruido
Om1 =
0.0589
0 0.005 0.01 0.015 0.02 0.025 0.03 0.035
-0.1
-0.08
-0.06
-0.04
-0.02
0
0.02
0.04
0.06
0.08
ruido(t)
-4 -3 -2 -1 0 1 2 3 4
0
0.2
0.4
0.6
0.8
1
x 10
-3
rads
DEP de Ruido
0 0.005 0.01 0.015 0.02 0.025 0.03 0.035
-0.5
-0.4
-0.3
-0.2
-0.1
0
0.1
0.2
0.3
0.4
0.5
Seal (entra al Receptor) ,atenuada y CON ruido
0 0.5 1 1.5 2 2.5 3
0
0.01
0.02
0.03
0.04
0.05
0.06
0.07
0.08
0 - Pot.de Seal - Pot. de Ruido - 0
La grfica est referida a la densidad espectral de ruido.
SNRe=1.0577
Luego en unidades logartmicas: 10*log(1.0577)= 0.5610[dB]
0 0.005 0.01 0.015 0.02 0.025 0.03 0.035
-0.5
-0.4
-0.3
-0.2
-0.1
0
0.1
0.2
0.3
0.4
0.5
Seal Recibida
-0.2 -0.1 0 0.1 0.2 0.3 0.4 0.5
0
0.5
1
1.5
2
2.5
3
3.5
4
rads
Espectro de la Seal Recibida
Doble Banda Lateral
-0.3 -0.2 -0.1 0 0.1 0.2
0
5
10
15
20
25
rads
DEP de la Seal en Banda Base
Om1 =
0.0589
0 0.005 0.01 0.015 0.02 0.025 0.03 0.035
-1
-0.8
-0.6
-0.4
-0.2
0
0.2
0.4
0.6
0.8
1
Seal Modulada en DBL
-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8
0
1
2
3
4
5
6
rads
Espectro de la Seal Modulada en DBL
0.6 0.65 0.7 0.75 0.8 0.85 0.9
0
1
2
3
4
5
6
rads
Espectro de la Seal Modulada en DBL
BLS
Om1 =
0.0589
Omc1 =
0.7854
suma =
0.8443
Matlab devuelve al usuario las Omegas central y de una de las bandas.
0 0.005 0.01 0.015 0.02 0.025 0.03 0.035
-0.4
-0.3
-0.2
-0.1
0
0.1
0.2
0.3
0.4
Seal (por el medio),atenuada y sin ruido
-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1
rads
Espectro de la Seal(por el medio),atenuada y sin ruido
0 0.005 0.01 0.015 0.02 0.025 0.03 0.035
-0.025
-0.02
-0.015
-0.01
-0.005
0
0.005
0.01
0.015
0.02
ruido(t)
-4 -3 -2 -1 0 1 2 3 4
0
1
2
3
4
5
6
7
x 10
-5
rads
DEP de Ruido
0 0.005 0.01 0.015 0.02 0.025 0.03 0.035
-0.5
-0.4
-0.3
-0.2
-0.1
0
0.1
0.2
0.3
0.4
0.5
Seal (entra al Receptor) ,atenuada y CON ruido
0 0.5 1 1.5 2 2.5 3
0
0.005
0.01
0.015
0.02
0.025
0.03
0.035
0.04
0 - Pot.de Seal - Pot. de Ruido - 0
SNRe=2.6667
Luego en unidades logartmicas: 10*log(2.6667)= 9.8083 [dB]
0 0.005 0.01 0.015 0.02 0.025 0.03 0.035
-0.5
-0.4
-0.3
-0.2
-0.1
0
0.1
0.2
0.3
0.4
Seal Modulada pasada por el medio y filtrada paso banda
-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1
rads
DEP de Seal Modulada,pasada por el medio y filtrada paso banda
0 0.005 0.01 0.015 0.02 0.025 0.03 0.035
-0.2
-0.15
-0.1
-0.05
0
0.05
0.1
0.15
Seal Recibida
-0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4 0.5
0.05
0.1
0.15
0.2
0.25
0.3
0.35
0.4
0.45
0.5
rads
Espectro de la Seal Recibida
2.2TRIANGULO
Banda Base
0 10 20 30 40 50 60 70
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1
Seal en banda base
-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
rads
DEP de la Seal en Banda Base
0 10 20 30 40 50 60 70
0
0.05
0.1
0.15
0.2
0.25
0.3
0.35
0.4
Seal (por el medio),atenuada y sin ruido
-1 -0.5 0 0.5 1
0.01
0.02
0.03
0.04
0.05
0.06
0.07
0.08
0.09
0.1
0.11
rads
Espectro de la Seal(por el medio),atenuada y sin ruido
0 10 20 30 40 50 60 70
-0.02
-0.015
-0.01
-0.005
0
0.005
0.01
0.015
0.02
ruido(t)
-4 -3 -2 -1 0 1 2 3 4
0
0.5
1
1.5
2
2.5
3
3.5
4
4.5
x 10
-5
rads
DEP de Ruido
0 10 20 30 40 50 60 70
-0.05
0
0.05
0.1
0.15
0.2
0.25
0.3
0.35
0.4
0.45
Seal (entra al Receptor) ,atenuada y CON ruido
0 0.5 1 1.5 2 2.5 3
0
0.005
0.01
0.015
0.02
0.025
0.03
0 - Pot.de Seal - Pot. de Ruido - 0
Donte matlab devuelve: SNRePRE=1.0893e+004
En Unidades logartmicas: 10*log(1.0893e+004)=92.9588[dB]
0 10 20 30 40 50 60 70
-0.05
0
0.05
0.1
0.15
0.2
0.25
0.3
Seal Recibida
-1.5 -1 -0.5 0 0.5 1 1.5
0
0.02
0.04
0.06
0.08
0.1
rads
Espectro de la Seal Recibida
Doble Banda Lateral
0 10 20 30 40 50 60 70
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1
Seal en Banda Base
-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
rads
DEP de la Seal en Banda Base
0 10 20 30 40 50 60 70
-0.6
-0.4
-0.2
0
0.2
0.4
0.6
0.8
1
Seal Modulada en DBL
-4 -3 -2 -1 0 1 2 3 4
0
0.02
0.04
0.06
0.08
0.1
0.12
0.14
0.16
0.18
rads
Espectro de la Seal Modulada en DBL
0 10 20 30 40 50 60 70
-0.3
-0.2
-0.1
0
0.1
0.2
0.3
0.4
Seal (por el medio),atenuada y sin ruido
-4 -3 -2 -1 0 1 2 3 4
0
0.005
0.01
0.015
0.02
0.025
0.03
rads
Espectro de la Seal(por el medio),atenuada y sin ruido
Om7 =
0.0299
Omc3 =0.7480
0 10 20 30 40 50 60 70
-0.02
-0.015
-0.01
-0.005
0
0.005
0.01
0.015
0.02
ruido(t)
-4 -3 -2 -1 0 1 2 3 4
0
1
2
3
4
5
6
x 10
-5
rads
DEP de Ruido
N0m =
5.0000e-005
0 10 20 30 40 50 60 70
-0.3
-0.2
-0.1
0
0.1
0.2
0.3
0.4
Seal (entra al Receptor) ,atenuada y CON ruido
0 0.5 1 1.5 2 2.5 3
0
0.002
0.004
0.006
0.008
0.01
0.012
0.014
0 - Pot.de Seal - Pot. de Ruido - 0
SNRsF =
255.4256
Tal que: 10*log(255.4256)= 55.4293[dB]
0 10 20 30 40 50 60 70
-0.2
-0.15
-0.1
-0.05
0
0.05
0.1
0.15
0.2
0.25
Seal Modulada pasada por el medio y filtrada paso banda
-4 -3 -2 -1 0 1 2 3 4
0
0.002
0.004
0.006
0.008
0.01
0.012
0.014
0.016
0.018
0.02
rads
DEP de Seal Modulada,pasada por el medio y filtrada paso banda
0 10 20 30 40 50 60 70
-0.02
0
0.02
0.04
0.06
0.08
0.1
0.12
Seal Recibida
-0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 1.2
1
2
3
4
5
6
7
8
9
10
x 10
-3
rads
Espectro de la Seal Recibida
Mayo
.2010.