基于MATLAB的通信系统的设计与仿真
2011届学士学位论文
基于MATLAB的通信系统的设计与仿真学院、专业物理与电子信息学院
电子信息工程
研究方向系统仿真
学生姓名
学号
指导教师姓名
指导教师职称教授
2011年4月29日
基于MATLAB的通信系统的设计与仿真
摘要通信是通过某种媒体进行的信息传递，目的是传输信息，通信系统是用以完成信息传输过程的技术系统的总称，作用是将信息从信源发送到一个或多个目的地。调制与解调在信息的传输过程中占据着重要的地位，是不可或缺的，因此研究系统的调制和解调过程就极为重要。MATLAB是集数值计算、图形绘制、图像处理及系统仿真等强大功能于一体的科学计算语言，它强大的矩阵运算和图形可视化的功能以及丰富的工具箱，为通信系统的调制和解调过程的分析提供了极大的方便。
本论文首先介绍了通信系统的概念，进而引出调制和解调，然后介绍了我们常用的几种调制和解调的方法。由于MATLAB具有的强大功能所以详细介绍了MATLAB通信系统工具箱，并给出了基于MATLAB的通信系统的调制与解调的实现，运用MATLAB仿真软件进行仿真。
关键词通信系统；调制与解调；MA TLAB
Simulation And Design Of Communication Systems Based
On MATLAB
Abstract Communication is through a media for transportation. Communication system which is used to complete the process of information transmission systems ,in general, is to send the information from the source to one or more destinations. Modulation and demodulation occupied an important position in the transmission of information which is essential, so the research about the modulation and demodulation process in the communication system is extremely important. MATLAB is a numerical computation, graphics rendering, image processing and system simulation and other powerful features in one of the scientific computing language, it is a powerful matrix calculation and graphical visualization features and a rich toolbox provides a great convenience for the communication system of modulation and demodulation process.
This paper introduces the concept of the communication system, and then leads to modulation and demodulation, and then introduced several of our commonly used method of modulation and demodulation. As the power of MATLAB so we introduced the communication system toolbox in the MATLAB. We gives several examples about the communication system based on MATLAB modulation and demodulation and use the software of MATLAB to simulate them.
Keywords Communication Systems;Modulation and demodulation; MATLAB
目录
1 引言 (1)
2 通信系统的调制与解调 (2)
2.1通信系统的概念 (2)
2.2调制和解调的概念 (4)
2.3幅度调制和相干解调的原理 (5)
3 MATLAB通信系统工具箱 (7)
4 基于MATLAB的模拟调制和解调实例 (8)
4.1用MATLAB分析双边带幅度调制（DSM-AM） (8)
4.2用MATLAB分析相干解调过程 (10)
总结 (12)
参考文献 (13)
致谢 (14)
1 引言
通信按照传统的理解就是信息的传输，在当今高度信息化得社会，信息和通信已经成为现代社会的“命脉”。信息作为一种资源，只有通过广泛的传播与交流，才能产生利用价值，促进社会成员之间的合作，推动社会生产力的发展，创造出巨大的经济效益。而通信作为传输信息的手段或方式，与传感技术、计算机技术相互融合，已经成为21世纪国际社会和世界经济发展的强大推动力，所以未来的通信对人们的生活方式和社会的发展将会产生更加重大和意义深远的影响。
在信息传输过程中，要求天线的尺寸要和信号的波长相比拟，信号才能有效的被辐射。对于语音信号来说，相应的天线尺寸要在几十公里以上，实际上不可能实现，所以需要经过调制将信号频谱搬移到较高的频率范围，如果不进行调制就把信号直接辐射出去，那么各电台所发出信号的频率就会相同。调制作用的实质就是使相同频率范围的信号分别依托于不同频率的载波上，接收机就可以分离出所需的频率信号，不致互相干扰。
有时信号过于复杂，人工计算其调制和解调过程较难实现，对其结果的分析又缺乏可视化的直观表现，影响了所得结果在实际生活中的应用，美国MathWorks 公司开发的MATLAB解决了这一问题。它应用于自动控制、数学计算、信号分析、信号处理等诸多领域,也是国内高校和研究部门进行许多科学研究的重要工具。MATLAB的出现给通信系统的分析提供了极大的方便。
本文主要包括三部分：第一部分介绍了通信系统的调制和解调，第二部分介绍了MATLAB通信工具箱，第三部分给出了基于MATLAB的通信系统的分析。
2 通信系统的调制与解调
2.1通信系统的概念
通信是为了传输信息，通信系统就是将信息从信源发送到一个或多个目的地，对于电通信来说，首先要把消息转变成电信号，然后经过发送设备，将信号送入信道，在接收端利用接受设备对接收信号作相应的处理后，送给信宿再转换为原来的消息，这一过程可利用图1所示的通信系统一般模型来概括。
（发送端） （接收端）
图 1 通信系统的一般模型
通信(Communication)传输的消息是多种多样的，可以分成两大类：一类称为连续消息；另一类称为离散消息。消息的传递是通过它的物理载体电信号来实现的，按信号参量的取值不同，可以把信号分为两类：模拟信号和数字信号。通常按照信道中传输的是模拟信号还是数字信号，相应的可以把通信系统分为模拟通信系统和数字通信系统。
1． 模拟通信系统
模拟通信系统是利用模拟信号来传递信息的通信系统，其模型如图2所示，其中包含两种重要变换。第一种变换是，在发送端把连续消息变换成原始电信号，在接收端进行相反的变换，这种变换由信源和信宿来完成，通常称为原始电信号为基带信号，基带的意思是指信号的频谱从零频附近开始。有些信道可以直接传输基带信号，而以自由空间作为信道的无线电传输却无法直接传输这些信号。因此，模拟通信系统中常常需要进行第二种变换：把基带信号变换成适合在信道中传输的信号，并在接收端进行反变换。完成这种变换和反变换的通常是调制器和解调器[1]。
图2 模拟通信系统模型型
信息源 发送设备 信道 接收设备 噪声源 受信者 受信者 模拟信息源 解调器 调制器 噪声源 信
道
2. 数字通信系统
数字通信系统（Digital Communication System,DCS ）是利用数字信号来传递信息的通信系统，如图3所示。数字通信所涉及的技术问题很多，其中主要有信源编码与译码、信道编码与译码、数字调制与解调、同步以及加密与解密等。
图3 数字通信系统模型
与模拟通信相比，数字通信具有许多优良的特性：
（1） 抗干扰能力强，且噪声不积累。数字通信系统中传输的是离散取值的
数字波形，接收端的目标不是精确的还原被传输的波形，而是从收到
噪声干扰的信号中判决出发送端所发送的是哪一个波形。
（2） 传输差错可控。在数字同喜系统中，可以通过信道编码技术进行检错
与纠错，降低误码率，提高传输质量。
（3） 便于用现代数字信号处理技术对数字信息进行处理、变换、存储。这
种数字处理的灵活性表现为可以将来自不同信源的信号综合到一起
传输。
（4） 易于集成，使通信设备微型化，重量轻。
（5） 易于加密处理，且保密性好。
数字通信的缺点是，一般需要较大的传输带宽，设备复杂。近年来，随着大规模集成电路的出现，数字系统的设备复杂程度和技术难度大大降低了，数字传输方式日益受到欢迎。数字处理的灵活性使得数字传输系统中传输的数字信息既可以来自计算机、电传机等数据终端的各种数字代码，也可以来自模拟信号经过数字化处理后的脉冲编码（PCM ）信号等。原理上，数字信息可以直接用数字代码序列表示和传输，但在世纪传输中，视系统的要求和信道情况，一般需要进行不同形式的编码，并且选用一组取值有限的离散波形来表示，这些取值离散的波形可以是未经调制的电信号，也可以是调制后的信号。未经调制的数字信号所占据的频谱是从零频或很低频率开始，称为数字基带信号。
信息源 信源编码 加密 信道编码 数字调制 信道 噪声源 数字解调 信道译码 解密 信源译码 受信者
2.2调制和解调的概念
调制：把信号转换成适合在信道中传输的形式的一种过程，广义的调制分为基带调制和带通调制（也称为载波调制）。在无线通信中和其他大多数场合，调制均指载波调制【2】。
载波调制，就是用调制信号去控制载波的参数的过程，使载波的某一个或某几个参数按照调制信号的规律而变化。调制信号是指来自信源的消息信号（基带信号），这些信号可以是模拟的，也可以是数字的。未受调制的周期性振荡信号称为载波，它可以是正弦波，也可以是非正弦波（如周期性脉冲序列）。载波调制后称为已调信号，它含有调制信号的全部特征。
解调：将已调信号中的调制信号恢复出来，是调制的逆过程。调制方式不同，解调方法也不一样。解调可以分为正弦波解调（有时也称为连续波解调）。正弦波解调还可再分为幅度解调、频率解调和相位解调，此外还有一些变种如单边带信号解调、残留边带信号解调等。同样，脉冲波解调也可分为脉冲幅度解调、脉冲相位解调、脉冲宽度解调和脉冲编码解调等。
调制方式有很多，根据调制信号是模拟信号还是数字信号，载波是连续波（通常是正弦波）还是脉冲序列，相应的调制方式有模拟连续波调制（简称模拟调制）、数字连续波调制（简称数字调制）、模拟脉冲调制和数字脉冲调制等。最重要和最常用的模拟调制方式是用正弦波作为载波的幅度调制和角度调制。常见的调幅、双边带、单边带和残留边带等调制就是幅度调制的几个典型实例。
解调的方法可以分为两类：相干解调和非相干解调（包络检波）。解调过程与采用何种解调方式有关，对于常规幅度调制，一般用包络检波进行解调，由于在这种解调方式中，接收机对载波频率和相应精度的了解是无关紧要的，所以解调过程相对简单。对于DSB-AM调制和SSB-AM调制，用相干解调的方法，它要求在接收机中有一个与载波同频同相的信号，接收机中产生所需要的正弦波振荡器，为本地振荡器。
数字通信系统中，发送信号某一参数的离散取值（如正弦波的振幅、频率、相位等）于所承载消息的数字信号（即离散符号）之间通常是一种一一对应关系。如同模拟信号的频带传输时一样，传输数字信号时也有三种基本的调制方式：振幅键控（ASK）、频移键控（FSK）和相移键控（PSK）。它们分别对应于利用载波（正弦波）的振幅、频率和相位来承载数字基带信号，可以是模拟线性调制和角度调制的特殊情况。
理论上数字调制与模拟调制在本质上没有什么不同，它们都属于正弦波调制。但是，数字调制是源信号为离散型的正弦波调制，而模拟调制则是源信号为连续型的正弦波调制，因而，数字调制具有由数字信号带来的一些特点。这些特点主要包括两个方面：第一，数字调制信号的产生，除把数字的调制信号当做模拟信号的特例而直接采用模拟调制方式产生数字调制信号外，还可以采用键控载波的方法。第二，对于数字调制信号的解调，为提高系统的抗噪声性能，通常采用与模拟调制系统中不同的解调方式。
数字信息由二进制和多进制之分，因此，数字调制可分为二进制调制和多进制调制。在二进制调制中，信号参量只有两种可能取值；而在多进制调制中，信号参量可能有M（M>2）中取值。
2.3 幅度调制和相干解调的原理
1）幅度调制（线性调制）的原理
幅度调制是由调制信号去控制高频载波的幅度，使之随调制信号作线性变化的过程，设正弦型载波为
c(t)=Acos(?ω+ct 0) (2.3-1)
式中：A 为载波幅度；c ω为载波角频率；?0为载波初始相位（以后可以假定为?0为0，而不失讨论的一般性）。
根据调制定义，幅度调制信号（已调信号）一般可表示成
Sm(t)=Am(t)cos (c ωt) (2.3-2)
式中：m(t)为基带调制信号。
设调制信号m(t)的频谱为)(ωM ,则由上式（2.3-2）不难得到已调信号Sm(t)的频谱)(ωS ；
)(ωS = )]()(c M c M ωωωω-++ (2.3-3)
由以上表示式可见，在波形上，幅度已调信号的幅度随基带信号的规律而呈正比地变化；在频谱结构上，它的频谱完全是基信号频谱在频域内的简单搬移（精确到常数因子）。由于这种搬移是线性的，因此，幅度调制通常又称为线性调制，这里的“线性”并不意味着已调信号与调制信号之间符合线性变换关系，实际上，
任何调制过程都是一种非线性的变换过程【3】。
标准调幅就是常规双边带调制，简称调幅（AM ）。假设调制信号m(t)的平均值为0，将其叠加一个直流偏量A0后与载波相乘（如图4所示），即可形成调幅信号。
m(t) Sm(t)
A0 ct ωcos
图 3 AM 调制模型
其时域表示式为
S AM (t)=[A0+m(t)] ct ωcos =A 0 ct ωcos +m(t) ct ωcos (2.3-4)
式中：A0为外加的直流分量；m(t)可以是确知信号，也可以是随机信号。 若m(t)是确知信号，则AM 信号的频谱为
S AM (ω)=∏A0[)()(c c ωωδωωδ-++]+1/2[)()(c M c M ωωωω-++] (2.3-5)
2) 相干解调
相干解调也叫同步检波，解调与调制的实质一样，均是频谱搬移。调制是把基带信号的频谱搬到了载波位置，这一过程可以通过一个相乘器与载波相乘来实现。解调则是调制的反过程，即把在载频位置的已调信号的频谱搬回到原始基带位置，因此同样可以用相乘器与载波相乘来实现，相干解调的一般模型如图4所示。
Sm(t) Sp(t) Sd(t)
c(t)= ct ωcos
图 4 相干解调器的一般模型
相干解调时，为了无失真的恢复原基带信号，接收端必须提供一个与接收的已调载波严格同步（同频同相）的本地载波（称为相干载波），它与接收的已调信号相乘后，经低通滤波器取出低频分量，即可得到原始的基带调制信号。 相干解调器适用于所有的线性调制信号的解调。送入解调器的已调信号的一般表达式为
Sm(t)=S I (t) ct ωcos +S Q (t) ct ωsin (2.3-6)
与同频同相的相干载波c(t)相乘后，得
S p (t)= Sm(t) ct ωcos =
1/2 S I (t)+1/2 S I (t) t ω2cos +1/2 S Q (t) ct ω2sin (2.3-7)
经低通滤波器（LPF ）后，得
Sd(t)=1/2 S I (t) （2.3-8）
其中，S I (t)是m(t)通过一全通滤波器H1（ω）后的结果。因此，Sd(t)就是解调输出，即
Sd(t)=1/2 S I (t) ∝ m(t)
由此可见，相干解调器适用于所有线性调制信号的解调，即对于AM 、DSB 、SSB 和VSB 都是适用的。只是AM 信号的解调结果中含有直流成分A0，这时在解调后加上一个简单的隔直流电容即可。
LPF
3 MATLAB通信系统工具箱
MATLAB通信系统工具箱主要应用于参数化模型，频谱分析和估计等。通信系统工具箱实际上是用MATLAB的基本语句变成的各种子程序集，用于解决某一方面的专门问题或实现某一类的新算法。MATLAB的通信系统工具箱可以任意增减，不同的工具箱可以给不同领域的用户提供丰富强大的功能。
MATLAB中的通信系统工具箱是目前比较成熟的通信系统仿真工具，MATLAB通信系统工具箱提供的函数主要用于处理通信及传输问题，为我们分析通信系统中每个过程提供了极大的帮助，使我们能够轻松直观的解决问题【4】。
常用的MATLAB函数主要有以下几项：
1、figure
功能：创建新的图形窗口(用于输出图形的窗口)。
格式：figure
说明：figure函数创建一个新的图形窗口, 并成为当前图形窗口, 所创建的图形窗口的序号是按同一MATLAB程序中创建的顺序号。
2、plot
功能：线型绘图函数。
格式：plot(x)
plot(x,y)
说明：plot(x)是一种最简单的调用方式, x是长度为n的数值向量。plot(x)的作用是在坐标系中顺序地用直接连接顶点{i,x(i),i=1,2,…,n}生成一条折线。当向量元素充分多时, 即可生成一条光滑的曲线。
3、subplot
功能：多坐标设置与定位当前坐标系。
格式：subplot(m,n,k)
说明：subp lot(m,n,k)将图形窗口分成m行n列m*n块子区域，按从上到下,从左到右的顺序,在第k块子区域定义一个坐标系, 使其成为当前坐标系, 随后的绘图函数将在该坐标系输出图形。
4、ezplot
功能：该命令用来绘制符号表达式的自变量和对应各函数值的二维曲线。
格式： ezplot(F, [xmin,xmax],fig)
说明：其中F是要画的符号函数；[xmin,xmax]是绘图的自变量范围，fip是窗口。
5、axis
功能：该命令用来控制坐标轴的特性。
格式： axis([xmin,xmax], [ymin,ymax])
说明：此为坐标范围,其中xmin＜xmax ，ymin＜ymax 。
7、angle
功能：求复数的相角。
格式：P=angle(Z)
说明：当Z为复数矩阵时，用来求矩阵Z中每个元素的相角，相角位于[-π,π]。
4 基于MATLAB 的模拟调制和解调实例
4.1 用MATLAB 分析双边带幅度调制（DSM-AM ）
我们可以利用 MATLAB 强大的符号运算功能来进行运算,再根据MATLAB 的可视化结果进行分析。
在DSB-AM 中，已调信号的时域表示为：
u(t)=m(t)c(t)=Ac m (t)
∏+)2cos(c fct φ 式中，m(t)是消息信号，c(t)=Ac ∏+)2cos(c fct φ为载波，fc 是载波的频率（单位：HZ ）, c φ是初始相位。为了讨论方便取初相c φ=0（以下类似）。
随u(t)作傅里叶变换，即可得到信号的频域表示：
U （f ）=Ac/2 )(fc f M -+Ac/2 )(fc f M +
传输带宽Bt 是消息信号带宽W 的两倍，即：Bt=2W 。
1 0＜t ＜t0/3
例1 某消息信号 m(t)= -2 t0/3＜t ＜2t0/3
0 其他 用信号m(t)以DSB-AM 方式调制载波c(t)=cos(2∏f c t ),所得到的已调制信号记为u(t).设t 。=0.15s,f c =250Hz 。试比较消息信号与已调信号，并绘制它们的频谱。 运用如下 MATLAB 程序【5】:
t 0=0.15; ％信号持续时间
t s =0.001; ％采样时间间隔
F c =250; ％载波频率
F S =1/t s; ％采样频率
df = 0.3; ％频率分辨率
t =[0:t s:t 0]; ％时间矢量
m=[ones(1,t 0/(3*t s)), -2*ones(1,t 0/(3*t s)),zeros(1,t 0/(3*t s)+1)]; ％定义信号序列
C =cos(2*pi*F c.*t); ％载波信号
u=m.*c; ％调制信号
[M,m,df1]=fft_seq(m,ts,df); ％傅里叶变换
M=M/F S;
[U,u,df1]=fft_seq(m,ts,df);
U=U/FS;
[C,c,df1]=fft_seq(c,ts,df);
f=[0:df1:df1*(length(m)-1)]-Fs/2; ％频率矢量
subplot(2,2,1);plot(t,m(1:length(t)));％未调制信号
title(‘未调制信号‘)；
subplot(2,2,2);plot(t,u(1:length(t))); ％已调制信号
title(‘已调制信号‘)；
subplot(2,2,3)；plot(f,abs(fftshift(M))); ％未调制信号频谱title(‘未调制信号频谱‘)；
subplot（2,2,4）；plot(f,abs(fftshift(U))); ％已调制信号频谱title(‘已调制信号频谱‘)；
傅里叶变换函数fft_seq,源代码如下：
function[M,m,df]=fft_seq(m,ts,df)
％[M,m,df]=fft_seq(m,ts,df)
％[M,m,df]=fft_seq(m,ts)
％M为输入序列m的傅里叶变换，ts为抽样间隔，输入df为频率分辨率％输出序列m按要求的频率分辨率df进行补零后的序列
％输出df为最终的频率分辨率
fs=1/ts;
if nargin==2,n1=0;
else, n1=fs/df;
end
n2=length(m);
n=2^(max(nextpow2(n1),nexpow2(n2)));
M=fft(m,n);
m=[m,zeros(n-n2)];
df=fs/n;
运行后得到的信号和调制信号如图所示：
4.2 用MATLAB 分析相干解调过程
1 0＜t＜t0/3
例某消息信号m(t)= -2 t0/3＜t＜2t0/3
0其他
用信号m(t)以DSB-AM方式调制载波c(t)=cos(2∏f c t),所得到的已调制信号记为
u(t).设t0=0.15s,f c=250Hz。得到DSM-AM的调制信号，现试对该调制信号进
行相干解调，并绘出消息信号的时域频域曲线。
MATLAB设计程序如下【6】：
t0=0.15; ％信号持续时间
t s=1/1500; ％采样时间间隔
F c=250; ％载波频率
F S=1/t s; ％采样频率
df= 0.3; ％频率分辨率
t=[0:t s:t0]; ％时间矢量
m=[ones(1,t0/(3*t s)), -2*ones(1,t0/(3*t s)),zeros(1,t0/(3*t s)+1)];
％定义信号序列
C=cos(2*pi*Fc.*t); ％载波信号
u=m.*c; ％调制信号
y=u.*c; ％混频
[M,m,df1]=fft_seq(m,ts,df);％傅里叶变换
M=M/Fs; ％缩放
[U,u,df1]=fft_seq(m,ts,df);
U=U/FS;
[Y,y,df1]=fft_seq(y,ts,df);
Y=Y/Fs;
f_cutoff=150; ％滤波器的截止频率
n_cutoff=floor(150/df1); ％设计滤波器
f=[0:df1:df1*(length(m)-1)]-Fs/2; ％频率矢量
H=zero(size(f));
H(1:n_cutoff)=2*ones(1,n_cutoff);
H(length(f)-n_cutoff+1:length(f))=2*ones(1,n_cutoff);
DEM=H.*Y; ％滤波器输出的频谱
dem=real(ifft(DEM))*Fs; ％滤波器的输出
subplot(2,2,1);plot(t,m(1:length(t)));％未调制信号
title(‘未调制信号‘)；
subplot(2,2,2);plot(t,dem(1:length(t)));％解调信号
title(‘解调信号‘)；
subplot(2,2,3)；plot(f,abs(fftshift(M))); ％未调制信号频谱title(‘未调制信号频谱‘)；
subplot（2,2,4）；plot(f,abs(fftshift(DEM)))；％解调信号的频谱title(‘已调制信号谱‘)【7】；
运行该程序后得到的信号和调制信号，以及信号调制前、后的频谱对比如图所示：
总结
本论文主要介绍了通信系统的含义，进而引出模拟调制和解调的概念。人工计算调制和解调过程比较困难，而MATLAB具有的强大功能可以很容易的解决这一问难题所以详细介绍了MATLAB通信系统工具箱；并给出了基于MATLAB的通信系统的设计和仿真的实现，运用MATLAB仿真软件进行仿真。通过对本次应用MATLAB语言对通信系统的分析，使我系统的理解了调制和解调的功能，对掌握信息的传输有了进一步提高。从上面的分析我们可以看出用MATLAB辅助分析信号频谱，具有编程简单、计算准确、绘图方便、结果直观等特点。用MATLAB的计算功能，不仅使大量的手工计算得以简化, 也使得离散系统的分析更为简便和高效。利用MATLAB的绘图功能，有利于分析结果的直观理解，也有利于深入掌握所学的内容。
参考文献
[1]樊昌信，曹丽娜. 通信原理（第六版）[M],北京：国防工业出版社，2009
[2]程相君，陈生潭信号与系统[M],西安：西安电子科技大学出版社,1990
[3]曹志刚，钱压生现代通信原理[M],北京：清华大学出版社，1992
[4]刘敏，魏玲MATLAB通信仿真与应用[M]，北京：国防工业出版社，2001
[5]唐向宏，岳恒立，郑学峰MATLAB及在电子信息类课程中的应用[M]，北京：
电子工业出版社，2006
[6]陈杰MATLAB宝典[M]，北京：电子工业出版社
[7]张德丰MATLAB/Simulink建模与仿真[M]，北京电子工业出版社
致谢
在我做论文期间，我的指导老师李素文老师给予我极大的帮助，正是在李素文老师的悉心指导和无微不至的关怀下完成的，在论文即将结束之际，我衷心的感谢我的指导老师李素文老师，谢谢她的帮助和支持！另外在论文期间，我周围的同学特别是我宿舍里的室友也给了我很大的帮助，再次对我的室友说一声谢谢，也对我周围的朋友、同学以及我的家人表示感激，谢谢你们的支持!
Matlab通信系统仿真实验报告
Matlab通信原理仿真 学号： 2142402 姓名：圣斌
实验一Matlab 基本语法与信号系统分析 一、实验目的： 1、掌握MATLAB的基本绘图方法； 2、实现绘制复指数信号的时域波形。 二、实验设备与软件环境： 1、实验设备：计算机 2、软件环境：MATLAB R2009a 三、实验内容： 1、MATLAB为用户提供了结果可视化功能，只要在命令行窗口输入相应的命令，结果就会用图形直接表示出来。 MATLAB程序如下： x = -pi::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’运行结果如下图： 2、上例中的图形使用的是默认的颜色和线型，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);
基于MATLAB的MIMO通信系统仿真(DOC)
目录 （一）基于MATLAB的MIMO通信系统仿真………………………… 一、基本原理……………………………………………………… 二、仿真…………………………………………………………… 三、仿真结果……………………………………………………… 四、仿真结果分析…………………………………………………（二）自选习题部分…………………………………………………（三）总结与体会……………………………………………………（四）参考文献…………………………………………………… 实训报告 （一）基于MATLAB的MIMO通信系统仿真 一、基本原理 二、仿真 三、仿真结果 四、仿真结果分析 OFDM技术通过将频率选择性多径衰落信道在频域内转换为平坦信道，减小了多径衰落的影响。OFDM技术如果要提高传输速率，则要增加带宽、发送功率、子载波数目，这对于频谱资源紧张的无线通信时不现实的。 MIMO能够在空间中产生独立并行信道同时传输多路数据流，即传输速率很高。这些增加的信道容量可以用来提高信息传输速率，也可以通过增加信息冗余来提高通信系统的传输可靠性。但是MIMO却不能够克服频率选择性深衰落。 所以OFDM和MIMO这一对互补的技术自然走到了一起，现在是3G，未来也是4G，以及新一代WLAN技术的核心。总之，是核心物理层技术之一。 1、MIMO系统理论：
核心思想：时间上空时信号处理同空间上分集结合。 时间上空时通过在发送端采用空时码实现： 空时分组、空时格码，分层空时码。 空间上分集通过增加空间上天线分布实现。此举可以把原来对用户来说是有害的无线电波多径传播转变为对用户有利。 2、MIMO 系统模型： 11h 12 h 21 h 22 h r n h 1r n h 21 R n h 2 R n h 1 n n R h 可以看到，MIMO 模型中有一个空时编码器，有多根天线，其系统模型和上述MIMO 系统理论一致。为什么说nt>nr ，因为一般来说，移动终端所支持的天线数目总是比基站端要少。 接收矢量为：y Hx n =+，即接收信号为信道衰落系数X 发射信号+接收端噪声 3、MIMO 系统容量分析： （附MIMO 系统容量分析程序） 香农公式的信道容量（即信息传送速率）为： 2log (1/)C B S N =+ 4、在MIMO 中计算信道容量分两种情况： 未知CSI 和已知CSI （CSI 即为信道状态信息），其公式推导较为复杂，推导结果为信道容量是信噪比与接收、发射天线的函数。 在推导已知CSI 中，常用的有waterfilling ，即著名的注水原理。但是，根据相关文献资料，通常情况下CSI 可以当做已知，因为发送，接收端会根据具体信道情况估算CSI 的相关参数。 在这里对注水原理做一个简单介绍：之所以成为注水原理是因为理想的注水原理是在噪声大的时候少分配功率，噪声小时多分配功率，最后噪声+功率=定值，这如果用图形来表示，则类似于给水池注水的时候，水池低的地方就多注水，也就是噪声小分配的功率就多，故称这种达到容量的功率分配方式叫做注水原理。通过给各个天线分配不同的发射功率，增加系统容量。核心思想就是上面所阐述的，信道条件好，则分配更多功率；信道条件差，则分配较少的功率。 在MIMO 的信道容量当中要注意几个问题：（下面说已知CSI 都是加入了估计CSI 的算法，并且采用了注水原理。） 1. 已知CSI 的情况下的信道容量要比发送端未知CSI 的情况下的信道容量高，这是 由于当发送端已知CSI 的时候，发送端可以优化发送信号的协方差矩阵。也就是
机械优化设计MATLAB程序文件
机械优化设计作业1.用二次插值法求函数()()()22 ?极小值，精度e=0.01。 t t =t 1- + 在MATLAB的M文件编辑器中编写的M文件，如下： f=inline('(t+1)*(t-2)^2','t') a=0;b=3;epsilon=0.01; t1=a;f1=f(t1); t3=b;f3=f(t3); t2=0.5*(t1+t3);f2=f(t2); c1=(f3-f1)/(t3-t1); c2=((f2-f1)/(t2-t1)-c1)/(t2-t3); t4=0.5*(t1+t3-c1/c2);f4=f(t4); k=0; while(abs(t4-t2)>=epsilon) if t2
运行结果如下： 迭代计算k= 7 极小点坐标t= 2 函数值f=0.0001 2.用黄金分割法求函数()32321+-=t t t ?的极小值，精度e=0.01。 在MATLAB 的M 文件编辑器中编写的M 文件，如下： f=inline('t^(2/3)-(t^2+1)^(1/3)','t'); a=0;b=3;epsilon=0.01; t1=b-0.618*(b-a);f1=f(t1); t2=a+0.618*(b-a);f2=f(t2); k=1; while abs(b-a)>=epsilon if f1 中南大学MATLAB程序设计实践学长有爱奉献，下载填上信息即可上交，没有下载券的自行百度。所需m文件照本文档做即可，即新建（FILE）→脚本（NEW-Sscript）→复制本文档代码→运行（会跳出保存界面，文件名默认不要修改，保存）→结果。第一题需要把数据文本文档和m文件放在一起。全部测试无误，放心使用。本文档针对做牛顿法求非线性函数题目的同学，当然第一题都一样，所有人都可以用。←记得删掉这段话 班级： ？ 学号： 姓名： 一、《MATLAB程序设计实践》Matlab基础 表示多晶体材料织构的三维取向分布函数（f＝f（φ1，φ，φ2））是一个非常复杂的函数，难以精确的用解析函数表达，通常采用离散 空间函数值来表示取向分布函数，是三维取向分布函数的一个实例。 由于数据量非常大，不便于分析，需要借助图形来分析。请你编写一 个matlab程序画出如下的几种图形来分析其取向分布特征： （1）用Slice函数给出其整体分布特征； " ~ （2）用pcolor或contour函数分别给出(φ2＝0, 5, 10, 15, 20, 25, 30, 35 … 90)切面上f分布情况(需要用到subplot函数)； (3) 用plot函数给出沿α取向线(φ1=0~90，φ＝45，φ2＝0)的f分布情况。 ( 备注：数据格式说明 解： （1）（ （2）将文件内的数据按照要求读取到矩阵f(phi1,phi,phi2)中，代码如 下： fid=fopen(''); for i=1:18 tline=fgetl(fid); end phi1=1;phi=1;phi2=1;line=0; f=zeros(19,19,19); [ while ~feof(fid) tline=fgetl(fid); data=str2num(tline); line=line+1;数据说明部分，与 作图无关此方向表示f随着 φ1从0,5,10,15, 20 …到90的变化而 变化 此方向表示f随着φ 从0,5,10,15, 20 … 到90的变化而变化 表示以下数据为φ2＝0的数据，即f（φ1，φ，0） MATLAB通信系统仿真心得体会 课程名称(中文) MATLAB通信系统仿真成绩姓名班级学号日期 学习MATLAB通信系统仿真心得体会 经过一学期的MATLAB通信系统仿真的学习，使我对通信原 理及仿真实践有了更深层次的理解。在学习过程当中，了解到了MATLAB的语言基础以及应用的界面环境，基本操作和语法，通信仿真工具箱的应用，simulink 仿真基础，信号系统分析等一系列的内容。我明白学好这门课程是非常的重要。 在学习当中，我首先明白了通信系统仿真的现实意义，系统模型是对实际系统的一种抽象，是对系统本质(或是系统的某种特性)的一种描述。模型可视为对真实世界中物体或过程的信息进行形式化的结果。模型具有与系统相似的特性，可以以各种形式给出我们所感兴趣的信息。知道了通信系统仿真的必要性，利用系统建模和软件仿真技术，我们几乎可以对所有的设计细节进行分层次的建模和评估。通过仿真技术和方法，我们可以有效地将数学分析模型和经验模型结合起来。利用系统仿真方法，可以迅速构建一个通信系统模型，提供一个便捷，高效和精确的评估平台。明白了MATLAB通信系统仿真课程重点就是系统仿真软件 Matlab / Simulink 在通信系统建模仿真和性能评估方面的应用原理，通信系统仿真的一般原理和方法。 MATLAB集成度高，使用方便，输入简捷，运算高效，内容丰富，并且很容易由用户自行扩展，与其它计算机语言相比， MATLAB有以下显著特点:1.MATLAB是一种解释性语言;2(变量的“多功能性”;3.运算符号的“多功能性”;4(人机界面适合科技人员;5(强大而简易的作图功能;6(智能化程度高;7(功能丰富，可扩展性强。在MATLAB的Communication Toolbox(通 信工具箱)中提供了许多仿真函数和模块，用于对通信系统进行仿真和分析。 基于MATLAB 的曲柄摇杆机构优化设计 1. 问题的提出 根据机械的用途和性能要求的不同，对连杆机构设计的要求是多种多样的，但这些设计要求可归纳为以下三种问题：（1）满足预定的运动规律要求；（2）满足预定的连杆位置要求；（3）满足预定的轨迹要求。在在第一个问题里按照期望函数设计的思想，要求曲柄摇杆机构的曲柄与摇杆转角之间按照()f φ?=（称为期望函数）的关系实现运动，由于机构的待定参数较少，故一般不能准确实现该期望函数，设实际的函数为()F φ?=（称为再现函数），而再现函数一般是与期望函数不一致的，因此在设计时应使机构再现函数()F φ?=尽可能逼近所要求的期望函数()f φ?=。这时需按机械优化设计方法来设计曲柄连杆，建立优化数学模型，研究并提出其优化求解算法，并应用于优化模型的求解，求解得到更优的设计参数。 2. 曲柄摇杆机构的设计 在图 1 所示的曲柄摇杆机构中，1l 、2l 、3l 、 4l 分别是曲柄AB 、连杆BC 、摇杆CD 和机架AD 的长度。这里规定0?为摇杆在右极限位置0φ时的曲柄起始位置角，它们由1l 、2l 、3l 和4l 确定。 图1 曲柄摇杆机构简图 设计时，可在给定最大和最小传动角的前提下，当曲柄从0?转到090??+时，要求摇杆的输出角最优地实现一个给定的运动规律()f ?。这里假设要求： ()()2 0023E f φ?φ??π ==+ - （1） 对于这样的设计问题，可以取机构的期望输出角()E f φ?=和实际输出角 ()F φ?=的平方误差之和作为目标函数，使得它的值达到最小。 2.1 设计变量的确定 决定机构尺寸的各杆长度1l 、2l 、3l 和4l ，以及当摇杆按已知运动规律开始运行时，曲柄所处的位置角0?应列为设计变量，即: []12340T x l l l l ?= (2) 考虑到机构的杆长按比例变化时，不会改变其运动规律，通常设定曲柄长度 1l =1.0，在这里可给定4l =5.0，其他杆长则按比例取为1l 的倍数。若取曲柄的初始 位置角为极位角，则?及相应的摇杆l 位置角φ均为杆长的函数，其关系式为: ()()()()222221243230124225arccos 210l l l l l l l l l l l l ?????++-+-+==????++???????? （3） ()()22222124323034325arccos 210l l l l l l l l l l ????? +--+--==???????????? （4） 因此，只有2l 、3l 为独立变量，则设计变量为[][]2312T T x l l x x ==。 2.2目标函数的建立 目标函数可根据已知的运动规律与机构实际运动规律之间的偏差最小为指标来建立，即： ()()2 1min m Ei i i f x φφ==-→∑ （5） 式中，Ei φ-期望输出角；m -输出角的等分数；i φ-实际输出角，由图 1 可知： ()()02i i i i i i i παβ?πφπαβπ?π--≤≤??=?-+≤≤?? (6) 式中，222222322132arccos arccos 22i i i i i r l l r x x rl r x α???? +-+-== ? ????? (7) 222241424arccos arccos 210i i i i i r l l r rl r β???? +-+== ? ????? (8) i r == (9) 2.3约束条件 实验一 Matlab 程序设计与M 文件 一、实验目的 1. 掌握Matlab 程序设计常用命令，如，循环、选择、暂停、显示输出、输入变量值等。 2. 掌握Matlab 的M 文件的创建于使用，包括脚本文件和函数文件。 3. 掌握常用的编程技巧。 二、实验内容 1. 输入一个百分制成绩，要求输出成绩等级A 、B 、C 、D 、E 。其中90分~100分为A ，80分~89分为B ，79分~79分为C ，60分~69分为D ，60分以下为E 。 要求： (1) 分别用if 语句和switch 语句实现。 (2) 输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出出错信息。 2. 硅谷公司员工的工资计算方法如下： (1) 工作时数超过120小时者，超过部分加发15%。 (2) 工作时数低于60小时者，扣发700元。 (3) 其余按每小时84元计发。 试编程按输入的工号和该号员工的工时数，计算应发工资。 3. 根据2 222211116123n π=++++，求π的近似值。当n 分别取100、1000、10000时，结果是多少？ 要求：分别用循环结构和向量运算（使用sum 函数）来实现。 4. 考虑以下迭代公式： 1n n a x b x +=+ 其中a 、b 为正的学数。 (1) 编写程序求迭代的结果，迭代的终止条件为|x n+1-x n |≤10-5，迭代初值x 0=1.0，迭代次数不超过500次。 (2) 如果迭代过程收敛于r ，那么r 的准确值是，当(a,b)的值 取(1,1)、(8,3)、(10,0.1)时，分别对迭代结果和准确值进行比较。 5. 已知 12312311021 323 n n n n f n f n f n f f f f n ---==??==??==??=-+>? 求f 1~f 100中： (1) 最大值、最小值、各数之和。 (2) 正数、零、负数的个数。 6. 若两个连续自然数的乘积减1是素数，则称这两个边疆自然数是亲密数对，该素数是亲密素数。例如，2×3-1=5，由于5是素数，所以2和3是亲密数，5是亲密素数。求[2,50]区间内： (1) 亲密数对的对数。 (2) 与上述亲密数对对应的所有亲密素数之和。 7. 设2411()(2)0.1(3)0.01 f x x x =+-+-+，编写一个MATLAB 函数文件fx.m ，使得调用f(x)时，x 可用矩阵代入，得出的f(x)为同阶矩阵。 8. 一物理系统可用下列方程组来表示： 111211 12 220cos sin 0sin 0cos 000sin 000cos 1a m m a m g m N m N m g θθθθθθ--??????????????????=??????-??????-?????? 从键盘输入m 1、m 2和θ的值，求a 1、a 2、N 1和N 2的值。其中g 取9.8，输入θ时以角度为单位。 要求：定义一个求解线性方程组AX=B 的函数文件，然后在命令文件中调用该函数文件。 实验一 2PSK调制数字通信系统 一实验题目 设计一个采用2PSK调制的数字通信系统 设计系统整体框图及数学模型； 产生离散二进制信源，进行信道编码（汉明码），产生BPSK信号； 加入信道噪声（高斯白噪声）； BPSK信号相干解调，信道解码； 系统性能分析（信号波形、频谱，白噪声的波形、频谱，信道编解 二实验基本原理 数字信号的传输方式分为基带传输和带通传输，在实际应用中，大多数信道具有带通特性而不能直接传输基带信号。为了使数字信号在带通信道中传输，必须使用数字基带信号对载波进行调制，以使信号与信道的特性相匹配。这种用数字基带信号控制载波，把数字基带信号变换为数字带通信号的过程称为数字调制。 数字调制技术的两种方法：①利用模拟调制的方法去实现数字式调制，即把数字调制看成是模拟调制的一个特例，把数字基带信号当做模拟信号的特殊情况处理；②利用数字信号的离散取值特点通过开关键控载波，从而实现数字调制。这种方法通常称为键控法，比如对载波的相位进行键控，便可获得相移键控（PSK）基本的调制方式。 图1 相应的信号波形的示例 1 0 1 调制原理 数字调相：如果两个频率相同的载波同时开始振荡，这两个频率同时达到正最大值，同时达到零值，同时达到负最大值，它们应处于"同相"状态；如果其中一个开始得迟了一点，就可能不相同了。如果一个达到正最大值时，另一个达到负最大值，则称为"反相"。一般把信号振荡一次（一周）作为360度。如果一个波比另一个波相差半个周期，我们说两个波的 相位差180度，也就是反相。当传输数字信号时，"1"码控制发0度相位，"0"码控制发180度相位。载波的初始相位就有了移动，也就带上了信息。 相移键控是利用载波的相位变化来传递数字信息，而振幅和频率保持不变。在2PSK中，通常用初始相位0和π分别表示二进制“1”和“0”。因此，2PSK信号的时域表达式为(t)=Acos t+) 其中，表示第n个符号的绝对相位： = 因此，上式可以改写为 图2 2PSK信号波形 解调原理 2PSK信号的解调方法是相干解调法。由于PSK信号本身就是利用相位传递信息的，所以在接收端必须利用信号的相位信息来解调信号。下图2-3中给出了一种2PSK信号相干接收设备的原理框图。图中经过带通滤波的信号在相乘器中与本地载波相乘，然后用低通滤波器滤除高频分量，在进行抽样判决。判决器是按极性来判决的。即正抽样值判为1，负抽样值判为0. 2PSK信号相干解调各点时间波形如图 3 所示. 当恢复的相干载波产生180°倒相时,解调出的数字基带信号将与发送的数字基带信号正好是相反,解调器输出数字基带信号全部出错. Matlab程序设计作业 姓名: 学号: 专业: ? MATLAB 程序设计》作业 1、考虑如下x-y 一组实验数据: x=[1, 2, 3, 4, 5, 6, 7, 8, 9, 10] y 二[1.2, 3, 4, 4, 5, 4.7, 5, 5.2, 6, 7.2] 分别绘出plot 的原始数据、一次拟合曲线和三次拟合曲线，给出 原始曲线 MATLAB 代码和运行结果。 7 6 5 4 3 2 2 3 4 5 6 7 8 9 10 7 6.5 6 5.5 5 4.5 4 3.5 3 2.5 10 一次拟合 三次拟合 x=[1, 2, 3, 4, 5, 6, 7, 8, 9, 10]; y=[1.2, 3, 4, 4, 5, 4.7, 5, 5.2, 6, 7.2]; figure; plot(x,y) p1=polyfit(x,y,1); y1=polyval(p1,x); figure; plot(x,y1) p2=polyfit(x,y,3); y2=polyval(p2,x); figure; plot(x,y2) 2、在［0, 3n区间，绘制y二Sin(x)曲线(要求消去负半波，即(n 2n)区间内的函数值置零)，求出曲线y 的平均值，以及y 的最大值及其最大值的位置。给出执行代码和运行结果。 x=0:pi/1000:3*pi; y=Sin(x); y1=(y>=0).*y; %消去负半波figure(1); plot(x,y1, 'b' ); a=mean(y1) %求出y1 的平均值 b=max(y1) %求出y1 的最大值b, 以及最大值在矩阵中的位置; d=x(find(y1==b)) >> ex1 a = 0.4243 b = 1 d = 1.5708 7.8540 >> 补充内容：模拟调制系统的MATLAB 仿真 1.抽样定理 为了用实验的手段对连续信号分析，需要先对信号进行抽样（时间上的离散化），把连续数据转变为离散数据分析。抽样（时间离散化）是模拟信号数字化的第一步。 Nyquist 抽样定律：要无失真地恢复出抽样前的信号，要求抽样频率要大于等于两倍基带信号带宽。 抽样定理建立了模拟信号和离散信号之间的关系，在Matlab 中对模拟信号的实验仿真都是通过先抽样，转变成离散信号，然后用该离散信号近似替代原来的模拟信号进行分析的。 【例1】用图形表示DSB 调制波形)4cos()2cos(t t y ππ= 及其包络线。 clf %%计算抽样时间间隔 fh=1;%%调制信号带宽(Hz) fs=100*fh;%%一般选取的抽样频率要远大于基带信号频率，即抽样时间间隔要尽可能短。 ts=1/fs; %%根据抽样时间间隔进行抽样，并计算出信号和包络 t=(0:ts:pi/2)';%抽样时间间隔要足够小，要满足抽样定理。 envelop=cos(2*pi*t);%%DSB 信号包络 y=cos(2*pi*t).*cos(4*pi*t);%已调信号 %画出已调信号包络线 plot(t,envelop,'r:','LineWidth',3); hold on plot(t,-envelop,'r:','LineWidth',3); %画出已调信号波形 plot(t,y,'b','LineWidth',3); axis([0,pi/2,-1,1])% hold off% xlabel('t'); %写出图例 【例2】用图形表示DSB 调制波形)6cos()2cos(t t y ππ= 及其包络线。 clf %%计算抽样时间间隔 fh=1;%%调制信号带宽(Hz) fs=100*fh;%抽样时间间隔要足够小，要满足抽样定理。 ts=1/fs; %%根据抽样时间间隔进行抽样 第三章1. (1)A=eye(3) (2)A=100+100*rand(5,6) (3)A=1+sqrt(0.2)*randn(10,50) (4)B=ones(size(A)) (5)A+30*eye(size(A)) (6)B=diag(diag(A)) 2. B=rot90(A) C=rot90(A,-1) 3. B=inv(A) ;A的逆矩阵 C=det(A) ;A的行列式的值 D=A*B E=B*A D=E 因此A与A-1是互逆的。 4. A=[4 2 -1;3 -1 2;12 3 0]; b=[2;10;8]; x=inv(A)*b x = -6.0000 26.6667 27.3333 5. (1) diag(A) ；主对角线元素 ans = 1 1 5 9 triu(A) ；上三角阵 ans = 1 -1 2 3 0 1 -4 2 0 0 5 2 0 0 0 9 tril(A) ；下三角阵 ans = 1 0 0 0 5 1 0 0 3 0 5 0 11 15 0 9 rank(A) ；秩 ans = 4 norm(A) ；范数 ans = 21.3005 cond(A) ；条件数 ans = 11.1739 trace(A) ；迹 ans = 16 （2）略 6. A=[1 1 0.5;1 1 0.25;0.5 0.25 2] A = 1.0000 1.0000 0.5000 1.0000 1.0000 0.2500 0.5000 0.2500 2.0000 [V,D]=eig(A) V = 0.7212 0.4443 0.5315 -0.6863 0.5621 0.4615 -0.0937 -0.6976 0.7103 D = -0.0166 0 0 0 1.4801 0 0 0 2.5365 附件2 《matlab程序设计与应用》实验指导书 山东建筑大学信息与电气工程学院 前言 一、实验目的 本课程是电气工程及其自动化、自动化、电力工程与管理专业本科生的 学科基础选修课，它在线性代数、信号分析和处理、控制系统设计和仿真等 方面有着广泛的应用。主要是学习MATLAB的语法规则、基本命令和使用环境，使学生掌握MATLAB的基本命令和基本程序设计方法，提高使用该语言 的应用能力，具有使用MATLAB语言编程和调试的能力，以便为后续多门课 程使用该语言奠定必要的基础。 通过上机实验，使学生掌握MATLAB在线帮助功能的使用、熟悉MATLAB运行环境和MATLAB语言的主要特点，掌握MA TLAB语言的基本 语法规则及基本操作命令的使用，学会M文件的建立和使用方法以及应用MATLAB实现二维和三维图形的绘制方法，具有使用MATLAB语言编程和 调试的能力。 二、实验前预习 每次实验前，学生须仔细阅读本实验指导书的相关内容，明确实验目的、要求；明确实验步骤、测试数据及需观察的现象；复习与实验内容有关的理论知识；预习仪器设备的使用方法、操作规程及注意事项；做好预习要求中提出的其它事项。 三、实验注意事项 1．实验开始前，应先检查本组电脑设备是否工作正常，matlab软件工作是否正常。 2．实验时每位同学应独立完成实验任务，避免抄袭。 3．实验后应及时将实验数据进行记录与存盘，避免因电脑故障或其它原因造成实验数据的丢失。 4．实验中严格遵循电脑操作规程，如电脑发生异常现象，应立即切断电源， 报告指导教师检查处理。 5．测量数据或观察现象要认真细致，实事求是。 6．未经许可，不得对电脑进行软件及硬件的更改操作。 7．实验结束后，实验记录交指导教师查看并认为无误后，方可关掉电脑离开。 8．爱护公物，发生仪器设备等损坏事故时，应及时报告指导教师，按有关实验管理规定处理。 9．自觉遵守学校和实验室管理的其它有关规定。 四、实验总结 每次实验后，应对实验进行总结，即实验数据进行整理，绘制波形和图表，分析实验现象，撰写实验报告。实验报告除写明实验名称、日期、实验者姓名、同组实验者姓名外，还包括： 1．实验目的； 2．实验内容 3．实验流程图； 4．程序命令清单； 5．运行结果； 6．实验的收获与体会； 7．回答每项实验的有关问答题。 t t t 机械优化设计作业 1.用二次插值法求函数?( )= ( +1)( - 2)2 极小值，精度 e=0.01。 在 MA TLAB 的 M 文件编辑器中编写的 M 文件，如下： f=inline('(t+1)*(t -2)^2','t') a=0;b=3;epsilon=0.01; t1=a;f1=f(t1); t3=b;f3=f(t3); t2=0.5*(t1+t3);f2=f(t2); c1=(f3-f1)/(t3-t1); c2=((f2-f1)/(t2-t1)-c1)/(t2-t3); t4=0.5*(t1+t3-c1/c2);f4=f(t4); k=0; while(abs(t4-t2)>=epsilon) if t2 3.用牛顿法、阻尼牛顿法及变尺度法求函数 的极小点。( ) ( ) ( )21121 22, xxxxxf -+-= 4 2 (1)在用牛顿法在 MATLAB 的 M 文件编辑器中编写的 M 文件，如下： function [x,fx,k]=niudunfa(x0) syms x1 x2 f=(x1-2)^4+(x1-2*x2)^2; fx=0; v=[x1,x2]; df=jacobian(f,v); df=df.'; G=jacobian(df,v); epson=1e -12; g1=subs(df,{x1,x2},{x0(1,1),x0(2,1)}); G1=subs(G,{x1,x2},{x0(1,1),x0(2,1)}); k=0; p=-G1\g1; x0=x0+p; while(norm(g1)>epson) p=-G1\g1; x0=x0+p; g1=subs(df,{x1,x2},{x0(1,1),x0(2,1)}); G1=subs(G,{x1,x2},{x0(1,1),x0(2,1)}); k=k+1; end x=x0; fx=subs(f,{x1,x2},{x(1,1),x(2,1)}); 运行结果如下： >> [x,fx,k]=niudunfa([1;1]) x =1.9999554476059523381489991377897 0.99997772380297616907449956889483 fx =0.0000000000000000039398907941382470301534502947647 k =23 (2)用阻尼牛顿法在 MA TLAB 的 M 文件编辑器中编写的 M 文件，如下： function [x,fx,k]=zuniniudunfa(x0)%阻尼牛顿法 syms x1 x2 f=(x1-2)^4+(x1-2*x2)^2; fx=0; v=[x1,x2]; df=jacobian(f,v); df=df.'; G=jacobian(df,v); epson=1e -12;%停机原则 创新实践报告 二 O 一四年十月十五日 实习报告 实习名称Matlab程序设计实训 专业班级**** 姓名*** 学号*** 成绩评定 电气与信息工程学院 和谐勤奋求是创新 实习考核和成绩评定办法 1．实习成绩的考核由指导教师根据实习表现、实习报告、实习成果、现场操作、设计、口试或笔试等几个方面，给出各项权重，综合评定。该实习考核教研室主任审核，主管院长审批备案。 2．成绩评定采用五级分制，即优、良、中、及格、不及格。 3．参加本次实习时间不足三分之二或旷课四天以上者，不得参加本次考核，按不及格处理。 4．实习结束一周内，指导教师提交实习成绩和实习总结。 5．实习过程考核和实习成绩在教师手册中有记载。 实习报告内容 实习报告内容、格式各专业根据实习类别（技能实习、认识实习、生产实习、毕业实习等）统一规范，经教研室主任审核、主管院长审批备案。 注：1. 实习任务书和实习指导书在实习前发给学生，实习任务书放置在实习报告封面后和正文目录前。 2. 为了节省纸张，保护环境，便于保管实习报告，统一采用A4纸，实习报告建议双面打印（正文采 用宋体五号字）或手写，右侧装订。 2016/2017学年第二学期 《Matlab程序设计实训》计划 指导教师：****** 班级：自动化1541、2班 实习地点：1教8楼机房 一、实习目的 实习的目的：了解并掌握MATLAB软件的各种应用开发环境；了解并掌握MATLAB软件强大的科学运算功能；了解并掌握利用MATLAB语言程序设计流程进行程序设计的方法；了解并掌握MATLAB的图形处理功能；了解MATLAB软件的交互式仿真功能；了解并掌握MATLAB软件在控制系统建模、分析以及设计中的应用。 实习的任务：通过实习可以使学生将MATLAB软件与控制系统理论及仿真相结合，加深对控制系统理论知识的理解，验证理论知识的正确性，提高软件的实际操作能力，可以有效地培养学生分析及解决实际问题的能力，不断提高自学能力，养成良好的科学态度以及实事求是、严谨踏实的学习风气，成为真正的应用型本科人才。通过实训达到边学边练、理论与实际操作相结合的学习目的。 二、实习内容 掌握MATLAB编程环境、基础知识；掌握矩阵生成及运算、数组运算；矩阵函数；掌握常用绘图命令及方法；掌握基本图形处理、二维图形处理方法；掌握可视化M文件程序设计及调试；掌握Simulink仿真工具的建立和基本仿真方法。 了解并掌握基于MATLAB的控制系统数学模型及分析方法；了解并掌握基于MATLAB的控制系统PID控制器的基本设计方法。综合应用MATLAB程序设计语言，根据设计任务要求，完成软件程序的编制并调试运行。 四、实习要求 实习期间，严格遵守课堂纪律，不许迟到、早退；实习期间，学习态度端正，认真上机练习；熟练掌握MATLAB软件的应用方法；完成实习报告一份，内容充实，写出实习的体会与收获；实习报告打印成文。 创新实践报告 报告题目: 基于matlab的通信系统仿真学院名称: 信息工程学院 姓名: 班级学号: 指导老师: 二O一四年十月十五日 一、引言 现代社会发展要求通信系统功能越来越强,性能越来越高,构成越来越复杂;另一方面,要求通信系统技术研究与产品开发缩短周期,降低成本,提高水平。这样尖锐对立的两个方面的要求,只有通过使用强大的计算机辅助分析设计技术与工具才能实现。在这种迫切的需求之下,MA TLAB应运而生。它使得通信系统仿真的设计与分析过程变得相对直观与便捷,由此也使得通信系统仿真技术得到了更快的发展。通信系统仿真贯穿着通信系统工程设计的全过程,对通信系统的发展起着举足轻重的作用。通信系统仿真具有广泛的适应性与极好的灵活性,有助于我们更好地研究通信系统性能。通信系统仿真的基本步骤如下图所示: 二、仿真分析与测试 (1)随机信号的生成 利用Matlab中自带的函数randsrc来产生0、1等概分布的随机信号。源代码如下所示: global N N=300; global p p=0、5; source=randsrc(1,N,[1,0;p,1-p]); (2)信道编译码 1、卷积码的原理 卷积码(convolutional code)就是由伊利亚斯(p 、Elias)发明的一种非分组码。在前向纠错系统中,卷积码在实际应用中的性能优于分组码,并且运算较简单。 卷积码在编码时将k 比特的信息段编成n 个比特的码组,监督码元不仅与当前的k 比特信息段有关,而且还同前面m=(N-1)个信息段有关。 通常将N 称为编码约束长度,将nN 称为编码约束长度。一般来说,卷积码中k 与n 的值就是比较小的整数。将卷积码记作(n,k,N)。卷积码的编码流程如下所示。 可以瞧出:输出的数据位V1,V2与寄存器D0,D1,D2,D3之间的关系。根据模2加运算特点可以得知奇数个1模2运算后结果仍就是1,偶数个1模2运算后结果就是0。 2、译码原理 卷积码译码方法主要有两类:代数译码与概率译码。代数译码主要根据码本身的代数特性进行译码,而信道的统计特性并没有考虑在内。目前,代数译码的主要代表就是大数逻辑解码。该译码方法对于约束长度较短的卷积码有较好的效果,并且设备较简单。概率译码,又称最大似然译码,就是基于信道的统计特性与卷积 码的特点进行计算。在现代通信系统中,维特比译码就是目前使用最广泛的概率 译码方法。 02 1V D D =⊕01232V D D D D =⊕⊕⊕matlab程序设计实践-牛顿法解非线性方程
MATLAB通信系统仿真心得体会
简述基于MATLAB的优化设计
1实验一 Matlab程序设计与M文件(1)
MATLAB 2psk通信系统仿真报告
matlab程序设计作业
MATLAB实现通信系统仿真实例
matlab程序设计与应用(第二版)第三章部分课后答案
matlab程序设计与应用实验指导书
机械优化设计MATLAB程序
(完整版)基于matlab的通信系统仿真毕业论文
报 告 题 目： 学 院 名 称： 姓 名：
基于 matlab 的通信系统仿真 信息工程学院 余盛泽
班 级 学 号： 指 导 老 师： 温 靖
目录
一、引言........................................................................................................................ 3 二、仿真分析与测试 ................................................................................................... 4
2.1 随机信号的生成 ............................................................................................................... 4 2.2 信道编译码 ........................................................................................................................ 4 2.2.1 卷积码的原理 ........................................................................................................ 4 2.2.2 译码原理 ................................................................................................................ 5 2.3 调制与解调 ....................................................................................................................... 5 2.3.1 BPSK 的调制原理 .................................................................................................. 5 2.3.2 BPSK 解调原理 ...................................................................................................... 6 2.3.3 QPSK 调制与解调 ................................................................................................. 7 2.4 信道 .................................................................................................................................... 8Matlab程序设计实训
基于matlab的通信系统仿真