Файл: Курсовая работа по дисциплине Теория электрической связи (наименование учебной дисциплины согласно учебному плану).docx

ylabel( 'G_I(w), G_B(w), G_C(w)' )

% импульсы Найквиста-------------------------------------------------------

syms t;

T = 1;

% для в = 0

X_0 = ( sin( pi*t/T ) )/( pi*t/T );

% для в = 0.25

X_0_25 = ( sin( pi*t/T ) * ( cos( pi*0.25*t/T ) ) ) / ( ( pi*t/T ) * ( 1 - ( 2*0.25*t )^2 /( T^2 ) ) );

% для в = 0.5

X_0_5 = ( sin( pi*t/T ) * ( cos( pi*0.5*t/T ) ) ) / ( ( pi*t/T ) * ( 1 - ( 2*0.5*t )^2 /( T^2 ) ) );

% для в = 1

X_1 = ( sin( pi*t/T ) * ( cos( pi*t/T ) ) ) / ( ( pi*t/T ) * ( 1 - ( 2*t )^2 /( T^2 ) ) );

fplot( X_0, 'g--', 'Linewidth', 2 )

hold on

fplot( X_0_25, 'b', 'Linewidth', 2 )

hold on

fplot( X_0_5, 'r', 'Linewidth', 2 )

hold on

fplot( X_1, 'k--', 'Linewidth', 2 )

grid

xlabel( 't' )

ylabel( 'x(t)' )

legend( '\beta = 0', '\beta = 0.25', '\beta = 0.5', '\beta = 1' )

% спектральная плотность импульса Найквиста--------------------------------

syms f;

T = 1;

b_0 = 0;

b_0_25 = 0.25;

b_0_5 = 0.5;

b_1 = 1;

S_0 = piecewise( 0 <= abs(f) <= ((1-b_0)/(2*T)), T, ((1-b_0)/(2*T)) <= abs(f) <= ((1+b_0)/(2*T)), T/2*(1+cos(pi*T/b_0*(abs(f)-((1-b_0)/(2*T))))), ( (1+b_0)/(2*T) ) < abs(f), 0 );

S_0_25 = piecewise( 0 <= abs(f) <= ((1-b_0_25)/(2*T)), T, ((1-b_0_25)/(2*T)) <= abs(f) <= ((1+b_0_25)/(2*T)), T/2*(1+cos(pi*T/b_0_25*(abs(f)-((1-b_0_25)/(2*T))))), ( (1+b_0_25)/(2*T) ) < abs(f), 0 );

S_0_5 = piecewise( 0 <= abs(f) <= ((1-b_0_5)/(2*T)), T, ((1-b_0_5)/(2*T)) <= abs(f) <= ((1+b_0_5)/(2*T)), T/2*(1+cos(pi*T/b_0_5*(abs(f)-((1-b_0_5)/(2*T))))), ( (1+b_0_5)/(2*T) ) < abs(f), 0 );

S_1 = piecewise( 0 <= abs(f) <= ((1-b_1)/(2*T)), T, ((1-b_1)/(2*T)) <= abs(f) <= ((1+b_1)/(2*T)), T/2*(1+cos(pi*T/b_1*(abs(f)-((1-b_1)/(2*T))))), ( (1+b_1)/(2*T) ) < abs(f), 0 );

fplot( S_0, 'g--', 'Linewidth', 2 )

hold on

fplot( S_0_25, 'b', 'Linewidth', 2 )

hold on

fplot( S_0_5, 'r', 'Linewidth', 2 )

hold on

fplot( S_1, 'k--', 'Linewidth', 2 )

xlim([ -2 2 ])

grid

xlabel( 'f' )

ylabel( 'S_x(f)' )

legend( '\beta = 0', '\beta = 0.25', '\beta = 0.5', '\beta = 1' )

% два графика ипульс Найквиста и импульс со спектр. плотн. под корнем------

syms w;

S = piecewise( -2*pi/T <= w <= 2*pi/T, T/2*( 1 + cos(w*T/2) ), abs(w) > 2*pi/T, 0 );

S_x1 = piecewise( -2*pi/T <= w <= 2*pi/T, sqrt(T)*cos(w*T/4), abs(w) > 2*pi/T, 0 );

fplot( S, 'r--', 'Linewidth', 2 )

hold on

fplot( S_x1, 'b', 'Linewidth', 2 )

grid

xlim([ -7.5 7.5 ])

xlabel( '\omega' )

ylabel( 'S(\omega)' )

legend( 'S_x(\omega) ', 'S_x_1(\omega)')

% импульс Найквиста c в = 1 и искомый ипульс x_1(t) от S_x1----------------

syms t;

T = 1;

X_1 = ( sin( pi*t/T ) * ( cos( pi*t/T ) ) ) / ( ( pi*t/T ) * ( 1 - ( 2*t )^2 /( T^2 ) ) );

X1 = 1/sqrt(T)*( (sin((2*pi/T)*(t-T/4)))/((2*pi/T)*(t-T/4)) + (sin((2*pi/T)*(t+T/4)))/((2*pi/T)*(t+T/4)) );

fplot( X_1, 'r--', 'Linewidth', 2 )

hold on

fplot( X1, 'b', 'Linewidth', 2 )

grid

xlabel( 't' )

ylabel( 'x(t)' )

legend( 'x(t) ', 'x_1(t)')

xlim([ -4 4 ])

% корреляционные функции случ проц Iф(t) Qф(t)-----------------------------

syms t;

T = 1;

X_1 = ( ( sin( pi*t/T ) * ( cos( pi*t/T ) ) ) / ( ( pi*t/T ) * ( 1 - ( 2*t )^2 /( T^2 ) ) ) )/1.27^2;

fplot( X_1, 'b', 'Linewidth', 2 );

grid

xlabel( '\tau' )

ylabel( 'B_I_ф(\tau) = B_Q_ф(\tau)' )

xlim([ -4 4 ])

% спектральные плотности случ проц Iф(t) Qф(t)-----------------------------

syms w;

T = (13.55*10^(-6))/2;

G = piecewise( abs( w ) <= ( 2*pi/T ), T/( 1.27^2 * 2 ) * ( 1 + cos( w*T/2 ) ), abs( w ) > ( 2*pi/T ), 0 );

fplot( G, 'b', 'Linewidth', 2 )

grid

xlabel( '\omega' )

ylabel( 'G_I_ф(\omega) = G_Q_ф(\omega)' )

xlim([ -9 9 ])

% корреляционная ф-ия B_S(t)-----------------------------------------------

---

syms t;

T_c = (13.55*10^(-6))/2;

w_c = 4*pi/( 13.55*10^(-6) );

B_S = ( ( sin( pi*t/T_c ) * ( cos( pi*t/T_c ) ) ) / ( ( pi*t/T_c ) * ( 1 - ( 2*t )^2 /( T_c^2 ) ) ) ) * cos(w_c*t)/ 1.27^2 ;

fplot( B_S, 'b', 'Linewidth', 2 );

grid

xlabel( '\tau' )

ylabel( 'B_S(\tau)' )

xlim([ -20*10^(-6) 20*10^(-6) ])

ylim([ -0.55 0.65 ])

% спектральная плотность мощности G_S(t)-----------------------------------

syms w;

T = (13.55*10^(-6));

w_c = 4*pi/T;

G = piecewise( w_c/2 <= w <= ( w_c + 2*pi/T ), T/( 1.27^2 * 2 * 2 ) * ( 1 + cos( ( w - w_c )*T/2 ) ), w > ( w_c - 2*pi/T ), 0, -( w_c + 2*pi/T ) <= w <= -w_c/2, T/( 1.27^2 * 2 * 2 ) * ( 1 + cos( ( w + w_c )*T/2 ) ), w < ( w_c + 2*pi/T ), 0 );

fplot( G, 'b', 'Linewidth', 2 );

grid;

xlabel( '\omega' );

ylabel( 'G_S(\omega)' );

xlim([ -2.1*10^6 2.1*10^6 ]);

ylim([0 4.2005*10^-6]);

---