t=(0:1/Fs:1); % ?%
n=1:N+1;
k=random('unif',0,n-1);
Phi=2*pi*k;
s2=(A*cos(2*pi*Fc*t)*exp(2*pi*(Fc-Fd)*t*i+Phi));
 
x=Signal2+randn(size(t));

s2=(A*cos(t*w1)*exp(t*w2*i+Phi));

clear all; 
close all;
clc;
 
G=6.67384e-20;
M=5.9726e24;
F=50; 
Fs=160;  %! 
Ts=1/Fs;
Fc=1.6e9;
Tc=1/Fc;
A=1; 
R=6371;
h=1400;
c=3e8; 
 
    Speed=(G*M/(R+h))^(1/2);
    Fd=Fs/(1-Speed/c);
    t=(0:1/Fs:1); 
N=c/F;
n=1:N+1;%1132356
k=random('norm',0,n-1); % 'norm' or 'unif'
Phi=2*pi*k;
w1=2*pi*Fc;
f=Fc-Fd;
w2=2*pi*(Fc-Fd);
s1=A*cos(w1*t);
Signal1=s1;
s2=A*cos(w1*t).*exp(i*w2*t+Phi);
Signal2=s2;
 
x=Signal2+randn(size(t));
Signal3=x;

n=1:N+1;%1132356
k=random('norm',0,n-1); % 'norm' or 'unif'

n=1:N+1;
k=random('unif',0,n-1);
Phi=2*pi*k;
s2=(A*cos(2*pi*Fc*t)*exp(2*pi*(Fc-Fd)*t*i+Phi));

clear all; 
close all;
clc;
 
G=6.67384e-20;
M=5.9726e24;
F=50; 
Fs=160;  %! 
Ts=1/Fs;
Fc=1.6e9;
Tc=1/Fc;
A=1; 
R=6371;
h=1400;
c=3e8; 
L=c*Ts;
 
    Speed=(G*M/(R+h))^(1/2);
    Fd=Fs/(1-Speed/c);
        
 
t=(0:1/Fs:1); 
 
k=random('norm',0,t); % 'norm' or 'unif'
 
Phi=2*pi*k;
w1=2*pi*Fc;
f=Fc-Fd;
w2=2*pi*(Fc-Fd);
 
 
s1=A*cos(w1*t);
Signal1=s1;
s2=A*cos(w1*t).*exp(i*w2*t+Phi);
Signal2=s2;
 
x=Signal2+randn(size(t));
Signal3=x;
 
 
NFFT = 2^nextpow2(length(s1)); 
FrequencyAxeFs=Fs/2*linspace(0,1,NFFT);
Spectr=(1/NFFT).*fft(Signal1, NFFT); 
SpectrSignalFs=abs(Spectr(1:NFFT));
PhaseSpectrFs=angle(fft(Signal3));
 
NFFT = 2^nextpow2(length(s1)); 
FrequencyAxeFd=Fd/2*linspace(0,1,NFFT);
Spectr=(1/NFFT).*fft(Signal3, NFFT); 
SpectrSignalFd=abs(Spectr(1:NFFT));
PhaseSpectrFd=angle(fft(Signal3));
 
 
plot(FrequencyAxeFs,SpectrSignalFs,'r') 
hold on
xlabel('Spectr of the signal');
ylabel('Frequency');
plot(SpectrSignalFd,FrequencyAxeFs);
xlabel('Spectr of the signal');
ylabel('Frequency');
 
 
figure
plot(SpectrSignalFs, fft(Signal3),'--or'), 
hold on
xlabel('PhaseSpectr of the signal');
ylabel('Frequency');
plot(FrequencyAxeFd, PhaseSpectrFd);
xlabel('PhaseSpectr of the signal');
ylabel('Frequency');

clear all;
close all;
clc;
 
G=6.67384e-11; % м^3/(кг*c^2)
M=5.9726e24; % кг
F = 50; % чатсота сигнала, Гц
Fs = 30*F; % чатсота дискретизации сигнала
Fc = 1.6e9; % несущая, Гц
Tc = 1/Fc; % период несущей
A = 1; % амплитуда
R = 6371e3; % метры
h = 1400e3; % метры
c = 3e8; % м/с
%L=c*Tc;
 
Speed=(G*M/(R+h))^(1/2); % скорость спутника
 
K = 1/(1-Speed/c); % коэф-т изменения частоты
Fcd = Fc*K;
Fsd = Fs*K;
 
 
t = 0:1/Fs:3/F; % масив времени
 
k = 0.03*randn(size(t));
Phi=2*pi*k;
w1=2*pi*F;
 
s1=A*cos(w1*t); % сигнал
Signal1=s1;
s2=A*cos(w1*t + Phi); % добавили скачки фазы
Signal2=s2;
x=Signal2 + 0.03*randn(size(t)); % добавили скачки амплитуды
Signal3=x;
 
% спектр исходного сигнала:
N = 2^nextpow2(length(s1));
fff = Fc + Fs/2*linspace(-1,1,N);
Spectr1 = (1/N)*fft(Signal1,N);
ASp1 = abs( fftshift(Spectr1) );
PhSp1 = angle( fftshift(Spectr1) );
 
% спектр сигнала с учетом эффекта:
ffd = Fcd + Fsd/2*linspace(-1,1,N);
SpectrD = (1/N).*fft(Signal3, N);
ASpD = abs( fftshift(SpectrD) );
PhSpD = angle( fftshift(SpectrD) );
 
plot(fff,ASp1,'r', ffd,ASpD,'b' );
title('Spectr of the signal');
xlabel('Frequency');
legend('исходный сигнал','c учетом эффекта')