通信原理Matlab仿真实验报告
Matlab在《现代通信原理与系统》实验中的应用[优秀范文五篇]
![Matlab在《现代通信原理与系统》实验中的应用[优秀范文五篇]](https://img.taocdn.com/s3/m/17cafdafd1d233d4b14e852458fb770bf78a3bd4.png)
Matlab在《现代通信原理与系统》实验中的应用[优秀范文五篇]第一篇:Matlab在《现代通信原理与系统》实验中的应用Matlab在《现代通信原理与系统》实验中的应用摘要:为了提高研究生教学质量,提高学生学习兴趣和学习热情,使学生更加透彻地理解所学知识,拓展学生向研究性发展的外延培养,训练学生创新能力的培养,开发了《现代通信原理与系统》课程相关的仿真演示实验。
教学实践中,通过Matlab仿真实验演示,有效地激发了学生学习的主动性和积极性,增强了学生的感性认识,提高了?n 程教学效果,提高了人才培养质量。
关键词:通信原理;Matlab;实验教学;系统仿真中图分类号:TN911 文献标志码:A 文章编号:1674-9324(2018)25-0267-03一、引言《现代通信原理与系统》课程是光纤通信、移动通信、卫星通信等等课程的重要基础,该门课程数学知识复杂,理论性内容较多,部分涉及非线性电子线路,比较抽象,缺乏直观性,学生难以想象,不好理解,相关实验也是验证性实验,学生对实验的感受不深,对设备的运行原理、运行情况了解不深,这对培养学生综合思维能力、创新能力没有起到任何作用。
为了提高学生学习兴趣和学习热情,使学生更加透彻地理解所学知识,拓展学生向研究性发展的外延培养,训练学生创新能力的培养,笔者通过该门课程典型实验仿真,动态演示,在课堂上形象生动展现波形,帮助学生深入了解课程内容,提高学习效率。
二、模拟调制实验仿真让载波的某个参量随模拟调制信号的变化而变化的方式叫作模拟调制,模拟调制有线性模拟调制与非线性模拟调制。
通过线性模拟调制与非线性模拟调制,利用Matlab仿真,加深学生对于调制、解调概念的理解,掌握线性调制与非线性调制的区别。
通俗地讲,线性模拟调制就是将调制信号“放”到了载波的振幅参量上,在频域发生频谱的搬移,经过解调,将调制信号从载波的振幅参量上“取”出来,恢复成原始的调制信号。
这样做的目的有三:第一方面,把低频信号变换成利于无线发送或在信道中传输的高频信号;第二方面,使得多路信号在一个信道中同时传输,实现信道多路复用;第三方面,可以改善传输系统的性能。
通信原理MATlAB简单仿真报告
![通信原理MATlAB简单仿真报告](https://img.taocdn.com/s3/m/bdafd7800b4e767f5bcfce7c.png)
MATLAB简单仿真实验一、实验目的:学会利用MATLAB软件进行简单的仿真。
通过实验提高学生实际动手能力和编程能力,为日后从事通信工作奠定良好的基础。
二、实验内容:(1)绘制函数y=xe-x在0≤x≤1时的曲线。
(2)将输入的一段二进制代码编成单极性不归零码和双极性不归零码。
(3)学习使用simulink进行仿真建模三、仿真和实验结果:(1)x=0:0.1:1 %定义自变量的采样点取值数组y=x.*exp(-x) %利用数组运算计算各自变量采样点上的函数值plot(x,y),xlabel('x'),ylabel('y'),title('y=x*exp(-x)') %绘图(2)(a)单极性不归零码程序function y=snrz(x)t0=300;t=0:1/t0:length(x);for i=1:length(x)if(x(i)==1) %如果信息位为1for j=1:t0y((i-1)*t0+j)=1;%该码元对应的点值取1endelsefor j=1:t0y((i-1)*t0+j)=0; %否则,取0endendendy=[y,x(i)];%为了画图,要将y序列加上最后一位M=max(y);m=min(y);subplot(2,1,1)plot(t,y);grid on;axis([0,i,m-0.1,M+0.1]);title('1 0 0 1 1 0 0 0 0 1 0 1');(b)双极性码程序(3)a)在MATLAB的命令窗运行指令simulink,或点击命令窗中的图标,便打开如图所示的SIMULINK模型库浏览器(simulink Library Browser)。
b)在库浏览器中直接点击左侧分类目录中的Source子库,便可以看到各种信源模块,如图所示。
c)点击工具条上的图标,打开一个名为untitled的空白模型窗口。
中国地质大学基于MATLAB通信仿真实验报告
![中国地质大学基于MATLAB通信仿真实验报告](https://img.taocdn.com/s3/m/c2a81101844769eae009ed8c.png)
中国地质大学(武汉)通信原理matlab仿真报告姓名: 张彪_班号: 075112 _学号:20111002253院系:_机电学院专业:_通信工程指导教师:_赵娟老师联系方式:预留指导时间:1、选题为3.周期信号m(t),周期为2,在区间[0,2]内定义为:m(t)=t , 1=<t=<2;m(t)=2-t , 1=<t<=2。
(1)以时间间隔为0.1 的步长对此信号进行8 电平的均匀量化,画出量化输出波形;(2)求各量化点的量化误差,画出量化误差波形;(3)通过计算误差信号的功率,求该系统的SNRq(以dB 计)。
2、解题原理与思路模拟信号的取值范围:a~b,量化电平数=M则均匀量化时的量化间隔为:(b-a)/M且量化区间的端点mi=a+i v i=0,1,…,M 若量化输出电平qi取量化间隔的中点,则qi=[mi+m (i-1)]/2 i=1,2…,M .我的解题思路很简单,大致分为三部分,第一部分构造出题中所给的分段函数;第二部分自定义均匀量化函数,并一并写入计算信噪比等功能;第三部分将第一部分带入自定义函数并求出量化误差波形,最后输出波形,完成题目要求。
思路框图3、实验代码及注释:(1)分段函数的构造并代入自定义函数及求量化误差。
clear all;x1=0:0.1:1;y1=x1;x2=1:0.1:2y2=-x2+2;y3(1:11)=y1;y3(11:21)=y2; %构造出分段函数x=0:0.1:2;plot(x,y3); %输出分段函数w=juny(y3,1,8); %代入自定义函数求值plot(x,y3,x,w); %输出原波形与量化波形title('原波形与量化波形')figure(2)z=abs(y3-w); %求量化误差plot(x,z) %输出量化误差波形title('量化误差波形')(2)自定义函数:量化函数的实现function h=juny(f,v,L)%f是量化的信号序列%v是量化电平最大值%L是量化电平数n=length(f);t=2*v/L; %纵坐标取间隔p=zeros(1,L+1); %确定范围及初始化for i=1:L+1,p(i)=-v+(i-1)*t; %确定原信号抽样值endfor i=1:n %求量化值if f(i)>=v %在此范围内的量化取值(大于量化电平最大值时)h(i)=v;endif f(i)<=-v %在此范围内的量化取值(小于量化电平最小值时)h(i)=-v;endflag=0; %设立标志for j=2:L/2+1 %对量化电平数即纵坐标取半if flag==0if f(i)<p(j) %原信号与抽样信号循环比较h(i)=(p(j-1)+p(j))/2; %确定量化电平flag=1;endendendfor j=L/2+2:L+1 %另一部分的量化(同上)if flag==0if f(i)<p(j)h(i)=(p(j-1)+p(j))/2;flag=1;endendendendnq=v^2/(3*L^2); %计算量化噪声功率nq %显示量化噪声功率Am=2;snr=(3/2)*(Am/v)^2*L^2; %计算信噪比snr %显示信噪比snr=10*log(snr)/log(10); %将信噪比改为dB表示snr4、运行结果:运算结果:nq =0.0052 snr =384 snr =25.8433dB5、分析:在实验过程中出现许多问题,比如分段函数的构造写了一下午,本来是很简单,可以直接是一个三角函数,调用即可,不过如果是其他形式的分段函数就不能了,所以过程中用if 语句等等去做,不过都没做好,最后想起来用的这种方法构造的。
通信仿真实验报告
![通信仿真实验报告](https://img.taocdn.com/s3/m/bd285a00f6ec4afe04a1b0717fd5360cbb1a8d68.png)
通信仿真实验报告一、实验目的本次实验旨在通过通信仿真软件对一些典型通信系统进行仿真实验,通过实验数据分析和结果对比,加深对通信系统原理的理解。
二、实验设备和软件1. 通信仿真软件:MATLAB/Simulink2.实验设备:个人电脑3.其他实验所需设备:无三、实验内容1.AM调制与解调2.FM调制与解调3.FSK调制与解调4.PSK调制与解调5.QPSK调制与解调四、实验步骤1.AM调制与解调实验首先,使用MATLAB/Simulink搭建AM调制系统。
将一个正弦信号作为载波信号,用一个矩形脉冲信号进行调制,调制结果通过一个图示仪表进行查看。
然后,再搭建相应的AM解调系统,将调制后的信号经过解调系统,恢复为原始的矩形脉冲信号。
通过调整调制信号的幅度、频率等参数,观察调制和解调系统的输入输出波形变化情况,分析调制和解调的效果。
2.FM调制与解调实验同样使用MATLAB/Simulink搭建FM调制系统,将一个正弦信号作为载波信号,用一个矩形脉冲信号进行调制。
调制结果通过一个图示仪表进行查看。
接着,搭建相应的FM解调系统,将调制后的信号经过解调系统,恢复为原始的矩形脉冲信号。
通过调整调制信号的幅度、频率调制指数等参数,观察调制和解调系统的输入输出波形变化情况,并进行分析比较。
3.FSK调制与解调实验使用MATLAB/Simulink搭建FSK调制系统,将两个正弦信号分别作为两种调制信号,用一个矩形脉冲信号进行调制。
调制结果通过一个图示仪表进行查看。
接着,搭建相应的FSK解调系统,将调制后的信号经过解调系统,恢复为原始的矩形脉冲信号。
通过调整调制信号的幅度、频率等参数,观察调制和解调系统的输入输出波形变化情况,并进行分析比较。
4.PSK调制与解调实验使用MATLAB/Simulink搭建PSK调制系统,将一个正弦信号作为载波信号,用一个矩形脉冲信号进行调制。
调制结果通过一个图示仪表进行查看。
接着,搭建相应的PSK解调系统,将调制后的信号经过解调系统,恢复为原始的矩形脉冲信号。
通信原理matlab实验
![通信原理matlab实验](https://img.taocdn.com/s3/m/c5a62a87d0d233d4b14e6947.png)
在无码间串扰和噪声的理想情况下,波形无失真,每个码元将重叠在一起,最终在示波器上看到的是迹线又细又清晰的“眼睛”,“眼”开启得最大。当有码间串扰时,波形失真,码元不完全重合,眼图的迹线就会不清晰,引起“眼”部分闭合。若再加上噪声的影响,则使眼图的线条变得模糊,“眼”开启得小了,因此,“眼”张开的大小表示了失真的程度,反映了码间串扰的强弱。由此可知,眼图能直观地表明码间串扰和噪声的影响,可评价一个基带传输系统性能的优劣。另外也可以用此图形对接收滤波器的特性加以调整,以减小码间串扰和改善系统的传输性能。
从以上分析可知,眼图可以定性反映码间串扰的大小和噪声的大小,眼图还可以用来指示接受滤波器的调整,以减少码间串扰,改善系统性能。
使用了随机数,产生的基带信号也是随机的,经过加白噪声,产生的眼图也是不一样的。选取的采样频率也影响着最后的博兴和眼图。一开始的时候,提示加噪声的部分有问题,通过调试,发现工作空间设置的不准确,相加会出现错误,将工作空间加大,语法没错了,但是眼图出不来。后来经过仔细观察,发现做卷积的时候把一个元素弄错了,经过修改,最终成功运行。
2015年6月25日
基于MATLAB的第一类部分响应系统的
基带信号及其眼图实现
摘要
所谓眼图,是指通过用示波器观察接收端的基带信号波形,从而估计和调整系统性能的一种方法。这种方法具体做法是:用一个示波器跨接在抽样判决器的输入端,然后调整示波器水平扫描周期,使其与接受码元的周期同步。此时可以从示波器显示的图形上,观察码间干扰和信道噪声等因素影响的情况,从而估计系统性能的优劣程度。因为在传输二进制信号波形时,示波器显示的图像很像人的眼睛,故名“眼图”。
Matlab通信系统仿真实验报告
![Matlab通信系统仿真实验报告](https://img.taocdn.com/s3/m/6bd1bc67f8c75fbfc67db252.png)
Matlab通信原理仿真
学号:2142402
姓名:圣斌
实验一
Matlab 基本语法与信号系统分析
一、实验目的:1、掌握MATLAB 的基本绘图方法;
2、实现绘制复指数信号的时域波形。
二、实验设备与软件环境:
1、实验设备:计算机
2、软件环境:MATLAB R2009a
三、实验内容:
1、MATLAB 为用户提供了结果可视化功能,只要在命令行窗口输入相应的命令,结果就会用图形直接表示出来。
MATLAB 程序如下:
x = -pi:0.1:pi;
y1 = sin(x);
y2 = cos(x);
%准备绘图数据figure(1);
%打开图形窗口subplot(2,1,1);
%确定第一幅图绘图窗口plot(x,y1);
%以x ,y1绘图title('plot(x,y1)');
%为第一幅图取名为’plot(x,y1)’grid on;
%为第一幅图绘制网格线subplot(2,1,2)
%确定第二幅图绘图窗口plot(x,y2);
%以x ,y2绘图xlabel('time'),ylabel('y') %第二幅图横坐标为’time ’,纵坐标为’y ’运行结果如下图:
-4-3-2-101234-1-0.5
0.5
1
plot(x,y1)
-4-3-2-10
1234-1-0.50
0.5
1
time y。
Matlab与通信仿真课程设计报告
![Matlab与通信仿真课程设计报告](https://img.taocdn.com/s3/m/b1e09d503169a4517623a33e.png)
Matlab与通信仿真课程设计报告班级:12通信(1)班姓名:诸葛媛学号:Xb12680129实验一S-函数&锁相环建模仿真一、实验目的1.了解S函数和锁相环的工作原理2.掌握用S函数建模过程,锁相环载波提取仿真二、实验内容1、用S函数编写Similink基本模块(1)信源模块实现一个正弦波信号源,要求其幅度、频率和初始相位参数可以由外部设置,并将这个信号源进行封装。
(2)信宿和信号显示模块实现一个示波器木块,要求能够设定示波器显示的扫描周期,并用这个示波器观察(1)的信源模块(3)信号传输模块实现调幅功能,输入用(1)信源模块,输出用(2)信宿模块;基带信号频率1KHz,幅度1V;载波频率10KHz,幅度5V实现一个压控正弦振荡器,输入电压u(t)的范围为[v1,v2]V,输出正弦波的中心频率为f0Hz,正弦波的瞬时频率f随控制电压线性变化,控制灵敏度为kHz/V。
输入用(1)信源模块,输出用(2)信宿模块2、锁相环载波提取的仿真(1)利用压控振荡器模块产生一个受10Hz正弦波控制的,中心频率为100Hz,频偏范围为50Hz到150Hz的振荡信号,并用示波器模块和频谱仪模块观察输出信号的波形和频谱。
(2)构建一个抑制载波的双边带调制解调系统。
载波频率为10KHz,被调信号为1KHz正弦波,试用平方环恢复载波并进行解调。
(3)构建一个抑制载波的双边带调制解调系统。
载波频率为10KHz,被调信号为1KHz正弦波,试用科斯塔斯环恢复载波并进行解调。
(4)设参考频率源的频率为100Hz,要求设计并仿真一个频率合成器,其输出频率为300Hz。
并说明模型设计上与实例3.26的主要区别三、实验结果分析1、用S函数编写Similink基本模块(3)为了使S函数中输入信号包含多个,需要将其输入变量u初始为制定维数或自适应维数,而在S函数模块外部采用Simulink基本库中的复用器(Mux)将3个信号复用在一根信号线上。
Matlab通信系统仿真实验报告材料
![Matlab通信系统仿真实验报告材料](https://img.taocdn.com/s3/m/4323791202d276a200292ea6.png)
Matlab通信原理仿真学号:2142402:圣斌实验一 Matlab 基本语法与信号系统分析一、 实验目的:1、掌握MATLAB 的基本绘图方法;2、实现绘制复指数信号的时域波形。
二、 实验设备与软件环境:1、实验设备:计算机2、软件环境:MATLAB R2009a三、 实验容:1、MATLAB 为用户提供了结果可视化功能,只要在命令行窗口输入相应的命令,结果就会用图形直接表示出来。
MATLAB 程序如下:x = -pi:0.1:pi; y1 = sin(x); y2 = cos(x); %准备绘图数据 figure(1); %打开图形窗口subplot(2,1,1); %确定第一幅图绘图窗口 plot(x,y1); %以x ,y1绘图title('plot(x,y1)'); %为第一幅图取名为’plot(x,y1)’ grid on; %为第一幅图绘制网格线 subplot(2,1,2) %确定第二幅图绘图窗口 plot(x,y2); %以x ,y2绘图xlabel('time'),ylabel('y')%第二幅图横坐标为’time ’,纵坐标为’y ’运行结果如下图:-1-0.500.51plot(x,y1)-1-0.500.51timey2、上例中的图形使用的是默认的颜色和线型,MATLAB中提供了多种颜色和线型,并且可以绘制出脉冲图、误差条形图等多种形式图:MATLAB程序如下:x=-pi:.1:pi;y1=sin (x);y2=cos (x);figure (1);%subplot (2,1,1);plot (x,y1);title ('plot (x,y1)');grid on%subplot (2,1,2);plot (x,y2);xlabel ('time');ylabel ('y')subplot(1,2,1),stem(x,y1,'r') %绘制红色的脉冲图subplot(1,2,2),stem(x,y1,'g') %绘制绿色的误差条形图运行结果如下图:3、一个复指数信号可以分解为实部和虚部两部分。
通信原理matlab实验报告
![通信原理matlab实验报告](https://img.taocdn.com/s3/m/628a219ad05abe23482fb4daa58da0116c171f1c.png)
通信原理matlab实验报告《通信原理matlab实验报告》在现代通信系统中,通信原理是至关重要的一部分。
为了更好地理解和应用通信原理,我们进行了一系列的实验,并在本报告中分享我们的实验结果和分析。
首先,我们使用了Matlab软件进行了频谱分析实验。
通过对信号的频谱进行分析,我们能够更好地了解信号的频率分布特性,从而为信号的传输和处理提供了重要的参考。
在实验中,我们使用了不同的信号类型,并通过Matlab的频谱分析工具对其进行了分析。
通过实验结果,我们发现不同类型的信号在频谱上呈现出不同的特征,这为我们在实际通信系统中的信号处理提供了重要的指导。
其次,我们进行了调制解调实验。
调制是将数字信号转换为模拟信号的过程,而解调则是将模拟信号转换为数字信号的过程。
在实验中,我们使用Matlab模拟了调制解调过程,并通过实验结果验证了调制解调的正确性。
通过这一实验,我们深入理解了调制解调的原理和过程,并为实际通信系统中的信号处理提供了重要的参考。
最后,我们进行了信道编码解码实验。
信道编码是为了提高通信系统的可靠性和抗干扰能力而进行的一种技术手段。
在实验中,我们使用Matlab对信道编码进行了模拟,并通过实验结果验证了信道编码的效果。
通过这一实验,我们更加深入地理解了信道编码的原理和作用,为实际通信系统中的信号处理提供了重要的参考。
综上所述,通过本次实验,我们更加深入地理解了通信原理的相关知识,并通过Matlab软件进行了实际操作,加深了对通信原理的理解和应用。
这些实验结果对我们今后在通信系统设计和应用中将起到重要的指导作用。
希望通过这份实验报告的分享,能够对通信原理的学习和应用有所帮助。
通信原理实验报告matlab
![通信原理实验报告matlab](https://img.taocdn.com/s3/m/88454198a48da0116c175f0e7cd184254b351bfc.png)
通信原理实验报告matlab《通信原理实验报告:MATLAB》摘要:本实验报告基于通信原理课程的实验要求,利用MATLAB软件进行了一系列的实验。
通过实验,我们深入了解了通信原理中的一些重要概念和技术,并通过MATLAB软件进行了模拟和分析。
本实验报告将详细介绍实验的目的、原理、实验步骤、实验结果和分析,以及对实验过程中遇到的问题和解决方法进行了总结和讨论。
1. 实验目的本实验旨在通过使用MATLAB软件进行通信原理相关的实验,加深对通信原理中的相关概念和技术的理解,并通过实际操作加强对课程知识的掌握和应用能力。
2. 实验原理在本实验中,我们将涉及到通信原理中的一些重要概念和技术,包括信号的调制与解调、信道编码、信道调制等内容。
通过MATLAB软件,我们可以对这些概念和技术进行模拟和分析,从而更好地理解其原理和应用。
3. 实验步骤本实验中,我们将根据实验要求,依次进行一系列的实验步骤,包括信号的调制与解调、信道编码、信道调制等内容。
通过MATLAB软件,我们将对这些实验步骤进行模拟和分析,得到实验结果。
4. 实验结果和分析在实验过程中,我们得到了一系列的实验结果,并进行了详细的分析。
通过对这些实验结果的分析,我们可以更好地理解通信原理中的相关概念和技术,并加深对课程知识的理解和掌握。
5. 实验总结和讨论在实验过程中,我们也遇到了一些问题,并通过一些方法进行了解决。
在本部分,我们将对实验过程中遇到的问题和解决方法进行总结和讨论,以便更好地应对类似的实验问题。
通过本次实验,我们加深了对通信原理中的相关概念和技术的理解,并通过MATLAB软件进行了模拟和分析,得到了一系列的实验结果。
这些实验结果将有助于我们更好地理解通信原理中的相关知识,并加强对课程知识的掌握和应用能力。
同时,本次实验也为我们今后的学习和研究提供了一定的参考和借鉴。
通信原理Matlab仿真实验报告
![通信原理Matlab仿真实验报告](https://img.taocdn.com/s3/m/2a33abe026fff705cc170a47.png)
通信原理Matlab仿真实验报告学号:姓名:实验一1、利用Matlab实现矩形信号串信号分解与合成⑴ Matlab程序代码:①矩形信号串信号分解与合成函数rectexpd():function rectexpd(T1,T0,m)%矩形信号串信号分解与合成%T1:矩信号区间为(-T1/2,T1/2)%T0:矩形矩信信号串周期%m:傅里叶级数展开项次数t1=-T1/2:0.01:T1/2;t2=T1/2:0.01:(T0-T1/2); t=[(t1-T0)';(t2-T0)';t1';t2';(t1+T0)'];n1=length(t1);n2=length(t2); %根据周期矩形信号函数周期,计算点数f=[ones(n1,1);zeros(n2,1);ones(n1,1);zeros(n2 ,1);ones(n1,1)]; %构造周期矩形信号串y=zeros(m+1,length(t));y(m+1,:)=f';figure(1);plot(t,y(m+1,:)); %绘制周期矩形信号串axis([-(T0+T1/2)-0.5,(T0+T1/2)+0.5,0,1.2]); set(gca,'XTick',[-T0,-T1/2,T1/2,T0]);set(gca,'XTickLabel',{'-T0','-T1/2','T1/2','T0'}) ;title('矩形信号串');grid on; a=T1/T0;pause; %绘制离散幅度谱freq=[-20:1:20];mag=abs(a*sinc(a*freq)); stem(freq,mag);x=a*ones(size(t));for k=1:m %循环显示谐波叠加图形pause;x=x+2*a*sinc(a*k)*cos(2*pi*t*k/T0);y(k,:)=x;plot(t,y(m+1,:));hold on;plot(t,y(k,:));hold off;grid on;axis([-(T0+T1/2)-0.5,[T0+T1/2]+0.5,-0.5,1.5]) ;title(strcat(num2str(k),'次谐波叠加')); xlabel('t');endpause;plot(t,y(1:m+1,:));grid on;axis([-T0/2,T0/2,-0.5,1.5]);title('各次谐波叠加');xlabel('t');②在命令窗口调用rectexpd()函数:T1=5;T0=10;m=5;rectexpd(T1,T0,m)⑵Matlab仿真结果:1次谐波叠加t3次谐波叠加tt2次谐波叠加t4次谐波叠加tt2、利用Matlab实现连续信号卷积运算⑴ Matlab程序代码:①连续信号卷积运算的通用函数sconv():function [f,k]=sconv(f1,f2,k1,k2,p) f=conv(f1,f2);f=f*p;k0=k1(1)+k2(1);k3=length(f1)+length(f2)-2;k=k0:p:k3*p; subplot(2,2,1); plot(k1,f1); title('f1(t)'); xlabel('t'); ylabel('f1(t)');subplot(2,2,2);plot(k2,f2);title('f2(t)'); xlabel('t'); ylabel('f2(t)'); subplot(2,2,3); plot(k,f); title('f(t)=f1(t)*f2(t)'); xlabel('t'); ylabel('f(t)');h=get(gca,'position'); h(3)=2.5*h(3);set(gca,'position',h);②在命令窗口调用sconv ()函数: p=0.01; k1=0:p:2; f1=exp(-k1); k2=0:p:3;f2=ones(1,length(k2)); [f,k]=sconv(f1,f2,k1,k2,p)⑵Matlab 仿真结果:0.511.5200.51f1(t)tf 1(t )01230.511.52f2(t)tf 2(t )00.51 1.522.533.544.550.51f(t)=f1(t)*f2(t)tf (t )实验二1、用Matlab模拟图形A律解码⑴Matlab程序代码:①建立函数ADecode()function y=ADecode(code,n) codesize=size(code);cr_len=codesize(1);cl_len=codesize(2);ca=zeros(1,cl_len-1);for i=1:cr_lenca=code(i,2:n);s=0;for j=1:n-1s=s+ca(j)*2^(n-1-j);enda=code(i,1);y(i)=s*((-1)^(a+1)); endy=y/(2^(n-1));A=87.6;A1=1+log(A); for j=1:length(y)if(y(j)>=0)if(y(i)<=1/A1)y(j)=y(j)*A1/A;elsey(j)=exp(y(j)*A1-1/A);endelsetemp=-y(j);if(temp<=1/A1)y(j)=-temp*A1/A;elsey(j)=-exp(temp*A1-1)/A;endendend②建立函数APCM()function code=APCM(x,n)xmax=max(abs(x));x=x/xmax;xlen=length(x);y=zeros(1,xlen);A=87.6;A1=1+log(A);for i=1:xlenif x(i)>=0if x(i)<=1/Ay(i)=(A*x(i))/A1;elsey(i)=(1+log(A*x(i)))/A1;endelsex1=-x(i);if x1<=1/Ay(i)=-(A*x1)/A1;elsey(i)=-(1+log(A*x1))/A1;endendendy1=y*(2^(n-1)-1);y1=round(y1);code=zeros(length(y1),n); c2=zeros(1,n-1);for i=1:length(y1)if(y1(i)>0)c1=1;elsec1=0;y1(i)=-y1(i);endfor j=1:n-1r=rem(y1(i),2);y1(i)=(y1(i)-r)/2;c2(j)=r;endc2=fliplr(c2);code(i,:)=[c1 c2];end③在新函数中调用前两个函数 t=0:0.01:1; x=sin(2*pi*t); code=APCM(x,7); y=ADecode(code,7); subplot(2,1,1) plot(t,x);title('原函数的图形'); subplot(2,1,2) plot(t,y);title('解码后函数的图形');⑵Matlab 仿真结果:00.10.20.30.40.50.60.70.80.91-1-0.500.51原函数的图形00.10.20.30.40.50.60.70.80.91-1-0.50.5解码后函数的图形2、用Matlab 模拟DSB 调制及解调过程 ⑴Matlab 程序代码: close all; clear all;dt=0.001;%采样时间间隔 fm=1; %信源最高频率 fc=10; %载波中心频率 N=4096;T=N*dt; t=0:dt:T-dt;mt=sqrt(2)*cos(2*pi*fm*t);%信源s_dsb=mt.*cos(2*pi*fc*t); %DSB-SC 双边带抑制载波调幅B=2*fm;figure(1); subplot(311) plot(t,s_dsb,'b-');hold on; %画出DSB 信号波形plot(t,mt,'r--');%画出m(t)信号波形 title('DSB 调制信号'); xlabel('t'); grid on;%DSB demodulation rt=s_dsb.*cos(2*pi*fc*t); rt=rt-mean(rt); [f,rf]=T2F(t,rt); [t,rt]=lpf(f,rf,B); subplot(312) plot(t,rt,'k-'); hold on;plot(t,mt/2,'r--');title('相干解调后的信号波形与输入信号的比较'); xlabel('t'); grid on; subplot(313)[f,sf]=T2F(t,s_dsb); %求调制信号的频谱psf=(abs(sf).^2)/T; %求调制信号的功率谱密度 plot(f,psf);axis([-2*fc 2*fc 0 max(psf)]); title('DSB 信号功率谱'); xlabel('f'); grid on;⑵Matlab 仿真结果:00.51 1.52 2.53 3.54 4.5-22DSB 调制信号t0.511.522.533.544.5-101相干解调后的信号波形与输入信号的比较t-20-15-10-5051015200.5DSB 信号功率谱f实验三1、用Matlab模拟双极性归零码⑴ Matlab程序代码:function y=drz(x)t0=300;x=[1 0 0 1 1 0 0 0 0 1 0 1];t=0:1/t0:length(x);for i=1:length(x);if(x(i)==1)for j=1:t0/2y(t0/2*(2*i-2)+j)=1;y(t0/2*(2*i-1)+j)=0;endelsefor j=1:t0/2y(t0/2*(2*i-2)+j)=-1;y(t0/2*(2*i-1)+j)=0;endendendy=[y,x(i)];M=max(y);m=min(y); subplot(211)plot(t,y);grid on;axis([0,i,m-0.1,M+0.1]);title('1 0 0 1 1 0 0 0 0 1 0 1 ');⑵Matlab仿真结果:2、用Matlab产生2FSK信号⑴ Matlab程序代码:①建立函数fskdigital()function fskdigital(s,f1,f2) t=0:2*pi/99:2*pi;m1=[];c1=[];b1=[];for n=1:length(s)if s(n)==0;m=ones(1,100);c=sin(f2*t);b=zeros(1,100) else s(n)==1;m=ones(1,100);c=sin(f1*t);b=ones(1,100)endm1=[m1 m]; c1=[c1 c]; b1=[b1 b]; endfsk=c1.*m1;subplot(211);plot(b1,'r')title('原始信号');axis([0 100*length(s) -0.1 1.1]);grid on;subplot(212);plot(fsk)title('2FSK信号');grid on;②在命令窗口调用函数fskdigital() >> s=[1 0 1 1 0 0 1 0]; f1=200; f2=100;fskdigital(s,f1,f2)⑵Matlab 仿真结果:10020030040050060070080000.20.40.60.81原始信号100200300400500600700800-1-0.500.512FSK 信号3、用Matlab 的simulink 模块模拟三角波的分解与还原 ⑴模块图⑵仿真结果:①三角波原图②用矩形波近似的冲击函数③分解后的三角波④经低通还原的三角波。
MATLAB通信系统仿真实验报告
![MATLAB通信系统仿真实验报告](https://img.taocdn.com/s3/m/86d4b8dabe23482fb4da4c89.png)
MATLAB通信系统仿真实验报告实验一、MATLAB的基本使用与数学运算目的:学习MATLAB的基本操作,实现简单的数学运算程序。
内容:1-1 要求在闭区间[0,2π]上产生具有10个等间距采样点的一维数组。
试用两种不同的指令实现。
运行代码:x=[0:2*pi/9:2*pi]运行结果:1-2 用M文件建立大矩阵xx=[ 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.91.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.92.1 2.2 2.3 2.4 2.5 2.6 2.7 2.8 2.93.1 3.2 3.3 3.4 3.5 3.6 3.7 3.8 3.9]代码:x=[ 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.91.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.92.1 2.2 2.3 2.4 2.5 2.6 2.7 2.8 2.93.1 3.2 3.3 3.4 3.5 3.6 3.7 3.8 3.9]m_mat运行结果:1-3已知A=[5,6;7,8],B=[9,10;11,12],试用MATLAB分别计算A+B,A*B,A.*B,A^3,A.^3,A/B,A\B.代码:A=[5 6;7 8] B=[9 10;11 12] x1=A+B X2=A-B X3=A*B X4=A.*B X5=A^3 X6=A.^3 X7=A/B X8=A\B运行结果:1-4任意建立矩阵A,然后找出在[10,20]区间的元素位置。
程序代码及运行结果:代码:A=[12 52 22 14 17;11 10 24 03 0;55 23 15 86 5 ] c=A>=10&A<=20运行结果:1-5 总结:实验过程中,因为对软件太过生疏遇到了些许困难,不过最后通过查书与同学交流都解决了。
例如第二题中,将文件保存在了D盘,而导致频频出错,最后发现必须保存在MATLAB文件之下才可以。
matlab 仿真实验报告
![matlab 仿真实验报告](https://img.taocdn.com/s3/m/15992fae112de2bd960590c69ec3d5bbfd0adabf.png)
matlab 仿真实验报告Matlab 仿真实验报告引言:在科学研究和工程应用中,仿真实验是一种非常重要的手段。
通过在计算机上建立数学模型和进行仿真实验,我们可以更好地理解和预测现实世界中的各种现象和问题。
Matlab作为一种强大的科学计算软件,被广泛应用于各个领域的仿真实验中。
本文将介绍我进行的一次基于Matlab的仿真实验,并对实验结果进行分析和讨论。
实验背景:在电子通信领域中,信号的传输和接收是一个重要的研究方向。
而在进行信号传输时,会受到各种信道的影响,如噪声、衰落等。
为了更好地理解信道的特性和优化信号传输方案,我进行了一次关于信道传输的仿真实验。
实验目的:本次实验的目的是通过Matlab仿真,研究不同信道条件下信号传输的性能,并对比分析不同传输方案的优劣。
实验步骤:1. 信道建模:首先,我需要建立信道的数学模型。
根据实际情况,我选择了常见的高斯信道模型作为仿真对象。
通过Matlab提供的函数,我可以很方便地生成高斯噪声,并将其加入到信号中。
2. 信号传输方案设计:接下来,我需要设计不同的信号传输方案。
在实验中,我选择了两种常见的调制方式:频移键控(FSK)和相移键控(PSK)。
通过调整不同的调制参数,我可以模拟不同的传输效果。
3. 信号传输仿真:在信道模型和传输方案设计完成后,我开始进行信号传输的仿真实验。
通过Matlab提供的信号处理函数,我可以很方便地生成调制后的信号,并将其传输到信道中。
4. 信号接收和解调:在信号传输完成后,我需要进行信号接收和解调。
通过Matlab提供的信号处理函数,我可以很方便地对接收到的信号进行解调,并还原出原始的信息信号。
5. 仿真结果分析:最后,我对仿真结果进行分析和讨论。
通过对比不同信道条件下的传输性能,我可以评估不同传输方案的优劣,并得出一些有价值的结论。
实验结果与讨论:通过对不同信道条件下的信号传输仿真实验,我得到了一些有价值的结果。
首先,我观察到在高斯噪声较大的信道条件下,PSK调制比FSK调制具有更好的抗干扰性能。
实验一 模拟通信的MATLAB仿真
![实验一 模拟通信的MATLAB仿真](https://img.taocdn.com/s3/m/facd2d224b73f242336c5fbf.png)
实验一 模拟通信的MATLAB 仿真姓名:左立刚 学号:031040522简要说明:实验报告注意包括AM ,DSB ,SSB ,VSB ,FM 五种调制与解调方式的实验原理,程序流程图,程序运行波形图,simulink 仿真模型及波形,心得体会,最后在附录中给出了m 语言的源程序代码。
一.实验原理1.幅度调制(AM )幅度调制(AM )是指用调制信号去控制高频载波的幅度,使其随调制信号呈线性变化的过程。
AM 信号的数学模型如图3-1所示。
图2-1 AM 信号的数学模型为了分析问题的方便,令δ=0,1.1 AM 信号的时域和频域表达式()t S AM=[A+m ()t ]cos t c ω (2-1)()t S AM =A 0π[()()ωωωωδC C ++-]+()()[]ωωωωc c M M ++-21(2-2)AM 信号的带宽2=BAMfH(2-3)式中,fH为调制信号的最高频率。
2.1.3 AM 信号的功率P AM 与调制效率ηAMPAM=()22220t m A+=PPmc+(2-4)式中,P C=2A为不携带信息的载波功率;()22t m P m=为携带信息的边带功率。
()()t t m A m PPAMCAM222+==η(2-5) AM 调制的优点是可用包络检波法解调,不需要本地同步载波信号,设备简单。
AM 调制的最大缺点是调制效率低。
2.2、双边带调制(DSB )如果将在AM 信号中载波抑制,只需在图3-1中将直流 A去掉,即可输出抑制载波双边带信号。
2.2.1 DSB 信号的时域和频域表达式()()t t m t c DSB S ωcos = (2-6)()()()[]ωωωωωCCDSB M M S ++-=21 (2-7)DSB 信号的带宽fBBHAMDSB2==(2-8)DSB 信号的功率及调制效率由于不再包含载波成分,因此,DSB 信号的功率就等于边带功率,是调制信号功率的一半,即()()t t m P S PC DSB DSB2221===(2-9) 显然,DSB 信号的调制效率为100%。
matlab与通信仿真实验报告
![matlab与通信仿真实验报告](https://img.taocdn.com/s3/m/8c4318bfc9d376eeaeaad1f34693daef5ef713e8.png)
matlab与通信仿真实验报告《Matlab与通信仿真实验报告》摘要:本实验报告通过使用Matlab软件进行通信仿真实验,对通信系统的性能进行了评估和分析。
首先介绍了通信系统的基本原理和模型,然后利用Matlab软件搭建了通信系统的仿真模型,并进行了实验验证。
通过实验结果的分析,得出了通信系统的性能指标,为通信系统的设计和优化提供了重要参考。
一、引言通信系统是现代信息社会中不可或缺的基础设施,它承载着各种类型的信息传输和交换。
通信系统的性能直接影响着信息传输的质量和效率,因此对通信系统的性能评估和分析具有重要意义。
Matlab软件是一种功能强大的科学计算软件,它提供了丰富的工具和函数库,可以用于通信系统的建模、仿真和分析。
本实验报告将利用Matlab软件进行通信系统的仿真实验,对通信系统的性能进行评估和分析。
二、通信系统的基本原理和模型通信系统由发送端、信道和接收端组成,发送端将信息转换成电信号发送出去,经过信道传输后,接收端将电信号转换成信息。
通信系统的性能评估主要包括信号传输质量、误码率、信噪比等指标。
在本实验中,我们将以常见的调制解调技术为例,建立通信系统的仿真模型。
三、Matlab软件在通信系统仿真中的应用Matlab软件提供了丰富的工具和函数库,可以用于通信系统的建模、仿真和分析。
在本实验中,我们将利用Matlab软件搭建通信系统的仿真模型,包括信号调制、信道传输、信号解调等过程。
通过Matlab软件的仿真实验,我们可以得到通信系统的性能指标,如误码率、信噪比等。
四、实验结果分析通过Matlab软件进行通信系统的仿真实验,我们得到了一系列实验结果。
通过对实验结果的分析,我们可以评估通信系统的性能,比如误码率随信噪比的变化规律、不同调制方式的性能比较等。
这些实验结果对于通信系统的设计和优化具有重要的参考价值。
五、结论本实验报告利用Matlab软件进行通信系统的仿真实验,对通信系统的性能进行了评估和分析。
用MatLab仿真通信原理系列实验
![用MatLab仿真通信原理系列实验](https://img.taocdn.com/s3/m/78cc1576c950ad02de80d4d8d15abe23482f03f8.png)
用MatLab仿真通信原理系列实验一、引言通信原理是现代通信领域的基础理论,通过对通信原理的研究和仿真实验可以更好地理解通信系统的工作原理和性能特点。
MatLab作为一种强大的数学计算软件,被广泛应用于通信原理的仿真实验中。
本文将以MatLab为工具,介绍通信原理系列实验的仿真步骤和结果。
二、实验一:调制与解调1. 实验目的通过MatLab仿真,了解调制与解调的基本原理,并观察不同调制方式下的信号特征。
2. 实验步骤(1)生成基带信号:使用MatLab生成一个基带信号,可以是正弦波、方波或任意复杂的波形。
(2)调制:选择一种调制方式,如调幅(AM)、调频(FM)或相移键控(PSK),将基带信号调制到载波上。
(3)观察调制后的信号:绘制调制后的信号波形和频谱图,观察信号的频谱特性。
(4)解调:对调制后的信号进行解调,还原出原始的基带信号。
(5)观察解调后的信号:绘制解调后的信号波形和频谱图,与原始基带信号进行对比。
3. 实验结果通过MatLab仿真,可以得到不同调制方式下的信号波形和频谱图,观察到调制后信号的频谱特性和解调后信号的还原效果。
可以进一步分析不同调制方式的优缺点,为通信系统设计提供参考。
三、实验二:信道编码与解码1. 实验目的通过MatLab仿真,了解信道编码和解码的基本原理,并观察不同编码方式下的误码率性能。
2. 实验步骤(1)选择一种信道编码方式,如卷积码、纠错码等。
(2)生成随机比特序列:使用MatLab生成一组随机的比特序列作为输入。
(3)编码:将输入比特序列进行编码,生成编码后的比特序列。
(4)引入信道:模拟信道传输过程,引入噪声和干扰。
(5)解码:对接收到的信号进行解码,还原出原始的比特序列。
(6)计算误码率:比较解码后的比特序列与原始比特序列的差异,计算误码率。
3. 实验结果通过MatLab仿真,可以得到不同编码方式下的误码率曲线,观察不同信道编码方式对信号传输性能的影响。
通信原理仿真实验报告
![通信原理仿真实验报告](https://img.taocdn.com/s3/m/801bd5c3951ea76e58fafab069dc5022abea464c.png)
通信原理仿真实验报告一、引言通信原理是现代社会中不可或缺的一部分,它涉及到信息的传输和交流。
为了更好地理解通信原理的工作原理和效果,我们进行了一次仿真实验。
本报告将详细介绍实验的目的、方法、结果和分析。
二、实验目的本次实验的目的是通过仿真实验,深入了解通信原理的基本原理和信号传输过程,掌握通信系统中常见的调制解调技术,并通过实验验证理论知识的正确性。
三、实验方法1. 实验平台:我们使用MATLAB软件进行仿真实验,该软件具有强大的信号处理和仿真功能,可以模拟真实的通信环境。
2. 实验步骤:a. 设计信号源:根据实验要求,我们设计了一种特定的信号源,包括信号的频率、幅度和相位等参数。
b. 调制过程:通过调制技术将信号源与载波信号进行合成,得到调制后的信号。
c. 信道传输:模拟信号在信道中的传输过程,包括信号的衰减、噪声的干扰等。
d. 解调过程:通过解调技术将接收到的信号还原为原始信号。
e. 信号分析:对解调后的信号进行频谱分析、时域分析等,以验证实验结果的准确性。
四、实验结果我们进行了多组实验,得到了一系列的实验结果。
以下是其中两组实验结果的示例:1. 实验一:调幅调制a. 信号源:频率为1kHz的正弦信号。
b. 载波信号:频率为10kHz的正弦信号。
c. 调制后的信号:将信号源与载波信号相乘,得到调制后的信号。
d. 信号分析:对调制后的信号进行频谱分析,得到频谱图。
e. 解调过程:通过解调技术,将接收到的信号还原为原始信号。
f. 结果分析:通过对比解调后的信号与原始信号,验证了调幅调制的正确性。
2. 实验二:频移键控调制a. 信号源:频率为1kHz的正弦信号。
b. 载波信号:频率为10kHz的正弦信号。
c. 调制后的信号:将信号源与载波信号相加,得到调制后的信号。
d. 信号分析:对调制后的信号进行频谱分析,得到频谱图。
e. 解调过程:通过解调技术,将接收到的信号还原为原始信号。
f. 结果分析:通过对比解调后的信号与原始信号,验证了频移键控调制的正确性。
通信原理Matlab仿真实验报告
![通信原理Matlab仿真实验报告](https://img.taocdn.com/s3/m/df107f8cc77da26924c5b00f.png)
实验二1、用Matlab模拟图形A律解码⑴Matlab程序代码:①建立函数ADecode()function y=ADecode(code,n) codesize=size(code);cr_len=codesize(1);cl_len=codesize(2);ca=zeros(1,cl_len-1);for i=1:cr_lenca=code(i,2:n);s=0;for j=1:n-1s=s+ca(j)*2^(n-1-j);enda=code(i,1);y(i)=s*((-1)^(a+1)); endy=y/(2^(n-1));A=87.6;A1=1+log(A);for j=1:length(y)if(y(j)>=0)if(y(i)<=1/A1)y(j)=y(j)*A1/A;elsey(j)=exp(y(j)*A1-1/A);endelsetemp=-y(j);if(temp<=1/A1)y(j)=-temp*A1/A;elsey(j)=-exp(temp*A1-1)/A;endendend②建立函数APCM()function code=APCM(x,n)xmax=max(abs(x));x=x/xmax;xlen=length(x);y=zeros(1,xlen);A=87.6;A1=1+log(A);for i=1:xlenif x(i)>=0if x(i)<=1/Ay(i)=(A*x(i))/A1;elsey(i)=(1+log(A*x(i)))/A1;endelsex1=-x(i);if x1<=1/Ay(i)=-(A*x1)/A1;elsey(i)=-(1+log(A*x1))/A1;endendendy1=y*(2^(n-1)-1);y1=round(y1);code=zeros(length(y1),n); c2=zeros(1,n-1);for i=1:length(y1)if(y1(i)>0)c1=1;elsec1=0;y1(i)=-y1(i);endfor j=1:n-1r=rem(y1(i),2);y1(i)=(y1(i)-r)/2;c2(j)=r;endc2=fliplr(c2);code(i,:)=[c1 c2];实验三1、用Matlab模拟双极性归零码⑴ Matlab程序代码:function y=drz(x)t0=300;x=[1 0 0 1 1 0 0 0 0 1 0 1];t=0:1/t0:length(x);for i=1:length(x);if(x(i)==1)for j=1:t0/2y(t0/2*(2*i-2)+j)=1;y(t0/2*(2*i-1)+j)=0;endelsefor j=1:t0/2y(t0/2*(2*i-2)+j)=-1;y(t0/2*(2*i-1)+j)=0;endendendy=[y,x(i)];M=max(y);m=min(y); subplot(211)plot(t,y);grid on;axis([0,i,m-0.1,M+0.1]);title('1 0 0 1 1 0 0 0 0 1 0 1 ');⑵Matlab仿真结果:2、用Matlab产生2FSK信号⑴ Matlab程序代码:①建立函数fskdigital()function fskdigital(s,f1,f2) t=0:2*pi/99:2*pi;m1=[];c1=[];b1=[];for n=1:length(s)if s(n)==0;m=ones(1,100);c=sin(f2*t);b=zeros(1,100) else s(n)==1;m=ones(1,100);c=sin(f1*t);b=ones(1,100)endm1=[m1 m]; c1=[c1 c]; b1=[b1 b]; endfsk=c1.*m1;subplot(211);plot(b1,'r')title('原始信号');axis([0 100*length(s) -0.1 1.1]);grid on;subplot(212);plot(fsk)title('2FSK信号');⑵仿真结果:①三角波原图②用矩形波近似的冲击函数③分解后的三角波④经低通还原的三角波。
用MatLab仿真通信原理系列实验
![用MatLab仿真通信原理系列实验](https://img.taocdn.com/s3/m/76237ff3ba0d4a7302763ad5.png)
实验三:低通信号的抽样定理…………………………….………….16
第四章小结……………………………………………...18
致谢………………………………………………………18
参考文献…………………………………………………18
摘要:本毕业设计用Matlab中的建模仿真工具Simulink对通信原理实验进行仿真。作为系列实验的第一部分,包括模拟信号的线性调制解调(AM、DSB、SSB)过程、扰码与解扰实验和低通信号的抽样定理实验。论文中讲述了Matlab的基础知识、Simulink仿真操作方法以及在通信系统中的应用,对被仿真实验的基本原理也进行了简要介绍。通过本设计对构造通信原理虚拟实验室这一课题进行了初步的探索。
MATLAB系统共有5个主要部分构成:
1)MATLAB语言;
2)MATLAB工作环境;
3)MATLAB数学函数库;
4)MATLAB图形处理系统;
5)MATLAB应用程序接口。
1.MATLAB语言
MATLAB语言是一种面向对象的高级语言,正如前面所述,它以矩阵作为最基本的数据结构。MATLAB语言由自己独特的数据结构、输入输出功能、流程控制语句和函数。MATLAB在工程计算方面具有其他高级语言无法比拟的优越性,它集计算、数据可视化、程序设计于一体,并能将数学问题和解决方案以用户熟悉的数学符号表示出来。
(1)实现了可视化建模,用户通过简单的鼠标操作就可建立起直观的系统模型,并进行仿真;
(2)实现了多工作环境间文件互用和数据交换,如Simulink与MATLAB、FORTRAN以及C的工作环境的信息交换都可以方便的实现;
(3)把理论研究和工程实现有机地结合在一起。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
通信原理Matlab仿真实验报告学号:姓名:实验一1、利用Matlab 实现矩形信号串信号分解与合成 ⑴ Matlab 程序代码:①矩形信号串信号分解与合成函数rectexpd (): function rectexpd(T1,T0,m) %矩形信号串信号分解与合成 %T1:矩信号区间为(-T1/2,T1/2) %T0:矩形矩信信号串周期 %m:傅里叶级数展开项次数 t1=-T1/2:0.01:T1/2;t2=T1/2:0.01:(T0-T1/2); t=[(t1-T0)';(t2-T0)';t1';t2';(t1+T0)']; n1=length(t1);n2=length(t2); %根据周期矩形信号函数周期,计算点数 f=[ones(n1,1);zeros(n2,1);ones(n1,1);zeros(n2,1);ones(n1,1)]; %构造周期矩形信号串 y=zeros(m+1,length(t));y(m+1,:)=f'; figure(1); plot(t,y(m+1,:)); %绘制周期矩形信号串 axis([-(T0+T1/2)-0.5,(T0+T1/2)+0.5,0,1.2]); set(gca,'XTick',[-T0,-T1/2,T1/2,T0]); set(gca,'XTickLabel',{'-T0','-T1/2','T1/2','T0'}); title('矩形信号串');grid on; a=T1/T0; pause; %绘制离散幅度谱 freq=[-20:1:20];mag=abs(a*sinc(a*freq)); stem(freq,mag); x=a*ones(size(t)); for k=1:m %循环显示谐波叠加图形 pause; x=x+2*a*sinc(a*k)*cos(2*pi*t*k/T0); y(k,:)=x; plot(t,y(m+1,:));hold on; plot(t,y(k,:));hold off;grid on; axis([-(T0+T1/2)-0.5,[T0+T1/2]+0.5,-0.5,1.5]); title(strcat(num2str(k),'次谐波叠加')); xlabel('t');end pause; plot(t,y(1:m+1,:));grid on; axis([-T0/2,T0/2,-0.5,1.5]); title('各次谐波叠加');xlabel('t'); ②在命令窗口调用rectexpd ()函数: T1=5; T0=10; m=5;rectexpd(T1,T0,m) ⑵Matlab 仿真结果:-T0-T1/2T1/2T00.20.40.60.81矩形信号串-20-15-10-551015200.050.10.150.20.250.30.350.40.450.5-10-50510-0.500.511.51次谐波叠加t-10-50510-0.500.511.53次谐波叠加t-10-50510-0.50.511.55次谐波叠加t-10-50510-0.50.511.52次谐波叠加t-10-50510-0.50.511.54次谐波叠加t-5-4-3-2-1012345-0.500.511.5各次谐波叠加t2、利用Matlab 实现连续信号卷积运算 ⑴ Matlab 程序代码:① 连续信号卷积运算的通用函数sconv (): function [f,k]=sconv(f1,f2,k1,k2,p) f=conv(f1,f2);f=f*p; k0=k1(1)+k2(1);k3=length(f1)+length(f2)-2; k=k0:p:k3*p;subplot(2,2,1); plot(k1,f1); title('f1(t)'); xlabel('t'); ylabel('f1(t)');subplot(2,2,2);plot(k2,f2);title('f2(t)'); xlabel('t'); ylabel('f2(t)'); subplot(2,2,3); plot(k,f); title('f(t)=f1(t)*f2(t)'); xlabel('t'); ylabel('f(t)');h=get(gca,'position'); h(3)=2.5*h(3);set(gca,'position',h);②在命令窗口调用sconv ()函数: p=0.01; k1=0:p:2; f1=exp(-k1); k2=0:p:3;f2=ones(1,length(k2)); [f,k]=sconv(f1,f2,k1,k2,p)⑵Matlab 仿真结果:0.511.5200.51f1(t)tf 1(t )01230.511.52f2(t)tf 2(t )00.51 1.522.533.544.550.51f(t)=f1(t)*f2(t)tf (t )实验二1、用Matlab模拟图形A律解码⑴Matlab程序代码:①建立函数ADecode()function y=ADecode(code,n) codesize=size(code);cr_len=codesize(1);cl_len=codesize(2);ca=zeros(1,cl_len-1);for i=1:cr_lenca=code(i,2:n);s=0;for j=1:n-1s=s+ca(j)*2^(n-1-j);enda=code(i,1);y(i)=s*((-1)^(a+1)); endy=y/(2^(n-1));A=87.6;A1=1+log(A); for j=1:length(y)if(y(j)>=0)if(y(i)<=1/A1)y(j)=y(j)*A1/A;elsey(j)=exp(y(j)*A1-1/A);endelsetemp=-y(j);if(temp<=1/A1)y(j)=-temp*A1/A;elsey(j)=-exp(temp*A1-1)/A;endendend②建立函数APCM()function code=APCM(x,n)xmax=max(abs(x));x=x/xmax;xlen=length(x);y=zeros(1,xlen);A=87.6;A1=1+log(A);for i=1:xlenif x(i)>=0if x(i)<=1/Ay(i)=(A*x(i))/A1;elsey(i)=(1+log(A*x(i)))/A1;endelsex1=-x(i);if x1<=1/Ay(i)=-(A*x1)/A1;elsey(i)=-(1+log(A*x1))/A1;endendendy1=y*(2^(n-1)-1);y1=round(y1);code=zeros(length(y1),n); c2=zeros(1,n-1);for i=1:length(y1)if(y1(i)>0)c1=1;elsec1=0;y1(i)=-y1(i);endfor j=1:n-1r=rem(y1(i),2);y1(i)=(y1(i)-r)/2;c2(j)=r;endc2=fliplr(c2);code(i,:)=[c1 c2];end③在新函数中调用前两个函数 t=0:0.01:1; x=sin(2*pi*t); code=APCM(x,7); y=ADecode(code,7); subplot(2,1,1) plot(t,x);title('原函数的图形'); subplot(2,1,2) plot(t,y);title('解码后函数的图形');⑵Matlab 仿真结果:00.10.20.30.40.50.60.70.80.91-1-0.500.51原函数的图形00.10.20.30.40.50.60.70.80.91-1-0.50.5解码后函数的图形2、用Matlab 模拟DSB 调制及解调过程 ⑴Matlab 程序代码: close all; clear all;dt=0.001;%采样时间间隔 fm=1; %信源最高频率 fc=10; %载波中心频率 N=4096;T=N*dt; t=0:dt:T-dt;mt=sqrt(2)*cos(2*pi*fm*t);%信源s_dsb=mt.*cos(2*pi*fc*t); %DSB-SC 双边带抑制载波调幅B=2*fm;figure(1); subplot(311) plot(t,s_dsb,'b-');hold on; %画出DSB 信号波形plot(t,mt,'r--');%画出m(t)信号波形 title('DSB 调制信号'); xlabel('t'); grid on;%DSB demodulation rt=s_dsb.*cos(2*pi*fc*t); rt=rt-mean(rt); [f,rf]=T2F(t,rt); [t,rt]=lpf(f,rf,B); subplot(312) plot(t,rt,'k-'); hold on;plot(t,mt/2,'r--');title('相干解调后的信号波形与输入信号的比较'); xlabel('t'); grid on; subplot(313)[f,sf]=T2F(t,s_dsb); %求调制信号的频谱psf=(abs(sf).^2)/T; %求调制信号的功率谱密度 plot(f,psf);axis([-2*fc 2*fc 0 max(psf)]); title('DSB 信号功率谱'); xlabel('f'); grid on;⑵Matlab 仿真结果:00.51 1.52 2.53 3.54 4.5-22DSB 调制信号t0.511.522.533.544.5-101相干解调后的信号波形与输入信号的比较t-20-15-10-5051015200.5DSB 信号功率谱f实验三1、用Matlab模拟双极性归零码⑴ Matlab程序代码:function y=drz(x)t0=300;x=[1 0 0 1 1 0 0 0 0 1 0 1];t=0:1/t0:length(x);for i=1:length(x);if(x(i)==1)for j=1:t0/2y(t0/2*(2*i-2)+j)=1;y(t0/2*(2*i-1)+j)=0;endelsefor j=1:t0/2y(t0/2*(2*i-2)+j)=-1;y(t0/2*(2*i-1)+j)=0;endendendy=[y,x(i)];M=max(y);m=min(y); subplot(211)plot(t,y);grid on;axis([0,i,m-0.1,M+0.1]);title('1 0 0 1 1 0 0 0 0 1 0 1 ');⑵Matlab仿真结果:2、用Matlab产生2FSK信号⑴ Matlab程序代码:①建立函数fskdigital()function fskdigital(s,f1,f2) t=0:2*pi/99:2*pi;m1=[];c1=[];b1=[];for n=1:length(s)if s(n)==0;m=ones(1,100);c=sin(f2*t);b=zeros(1,100) else s(n)==1;m=ones(1,100);c=sin(f1*t);b=ones(1,100)endm1=[m1 m]; c1=[c1 c]; b1=[b1 b]; endfsk=c1.*m1;subplot(211);plot(b1,'r')title('原始信号');axis([0 100*length(s) -0.1 1.1]);grid on;subplot(212);plot(fsk)title('2FSK信号');grid on;②在命令窗口调用函数fskdigital() >> s=[1 0 1 1 0 0 1 0]; f1=200; f2=100;fskdigital(s,f1,f2)⑵Matlab 仿真结果:10020030040050060070080000.20.40.60.81原始信号100200300400500600700800-1-0.500.512FSK 信号3、用Matlab 的simulink 模块模拟三角波的分解与还原 ⑴模块图⑵仿真结果:①三角波原图②用矩形波近似的冲击函数③分解后的三角波④经低通还原的三角波。