%function?cpfsk_signal?=?cpfsk（signal_num）
%signal_number為要產生的cpfsk信號的點數%
close?all
clear?all
clc

fc?=?1900;??%載波頻率
T?=?1/3600;?%碼元頻率為3600
fs?=?（1/T）*24;???%一個碼元周期采樣24點
ts?=?1/fs;
h?=?2/9;????%調制系數

%產生-1，1的隨機信號%
k_tmp?=?1;
In?=?2.^k_tmp?-?1;
In?=?[In-In];
pulse_num?=?100;??%-1，1的總個數為2×pulse_num個
In?=?repmat（In?1?pulse_num）;
Index?=?randperm（length（In））;
In?=?In（Index）;

phase0?=?0;?%載波相位，在相干解調情況下可設為0
cpfsk_signal?=?[];??%cpfsk信號
add_phase?=?[];???%附加相位
phase1?=?0;
time?=?ts:ts:T;?%一個碼元周期內采樣點數

for?i?=?1:length（In）
????q?=?time/（2*T）;
????if?i?~=?1
????????theta_k?=?pi*h*sum（In（1:i?-?1））;
????else
????????theta_k?=?0;
????end
????