MATLAB 产生带宽信号

clc
clear all
close all
%%%  产生中心频率150Mhz ，带宽50Mhz, 采样率252Mhz的带宽信号
Tp    = 65.015e-6;            % 发射时宽. 控制点数.比如要输出的点数是L， Tp = L /Fs .
                              % 16384 /  252e6 = 65.015e-6                              
B     = 50e6;                 % 发射带宽
K     = B/Tp;                 % 发射调频率
Fs    = 252e6;             % 发射信号采样率

t     = [-Tp/2:1/Fs:Tp/2];    % 发射采样时刻
 sig   = exp(i*pi*K*t.^2);     % 发射信号IQ两路
 
figure;
real_sig_tmp = real (sig);

real_sig = floor (16384 * real_sig_tmp);
L = length (t);
figure;
subplot(211),plot([-Fs/2:Fs/L:Fs/2-Fs/L],db(fftshift(fft(real_sig))));
title('发射频谱')

subplot(212), plot(real_sig);
title('发射信号I部')



fs1=252e6;%采样频率
f1=102e6;%信号频率Hz ,输出中心频率150M ，需要镜像输出，252 - 150 = 102
N=16384;%采样点数
t1=(0:N-1)/fs1;%采样时间s，fs的值越大，出来的波形失真越小
sin1=sin(2*pi*f1*t1);%信号采样值

real_sig_duc = real_sig.*sin1(1,:);               %上变频输出
figure;
subplot(211),plot([-Fs/2:Fs/L:Fs/2-Fs/L],db(fftshift(fft(real_sig_duc))));
title('发射频谱')

subplot(212), plot(real_sig_duc);
title('发射信号I部')


fid = fopen('real_sig_duc.txt','w');
fprintf (fid,'%0.0f,\r',real_sig_duc);  %%0.14f 保留小数点后几位，存小数的时候去掉科学计数法，很有用
fclose (fid);