Xlabel('Timpul in secunde'); hold off;

Xd1_f = fft(xd1_f);

in = 0:0.1:N-1;

Xdin1_f = zeros(1,length(in));

for l = 1:length(in)

Xdin1_f(l) = sum(xd1_f.*exp(-j*2*pi*in(l)*(0:N-1)/N));

end;

subplot(212), stem((0:N_2-1)*fs/N, abs(Xd1_f(1:N_2)/N));

hold on;

plot(in(in1)*fs/N, abs(Xdin1_f(in1)/N));

title(['Spectrul de amplitudine a semnalului xd1 cu efect ',text_1,' (f1 = ',... num2str(f1), 'Hz)']);

Xlabel('Frecventa in Hz');

ylabel('Modulul (DFT/N)');

2.6 Сформировать сигнал в виде прямоугольного одиночного импульса с амплитудой и шириной . Зададим шаг дисретизации , а длительность анализа .

Ts=0.01;T=1;

A=0.8;

w=0.5;N=T/Ts;t=0:Ts:T;

y=A*rectpuls(t,w);

plot(t(1:100),y(1:100)); grid

Применим к вектору процедуру и построим график зависимости модуля результата от частоты.

x=fft(y)/N;df=1/T;Fm=1/Ts;

f=0:df:Fm;a=abs(x);

stem(f,a);grid

График спектра на полученном рисунке представлен в интервале частот от 0 до . Для построения графика модуля спектра в привычной форме, в диапазоне частот от до используем команды

xp=fftshift(x);f1=-Fm/2:df:Fm/2;

a=abs(xp);stem(f1,a),grid;xlabel('Frecventa, Hz');

ylabel('Modulul')

Построим графики действительной и мнимой частей Фурье – изображения прямоугольного импульса:

Re=real(xp);Im=imag(xp);

plot(f1,Re,f1,Im),grid

Xlabel('Frecventa, Hz')

ylabel('Partea reala si imaginara')

Выводы: В лабораторной работе №2 были изучены методы спектрального анализа и синтеза сигналов а так же спектральный анализ сигналов с помощью системы MATLAB. Были выполнены все требования по расчетам и через программу были выявлены графики зависимостей. MATLAB является системой математического моделирования различных процессов и систем, включая, в том числе и хорошо развитые средства спектрального анализа сигналов.