软件能迅速测试设想构想，综合评测系统性能，快速设计更好方案来确保更高技术要求。同时MATLAB也是国家教委重点提倡的一种计算工具。
MATLAB主要由C语言编写而成，采用LAPACK 为底层支持软件包。 MATLAB的编程非常简单，它有着比其他任何计算机高级语言更高的编程效率、更好的代码可读性和移植性，以致被誉为“第四代”计算机语言，MATLAB是所有MATHWORKS公司产品的数值分析和图形基础环境。此外MATLAB 还拥有强大的2D和3D甚至动态图形的绘制功能，这样用户可以更直观、更迅速的进行多种算法的比较，从中找出最好的方案。
从通信系统分析与设计、滤波器设计、信号处理、小波分析、神经网络到控制系统、模糊控制等方面来看，MATLAB提供了大量的面向专业领域的工具箱。通过工具箱，以往需要复杂编程的算法开发任务往往只需一个函数就能实现，而且工具箱是开放的可扩展集，用户可以查看或修改其中的算法，甚至开发自己的算法。
目前， MATLAB已经广泛地应用于工程设计的各个领域，如电子、通信等领域；它已成为国际上最流行的计算机仿真软件设计工具。现在的MATLAB不再仅仅是一个矩阵实验室，而是一种实用的、功能强大的、不断更新的高级计算机编程语言。
现在从电子通信、自动控制图形分析处理到航天工业、汽车工业，甚至是财务工程。MATLAB都凭借其强大的功能获得了极大的用武之地。广大学生可以使用MATLAB来帮助进行信号处理、通信原理、线性系统、自动控制等课程的学习；科研工作者可以使用MATLAB进行理论研究和算法开发；工程师可以使用MATLAB进行系统级的设计与仿真。
3.2 8PSK调制部分
本设计采用相位选择法进行8PSK调制。其三位二进制序列与对应控制相位对应如下表1所示：
表1 三位二进制序列与对应控制相位
1 1 1------π/8 0 1 1------7π/8 1 0 0------13π/8 1 1 0------3π/8 0 0 1------9π/8 1 0 1------15π/8 0 1 0------5π/8 0 0 0------11π/8
6
对应上表，首先生成对应相位的同频载波，再根据输入二进制序列对每连续三位进行判决，输出相应相位的载波即可。其MATLAB程序如下： f=150
%抽样频率
g=(sign(rand(1,150)-0.5)+1)/2 %产生二进制序列 sn=randn(1,50*length(g)); dt=2*pi/149; t=0:dt:2*pi; si=[];co=[]; sit=[];sqt=[];
%si为正交分量，co为同相分量
%sit为同相分量幅度，sqt为正交分量幅度
%产生加性高斯白噪声
sb2=[]; %输入二进制序列
%8PSK调制过程
for n=1:3:length(g); %一次取3个二进制数 if g(n)==0 && g(n+1)==0 && g(n+2)==0 ±b2b3=000时正交分量和同相分量的幅值 it=-0.383*ones(1,150); qt=-0.924*ones(1,150);
b2=[zeros(1,50) zeros(1,50) zeros(1,50)] elseif g(n)==0 && g(n+1)==0 && g(n+2)==1 it=-0.924*ones(1,150); qt=-0.383*ones(1,150);
b2=[zeros(1,50) zeros(1,50) ones(1,50)] elseif g(n)==1 && g(n+1)==0 && g(n+2)==0 it=0.383*ones(1,150); qt=-0.924*ones(1,150);
b2=[ones(1,50) zeros(1,50) zeros(1,50) ] elseif g(n)==1 && g(n+1)==0 && g(n+2)==1 it=0.924*ones(1,150); qt=-0.383*ones(1,150);
b2=[ones(1,50) zeros(1,50) ones(1,50) ] elseif g(n)==0 && g(n+1)==1 && g(n+2)==0
±b2b3=001时
7
it=-0.383*ones(1,150); qt=0.924*ones(1,150);
b2=[zeros(1,50) ones(1,50) zeros(1,50) ] elseif g(n)==0 && g(n+1)==1 && g(n+2)==1 it=-0.924*ones(1,150); qt=0.383*ones(1,150);
b2=[zeros(1,50) ones(1,50) ones(1,50) ] elseif g(n)==1 && g(n+1)==1 && g(n+2)==1 it=0.924*ones(1,150); qt=0.383*ones(1,150);
b2=[ones(1,50) ones(1,50) ones(1,50) ] elseif g(n)==1 && g(n+1)==1 && g(n+2)==0 it=0.383*ones(1,150); qt=0.924*ones(1,150);
b2=[ones(1,50) ones(1,50) zeros(1,50) ] end
sb2=[sb2 b2];
c=cos(f*t); s=sin(f*t); sit=[sit it]; sqt=[sqt qt]; co=[co c]; si=[si s]; end
psk=sit.*co+sqt.*si;
%调制后的8psk信号
仿真结果如下图3.1和3.2所示：
8
图3.1 二进制序列
图3.2 8PSK调制
9
相关推荐：
loading