close all;
clear all;

A=10;
f=100;
fi=-40;

rad=2*pi*fi/360;
fvz=10000;

osx=linspace(0,1/f,fvz/f);
spektrum=[0 zeros(1,f-1) exp(1i*rad) zeros(1,fvz-2*f-1) exp(-1i*rad) zeros(1,f-1)];
cas=A*fvz/2*ifft(spektrum);

figure;
plot(osx,cas(1:fvz/f));