时域脉冲通信采用高斯脉冲且使用PAM调制的Matlab简易演示仿真

环境

matlab 2016a

指标

1 将声音信号转为二进制码
2 PAM调制
3 采用高斯脉冲

流程

在这里插入图片描述

代码

[OriginVoice,fs]=audioread('voice.m4a') ; 
OriginVoice=OriginVoice(:,2);
N=length(OriginVoice); % 计算信号x的长度
t=0:1/fs:(N-1)/fs; % 计算时间范围，样本数除以采样频率
OriginVoice=OriginVoice';
% 
% % sound(y,fs);
figure(1);
plot(t,OriginVoice)
OriginVoice_max=max(OriginVoice);
OriginVoice1=OriginVoice+OriginVoice_max;
OriginVoice1=OriginVoice1*100000;
OriginVoiceBit= dec2bin(OriginVoice1);		%将十进制数据转为字符
AA1 = boolean(OriginVoiceBit-'0');	%将字符转成逻辑量
AA1 = double(AA1);			%强制转为double型
AA = AA1.';
AA = AA(:);
OriginVoiceBitSend = AA.';



t=-5:0.1:5;
Lenth_t=length(t);
tao=1;
t0=0;

Info=OriginVoiceBitSend;
InfoBitNum=length(OriginVoiceBitSend);
MSign=zeros(InfoBitNum,Lenth_t);
% y0=zeros(Lenth_t,1);
y0=(-0.5).*(t.^2-tao.^2).*exp(-1.*((t-t0)/tao).^2)./(((2*pi).^(1/2)).*(tao.^5));
%%%一阶高斯脉冲ppt
%y=(t-t0).*exp(-4*pi.*((t-t0).^2)/(tao.*tao))./tao; 
%%%二阶高斯脉冲book
%y=(-t).*exp(-1.*((t-t0)/tao).^2)./(((2*pi).^(1/2)).*(tao.^3));
%%%二阶高斯脉冲
y=(-1).*(t.^2-tao.^2).*exp(-1.*((t-t0)/tao).^2)./(((2*pi).^(1/2)).*(tao.^5));
% plot(t,y)
% grid on
% axis([-10,10,-0.5,0.5])
% xlabel('t')
% ylabel('y')


for i=1:InfoBitNum
    if Info(i)==1
        MSign(i,:)=repmat(y,1);
    else 
        MSign(i,:)=repmat(y0,1);
    end
end

MSignDisplay=MSign.';
MSignDisplay=MSignDisplay(:);
MSignDisplay=MSignDisplay.';
SendSignal=MSignDisplay;
SendSignal=awgn(SendSignal,20,'measured');

MSignDisplay=MSignDisplay(1:101*16*2);

MSignDisplay = awgn(MSignDisplay,20,'measured');

tDisplay=linspace(0,100*16*2,Lenth_t*16*2);
figure(2);
plot(tDisplay,MSignDisplay)
grid on
axis([-1,100*16*2,-0.5,0.5])
xlabel('t')
ylabel('MSign1')

%%解调
sumsign=zeros((length(SendSignal)/101),1);
verdict=zeros((length(SendSignal)/101),1);
for i=1:(length(SendSignal)/101-1)
    muli=(SendSignal((i-1)*101+1:(i-1)*101+101)).*y;
    sumsign(i,:)=trapz(muli,t);
    if sumsign(i)<=-0.1
        verdict(i)=1;
    else 
        verdict(i)=0;
    end
end

verdict=verdict';
RcvInfoBit=verdict;
RcvInfoBitNum=length(RcvInfoBit);
RcvVoice=zeros(RcvInfoBitNum/16,1);
for i=1:(RcvInfoBitNum/16-1)
    TempChar=mat2str(RcvInfoBit((i-1)*16+1:(i-1)*16+17));
    TempChar=TempChar(2:32);
    TempNum=bin2dec(TempChar);
    TempNum=double(TempNum);
    TempNum=TempNum/100000;
    TempNum=TempNum-0.3124;
    RcvVoice(i)=TempNum;
end
RcvVoice=RcvVoice';


RcvN=length(RcvVoice); % 计算信号x的长度
t=0:1/fs:(RcvN-1)/fs; % 计算时间范围，样本数除以采样频率
sound(RcvVoice,fs);
figure(3);
plot(t,RcvVoice);