fangzhen.m

基于MATLAB的线性调频信号产生以及频谱分析

clear all;
%定义参数
B=10e6;       % B ——带宽
Fc=10e9;     % Fc——载频
Fs=10*B;     % Fs——采样率?
T2=10e-6;   % T1——脉冲宽度
Tr=375e-6;   % Tr——PRT
C=3e8;	

%产生线性调频脉冲信号
mu=B/T2;
NT2=[Fs*T2];               %采样点数
t=[-NT2/2:NT2/2-1]*T2/NT2; %取一个脉冲
fai=Fc+mu*t.*t/2;          %相位 
singal=exp(2*pi*j*fai);    %LFM信号，幅度为1
figure(1);
subplot(2,1,1)
plot(t*1e6,real(singal));
title('近距离时线性调频信号');
xlabel('延时 - us');
ylabel('归一化幅度');
grid on;

%线性调频信号频谱
singalfft=fftshift(fft(singal));
Nf=length(singalfft);
f=[-Nf/2:Nf/2-1]*Fs/Nf;
subplot(2,1,2)
plot(f*1e-6,abs(singalfft)/max(abs(singalfft)));
title('近距离时的信号频谱图');
xlabel('频率f/MHz');
ylabel('归一化幅度');
grid on;

%模糊函数
Taw=100;
i=0;
for t=[-Taw:Taw-1]*T2/Taw
    i=i+1;
    j=0;
    for Fd=[-Taw:Taw-1]*B/Taw;
        j=j+1;
        val1=1-abs(t)/T2;
        val2=pi*T2*val1;
        val3=Fd+mu.*t+eps;
        val4=val2*val3;
        ambg(j,i)=abs(val1*(sin(val4)/(val4))).^2;
    end
end
t=[-Taw:Taw-1]*T2/Taw;
Fd=[-Taw:Taw-1]*B/Taw;
figure(2);
mesh(t*1e6,Fd*1e-6,ambg);
title('近距离时的模糊函数图');
xlabel('延时-us');
ylabel('多普普勒频移fd/MHz');
zlabel('归一化幅度');
grid on;
%绘图
%模糊函数等高线图
figure(3);
contour(t*1e6,Fd*1e-6,ambg);
title('近距离时模糊函数等高线图');
xlabel('延时-us');
ylabel('多普普勒频移fd/MHz');
grid on;

%自相关函数
figure(4);
subplot(2,1,1);
plot(t*1e6,abs(ambg(:,Taw+1)));
title('近距离时自相关函数');
xlabel('延时-us');
ylabel('归一化幅度');
grid on;

%多普勒敏感度
subplot(2,1,2)
plot(Fd*1e-6,abs(ambg(Taw+1,:)));
title('近距离时多普勒敏感度');
xlabel('多普勒频移-MHz');
ylabel('归一化幅度');
grid on;