cpfsk信号的产生

%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