Еремеева спвтлаб1

Задание 1.

clear; clc;

syms t;

st = 4* cos(pi* t+pi/4);

a0=(1/(2*pi))* int(st,t,-pi,pi);

FourierSum=a0;

double(a0)

for n=1:5

a=double((1/pi)* int(st* cos(n* t),t,-pi,pi));

b=double((1/pi)* int(st* sin(n* t),t,-pi,pi));

result=sprintf('%f * cos(%d*t) + %f * sin(%d*t)',a,n,b,n);

disp(result)

end

ans = -0.1233

0.274437 * cos(1*t) + -0.087356 * sin(1*t)

-0.414704 * cos(2*t) + 0.264009 * sin(2*t)

2.799148 * cos(3*t) + -2.672989 * sin(3*t)

0.397063 * cos(4*t) + -0.505556 * sin(4*t)

-0.160878 * cos(5*t) + 0.256046 * sin(5*t)

clear, clc, clf

N = 2; % длина последовательности

z = [1,1,1,1,-1,-1,-1,-1,1,1,1,1,-1,-1,-1,-1,1,1,1,1,] % случайные данные

fc = 4; % часота несущего колебания

fz = 1; % частота символов

fd = 10*fc; % чатоста дискретизации

td = 1/fd; % период дискретизации

t = 0:td:N/fz-td; % массив времени

% Расширяем сигнал.

% Вариант №1:

M = fd/fz; % количество отсчетов на бит

Z = repmat(z,M,1); % повторяем последовательность по строкам

Z = Z(:)'; % вытягиваем по столбцам

% Вариант №2:

Z = z( fix(t*fz)+1 );

S = sin(2*pi*fc*t + Z*pi); % сигнал с БФМ

plot(t,Z,'linew',2)

hold on

grid on

plot(t,S,'r','linew',2)

z = 1 1 1 1 -1 -1 -1 -1 1 1 1 1 -1 -1 -1 -1 1 1 1 1