|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
if offtime<600
error('仿真时间必须大于600s,请重新输入');
end
o=zeros(1,offtime);
m=zeros(1,offtime);
n=zeros(1,offtime);
O=zeros(1,ceil(offtime/Ts));
M=zeros(1,ceil(offtime/Ts));
N=zeros(1,ceil(offtime/Ts));
%起始点坐标
o0=0;
m0=0;
n0=0;
a=0.05;c=0.2;
n(1)=1;
w=c*randn(1,offtime);
for k=1:offtime
o(k+1)=exp(m(k)*Ts);
m(k+1)=exp(n(k)*Ts);
n(k+1)=exp(-a*Ts)*n(k)+w(k);
end
% 得到雷达的观测数据
for n=0:Ts:offtime
N(n/Ts+1)=n(n+1);
O(n/Ts+1)=o(n+1);
M(n/Ts+1)=m(n+1);
end
%显示真实轨迹
plot(N,M,O,'LineWidth',2),grid on;
('目标真实航迹');
[ 本帖最后由 80522957 于 2007-5-28 19:02 编辑 ] |
|