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個文件信息