【通信】基于matlab数字基带传输系统【含Matlab源码 885期】

一、获取代码方式

获取代码方式1：
完整代码已上传我的资源：【通信】基于matlab数字基带传输系统【含Matlab源码 885期】

获取代码方式2：
通过订阅紫极神光博客付费专栏，凭支付凭证，私信博主，可获得此代码。

备注：
订阅紫极神光博客付费专栏，可免费获得1份代码（有效期为订阅日起，三天内有效）；

二、数字基带传输系统简介

二进制比特序列时间离散，取值离散，而实际的物理信道只能承载时间上连续的信号，为使二进制比特序列能在实际的物理信道传输，要用连续的波形表示比特0和比特1。
那么如何用连续的波形表示比特0和比特1？
要经过符号映射和波形成形这两个步骤：

1 符号映射
符号映射：符号映射是将抽象的比特0，比特1映射为实际的物理量，比如电压等。

1.1 符号映射的分类

1.1.1 一个比特映射为一个符号：
1）单极性码序列：比特0->0电平 比特1->+A电平， 如TTL电平
2）双极性码序列：比特0->-A电平 比特1->+A电平 如RS232
3) 三电平映射：0,1->0,-A,+A电平 如AMI码

1.1.2 多个比特映射为一个符号
4）多个比特映射为一个符号：k个比特为映射成一个符号：

1.2 波形成形

1.3 s(t)的功率谱（数字基带信号的由来）

2 数字传输系统的分类
2.1 基于传输信号的频率高低分为基带和频带传输系统。




矩形脉冲成形的数字信号只适合于在非带限信道中传输，带限脉冲成形的信号适合于在带限信道上传输。

三、部分源代码

fs=1;
Tc=4;
N=35;
a=1;
n=(-(N-1)/2):((N-1)/2);
[h]=umsf(N,Tc,a);
%发送信号波形
L=240;
A=4;
T=1;
[an]=Binary_sources(L);
[dnT]=Send_signal(A,T,L,an);
%发送滤波器输出波形
[xnT]=msf_signal(dnT,h);
%噪声为0 接收滤波器为直通型 所以
figure;
rnT=xnT;
stem(rnT(1:100),'.');
xlabel('n');
ylabel('r(nT)');
title('A=3接收滤波器(非匹配)的输出信号波形(a=0.3)');
%抽样判决点信号生成
[bn]=samples_signal(rnT,A,L,N);
%判决子函数
[cn]=judgement(bn);
%画眼图子函数
figure
[k]=eye_pattern(rnT,A);
title('二进制非匹配模式接收滤波器输出信号波形的眼图（a=0.3）')
figure;
scatterplot(bn);
%误比特率计算子函数
[Pe,p]=Pe_calculate(an,cn,L);
function [an]=Binary_sources(L)
  an=rand(1,L)%产生L个1以内的随机数
  for i=1:L
      if an(i)>0.5
          an(i)=1;
      else if an(i)<0.5
              an(i)=0;
          end
      end
  end
end
function [x]=idft(Xk,k,N)
n=(-(N-1)/2):((N-1)/2);
L=length(k); 
for a=1:N
    sum=0;
    for b=1:L
        sum=sum+(1/N)*Xk(b)*exp(j*2*pi*n(b)*k(a)/N);
        x(a)=sum;% 当a=1时即X的第一个点=w1*n的所有点
    end
end
function [cn]=judgement(bn)
for i=1:length(bn)
    if bn(i)>1/2
        cn(i)=1;
    else if bn(i)<1/2
            cn(i)=0;
        end
    end
end
end
function [h]=msf(N,Tc,a)
t=(-(N-1)/2):((N-1)/2);
n=-(N-1)/2:(N-1)/2;
k=n;
h=zeros(1,N);
for i=1:N;
    if(abs(t(i))==0)
        h(i)=1;
    elseif((1-4*a*a*t(i)*t(i)/Tc/Tc)==0)
        h(i)=Tc*sin(pi*t(i)/Tc)/t(i)/4;
    else
        h(i)=sin(pi*t(i)/Tc)*cos(a*pi*t(i)/Tc)/...
        (1-4*a*a*t(i)*t(i)/Tc/Tc)/(pi*t(i)/Tc);
    end
end
wd=(blackman(N))';
%h=h.*wd;
% [H]=dft(h,n,N);
H = h*exp(-1j*pi/N*k'*n);
H=sqrt(abs(H));
h = 1/N*H*exp(-1j*pi/N*k'*n);
% [h]=idft(H,k,N);
h=h.*wd;
h=real(h);
end
        

  

 

  
1
  
2
  
3
  
4
  
5
  
6
  
7
  
8
  
9
  
10
  
11
  
12
  
13
  
14
  
15
  
16
  
17
  
18
  
19
  
20
  
21
  
22
  
23
  
24
  
25
  
26
  
27
  
28
  
29
  
30
  
31
  
32
  
33
  
34
  
35
  
36
  
37
  
38
  
39
  
40
  
41
  
42
  
43
  
44
  
45
  
46
  
47
  
48
  
49
  
50
  
51
  
52
  
53
  
54
  
55
  
56
  
57
  
58
  
59
  
60
  
61
  
62
  
63
  
64
  
65
  
66
  
67
  
68
  
69
  
70
  
71
  
72
  
73
  
74
  
75
  
76
  
77
  
78
  
79
  
80
  
81
  
82
  
83
  
84
  
85
  
86
  
87
  
88
  
89
  
90
 

四、运行结果

五、matlab版本及参考文献

1 matlab版本
2014a

2 参考文献
[1] 沈再阳.精通MATLAB信号处理[M].清华大学出版社，2015.
[2]高宝建,彭进业,王琳,潘建寿.信号与系统——使用MATLAB分析与实现[M].清华大学出版社，2020.
[3]王文光,魏少明,任欣.信号处理与系统分析的MATLAB实现[M].电子工业出版社，2018.

文章来源: qq912100926.blog.csdn.net，作者：海神之光，版权归原作者所有，如需转载，请联系作者。

原文链接：qq912100926.blog.csdn.net/article/details/115979159