MATLAB大作业

发布 2020-02-25 07:47:28 阅读 5500

安阳工学院。

专业课课程大作业。

课程名称(中文) matlab编程应用

成绩。姓名黄红伟。

班级通信工程。

学号201002030060

日期2023年12月23日。

数字基带传输系统的**实现。

摘要]:matlab 语言是一种数据分析和处理功能十分强大的计算机应用软件 ,它可以将声音文件变换为离散的数据文件 , 然后利用其强大的矩阵运算能力处理数据 ,如数字滤波 、傅里叶变换 、时域和频域分析 、声音回放以及各种图的呈现等 , 信号处理是matlab 重要应用的领域之一 。数字通信系统有两种主要的通信模式:

数字频带传输通信系统,数字基带传输通信系统。数字基带传输系统指不经载波调制而直接传输数字基带信号的系统,常用于传输距离不太远的情况下。本次大作业通过分析数字基带信号传输的特性,运用数值**的方法,对数字基带传输系统作了模拟。

关键词]:matlab;数字基带传输系统。

1.背景介绍。

20世纪60年代出现了数字传输技术,它采用了数字信号来传递信息,从此通信进入了数字化时代。目前,通信网已基本实现数字化,在我国公众通信网中传输的信号主要是数字信号。数字通信技术的应用越来越广泛,例如数字移动通信、数字卫星通信、数字电视广播、数字光纤通信、数字微波通信、数字**通信、多**通信等等。

数字通信系统主要的两种通信模式:数字频带传输通信系统,数字基带传输通信系统。数字基带信号指未经调制的数字信号,它所占据的频谱是从零频或很低频率开始的。

数字基带传输系统指不经载波调制而直接传输数字基带信号的系统,常用于传输距离不太远的情况下。

研究数字基带传输系统的原因:实际中,基带传输不如频带传输应用广泛,但对基带传输的研究仍有非常重要的意义。这是因为:

第一,数字基带系统在近程数据通信系统中广泛采用;第二,数字基带系统的许多问题也是频带传输系统必须考虑的问题;第三,随着数字通信技术的发展,基带传输这种方式也有迅速发展的趋势,它不仅用于低速数据传输,而且还用于高速数据传输;第四,在理论上,任何一个线性调制的频带传输系统,总是可以有一个等效的基带载波调制系统所替代。因此,很有必要对基带传输系统进行综合系统的分析。

2.**理论基础。

matlab工具有很强的**能力可以**nr码、nrz码、ami码、pam码等各种编码及分析其功率谱。同时可以**通信系统的应用。这里先简单介绍下**的理论基础。

2.1.时域抽样和频域抽样。

若信号函数 s(t)为定义在时间区间(-,上的连续函数,按区间[-t/2,t/2]截短为st(t),再对st(t)按时间间隔t均匀取样,取样点数:n(t)=t/t;**时用这个样值函数来表示s(t)。如果信号的最高频率为fh,那么必须有fh〈〈0.

5/t,设为bs,bs=0.5/t 称为**系统的系统带宽。

若信号s(t)的频谱s(f)为定义在频率区间 (-上的连续函数,按区间[-bs,bs]截断s(f),然后f均匀取样。n(t)=2 bs /f=n;f=1/t如果**程序设定的频域采样间隔是f,那么就不能**截断时间1/f的信号。2.

2.频域分析。

函数t2f的功能是作傅立叶变换格式:x=t2f(x)函数f2t的功能是作傅立叶反变换格式:x=t2f(x)其中x是时域信号x(t)的截短抽样取得的取样值矢量。

x是傅立叶x(f)的取样值矢量。

2.2.取样判决和误码率的测量。

对于基带系统的研究,误码率是一个非常重要的概念,也是评价系统好坏的重要参数。在用matlab**系统中,在模拟实际的条件下,达到理想的误码率是我们的目标。

设矢量a表示发送的码序列,矢量y表示在判决观测到的叠加有噪声的接收信号。接收端在每隔一个码元间隔ts取一个点作为判决量。所有取样结果构成一个矢量;若取样时刻无偏差则b=y(1:

l:n),若取样时刻有k*△t的恒定偏差,则b=y(1+k:l:

n)。若判决门限为vth,则判决结果为c=sign(b-vth)双极性或c=sign(b-vth)+1)/2 单极性。若在一次观察中发送的码元是长度为m的矢量a,,对应的判决结果是c。

误码数是a与c中不相同的符号数,即n_err=length(find(a~=c)),于是误码率为pe=n_err/m。

3.基带传输系统的设计及**。

通信系统**原理:matlab用编程的方式在软件上对基带系统进行**。考虑到实际中有噪声的存在,以及对滤波器性能的要求,我们对基带传输函数选择升余弦滚降函数。

并且实际中,我们对系统有较低的误码率,并且系统在设计过程中易于完成,且成本低,我们采取匹配滤波器,在最大信噪比接收条件下对基带最佳化。并假设信道特性理想即,且我们忽略匹配器的延时。大大的减少了编程的复杂程度,但系统的误码率却相差很小。

为了分析**的可行性,并比较与理论上的研究值的差别,**时我们在绘出实测曲线的同时还考虑理论曲线以作以比较。

如图1就是一个以上所设计的基带传输系统。由于匹配器对于输入信号的要求,码序列我们采用随机的冲击序列信号。

在基带信道是理想限带情况下,适当地设计发送滤波器及接收滤波器,使得在接收端抽样时刻的码间干扰为零,则系统的合成传递函数必须满足以下两个条件3-1)

式中,,分别是发送滤波器,信道,接收滤波器的相频特性,是一时间延迟,,其中,,分别是发送滤波器,信道,接收滤波器引入的时延,w为升余弦滤波器的截止频率(此w值取决于符号速率及升余弦滚降因子的值)。

在接收端抽样时刻无码间干扰条件下,引起误码的是加性噪声,此时,最佳接收的滤波器应匹配于所接收的确定信号,使接收端抽样时刻的信噪比最大。设限带信道是理想低通特性,并设信道不引入时延()

则接收到的确定信号的频谱仅取决于发送滤波器的的特性,所以接收滤波器的应与发送滤波器共轭匹配,这样,在理想限带信道情况下,既要使接收端抽样时刻的抽样值无码间干扰,又要使得在抽样时刻抽样值的信噪比最大,则

综上所述,数字pam信号通过限带信道、并受到加性噪声干扰的情况下,在限带信道是理想低通条件下的最佳基带传输的发送及接收滤波器的设计是:总的收发系统的传递函数要符合无码间干扰基带传输的升余弦特性;且又要考虑在抽样时刻信噪比最大的收、发滤波共轭匹配的条件。在综合考虑这两方面因素的基础上,在设计发送及接受滤波器时,要使发送及接收滤波器的传递函数的模值分别是近似于升余弦的平方根频谱,其相移是线性的,时延是用来确保滤波器的物理可实现性。

编写的程序见附录主程序。

运行程序后,图2就是系统的**图。红色的代表实际曲线,蓝色的代表理论曲线。

图2基带接收系统。

从图2中可以看出,实测的曲线非常接近理论曲线。因此,系统的设计和**达到了本实验的要求。

参考文献。1] 南利平。 通信原理简明教程[m]. 清华大学出版社,2006

2] 朱光荣。 综合业务数字通信系统与通信网[j]. 电子科技文摘, 2006,5:61

3] 王旭,曹旭超。 移动通信系统的发展与未来[j]. 科技信息, 2007,18:41~65

4] 庞沁华,吴伟陵。 通信原理(合定本)[m]. 北京邮电大学出版社, 2005

5] 孙祥。 matlab7.0基础教程[m]. 北京: 清华大学出版社, 2005

6] 魏玲,刘敏。matlab通信**与应用[m]北京:国防工业出版社, 2001

7] 黄葆华,杨晓静。 通信原理[m]. 西安电子科技大学出版社, 2007

8] 樊昌信。 通信原理教程[m]. 北京:电子工业出版社, 2005

9] 申普兵。 数据通信技术[m]. 北京:国防工业出版社, 2006

10] 黎洪松。 数字通信原理[m]. 西安电子科技大学出版社, 2005

11] 王兴亮。 数字通信原理与技术[m]. 西安电子科技大学出版社, 2002

附录:数字基带传输系统的源程序。

global dt t f df n t

close all

clear eb_n0 pe

k=input('取样点数=2^k, k取13左右');

if isempty(k), k=13; end

z=input('每个信号取样点数=2^z, zif isempty(z), z=3; end

aa=input('滚降系数=[0.5]')

if aa==[aa=0.5;end

n=2^k;

l=2^z;m=n/l;

rb=2码速率是2mb/s

ts=1/rb;%码元间隔。

dt=ts/l;%时域采样间隔。

df=1/(n*dt);%频域采样间隔。

t=n*dt;%截短时间。

bs=n*df/2;%系统带宽。

f=[-bs+df/2:df:bs];%频域横坐标。

t=[-t/2+dt/2:dt:t/2];%时域横坐标。

g=sin(pi*t/ts).*cos(pi*t*aa/ts)./pi*t/ts.*(1-4*t.^2*aa^2/ts^2)];升余弦脉冲波形。

gg=t2f(g);gg=abs(gg);%升余弦脉冲的傅式变换。

gt=sqrt(gg);gr=gt;%最佳系统的发送接收滤波器的傅式变换。

for l1=1:20;

eb_n0(l1)=(l1-1);%eb/n0 in db

eb_n0(l1)=10^(eb_n0(l1)/10);

eb=1;n0=eb/eb_n0(l1);%信道的噪声谱密度。

sita=n0*bs;%信道中噪声功率。

n_err=0;%误码计数。

for l2=1:20;

b=sign(randn(1,m));

s=zeros(1,n);%产生冲激序列。

s(l/2:l:n)=b/dt;

ss=t2f(s);

s=ss.*gg;%信道的傅式变化。

a=f2t(s);a=real(a);%不加噪声的输出。

n_ch=sqrt(sita)*randn(size(t));信道噪声。

matlab大作业

matlab语言及应用。周一晚作业廖育州07010559 第一组程序。绘制选用潍柴发动机wd615.50时汽车行驶性能曲线图,ft ua图clear,clc 清除内存变量,清屏 n为发动机的转速。n 1000,1100,1200,1300,1400,1500,1600,1700,1800,1900,...

MATLAB大作业

实验七图像的形态 算。一 实验目的。1.掌握图像的奇异值分解与合成。2.能够利用模板操作实现图像的腐蚀运算3.能够利用模板操作实现图像的膨胀运算4.掌握图像的开 闭运算。二 实验内容。1.读下面的程序,分析该程序要完成什么样的功能。a 1 2 1 2 3 2 1 2 1 u d v svd a nx...

MATLAB大作业

1.设计一个程序,要求用户用键盘输入10次,输入内容可以是数字 英文 汉字 或者其他文字或者符号。如果输入的是数字,则把存入一个数字数组当中 如果输入的是英文的字符,则存入一个英文字符数组当中 如果是其他字符,则存入一个字符串数组中。当10次内容都输入完毕后,先输出英文数组中的内容,然后输出字符串数...