该楼层疑似违规已被系统折叠 隐藏此楼查看此楼

小女子 路遇强悍的毕业设计 来到贴吧 求大神们 出手相助 感激不尽！！

在此给出一点程序 看可不可以帮助到大家

线性调频信号的产生程序

T=10e-6; %pulse duration10us

B=30e6; %chirp frequency modulation bandwidth 30MHz

K=B/T; %chirp slope

Fs=2*B;Ts=1/Fs; %sampling frequency and sample spacing

N=T/Ts;

t=linspace(-T/2,T/2,N);

St=exp(j*pi*K*t.^2); %generate chirp signal

subplot(211)

plot(t*1e6,real(St));

xlabel('Time in u sec');

title('Real part of chirp signal');

grid on;axis tight;

subplot(212)

freq=linspace(-Fs/2,Fs/2,N);

plot(freq*1e-6,fftshift(abs(fft(St))));

xlabel('Frequency in MHz');

title('Magnitude spectrum of chirp signal');

grid on;axis tight;

线性调频信号的压缩

T=10e-6; %pulse duration10us

B=30e6; %chirp frequency modulation bandwidth 30MHz

K=B/T; %chirp slope

Fs=10*B;Ts=1/Fs; %sampling frequency and sample spacing

N=T/Ts;

t=linspace(-T/2,T/2,N);

St=exp(j*pi*K*t.^2); %chirp signal

Ht=exp(-j*pi*K*t.^2); %matched filter

Sot=conv(St,Ht); %chirp signal after matched filter

subplot(211)

L=2*N-1;

t1=linspace(-T,T,L);

Z=abs(Sot);Z=Z/max(Z); %normalize

Z=20*log10(Z+1e-6);

Z1=abs(sinc(B.*t1)); %sinc function

Z1=20*log10(Z1+1e-6);

t1=t1*B;

plot(t1,Z,t1,Z1,'r.');

axis([-15,15,-50,inf]);grid on;

legend('emulational','sinc');

xlabel('Time in sec \times\itB');

ylabel('Amplitude,dB');

title('Chirp signal after matched filter');

subplot(212) %zoom

N0=3*Fs/B;

t2=-N0*Ts:Ts:N0*Ts;

t2=B*t2;

plot(t2,Z(N-N0:N+N0),t2,Z1(N-N0:N+N0),'r.');

axis([-inf,inf,-50,inf]);grid on;

set(gca,'Ytick',[-13.4,-4,0],'Xtick',[-3,-2,-1,-0.5,0,0.5,1,2,3]);

xlabel('Time in sec \times\itB');

ylabel('Amplitude,dB');

title('Chirp signal after matched filter (Zoom)');

雷达的目标跟踪源程序

clc;

clear;

T=2;%雷达扫描周期

num=100;%滤波次数

%*******************************产生真实轨迹**********************************

N=800/T;

x=zeros(N,1);y=zeros(N,1);

vx=zeros(N,1);vy=zeros(N,1);

x(1)=-2000;

y(1)=1000;

vx=15;vy=0;

ax=0;ay=0;

var=100;

%产生真实轨迹

for i=1:N-1

x(i+1)=x(i)+vx*T+0.5*ax*T^2;

y(i+1)=y(i)+vy*T+0.5*ay*T^2;

end

nx=zeros(N,1);ny=zeros(N,1);

nx=100*randn(N,1);

ny=100*randn(N,1);

zx=x+nx;zy=y+ny;

%滤波50次

for m=1:num

z=2:1;

xks(1)=zx(1);

yks(1)=zy(1);

xks(2)=zx(2);

yks(2)=zy(2);

o=4:4;g=4:2;h=2:4;q=2:2;xk=4:1;perr=4:4;

o=[1,T,0,0;0,1,0,0;0,0,1,T;0,0,0,1];

h=[1 0 0 0;0 0 1 0];

g=[T/2,0;T/2,0;0,T/2;0,T/2];

q=[10000 0;0 10000];

perr=[var^2 var^2/T 0 0

var*var/T 2*var^2/(T^2) 0 0

0 0 var^2 var^2/T

0 0 var^2/T 2*var^2/(T^2)];

vx=(zx(2)-zx(1))/2;

vy=(zy(2)-zy(1))/2;

xk=[zx(1);vx;zy(1);vy];

%卡尔曼滤波开始

for r=3:N;

z=[zx(r);zy(r)];

xk1=o*xk;

perr1=o*perr*o';

k=perr1*h'*inv(h*perr1*h'+q);

xk=xk1+k*(z-h*xk1);

perr=(eye(4)-k*h)*perr1;

xks(r)=xk(1,1);

yks(r)=xk(3,1);

vkxs(r)=xk(2,1);

vkys(r)=xk(4,1);

xk1s(r)=xk1(1,1);

yk1s(r)=xk1(3,1);

perr11(r)=perr(1,1);

perr22(r)=perr(2,2);

rex(m,r)=xks(r);

rey(m,r)=yks(r);

end %结束一次滤波

end

ex=0;ey=0;

eqx=0;eqy=0;

ey1=0;

ex1=N:1;ey1=N:1;

%计算滤波的均值

%计算滤波误差的均值

for i=1:N

for j=1:num

ex=ex+x(i)-rex(j,i);

ey=ey+y(i)-rey(j,i);

end

ex1(i)=ex/num;

ey1(i)=ey/num;

ex=0;eqx=0;ey=0;eqy-0;

end

%绘图

figure(1);

plot(x,y,'k-',zx,zy,'b:',xks,yks,'r-.');

legend('真实轨迹','观测样本','估计轨迹')'

figure(2);

plot(ey1);

legend('x方向平均误差');