%function y = PhaseMod(OFDM_signal,number_of_bits)
%Input:OFDM signal,number of bits of OFDM signal,hp:mudulation of PM system 
%      ,Ac:Amplitude of OFDM signal,fc:Phase modulation carrier frequency
%Output:CE-OFDM-PM signal

function [CE_OFDM_PM_signal]=Phase_MOD(OFDM_signal,bits_per_ch,FFT_size,ts)
%代剛ノ計沮
%bits_per_ch=54;
%ts=.1
%FFT_size=16

hp=0.8
Ac=sqrt(2/FFT_size)
fc=1e+006
CE_OFDM_PM_signal=[]

for t=1:bits_per_ch*2
    CE_OFDM_PM_signal(t)=Ac*cos(2*pi*fc*((t+9)*ts)+hp*OFDM_signal(t))
end