Контрольная работа по "Программированию"

1. Построим  временные диаграммы, рассчитаем и построить в масштабе амплитудные

спектры :

1.  Листинг программы на языке Matlab:

>> t = 1:0.015:10; U0=15; T=3; phi=0;

>> s=U0*sin(2*pi*T*t+phi);

>> plot(t, square(t,20))

Временная диаграмма имеет вид:

[pic 1]

Рисунок 1

Рассчитаем и построим амплитудные спектры сигнала:

% Файл - сценарий для построения графиков спектров

 T = 3; % период сигнала

w0 = 2*pi/T; % основная частота

n1=1:10; % положительные номера гармоник

cn=-4*j./n1/pi.*sin(pi*n1/6).*sin(n1*pi/2).*exp(-j*n1*pi/3);

n2=-15:-1;1 % отрицательные номера гармоник

c_n=-4*j./n2/pi.*sin(pi*n2/6).*sin(n2*pi/2).*exp(-j*n2*pi/3);

Xn=[c_n 0 cn]; % вектор коэффициентов Фурье для n>0, n=0, n<0

n=-10:10; % вектор номеров гармоник

figure(1), subplot(121),stem(n*w0,abs(Xn))

Амплитудный спектр имеет вид:

[pic 2]

Рисунок 2

2. Листинг программы на языке Matlab:

>> t = 1:0.0075:10; U0=15; T=3; phi=0;

>> s=U0*sin(2*pi*T*t+phi);

>> plot(t, square(t,20))

Временная диаграмма имеет вид:

[pic 3]

Рисунок 3

Рассчитаем и построим амплитудные спектры сигнала:

% Файл - сценарий для построения графиков спектров

 T = 3; % период сигнала

w0 = 2*pi/T; % основная частота

n1=1:10; % положительные номера гармоник

cn=-4*j./n1/pi.*sin(pi*n1/6).*sin(n1*pi/2).*exp(-j*n1*pi/3);

n2=-15:-1;1 % отрицательные номера гармоник

c_n=-4*j./n2/pi.*sin(pi*n2/6).*sin(n2*pi/2).*exp(-j*n2*pi/3);

Xn=[c_n 0 cn]; % вектор коэффициентов Фурье для n>0, n=0, n<0

n=-10:10; % вектор номеров гармоник

figure(1), subplot(121),stem(n*w0,abs(Xn))

[pic 4]

Рисунок 4

3.  Листинг программы на языке Matlab:

>> t = 1:0.015:10; U0=15; T=6; phi=0;

>> s=U0*sin(2*pi*T*t+phi);

>> plot(t, square(t,20))

[pic 5]

Рисунок 5

Рассчитаем и построим амплитудные спектры сигнала:

% Файл - сценарий для построения графиков спектров

 T = 6; % период сигнала

w0 = 2*pi/T; % основная частота

n1=1:10; % положительные номера гармоник

cn=-4*j./n1/pi.*sin(pi*n1/6).*sin(n1*pi/2).*exp(-j*n1*pi/3);

n2=-15:-1;1 % отрицательные номера гармоник

c_n=-4*j./n2/pi.*sin(pi*n2/6).*sin(n2*pi/2).*exp(-j*n2*pi/3);

Xn=[c_n 0 cn]; % вектор коэффициентов Фурье для n>0, n=0, n<0

n=-15:15; % вектор номеров гармоник

figure(1), subplot(121),stem(n*w0,abs(Xn))

[pic 6]

Рисунок 6

4. Листинг программы на языке Matlab:

>> t = 1:0.015:10; U0=15; T=1; phi=0;

>> s=U0*sin(2*pi*T*t+phi);

>> plot(t, square(t,20))

[pic 7]

Рисунок 7

Рассчитаем и построим амплитудные спектры сигнала:

% Файл - сценарий для построения графиков спектров

T = 1; % период сигнала

w0 = 2*pi/T; % основная частота

n1=1:10; % положительные номера гармоник

cn=-4*j./n1/pi.*sin(pi*n1/6).*sin(n1*pi/2).*exp(-j*n1*pi/3);

n2=-15:-1;1 % отрицательные номера гармоник

c_n=-4*j./n2/pi.*sin(pi*n2/6).*sin(n2*pi/2).*exp(-j*n2*pi/3);

Xn=[c_n 0 cn]; % вектор коэффициентов Фурье для n>0, n=0, n<0

n=-15:15; % вектор номеров гармоник

figure(1), subplot(121),stem(n*w0,abs(Xn))