现代通信原理实验---模拟调制的MATLAB实现

画出频谱、功率谱密度图。

dt=0.001;

fmax=1;

fc=10;

T=5;

N=T/dt;

t=[0:N-1]*dt;

mt=sqrt(2)*cos(2*pi*fmax*t);

A=0;

s_ssb=real(hilbert(mt).*exp(j*2*pi*fc*t));

[f,Xf]=FFT_SHIFT(t,s_ssb);

PSD=(abs(Xf).^2)/T;

figure(1)

subplot(211);

plot(t,s_ssb);hold on ;

title('SSB 调制信号');

subplot(212);

plot(f,PSD);

axis([-2*fc 2*fc 0 1.5*max(PSD)]);

title('SSB 信号功率谱');

xlabel('f');

xlabel('f');

00.20.40.60.81 1.2 1.4 1.6 1.82

-2-1

1

2

SSB 调制信号

-20-15-10-50

51015200

1

2

3

SSB 信号功率谱

f

画出频谱、功率谱密度图。

dt=0.001; %时间采样频谱

fmax=1; %信源最高频谱

fc=10; %载波中心频率

T=5; %信号时长

N=T/dt;

t=[0:N-1]*dt;

mt=sqrt(2)*cos(2*pi*fmax*t); %信源

A=0;

s_dsb=mt.*cos(2*pi*fc*t);

[f,Xf]=FFT_SHIFT(t,s_dsb); %调制信号频谱

PSD=(abs(Xf).^2)/T; %调制信号功率谱密度

figure(1)

subplot(211);

plot(t,s_dsb);hold on; %画出DSB 信号波形

plot(t,A+mt,'r--'); %表示DSB 包络

plot(t,-A-mt,'r--');

title('DSB 调制信号及其包络');

xlabel('t');

subplot(212); %画出功率谱图形

plot(f,PSD);

axis([-2*fc 2*fc 0 1.5*max(PSD)]);

title('DSB 信号功率谱');

xlabel('f');

xlabel('f');

00.51 1.52 2.53

3.54

4.55

-2-1

1

2

DSB 调制信号

及其包络

t

-20-15-10-50

51015200

0.2

0.4

0.6

0.8

DSB 信号功率谱f

画出频谱、功率谱密度图。

程序：dt=0.001

T=5

N=T/dt

T=[0:N-1]*dt

mt=sqrt(2)*cos(2*pi*fmax*t)

s_am=10-cos(10*t(2*pi)*(5/0.3)*sin(0.3*t)) [f,Xf]=fft_shift(t,s_am)

PSD=(abs(Xf).^2)/T

figure(1)

Subplot(211)

Plot(t,s_am)

hold on

Plot(t,mt,’r--’)

title(‘FM调制信号及包络’)

xlabel(‘t’)

subplot(212)

plot(f,PSD)

axis([-2*fc 2*fc 0 1.5*max(PSD)])

tltle(‘FM信号功率谱’)

xlabel(‘f’)

基于Matlab的FM仿真实现

摘要 本次设计主要是以Matlab为基础平台，对FM信号进行仿真。介绍了FM信号，及其调制和解调的基本原理，并设计M文件，分析在混入噪声环境下的波形失真，以及分析FM的抗噪声性能。本设计的主要目的是对Matlab的熟悉和对模拟通信理论的更深化理解。 关键词：Matlab；FM；噪声

前言 (2) 1 设计基础 (3) 1.1 Matlab及M文件的简介 (3) 1.2模拟调制概述 (4) 1.2.1模拟调制系统各个环节分析 (5) 1.2.2 模拟调制的意义 (6) 2 FM基本原理与实现 (7) 2.1 FM的基本原理 (7) 2.1.1调制 (7) 2.1.2解调 (8) 2.2 FM的实现 (8) 2.2.1 FM调制的实现 (8) 2.2.2 FM解调的实现 (9) 2.3 调频系统的抗噪声性能 (10) 2.3.1 高斯白噪声信道特性 (10) 3 FM的仿真实现与分析 (14) 3.1 未加噪声的FM解调实现 (14) 3.2 叠加噪声时的 FM解调 (16) 总结 (20) 致谢 (21) 参考文献 (22) 附录 (23)

通信按照传统的理解就是信息的传输。在当今高度信息化的社会，信息和通信已成为现代社会的命脉。信息作为一种资源，只有通过广泛传播与交流，才能产生利用价值，促进社会成员之间的合作，推动社会生产力的发展，创造出巨大的经济效益。而通信作为传输信息的手段或方式，与传感技术、计算机技术相融合，已成为21世纪国际社会和世界经济发展的强大动力。可以预见，未来的通信对人们的生活方式和社会的发展将会产生更加重大和意义深远的影响。 在通信系统中，从消息变换过来的原始信号所占的有效频带往往具有频率较低的频谱分量(例如语音信号)，如果将这种信号直接在信道中进行传输，则会严重影响信息传送的有效性和可靠性，因此这种信号在许多信道中均是不适宜直接进行传输的。在通信系统的发射端通常需要有调制过程，将调制信号的频谱搬移到所希望的位置上，使之转换成适于信道传输或便于信道多路复用的已调信号；而在接收端则需要有解调过程，以恢复原来有用的信号。调制解调方式常常决定了一个通信系统的性能。随着数字化波形测量技术和计算机技术的发展，可以使用数字化方法实现调制与解调过程。 调制在通信系统中具有重要的作用。通过调制，不仅可以进行频谱搬移，把调制信号的频谱搬移到所希望的位置上，从而将调制信号转换成适合于信道传输或便于信道多路复用的已调信号，而且它对系统的传输有效性和传输可靠性有着很大的影响。调制方式往往决定了一个通信系统的性能。调制技术是指把基带信号变换成传输信号的技术。基带信号是原始的电信号，一般是指基本的信号波形，在数字通信中则指相应的电脉冲。在无线遥测遥控系统和无线电技术中调制就是用基带信号控制高频载波的参数（振幅、频率和相位），使这些参数随基带信号变化。用来控制高频载波参数的基带信号称为调制信号。未调制的高频电振荡称为载波（可以是正弦波，也可以是非正弦波，如方波、脉冲序列等）。被调制信号调制过的高频电振荡称为已调波或已调信号。已调信号通过信道传送到接收端，在接收端经解调后恢复成原始基带信号。

基于MATLAB的模拟调制系统仿真与测试(AM调制)

闽江学院 《通信原理设计报告》 题目：基于MATLAB的模拟调制系统仿真与测试学院：计算机科学系 专业：12通信工程 组长：曾锴（3121102220） 组员：薛兰兰（3121102236） 项施旭（3121102222） 施敏（3121102121） 杨帆（3121102106） 冯铭坚（3121102230） 叶少群（3121102203） 张浩（3121102226） 指导教师：余根坚 日期：2014年12月29日——2015年1月4日

摘要在通信技术的发展中，通信系统的仿真是一个重点技术，通过调制能够将信号转化成适用于无线信道传输的信号。 在模拟调制系统中最常用最重要的调制方式是用正弦波作为载波的幅度调制和角度调制。在幅度调制中，文中以调幅、双边带和单边带调制为研究对象，从原理等方面阐述并进行仿真分析；在角度调制中，以常用的调频和调相为研究对象，说明其调制原理，并进行仿真分析。利用MATLAB下的Simulink工具箱对模拟调制系统进行仿真，并对仿真结果进行时域及频域分析，比较各个调制方式的优缺点，从而更深入地掌握模拟调制系统的相关知识，通过研究发现调制方式的选取通常决定了一个通信系统的性能。 关键词模拟调制；仿真；Simulink 目录 第一章绪论 (1) 1.1 引言 (1) 1.2 关键技术 (1) 1.3 研究目的及意义 (2) 1.4 本文工作及内容安排 (2) 第二章模拟调制原理 (3) 2.1 幅度调制原理 (3) 2.1.1 AM调制 (4) 第三章基于Simulink的模拟调制系统仿真与分析 (6) 3.1 Simulink工具箱简介 (6) 3.2 幅度调制解调仿真与分析 (8) 3.2.1 AM调制解调仿真及分析 (8) 第四章总结 (12) 4.1 代码 (13) 4.2 总结 (14)

matlab模拟

1.1.1模拟技术 在许多数学方法中，一般都要用到解析论证和数值计算的技巧。但是，许多现实的系统是很复杂的，其中的随机性因素往往难以用数学公式表示出来，也就很难使用数学推导或数值计算机的手段来分析系统、预测系统的性能。因此，产生了对系统进行模拟的技术。 在使用计算机对系统进行模拟之前，一般要清楚模拟的一般步骤和方法。后面将从较小的模拟实例，对模拟技术进行简要的介绍。 模拟过程的一般过程为： （1）分析问题，收集资料。需要搞清楚问题要达到的目标，根据问题的性质收集有关随机性因素的资料。这里用得较多的知识为概率统计方面。 在这个阶段，还应当估计一下待建立的模拟系统的规模和条件，说明 哪些是可以控制的变量，哪些是不可控制的变量。 （2）建立模拟模型，编制模拟程序。按照一般的建模方法，对问题进行适当的假设。也就是说，模拟模型未必要将被模拟系统的每个细节全部 考虑。模拟模型的优劣将通过与实际系统有关资料的比较来评价。如 果一个“粗糙”的模拟模型已经比较符合实际系统的情况，也就没有 必要建立费时、复杂的模型。当然，如果开始建立的模型比较简单， 与实际系统相差较大，那么可以在建立了简单模型后，逐步加入一些 原先没有考虑的因素，直到模型达到预定的要求为止。编写模拟程序 之前，要现画出程序框图或写出算法步骤。然后选择合适的计算机语 言，编写模拟程序。模拟实现的工具较多，如数学软件类：Matlab、 Mathematica、Maple、MathCAD等，还有其它的高级语言工具，如： Visual C++、C++ Builder、Delphi、Borland C++3.1等。 （3）运行模拟程序，计算结果。为了减小模拟结果的随机性偏差，一般要多次运行模拟程序，还有就是增加模拟模型的时段次数。 （4）分析模拟结果，并检验。模拟结果一般说来反映的是统计特性，结果的合理性、有效性，都需要结合实际的系统来分析，检验。以便提出 合理的对策、方案。 以上步骤是一个反复的过程，在时间和步骤上是彼此交错的。比如模型的修改和改进，都需要重新编写和改动模拟程序。模拟结果的不合理，则要求检查模型，并修改模拟程序。 1.1.2模拟时间 利用计算机进行模拟时，有两种控制模拟时间的方法。一种是固定时间增量法，另一种是可变时间增量法，又叫面向事件法。 固定时间增量法，是选用一段合适的时间作单位，然后每隔一个单位时间就计算一次有关参数的值，到达预定的模拟时间后，模拟程序结束。在编写这种程序时，一般可以建立一个“模拟时钟”变量。程序的主体框架一般时个大的循环，循环变量，则为模拟时间；在每个循环体内，就是对每个时段作处理。例如，有些排队论模型，可能就是以每隔一段时间（一天或者一个月）进行处理。 采用可变时间增量法编写的模拟程序也有一个“模拟时钟”变量，但它是在一个事件发生时，“模拟时钟”才向前推进。需要注意的是，该模拟方法每一步

基于MATLAB的模拟信号频率调制(FM)与解调分析

课程设计任务书 学生姓名：杨刚专业班级：电信1302 指导教师：工作单位：武汉理工大学 题目：信号分析处理课程设计 －基于MATLAB的模拟信号频率调制（FM）与解调分析 初始条件： 1.Matlab6.5以上版本软件； 2.先修课程：通信原理等； 要求完成的主要任务：（包括课程设计工作量及其技术要求，以及说明书撰写等具体要求） 1、利用MATLAB中的simulink工具箱中的模块进行模拟频率（FM）调制与解调，观 察波形变化 2、画出程序设计框图，编写程序代码，上机运行调试程序，记录实验结果（含计算结 果和图表等），并对实验结果进行分析和总结； 3、课程设计说明书按学校统一规范来撰写，具体包括： ⑴目录；⑵理论分析； ⑶程序设计；⑷程序运行结果及图表分析和总结； ⑸课程设计的心得体会（至少800字，必须手写。）； ⑹参考文献（不少于5篇）。 时间安排： 周一、周二查阅资料，了解设计内容； 周三、周四程序设计，上机调试程序； 周五、整理实验结果，撰写课程设计说明书。 指导教师签名： 2013 年 7月 2 日 系主任（或责任教师）签名： 2013年 7月 2日

目录 1 Simulink简介 (1) 1.1 Matlab简介······················································错误！未定义书签。 1.2 Simulink介绍 ···················································错误！未定义书签。 2 原理分析 ·····························································错误！未定义书签。 2.1通信系统 ·························································错误！未定义书签。 2.1.1通信系统的一般模型 ···································错误！未定义书签。 2.1.2 模拟通信系统 (3) 2.2 FM调制与解调原理···········································错误！未定义书签。 3 基于Matlab方案设计 (6) 3.1 Matlab代码 (6) 3.2 Matlab仿真 (8) 4 基于Simulink方案设计 (12) 4.1 使用Simulink建模和仿真的过程 (12) 4.1.1 Simulink模块库简介 (12) 4.1.2 调制解调模块库简介 (13) 4.2 FM调制与解调电路及仿真 (14) 4.3 仿真结果分析 (17) 5 心得体会 ·····························································错误！未定义书签。 6 参考文献 (20) 本科生课程设计评定表

MATLAB编程用两种方法模拟光学实验

MATLAB编程用两种方法模拟光学实验 摘要： 利用MATLAB软件编程实现了用衍射积分的方法对单缝衍射、杨氏双缝干涉、黑白 光栅衍射的计算机模拟；以及用傅立叶变换方法对简单孔径衍射、黑白光栅及正弦光栅夫 琅和费衍射的模拟。 关键词： MATLAB；衍射积分；傅立叶变换；计算机模拟 引言： 美国Mathworks公司推出的MA TLAB，是一种集数值计算、符号预算、可视化建模、 仿真和图形处理等多种功能于一体的优秀图形化软件。本文介绍了通过MA TLAB软件编 程实现用衍射积分和傅立叶变换实现夫琅和费衍射计算机模拟的方法。 计算机模拟为衍射实验的验证提供一条简捷、直观的途径。从而加深了对物理原理、 概念和图像的理解。 正文： 大学教学课程中引入计算机模拟技术正日益受到重视，与Basic、C和Fortran相比，用MA TLAB软件做光学试验的模拟，只需要用数学方式表达和描述，省去了大量繁琐的编 程过程。下面来介绍利用MATLAB进行光学模拟的两种方法。 （一）衍射积分方法： 该方法首先是由衍射积分算出接收屏上的光强分布，然后根据该分布调制色彩作图，从而得到衍射图案。 1．单缝衍射。 把单缝看作是np个分立的相干光源，屏幕上任意一点复振幅为np个光源照射结果 的合成，对每个光源，光程差Δ=ypsinΦ,sinΦ=ys/D,光强I=I0(Σcosα)2+(Σsinα)2，其中α=2Δ/λ=πypys/λD 编写程序如下，得到图1 lam=500e-9; a=1e-3;D=1; ym=3*lam*D/a; ny=51; ys=linspace(-ym,ym,ny); np=51; yp=linspace(0,a,np); for i=1:ny sinphi=ys(i)/D; alpha=2*pi*yp*sinphi/lam; 图1 单缝衍射的光强分布 sumcos=sum(cos(alpha)); sumsin=sum(sin(alpha)); B(i,:)=(sumcos^2+sumsin^2)/np^2; end N=255; Br=(B/max(B))*N; subplot(1,2,1)

基于MATLAB模拟调制系统的仿真设计

1 线性模拟调制 1.1模拟调制原理 模拟调制是指用来自信源的基带模拟信号去调制某个载波，而载波是一个确知的周期性波形。模拟调制可分为线性调制和非线性调制，本文主要研究线性调制。 线性调制的原理模型如图1.1所示。设c(t)=Acos2t f o π,调制信号为m(t),已调信号为s(t)。 图1.1 线性调制的远离模型 调制信号m(t)和载波在乘法器中相乘的结果为：t A t m t s w o cos )()('=，然后通过一个传输函数为H(f)的带通滤波器，得出已调信号为。 从图1.1中可得已调信号的时域和频域表达式为： （1-1） 式（1-1）中，M(f)为调制信号m(t)的频谱。 由于调制信号m(t)和乘法器输出信号之间是线性关系，所以成为线性调制。带通滤波器H(f)可以有不同的设计，从而得到不同的调制种类。 1.2双边带调制DSB 的基本原理 在幅度调制的一般模型中，若假设滤波器为全通网络，调制信号m(t)中无直流分量，则输出的已调信号就是无载波分量的双边带调制信号，或称抑制载波双边带（DSB ）调制信号，简称双边带（DSB ）信号。 设正弦型载波c(t)=Acos( t) ，式中：A 为载波幅度， 为载波角频率。 根据调制定义，幅度调制信号（已调信号）一般可表示为： (t)=Am(t)cos(t) （1-2） ?? ???-++==) ()]()([21)()(*]cos )([)(f H f f M f f M f s t h t t m t s o o o w m(t) H(t) A os t w o c s(t) )(' t s

其中，m(t)为基带调制信号。 设调制信号m(t)的频谱为M()，则由公式2-2不难得到已调信号 (t)的频谱： )]()([2 )(c c m M M A s ωωωωω-++= （1-3） 由以上表示式可见，在波形上，幅度已调信号随基带信号的规律呈正比地变化；在频谱结构上，它的频谱完全是基带信号频谱在频域的简单搬移。 标准振幅就是常规双边带调制，简称调幅（AM ）。假设调制信号m(t)的平均值为0，将其叠加一个直流偏量 后与载波相乘，即可形成调幅信号。其时域表达式为: )cos())(()(0t t m t c AM A s ω+= （1-4） 式中： 为外加的直流分量；m(t)可以是确知信号，也可以是随机信号。 若为确知信号，则AM 信号的频谱为: （1-5） AM 信号的频谱由载频分量、上边带、下边带三部分组成。AM 信号的总功率包括载波功率和边带功率两部分。只有边带功率才与调制信号有关，也就是说，载波分量并不携带信息。因此，AM 信号的功率利用率比较低。 AM 调制器模型如下图所示。 图1.2 AM 调制器模型 AM 信号的时域和频域表达式分别为 （1-6） （1-7） 式中，A o 为外加的直流分量；m(t)可以是确知信号也可以是随机信号，但通常认为其平均值为0，即0)(=t m — 。 由频谱可以看出，AM 信号的频谱由载波分量、上边带、下边带三部分组成。上边带的频谱结构与原调制信号的频谱结构相同，下边带是上边带的镜像。因此，AM 信号是带有载波 分量的双边带信号，他的带宽是基带信号带宽 的2倍，即 ) (cos )()(cos ) (cos )]([)(t w c t m t w c A t w c t m A o t s o AM +=+=)]()([2 1)]()([)(w c w M w c w M w c w w c w A o t s AM -+++-++=δδπ)] ()([2 1)]()([)(0 ω ω ω ω ωωωδωδπωc c c c m M M A s -+++-++=f H

基于MATLAB的各类混沌系统的计算机模拟(教学版)

基于MATLAB 的各类混沌系统的计算机模拟 ―――《混沌实验教学平台的设计与实现》初期报告 物电05级1A 班 张丹伟 20050003101 摘要：本文利用数学软件MATLAB 对Lorenz 系统等六个重要的混沌模型进行数值计算，同 时模拟出各类混沌系统的独特性质，如混沌吸引子，倍周期，初值敏感性，相图，分岔图等。通过观察和分析上述特性，加深了我们对混沌现象的理解。 关键词：混沌； 微分方程； MA TLAB ； 引言． 混沌探秘 混沌是非线性系统所独有且广泛存在的一种非周期运动形式, 其覆盖面涉及到自然科学和社会科学的几乎每一个分支。1972年12月29日，美国麻省理工学院教授、混沌学开创人之一E.N.洛伦兹在美国科学发展学会第139次会议上发表了题为《蝴蝶效应》的论文，提出一个貌似荒谬的论断：在巴西一只蝴蝶翅膀的拍打能在美国得克萨斯州产生一个龙卷风，并由此提出了天气的不可准确预报性。为什么会出现这种情况呢？这是混沌在作怪！ “混沌”译自英语中“chaos”一词，原意是混乱、无序，在现代非线性理论中，混沌则是泛指在确定体系中出现的貌似无规则的、类随机的运动。 混沌现象是普遍的，就在我们身边，是与我们关系最密切的现象，我们就生活在混沌的海洋中。一支燃着的香烟，在平稳的气流中缓缓升起一缕青烟，突然卷成一团团剧烈搅动的烟雾，向四方飘散；打开水龙头，先是平稳的层流，然后水花四溅，流动变的不规则，这就是湍流；一个风和日丽的夏天，突然风起云涌，来了一场暴风雨。一面旗帜在风中飘扬，一片秋叶从树上落下，它们都在做混沌运动。可见混沌始终围绕在我们的周围，一直与人类为伴。 一． 混沌的基本概念 1. 混沌: 目前尚无通用的严格的定义, 一般认为,将不是由随机性外因引起的, 而是由确定性方程(内因)直接得到的具有随机性的运动状态称为混沌。 2. 相空间: 在连续动力系统中, 用一组一阶微分方程描述运动, 以状态变量(或状态向量)为坐标轴的空间构成系统的相空间。系统的一个状态用相空间的一个点表示, 通过该点有唯一的一条积分曲线。 3. 混沌运动: 是确定性系统中局限于有限相空间的高度不稳定的运动。所谓轨道高度不稳定, 是指近邻的轨道随时间的发展会指数地分离。由于这种不稳定性, 系统的长时间行为会显示出某种混乱性。 4. 分形和分维: 分形是 n 维空间一个点集的一种几何性质, 该点集具有无限精细的结构, 在任何尺度下都有自相似部分和整体相似性质, 具有小于所在空间维数 n 的非整数维数。分维就是用非整数维——分数维来定量地描述分形的基本性质。 5. 不动点: 又称平衡点、定态。不动点是系统状态变量所取的一组值, 对于这些值系统不随时间变化。在连续动力学系统中, 相空间中有一个点0x , 若满足当 t →∞时, 轨迹0()x t x →, 则称0x 为不动点。

实验一 模拟通信的MATLAB仿真

实验一 模拟通信的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 0 +m ()t ]cos t c ω （2-1） ()t S AM =A 0 π[()()ωωωωδC C ++-]+()()[]ωωωωc c M M ++-2 1 （2-2）

AM 信号的带宽 2 =B AM f H （2-3） 式中， f H 为调制信号的最高频率。 2.1.3 AM 信号的功率P AM 与调制效率 η AM P AM =()222 2 t m A +=P P m c + （2-4） 式中，P C =2 A 为不携带信息的载波功率；()2 2 t m P m =为携带信息的边带 功率。 ()() t t m A m P P AM C AM 2 2 2+= = η （2-5） AM 调制的优点是可用包络检波法解调，不需要本地同步载波信号，设备简单。AM 调制的最大缺点是调制效率低。 2.2、双边带调制（DSB ） 如果将在AM 信号中载波抑制，只需在图3-1中将直流 A 0 去掉，即可输出 抑制载波双边带信号。 2.2.1 DSB 信号的时域和频域表达式 ()()t t m t c DSB S ωcos = （2-6） ()()()[]ωωωωωC C DSB M M S ++-=2 1 （2-7） DSB 信号的带宽 f B B H AM DSB 2 == （2-8）

基于MATLAB的数字模拟仿真..

基于MATLAB的数字模拟仿真 摘要：本文阐述了计算机模拟仿真在解决实际问题时的重要性，并较为系统的介绍了使用计算机仿真的原理及方法。对于计算机模拟仿真的三大类方法：蒙特卡罗法、连续系统模拟和离散事件系统模拟，在本文中均给出了与之对应的实例及基于MATLAB模拟仿真的相关程序，并通过实例深入的分析了计算机模拟解决实际问题的优势及不足。 关键词：计算机模拟；仿真原理；数学模型；蒙特卡罗法；连续系统模拟；离散事件系统模拟 在实际问题中，我们通常会面对一些带随机因素的复杂系统，用分析方法建模常常需要作许多简化假设，这样进行处理过后的模型与我们面临的实际问题可能相差很远，以致求解得到答案根本无法应用，这时，计算机模拟几乎成为唯一的选择。本文通过对计算机模拟仿真进行系统地介绍，寻求利用模拟仿真来解决问题的一般方法，并深入探讨了这些方法的长处和不足。我们定义一些具有特定的功能、相互之间以一定的规律联系的对象所组成的总体为一个系统，模拟就是利用物理的、数学的模型以系统为问题解决对象，来类比、模仿现实系统及其演变过程，以寻求过程规律的一种方法。模拟的基本思想是建立一个实验的模型，这个模型包含所研究系统的主要特点，这样做的目的就是通过对这个实验模型的运行，获得所要研究系统的必要信息。另外，系统的运行离不开算法，仿真算法是将系统模型转换成仿真模型的一类算法，在数字仿真模型中起核心和关键作用。 1、所谓计算机仿真 计算机仿真是利用计算机对一个实际系统的结构和行为进行动态演示，以评价或预测该系统的行为效果。它是解决较复杂的实际问题的一条有效途径。针对一个确定的系统，根据运行的相似原理，利用计算机来逼真模仿研究对象（研究对象可以是真实的系统，也可以是设想中的系统），计算机仿真是将研究对象进行数学描述，建模编程，且在计算机中运行实现。 对比于物理模拟通常花费较大、周期较长，且在物理模型上改变系统结构和系数都较困难的诸多缺陷，计算机模拟不怕破坏、易修改、可重用，有更强的系统适应能力。但是计算机模拟也有缺陷，比如受限于系统建模技术，即系统数学模型不易建立、程序调试复杂等。 计算机仿真可以用于研制产品或设计系统的全过程中，包括方案论证、技术指标确定、设计分析、生产制造、试验测试、维护训练、故障处理等各个阶段。 2、计算机仿真的目的 对于一个系统，是否选择进行计算机模拟的问题，基于判断计算机模拟与非计算机模拟方法孰优孰劣的问题。归纳以下运用计算机模拟的情况： （1）在一个实际系统还没有建立起来之前，要对系统的行为或结果进行分析研究时，计算机仿真是一种行之有效的方法。 （2）在有些真实系统上做实验会影响系统的正常运行，这时进行计算机模拟就是为了避免给实际系统带来不必要的损失。如在生产中任意改变工艺参数可能会导致废品，在经济活动中随意将一个决策付诸行动可能会引起经济混乱。 （3）当人是系统的一部分时，他的行为往往会影响实验的效果，这时运用系统进行仿真研究，就是为了排除人的主观因素的影响。

MATLAB仿真 BPSK调制

matlab BPSK 调制与解调 1、调制 clear all; g=[1 0 1 0 1 0 0 1];%基带信号 f=100; %载波频率 t=0:2*pi/99:2*pi; cp=[];sp=[]; mod=[];mod1=[];bit=[]; for n=1:length(g); if g(n)==0; die=-ones(1,100); %Modulante se=zeros(1,100); % else g(n)==1; die=ones(1,100); %Modulante se=ones(1,100); % end c=sin(f*t); cp=[cp die]; mod=[mod c]; bit=[bit se]; end bpsk=cp.*mod; subplot(2,1,1);plot(bit,'LineWidth',1.5);grid on; title('Binary Signal'); axis([0 100*length(g) -2.5 2.5]); subplot(2,1,2);plot(bpsk,'LineWidth',1.5);grid on; title('ASK modulation'); axis([0 100*length(g) -2.5 2.5]); 2、调制解调加噪声 clc; close all; clear; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % % 假定：

% 2倍载波频率采样的bpsk信号 % 调制速率为在波频率的 N/2m % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% m=128; N=512; n=1:1:N; N0=0.5*randn(1,N) %噪声 h0=zeros(1,N); % 30阶低通滤波器 h0 f = [0 0.3 0.3 1]; w0 = [1 1 0 0]; b = fir2(30,f,w0); [h,w] = freqz(b,1,N/2); h0(1,1:N/2)=abs(h'); for i=1:N/2 h0(1,N-i+1)=h0(1,i); end; %%%%%%%%% 随机序列 a=rand(1,m); for i=1:m if(a(1,i)>0.5) a(1,i)=1; else a(1,i)=-1; end; end; %%% 生成BPSK信号 bpsk_m=zeros(1,N); j=1;k=1; for i=1:N if(j==(N/m+1)) j=1; k=k+1; end; % 0.05*pi 为初始相位，可以任意改变 bpsk_m(1,i)=a(1,k)*sin(2*pi*0.5*i+0.05*pi)+a(1,k)*cos(2*pi*0.5*i+ 0.05*pi); j=j+1; end; bpsk_m=bpsk_m+N0;% 信号加噪声，模拟过信道 % 接收处理用正交本振与信号相乘，变频 bpsk_m1=bpsk_m.*sin(2*pi*0.5*n); bpsk_m2=bpsk_m.*cos(2*pi*0.5*n); %滤波 tempx=fft(bpsk_m1);

matlabFM调制仿真

Matlab FM调制仿真

目录 引言.................................................................................. 一．课程设计的目的与要求 .............................................. 1.1课程设计的目的.................................................... 1.2课程设计的要求.................................................... 二．FM调制解调系统设计............................................... 2.1FM调制模型的建立............................................. 2.2调制过程分析........................................................ 2.3FM解调模型的建立............................................. 2.4解调过程分析........................................................ 2.5高斯白噪声信道特性 ............................................ 2.6调频系统的抗噪声性能分析 ................................ 三．仿真实现...................................................................... 3.1MATLAB源代码.................................................. 3.2仿真结果................................................................ 四．心得体会...................................................................... 五．参考文献...................................................................... 引言 本课程设计用于实现DSB信号的调制解调过程。信号的调制与解调在通信系统中具有重要的作用。调制过程是一个频谱搬移的过程，它是将低频信号的频谱搬移到载频位

我的基于MATLAB仿真的数字调制与解调设计

摘要：设计了二进制振幅键控(2ASK)、二进制移频键控(2FSK) 、二进制移相键控(2PSK)调制解调系统的工作流程图,并得用了MATLAB软件对该系统的动态进行了模拟仿真,得用仿真的结果,从而衡量数字信号的传输质量。（仿宋、小五号） 关键词：调制解调、2ASK、2FSK、2PSK、2DPSK、MATLAB（宋体、小五号） ABSTRACT(四号加粗居中放置)： The work stream diagrams of 2ASK、2FSK、2PSK are designed .MA TLAB softwave is used to simulate the modem system by the scatter diagrams and wave diagrams, then the transmit quality of digital signal can be measured.（小五号） Key word：Amodulate and ademodulate 、2ASK、2FSK、2PSK、2DPSK、MATLAB（小五号） （正文：宋体、五号 一级标题：黑体、四号，小标题上下空一行。） 一、数字调制解调相关原理 在通信系统中,信道的频段往往是很有限的,而原始的通信信号的频段与信道要求的频段是不匹配的,这就要求将原始信号进行调制再进行发送.相应的在接收端对调制的信号进行解调,恢复原始的信号,而且调制解调还可以在一定程度上抑制噪声对通信信号的干扰。 调制解调技术按照通信信号是模拟的还是数字的可分为模拟调制解调和数字调制解调。数字调制的基本方式可以归结为3类：振幅键控（ASK）、频移键控（FSK）和相移键控（PSK）。此外还有这3类的混合方式。 对于数字调制信号，为了提高系统的抗噪声性能，衡量系统性能的指标是误码率。1.1二进制振幅键控(2ASK) 振幅键控是正弦载波的幅度随数字基带信号而变化的数字调制。当数字基带信号为二进制时，则为二进制振幅键控。设发送的二进制符号序列由0,1序列组成，发送0符号的概率为P，发送1符号的概率为1-P，且相互独立.该二进制符号序列可表示为： 其中： Ts是二进制基带信号时间间隔，g(t)是持续时间为Ts的矩形脉冲， 为单极性不归零脉冲序列，则根据幅度调制的原理，一个二进制的振幅键控信号可以表示成一个单极性矩形脉冲序列与一个正弦型载波的相乘，即 2ASK信号的时间波形如果是通断方式，就称为通断键控信号(OOK信号)。 二进制振幅键控信号的产生可以采用数字键控的方法实现也可以采用模拟相乘的方法实现。2ASK信号与模拟调制中的AM信号类似。所以，对2ASK信号也能够采用非相干解调(包络检波法)和相干解调(同步检测法)，其相应原理方框图如图1.1所示。

利用MATLAB仿真模拟调制系统

利用MATLAB仿真模拟调制系统 MATLAB的名称源自Matrix Laboratory，专门以矩阵形式处理数据，是目前国际上流行的进行科学研究、工程计算的软件，广泛地应用于科学计算、控制系统、信息处理等领域的分析、仿真和设计工作中。MATLAB的出现使得通信系统的仿真能够用计算机模拟实现，只需要输入不同的参数就能得到不同情况下的系统性能，而且在结构的观测和数据的存储方面也比传统的方式有优势，因而MATLAB在通信仿真领域得到越来越多的应用。 本文中，我们对模拟调制系统、数字带通传输系统等列举了一些MATLAB仿真的实例，作为大家学习MATLAB的参考资料，让读者学会处理具体问题的建模编程方法，逐渐掌握MATLAB的通信系统仿真。 由本章的学习我们知道，各种信源所产生的基带信号并不能在大多数信道内直接传输，而是需要经调制后再送到信道中去。在接受端就必须通过相反的过程，即解调。本章中，我们以常规双边带调幅AM系统为例仿真模拟通信系统的各个过程。 我们假定信号频率为10Hz，载波频率为50Hz，采样率为1000Hz，信噪比SNR等于3。要求利用MATLAB软件仿真AM调制每一点的波形，包括信息信号、AM信号、载波信号、已调信号、通过带通滤波器后的信号，解调后的信号；并仿真AM信号频谱、已调信号频谱与解调信号频谱。 MATLAB程序如下：

% 标准调幅AM调制 a0=2;f0=10;fc=50;snr=3; fs=1000; % 变量定义 t=[-50:0.001:50]; am1=cos(2*pi*f0*t); % 产生信号频率为f0的基带信号 am=a0+am1; % 产生AM信号 c_am=cos(2*pi*fc*t); % 产生频率为fc的载波 AM_mod=am.*c_am; % 产生调制信号 am_f=fft(am); % AM频域 AM_modf=fft(AM_mod); y=awgn(AM_mod,snr); % 叠加噪声 figure(1); hold on; subplot(2,2,1); plot(t,am1); axis([0 0.4 -2 2]); title('基带信号波形'); % 绘图subplot(2,2,2); plot(t,am); axis([0 0.4 -2 6]); title('AM信号波形'); subplot(2,2,3); plot(t,c_am); axis([0 0.4 -2 2]); title('载波信号波形'); subplot(2,2,4); plot(t,AM_mod); axis([0 0.4 -8 8]); title('已调信号波形'); hold off; figure(2); hold on; subplot(2,2,1); plot(t,AM_mod); axis([0 0.4 -8 8]); title('已调信号波形'); subplot(2,2,2); plot(t,y); axis([0 0.4 -8 8]); title('叠加噪声后的信号波形');; a=[35,65];b=[30,70]; Wp=a/(fs/2);Ws=b/(fs/2);Rp=3; Rs=15; [N,Wn]= Buttord(Wp,Ws,Rp,Rs) ; % 计算巴特沃斯数字滤波器的阶数和 3db截止频率 [B,A]=Butter(N,Wn,'bandpass'); % 计算巴特沃斯模拟滤波器系统函数的分子、分母多项式系数向量 sig_bandpass=filtfilt(B,A,y); % 带通滤波后信号 subplot(2,2,3); plot(t,sig_bandpass); axis([0 0.4 -8 8]); title('经带通滤波后信号波形'); hold off; AM_dem=sig_bandpass.*c_am; Wp=15/(fs/2);Ws=40/(fs/2);Rp=3; Rs=20; [N,Wn]= Buttord(Wp,Ws,Rp,Rs) ; % 同上 [B,A]=Butter(N,Wn,'low'); AM_demod=filtfilt(B,A,AM_dem) % 低通滤波后信号 AM_demodf=fft(AM_demod); subplot(2,2,4); plot(t,AM_demod); axis([0 0.4 0 2]); title('解调信号波形'); hold off; f=(0:100000)*fs/100001-fs/2; figure(3); hold on;

基于MATLAB的交通流计算机模拟

基于MATLAB的交通流计算机模拟 摘要： 设计标准和各类出行，环境及社会的发展预测，对六车道的桥面的的交通流进行预测模拟，并实现其可视化，直观地了解未来桥面的车流模式，为评估和修正设计方案提供依据。 本作品利用MATLAB软件方便的技术方法来实现交通流的模拟和可视化，具有较强的可读性和可控制性。 1 引言 计算机模拟技术作为一门独立的学科始于20世纪40年代。70年代以来，随着系统科学与计算机科学技术的发展，模拟技术得到了迅猛的发展，已经广泛地应用几乎所有的学科。在交通运输系统的规划、设计、运营分析等方面的应用更是得到了长足发展，并在交通运输工程学科中形成了交通模拟这一崭新的领域。 交通模拟技术在分析、评价公路运输系统及其构成单元中起主要角色。它们通常与其他诸如供给-需求分析、通行能力分析、交通流模拟、跟车理论、波动理论等分析方法相结合来构造复杂的公路交通子系统，或一些子系统经过相互作用而组成的大系统的模拟框架。这些子系统可以是单个的信号交叉口、无信号交叉口、居民区或商业中心区的交通密集路网、线控或网控信号系统、高速公路、乡村双车道公路或多车道公路系统。 到现在为止，可以说交通问题的研究已经有三种方法——经验实测方法、理论分析方法、计算机模拟方法。最常用的方法是经验实测法。实测法的最大优点是基本数据都来源于实际现场，有限大的可信度，不需要什么假设条件。但是，其弱点是对于个别因素的影响情况很难确定。理论分析法，总是要采取一些基本

假设，这些假设受理论研究者水平的限制有些可能不正确，必定或多或少地与实际有些偏差。其优点是对于个别因素的影响有明确的数量关系表示。计算机模拟则间有以上两种方法的优点，由于计算机模拟模型是理论推演，抽象出来的，而一些基本数据则是来自现场实测，而且利用计算机模拟方法能产生很多像实测法那样得到的交通数据。 MATLAB 是一种用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，拥有友好的工作平台和编程环境，简单易用的程序语言，强大的科学计算机数据处理能力，出色的图形处理功能，应用广泛的模块集合工具箱，实用的程序接口和发布平台等诸多优点，方便交通流模拟的实现和可视化，大大减轻程序语言的复杂程度。 2 六车道交通流的运行模型 建立交通流模型的根本目的是要以足够的精度来再现客观的交通现象，在进行交通流的微观模拟的过程中，根据交通流分布的一般规律、交通调查和基本经验为基础，全面构造车辆的到达，车流分布，车辆类型，车速，在路段上自由行驶，跟驰等。这里的交通模型主要包括两部分：一是车辆的产生模型；二是车辆的行驶模型。 车辆产生模型就是车辆的输入部分，作用与被模拟路段的起始断面上，它依依靠随机技术产生符合给定参数的泊松分布，向系统提供初值，并以经验概率分布于六个不同的车道。车辆的行驶模型即是反应车辆在路段上行驶状态变化的模型，本文根据交通调查得出的一般车流车辆类型的分布将车辆分为十一种类型，前六种为固定车重的车辆，以各自不等的概率出现在模拟路段的车道上；后六种车是变载的车辆，将其分为空车和满载两种情况（不计车辆半载情况），并以不同的概率出现。车速根据交通调查所得的各类型的车辆的平均车速，并允许其随机产生上下数值为μ的震荡，这里我们取μ为0.1。

油藏数值模拟显式差分MATLAB源程序

%显式求解方法 %t为投产后某一时刻，单位：天; %d:迭代时间； %Pwf1:W1井底流压; %Q2:W2井产油量; function [P1,d,Pwf1,Q2]=explict(t) %油藏参数 Pini=20; u=5e-3; C=2e-4; Q1=30; Pwf2=15; dx=200; dy=200; dt=1.7; n=t*24/dt;%迭代时间步数 re=0.208*dx; rw=0.1; %渗透率 K=[0 259 222 200 190 180 185 0 0 0 0;259 259 222 200 190 180 185 185 0 0 0;310 310 240 235 228 210 195 195 0 0 0;330 330 290 270 250 230 205 197.5 180 185 0;350 350 300 280 259 222 200 190 180 185 185;340 340 320 290 310 240 235 228 210 195 195;355 355 335 315 310 290 270 250 230 205 205;0 0 0 0 325 300 280 240 210 215 215;0 0 0 0 340 320 290 260 235 225 225;0 0 0 0 355 335 315 295 275 255 0]; %厚度 H=K/50; %孔隙度 Fai=(K.*0.02+15)/100; %原始地层压力 P=Pini*ones(10,11); %P(8:10,1:4)=0; %P(1:3,9:11)=0; %P(1,1)=0; %P(1,8)=0; %P(10,11)=0; %P(4,11)=0; %生产后某一时刻地层压力 P1=Pini*ones(10,11); %P1(8:10,1:4)=0; %P1(1:3,9:11)=0; %P1(1,1)=0; %P1(1,8)=0; %P1(10,11)=0; %P1(4,11)=0;