G=6.67384e-20;
M=5.9726e24;
f0=1.6e9;
R=6371;
c=3e8;
h=300:100:1500;
a=0:0.1309:3.14;
 
 
for i = 1:length(h)
    V=(G*M/(R+h(i)))^(1/2);
    for j = 1:length(a)
        t(j)=(a(j)*(R+h(i)))/V;
        fd(j,i)=((V*R*f0*cos(a(j))))/((R+h(i))*c);
   
           end;
end;
plot(t,fd)
legend(strcat('h=',num2str(h')))
xlabel('t,s')
ylabel('fd,kHz')

clear, clc
 
G=6.67384e-11;
M=5.9726e24;
f0=1.6e9;
R=6371;
c=3e8;
h=300:100:1500;
a=0:pi/500:pi;
t = zeros(length(a),length(h));
Fd = zeros(size(t));
for i = 1:length(h)
    V = sqrt(G*M/(R+h(i)));
    t(:,i) = a*(R+h(i))/V;
    Fd(:,i)= V*R*f0*cos(a) / ((R+h(i))*c);
end;
plot(t,Fd)
legend(strcat('h=',num2str(h')))
xlabel('t')
ylabel('fd')
 
dFdt = diff(Fd) ./ diff(t); % производная
tt = cumsum(diff(t)); % вспомогательный маив времени
 
figure
plot(tt,dFdt)
legend(strcat('h=',num2str(h')))
xlabel('t')
ylabel('dFdt')

clear all; 
close all;
clc;
format longe;
f=50; 
Fs=30*f; 
Ts=1/Fs; 
Fc=1.6e9;
A=1; 
 
alfa=0; 
c=3e8; 
 
Speed=c/f;
fd=((Speed*(Fc+f)*cos(alfa(j))))/c;
  
t=(0:1/Fs:1);
N = length(t);
for i=1:4;
    fi=(2*i-1)*(pi/4);
    s=A*cos(2*pi*(Fc+f)*t+fi);
    FF=fft(s);
    s2=ifft(FF);
    s3=real(s2)*exp(j*2*pi*(Fc+f)*t);
 
x=s3+randn(1,length(t));
end;
Spectr=fft(x); 
SpectrSignal=abs(Spectr);
plot(SpectrSignal)

clear all; 
close all;
clc;
format longe;
f=50; 
Fs=30*f; 
Ts=1/Fs; 
Fc=1.6e9;
A=1; 
 
alfa=0; 
c=3e8; 
 
Speed=c/f;
fd=((Speed*(Fc+f)*cos(alfa(j))))/c;
  
t=(0:1/Fs:1);
N = length(t);
for i=1:4;
    fi=(2*i-1)*(pi/4);
    s=A*cos(2*pi*(Fc+f)*t+fi);
    FF=fft(s);
    s2=ifft(FF);
    s3=real(s2)*exp(j*2*pi*(Fc+f)*t);
 
x=s3+randn(1,length(t));
end;
Spectr=fft(x); 
SpectrSignal=abs(Spectr);
plot(SpectrSignal)

clear all; 
close all;
clc;
format longe;
f=50; %частота сигнала%
Fs=30*f; %частота дискретизации%
Ts=1/Fs; %период дискретизации%
Fc=1.6e9;%частота несущая%
A=1; %амплитуда сигнала%
 
alfa=0; %угол%
c=3e8; %скорость света%
 
Speed=c/f;
fd=((Speed*(Fc+f)*cos(alfa(j))))/c;%доплеровский сдвиг %
  
t=(0:1/Fs:1);%время сигнала%
N = length(t); %длина сигнала%
for i=1:4;
    fi=(2*i-1)*(pi/4);% фаза%
    s1(j,i)=A*cos(2*pi*(Fc-Fd(j,i)+F)); %сигнал%
    s2(j,i)=(A*cos(2*pi*(Fc-Fd(j,i)+F)))*exp(2*pi*(Fc))%умноение на комплекксную составляющую%
    x(j,i)=s2+randn(size(t))%сигнал шум%
    %фурье преобразование и амплитудный спектр%
    NFFT = 2^nextpow2(L);
FrequencyAxe=Fs*linspace(0,1,NFFT);
Spectr=(1/NFFT).*fft(Signal3, NFFT); 
SpectrSignal=abs(Spectr(1:NFFT))
 
PhaseSpectr=angle*fft(Spectr,Fftl);%фазовый спектр%
end;
%графики спектров %
plot(SpectrSignal,w);
figure
plot(PhaseSpectr,w);

clear, clc
 
 
fs = 100; % частота дискретизации
ts = 1/fs; % период дискретизации
F = 1;
 
t = 0:ts:2;
s = sin(2*pi*F*t);
 
plot(t,s), hold on
 
v = 0.1; % скорость источника
c = 3; % скорось волны
 
fd = fs/(1-v/c) % частота дискретизации после эффекта
tsd = 1/fd; % новый период дискретизации
 
t1 = 0: tsd: tsd*(length(s)-1); % новый массив времени
 
plot(t1,s, 'r') % стоит заметть, что отсчеты сигнала остались без изменения
legend('исходный сигнал','эффект Доплера')

clear, clc
 
 
fs = 10; % частота дискретизации
ts = 1/fs; % период дискретизации
F = 1;
 
t = 0:ts:100;
s = sin(2*pi*F*t);
 
NFFT = 2^nextpow2(length(s)); 
N2 = NFFT/2+1;
 
S = fft(s,NFFT)/length(s);
 
fff = fs/2*linspace(0,1,N2);
Sp = abs(S(1:N2));
 
plot(fff, Sp), hold on
 
v = 0.5; % скорость источника
c = 3; % скорось волны
fd = fs/(1-v/c); % частота дискретизации после эффекта
ffd = fd/2*linspace(0,1,N2);;
 
plot(ffd, Sp,'r')
xlim([0, 5*F])
legend('Исходный сигнал','Допплер')