这是我做的关于伪码捕获的仿真其中用了dmod 函数进行调制和解调  现在要在信道中加频偏  请问怎么加  不要信道衰减  不要延时  只要频偏

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%

%  直接序列扩频信号的产生与串行捕获模拟

%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

clear;

clc;

fPN =1023;                      % 伪码速率

fCarrier=fPN*4;

fSample =fCarrier*5 ;         % 载波的采样点频率

Ts=1/fSample;                 % 采样后符号持续时间

TxTemp1= fSample/ fPN;                      % 滤波器阶数

Delay=8;                     % 滤波器延迟

rolloff=0.25;

M=2;

chip = 1023;                    % 伪码长度，GOLD码10230位

dataLength = 10;               % 待发送数据位数

SNR=0;                      % 信噪比为3.5dB

FD=10000;                        % 设置多普勒频移

N=1000;                       % 蒙特卡罗仿真次数

num_of_det=0;                 % 检测次数

% 初始化随机数，每次不同

sourceData = randint(1, dataLength); % 0,1的数据源

sendDataLength = dataLength * chip; % 准备调制发送时的数据长度，目前一样，没有在数据前面加噪音

figure(1);

stem(sourceData);

title('数据源');

% 生成GOLD码

pnSequence =GoldGenerator();% GOLD码31位

pn_len=length(pnSequence);

figure(2);

stem(pnSequence(1:50));

title('10230位gold码');

% 扩频

spreadData = zeros(1, dataLength*chip);

temp = ones(1, chip);

for i = 1:dataLength

temp = ones(1, chip)*sourceData(i);

% 相加模2，异或

spreadData(((i-1)*chip+1):i*chip) = xor(temp, pnSequence);

end

figure(3);

stem(spreadData(1:50));

title('扩频后的数据');

[sendedMData,t] = dmod(spreadData, fCarrier, fPN, fSample, 'psk', 2);

sendedMDataLength = length(sendedMData);

figure(4);

plot(t(1:50),sendedMData(1:50));

title('经BPSK调制后的数据');

% 解调成数字信号

recv_da= ddemod(sendedMData, fCarrier, fPN, fSample, 'psk', 4);

recvDDataLen = length(recv_da);       % 捕获到的数据长度

figure(6);

stem(recv_da(1:50));

title('解调后的信号');

recvDData=2*recv_da-1;

RxGoldSequence=pnSequence;

RxGoldSequence1=2*pnSequence-1;

despreadData = zeros(1, dataLength);

RxDelayDelta = 1;

RxGoldSearchTimes = ceil(pn_len / RxDelayDelta);

RxSignalAmp=zeros(1,RxGoldSearchTimes);

amp=zeros(1,dataLength);

for i=1:RxGoldSearchTimes

amp=zeros(1,dataLength);

RxGoldSequence2=circshift(RxGoldSequence1',(i-1)*RxDelayDelta)';

for j=1:dataLength

recvDData1=recvDData((j-1)*chip+1:j*chip);

amp(j)=abs(sum(recvDData1.*RxGoldSequence2));

end

RxSignalAmp(i)=sum(amp);

end

RxSignalResAverage = mean(RxSignalAmp);

RxSignalResDeviation = std(RxSignalAmp);

RxSignalResMax=max(RxSignalAmp);

RxCatchThreshold = RxSignalResAverage + 3* RxSignalResDeviation;

figure(6);

stem(RxSignalAmp);

title('相关峰值');

%判断是否捕获

%设置门限

for i=1:RxGoldSearchTimes

if(RxSignalAmp(i)>=RxCatchThreshold)

%捕获成功

display('Success!!!');

c=i;

num_of_det=num_of_det+1;

break

else

display('fail');

end

end

Pn_sequence=circshift(RxGoldSequence',c*RxDelayDelta-1)';

for i=1:dataLength

De_spr(pn_len*(i-1)+1:pn_len*i)=xor(recv_da(pn_len*(i-1)+1:pn_len*i),Pn_sequence);

end

Rec_Data=downsample(De_spr,1023);

figure(7);

stem(Rec_Data);

title('恢复后的数据');