Tassos icon

Matlab labs

Tassos | PRO | 01/14/15 06:09:41 PM UTC | 0 ⭐ | 1925 👁️ | Never ⏰ | []
MatLab |

5.93 KB

|

None

|

0 👍

/

0 👎

=======================================================================================================
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 -> άπειρο
 
;;;

Comments