P1-P4相位编码信号生成

包括了雷达信号P1，P2，P3，P4的生成和时频图的绘制。MATLAB版本

clc;
clear;
NewFolder_Jpg=[pwd‘\P1_CWD‘‘\‘];
mkdir（NewFolder_Jpg）;%新建文件夹
%%%%%%%%%利用Frank码对相位进行编码
% N=8;
% %Creating the phase matrix
%        for k=1:N
%         for n=1:N
%             phi（kn）= （-pi/N）*（N-（2*k-1））*（（k-1）*N+（n-1））;
%         end
%        end
%        
% %Forming matrix vector    
% nn=0;
%     for n=1:N
%         for k=1:N
%             nn=nn+1;
%             phi2（nn）=phi（nk）;
%         end
%     end
p1_4=load（‘P1_4.mat‘）;
p1_5=load（‘P1_5.mat‘）;
p1_6=load（‘P1_6.mat‘）;
p1_7=load（‘P1_7.mat‘）;
p1_8=load（‘P1_8.mat‘）;
P1_sum=table（）;
h=[];
for i=1:100
position=randi（511）;
switch position
    case 1
        P1_select=p1_4.phi2;
    case 2
        P1_select=p1_5.phi2;
    case 3
        P1_select=p1_6.phi2;
    case 4
        P1_select=p1_7.phi2;
    case 5
        P1_select=p1_8.phi2;
end
num_p1=randi（411）+1;
XHCD=num_p1*numel（P1_select）;%序列时间长度是一定的
MYSLpsk=100;
% fcpsk=（（randi（71）*1e5+randi（91）*1e4）/1.2）/1000; %载波频率
fy=2e3;%采样频率
fcpsk=randi（511）+3;%[1/8~1/4]
fcpsk=fy*（1/fcpsk）;

% n=1:XHCD;
% a=1:numel（frank_select）;
% p1 = load（‘p1.mat‘）;
% code13=[1111100110101];%13位barker码
% seta=pi/numel（a）;%分配相位
% a=randi（MpskXHCD1）-1;%编码后的随机序列
% sh=code13.*seta;
% sh=（a-1）.^2*seta-pi*（a-1）;
% sh=（a-1）.^2*seta;
% frank5 = load（‘frank8.dat‘）;
sh=P1_select;
for j=1:num_p1
    h=[hsh];
end
clear t;
clear M;
clear m;
t =1/fy:1/fy:（XHCD*（1/MYSLpsk））;%码元长度*每个码元的持续时间=总时间
m=zeros（1length（t））;
time_duration=round（length（t）/XHCD）;
for j=1:XHCD
       m（（j-1）*time_duration+1:（j）*time_duration）=h（j）;
end
m=m（1:length（t））;

spsk=exp（1i*（2*pi*（fcpsk）*t+m））;%调制后的信号mpsk其中每个码元的相位叠加不一样
st=round（t*fy）;
[tfrtf] = tfrcw（spsk‘st512）;
frame=imagesc（tfabs（tfr））;
set （gca‘position‘[0011]）
name=[NewFolder_Jpgnum2str（i）‘.jpg‘];
saveas（ frame name）;%生成原始图片数据
% n=numel（spsk）;
% fre=（fy/2）/（n/2）:（fy/2）/（n/2）:（fy/2）;
% [IF Te tfr] = SET_Y（spsk‘40）;
% frame=imagesc（tfreabs（Te））;
% set （gca‘position‘[0011]）
% name=[NewFolder_Jpgnum2str（i）‘.jpg‘];
% saveas（ frame name）;%生成原始图片数据
%  spsk=real（hilbert（spsk））;
% spsk=spsk（1:fy/4）;
% xn1=（spsk）;
% xr=real（xn1）;
% xi=imag（xn1）;
% xh1=unwrap（atan2（xixr））;
% fs=fy;
% xhd1=fs*diff（xh1）/（2*pi）;
% % subplot（224）;
% figure;
% plot（xhd1）;
% xlabel（‘采样点数n‘）;ylabel（‘瞬时频率‘）;
% P1_sum.（i）={spsk};
disp（i）
end
% save（‘P1‘ ‘P1_sum‘‘-v7.3‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2017-12-02 00:47  P1-P4\
     文件        2551  2017-11-30 22:06  P1-P4\P1.m
     文件         271  2017-11-29 09:37  P1-P4\P1_4.mat
     文件         299  2017-11-29 09:37  P1-P4\P1_5.mat
     文件         414  2017-11-29 09:37  P1-P4\P1_6.mat
     文件         426  2017-11-29 09:37  P1-P4\P1_7.mat
     文件         606  2017-11-29 09:38  P1-P4\P1_8.mat
     目录           0  2017-12-02 00:47  P1-P4\P1_CWD\
     文件       21342  2017-11-30 22:07  P1-P4\P1_CWD\1.jpg
     文件       20874  2017-11-30 22:07  P1-P4\P1_CWD\2.jpg
     文件       20969  2017-11-30 22:07  P1-P4\P1_CWD\3.jpg
     文件       16785  2017-12-02 00:18  P1-P4\P1_CWD\P1.jpg
     文件        2561  2017-11-30 21:42  P1-P4\P2.m
     文件         226  2017-11-29 09:53  P1-P4\P2_4.mat
     文件         234  2017-11-29 09:53  P1-P4\P2_5.mat
     文件         296  2017-11-29 09:54  P1-P4\P2_6.mat
     文件         295  2017-11-29 09:54  P1-P4\P2_7.mat
     文件         376  2017-11-29 09:54  P1-P4\P2_8.mat
     目录           0  2017-12-02 00:48  P1-P4\P2_cwd\
     文件       21306  2017-11-30 21:42  P1-P4\P2_cwd\1.jpg
     文件       26156  2017-11-30 21:46  P1-P4\P2_cwd\100.jpg
     文件       21306  2017-11-30 21:45  P1-P4\P2_cwd\98.jpg
     文件       22196  2017-11-30 21:46  P1-P4\P2_cwd\99.jpg
     文件       24638  2017-12-02 00:19  P1-P4\P2_cwd\P2.jpg
     文件        1823  2017-11-30 22:27  P1-P4\P3.m
     目录           0  2017-12-02 00:48  P1-P4\P3_cwd\
     文件       20552  2017-11-30 22:27  P1-P4\P3_cwd\1.jpg
     文件       19936  2017-11-30 22:32  P1-P4\P3_cwd\100.jpg
     文件       20124  2017-11-30 22:32  P1-P4\P3_cwd\97.jpg
     文件       20170  2017-11-30 22:32  P1-P4\P3_cwd\98.jpg
     文件       19726  2017-11-30 22:32  P1-P4\P3_cwd\99.jpg
............此处省略9个文件信息