A=10;   %Амплитуда сигнала%
f=50;   %Частота сигнала%
phi=pi/4; %Начальная фаза сигнала%
t=0:0.001:0.04; %Задание временной оси%
y=A*sin(2*pi*f*t-phi);
%Построение графика%
plot(t,y,'ro-'), grid
xlabel('t, c  \rightarrow')
ylabel('y(t)  \rightarrow')
title('График гармоники')