Teoria transmisiunii informației

Scopul lucrării: Analiza şi sinteza semnalelor periodice şi aperiodice.

Mersul lucrării:

1.Să se efectueze calculul spectrului discret a succesiunii periodice de impulsuri dreptunghiulare cu ajutorul programului fourier2.

>> delta_tau = 0.1; % durata relativa a impulsurilor

>> imax = 50; % numarul armonicilor calculate

>> flag = 1;

>> d_ik = imax/200; % pasul de esantionare

>> i = 0:imax; % diapazonul de armonici

>> ik = 0:d_ik:imax; % formarea axei X - a numarului de esantioane

>> Aik = 2*(delta_tau)*(sin(ik*pi*delta_tau+eps)./(ik*pi*delta_tau+eps)); % Amplitudinea semnalului

>> phik = -(pi/2)*(1-sign(Aik)); % Faza semnalului

>> Ai = 2*(delta_tau)*(sin(i*pi*delta_tau+eps)./(i*pi*delta_tau+eps)); % Spectrul de amplitudine

>> phi = -(pi/2)*(1-sign(Ai)); % Spectrul de faza

>> Aik = abs(Aik);

>> Aik(1) = Aik(1)/2;

>> Ai = abs(Ai);

>> Ai(1) = Ai(1)/2;

>> if flag == 1

>> figure(1); clf;

>> subplot(211), plot(ik, Aik); hold on;

>> stem(i, Ai); grid;

>> title('Spectrul de amplitudine'); xlabel('Numarul de armonici')

>> hold off

>> subplot(212), plot(ik, phik); hold on;

>> stem(i, phi); grid;

>> title('Spectru de faza'); ylabel('rad');

>> xlabel(['Numarul de armonici (relative Pulsdauer =',

num2str(delta_tau),' )']);

>> hold off;

>> end; Pentru τ0 = 0.1, N= 50.

Pentru τ0 = 0.2, N= 50. Pentru τ0 = 0.5, N= 50.

Rezultatele calcului spectrelor pentru trei valori 0: 0.1; 0.05; 0.01, obținute cu ajutorul programului fourier3.

>> imax = 50; % Numarul maxim de armonic

>> figure(1); clf;

>> delta_tau = 0.1;

>> d_ik = imax/200; % pasul de esantionare

>> i = 0:imax; % diapazonul de armonici

>> ik = 0:d_ik:imax; % formarea axei X - a numarului de esantioane

>> Aik = 2*(delta_tau)*(sin(ik*pi*delta_tau+eps)./(ik*pi*delta_tau+eps)); % >> Amplitudinea semnalului

>> phik = -(pi/2)*(1-sign(Aik)); % Faza semnalului

>> Ai = 2*(delta_tau)*(sin(i*pi*delta_tau+eps)./(i*pi*delta_tau+eps)); % Spectrul de amplitudine

>> phi = -(pi/2)*(1-sign(Ai)); % Spectrul de faza

>> Aik = abs(Aik);

>> Aik(1) = Aik(1)/2;

>> Ai = abs(Ai);

>> Ai(1) = Ai(1)/2;

>> subplot(311), stem(i, Ai); hold on;

>> plot(ik, Aik);

>> title(['Spectrul de amplitudine (delta-tau = ', num2str(delta_tau),' )']);

>> hold off;

>> p = get(gca,'Position');

>> set(gca,'Position',[p(1),p(2),p(3),p(4)*0.9]);