程序简介：用来产生高动态的信号的多普勒频移（用GPS的频点1575.42MHZ）
仿真时间为8s，时间间隔为1ms。加加速度为100g，加速度为25g。两段仿真
时间为3s~3.5s和5.5s~6s。
function [fed]=freq_high_Dynamics
t=8; %·ÂÕæʱ¼ä(s)
g=9.8;
n=8000;
a1=zeros(1,n+1);%¼Ó¼ÓËÙ¶È
a1(1,3001:3500)=100*g*ones(1,500);
a1(1,5501:6000)=-100*g*ones(1,500);
a=zeros(1,n+1);%¼ÓËÙ¶È
a(1,1)=-25*g;
for i=1:8000
a(i+1)=a1(i)/1000+a(i);
end
v=zeros(1,n+1);%½ÓÊÕ»úËÙ¶È
v(1)=-100;
for i=1:8000
v(i+1)=a(i)/1000+v(i);
end
fed=zeros(1,n+1);%¶àÆÕÀÕƵÂÊÆ«ÒÆ
fed=v*1575.42e+6/3e+8;
figure
plot(fed)
title('¶àÆÕÀÕƵÒÆ');
xlabel('ʱ¼ä/s');
ylabel('¶àÆÕÀÕƵÒÆ/Hz');
grid on;
figure
p1=2*pi*fed.*[0:0.001:8];%²»Í¬Ê±¿Ì¶àÆÕÀÕÏàλƫÒÆ
plot(p1)
title('²»Í¬Ê±¿Ì¶àÆÕÀÕÏàλƫÒÆ');
xlabel('ʱ¼ä/s');
ylabel('¶àÆÕÀÕÏàλƫÒÆ');
grid on;
p=zeros(1,8001);
for i=1:8001
p(i+1)=p(i)+2*pi*fed(i)*0.001;%¶àÆÕÀÕÏàλƵÒÆÀÛ¼ÓÖµ
end
figure
plot(p);
title('¶àÆÕÀÕÏàλƫÒÆÀÛ¼ÓÖµ');
xlabel('ʱ¼ä/s');
ylabel('ÏàλƫÒÆÀÛ¼ÓÖµ');
