xxtpcaiyang.m

来自「功率谱分析是信号处理的重要内容」· M 代码 · 共 24 行

clear all;
clc;
%%%%%%%%%%%%%%%%%%%
%设置系统参数
%%%%%%%%%%%%%%%%%%%
Tp=300e-7;            %脉冲宽度
B=2e6;                %带宽
kr=B/Tp;              %调频率
fs=10e6;              %采样频率
R0=180e3;
theta=0;
c=3e8;
t=-Tp/2:1/fs:Tp/2;
echo=exp(j*0.5*kr*2*pi*t.^2);
plot(t,real(echo));
Fs=10e6;Ts=1/Fs;
n=-150:150;
echo=exp(j*0.5*kr*2*pi*(n*Ts).^2);
figure,stem(n,real(echo));
k=0:1024;
S=echo*exp(-j*2*pi/1024).^(n'*k);
w=2*pi*k*4e4/1024;
figure,plot(w,fftshift(abs(S)));