Matlab labs

=======================================================================================================
13-1-15
=======================================================================================================

amps ( fdma/fdd )

synoliko euros : 12.5 Mhz

sixnotita diaxorismou : 10 Khz

euros ana kanali : 30 Khz

Na brite : ton ari8mo xriston & ke tin apodosi.
Lisi :


>> clear all
>> bt = 12 * 5 * 10 ^ 6 ;
>> bguard = 10 * 10 ^ 3;
>> bc = 30 * 10 ^ 3 ;
>> n = ( bt - 2 * bguard ) / bc ;
>> bt = 12.5 * 10 ^ 6 ;
>> na = n * bc / bt

na =

    4.7984

>>  n = ( bt - 2 * bguard ) / bc ;
>> na = n * bc / bt

na =

    0.9984

>>



=========================================================
Sistima GSM ( TDMA )
exi plesio pou apotelite apo 8 slots.
ena to ka8e slot exi 156,25 bits
ke ri8mo metadosis 270,833 kbps
na ipologistoun :
i diarkia tou bit ( Tb )
i diarkia tou slot ( Tslot )
i diarkia tou plesiou ( Tframe )
meta apo poso xrono 8a epanametadosi enas xristis ?
tb = 1 / R
 lisis :

>> R = 270.833 * 10 ^ 3; // to 10 ^ 3 epidi prepi na to metatre4oume
??? R = 270.833 * 10 ^ 3; // to 10 ^ 3 epidi prepi na to metatre4oume
                          |
Error: Unexpected MATLAB operator.

>> R = 270.833 * 10 ^ 3;
>> tb = 1/R

tb =

  3.6923e-006

>> tslot = tb * 156.25;
>> tframe = 8 * tslot

tframe =

    0.0046



=======================================================================================================
16-12-14
=======================================================================================================

( κάτι δεν έχω γράψει πολύ καλά σε αυτό το εργαστήριο )

Τύποι :
L = (4pi*d/λ)^2 = (4pi*d*f/c)^2
Ldb = 10 * log10 (L)


Για :
d = 1m - 10Km
f = 5 GHz
L = ;
Ldb = ;

Matlab :

f = 5 * 10 ^ 9;     % Επειδή είναι 5 GHz.
d = [1:10:10000];
c = 3 * 10 ^ 8 ;    % Ταχύτητα φωτός γνωστή πάντα.
L = ( (4 * pi * d * f ) /c ).^2;
Ldb = 10 * log10 ( L );
Subplot (2,1,1);
plot (d,L);
Subplot (2,1,2);
plot(d, Ldb);


-------------------------

Τύποι εκπομπής Friis :
Pr : Δέκτης.
Pt : Πομπός.
Aet : Ενεργός περιοχή κεραίας.
Pr = Pt * (AeT * AeR / d^2 * λ^2 );


Επειδή η κεραία μας λέει πως είναι παρεμβολική :
Ae = 0.56 A
   = 0.56 pi * r^2
   = 0.56 pi

Διάμετρος : 2m = 2r
d = 10 km
f = 8 MHz
Pt = 10 - 1000
Pr = ;

Matlab :

f = 8 * 10 ^ 6;     % Επειδή είναι σε MHz.
d = 10000;
Ae = 0.56 * pi;
Pt = [10 : 100];
c = 3 * 10^8;       % Ταχύτητα φωτός ( πάντα γνωστή ).
l = c/f;        % Τύπος : λ = c/f .
Pr = Pt * ( (Ae^2) / (d^2 * l^2) );
plot(pr);
plot(pt);


=======================================================================================================
9-12-14
=======================================================================================================

πάντα για τις παραβολικές είναι : 0,56

>> k = 0.56

k =

    0.5600

>> d =41000/(32*14)

d =

   91.5179

>> G = d*k

G =

   51.2500

>> gdb = 10*log10(G)

gdb =

   17.0969





----------------------------------------------------------------------
SOS
Η απώλεια στο ελεύθερο διάστημα ( fresnel)

L = ( (4πd)/λ )^2
d = αποστάτη
λ = μήκος κύματος

λ = c/f
Ldb = 10 log10(L)

από ο πάνω ( SOS!! ) να αποδείξω αυτό :

Ldb = 32.44+20log10(d) + 20log10(f)
f : MHz
d : km

----------------------------------------------------------------------

f = 12 GHz
d = 10 Km
να βρούμε το L , Ldb

λ = c/f

===

>> f = 12*10^9

f =

  1.2000e+010

>> c = 3*10^8

c =

   300000000

>> l = c/f

l =

    0.0250

>> L = ( (4*pi*10)/l )^2

L =

  2.5266e+007

>> Ldb = 10*log10(L)

Ldb =

   74.0254


=======================================================================================================
18-11-2014
=======================================================================================================
Θεωρητικό :
C : Χωρητικότητα.
B : Εύρος ζώνης.
S : Ισχύς σήματος.
N : Ισχύς θορύβου.

C = B * log2 (1 + (S / No*b) )      % (bps)

S/No = 25 Db

B : 1 Hz...10 Khz


C  ^
   |
   |
   |
   |
   ___________________>
   1           b^10^4


Matlab :

>> s_no_db = 25;

>> s_no = 10^(s_no_db/10);  % Κάνει το db καθαρό αριθμό

>> b = [0:10,11:5:1000,1010:10:5000,5000:100:10000];

>> c = b.*log2(1+s_no./b);

>> semilogx(b,c);

>> xlabel('b')

>> ylabel('capacity (c)')

>> title('Channel capacity (c): b vs c')


Τι συμβαίνει με την χωρητικότητα (C) όταν το εύρος τέμνει στο 0 ή στο άπειρο ;
Δηλαδή :
C = ;
S/No ->
S/No -> άπειρο

______________

C = ;
B -> 0
B -> άπειρο


Άλλο που μοιάζει μάλλον από την ίδια θεωρία :

Θεωρητικό :
C : Χωρητικότητα.
B : Εύρος ζώνης.
S : Ισχύς σήματος.
N : Ισχύς θορύβου.

N = No * B

C = B * log2 (1 + (S / No*b) )      % (bps)

S/No : -20 .. έως .. 30 Db

B : 3 Khz

C= ;

Matlab :

S_No_dB = [-20:0.1:30];
S_No = 10 .^ ( S_No_dB / 10 );
B = 3000;

C = B * log2 ( 1 + (S_No / B ) );   % !!~SOS~!!

semilogx(S_No,C);
title('Channel Capacity(C) : S/No vs C ');
xlabel('S/No');
ylabel('Capacity (C)');

b1 = 4 KHz   -> C1
b2 = 30 KHz  -> C2

( εδώ έχω μια γραφική αναπαράσταση στο χαρτί με 3 καμπύλες B, B1 , B2 )

semilogx(S_No, C, S_No, C1, S_No, C2);

Τι συμβαίνει με την χωρητικότητα όταν :

C  -> S/No -> 0
C  -> S/No -> άπειρο

;;;