1引言
1.1设计目的
本报告依照传统模拟调制的规范，给出了AM调制解调的具体流程，重点研究了系统中各阶段信号时域和频谱波形以及频谱的搬移变化，为AM调制解调系统信号波形的进一步深入研究做了基础。
1.2术语定义
本文档使用以下关键术语和简略语。
1.3参考资料
[1]通信原理（第六版）樊昌信曹丽娜编著国防工业出版社2007年1月
[2]
1.4文档组织
报告第二部分给出了AM调制解调的基本原理；第三部分给出了系统在MATLAB里面的程序调试及仿真；第四部分给出了各仿真模块输出时域和频域波形，并对比发射信号和接收信号的时域波形；第五部分对报告进行了总结。
2 AM 调制解调
图1 AM 调制解调系统框图
图1显示给出了用于AM 调制解调的系统框图。从图中可知，发送端信源信号()m t 经AM 调制器的调制后搭载高频载波发送出去，发送信道中经历加性高斯白噪声的干扰。接收端信号经历AM 解调器的解调输出，最终得到信宿信号^
()m t 。
2.1 AM 调制
A o
cos(?c t)
图2 AM 调制模型
图2显示给出了AM 调制的原理模型。从图中可知发送信号()m t 和直流分量0A 叠加后乘以高频载波cos()c t ω后即可形成AM 调制信号。
具体时域表波形为：
[]00()()cos ()
cos ()()cos ()
AM c c c s t A m t t A t m t t ωωω=+=+ (1)
对应的频谱波形为：
[][]01
()()()()()22
AM c c c c A S f M f f M f f f f f f δδ=
++-+++- (2)
2.2 AM 解调
c
图3 AM 解调模型
图3显示给出了AM 解调的数学模型。由上图可知，解调端信道输出信号()m s t 乘以跟发送端同频同相的高频载波cos()c t ω后，经低通滤波器提取低频分量，即可得到原始的基带调制信号。具体理论推导如下：
送入解调器的AM 的表达为
[]0()()cos()m c s t A m t t ω=+? (3)
与同频同相的相干载波相()=cos()c c t t ω乘后得
[][][]2000()()cos ()
11
()()cos(2)22
p c c s t A m t t A m t A m t t ωω=+?=?++?+? (4) 经历低通滤波器滤除高频信号后得
[]01
()()2
d s t A m t =
+ (5) 再经过隔直流电容后
1
()()2
s t m t = (6)
3 基于MATLAB的AM仿真
本仿真利用MATLAB自带函数仿真了AM系统各阶段信号的时域和频域波形，并对比发射信号和解调恢复信号的波形。
3.1仿真基本参数
系统仿真前定义采样间隔ts,载波频率fc,噪声系数sigma,和时间长度t。具体如下：%----------------------------------------------------------------
ts=1.e-4;
t=-0.04:ts:0.04 ;
fc=500;
sigma=0.3 ;
3.2生成调制信号
仿真中采用时间长度[-0.02 0.02]高度为1的三角脉冲波形作为调制信号。
%---------------------------------------------------------------- % 生成调制信号
m_sig=tripuls(t+0.01,0.02)-tripuls(t-0.01,0.02)；%调整信号m_sig由两个三角脉冲波形左右平移后叠加生成。
3.3 AM调制器
在MA TLAB中使用如下的语句对发射信号进行调制：
%---------------------------------------------------------------- % 调制过程
s_am=(1+m_sig).*cos(2*pi*fc*t)+sigma*randn(size(t))；
%AM信号由调制信号m_sig叠加直流分量后乘以载波并混以噪声形成。
cos(2*pi*fc*t)；% 载波信号，fc为载波频率。
randn(size(t));%代表服从正态分布的噪声。
3.4相干解调器
乘以接收信号，再经相干解调依靠接收端提供跟发送端同频同相的高频载波cos()c t
低通滤波器提取低频分量来恢复出原调制信号。具体程序实现如下：
%---------------------------------------------------------------- % 产生本地接收载波
s_carr =cos(2*pi*fc*t)；
% 同步解调
s_dem=s_am.*s_carr；
% 定义lfft变量
Lfft=length(t)；%取时间区间的长度并赋值给Lfft。
Lfft=2^ceil(log2(Lfft)+1)；%将Lfft重新转化成为更大的偶数Lfft。
ceil；%向正无穷大的方向取整
% 绘制解调后信号频谱
S_dem=fftshift(fft(s_dem,Lfft)/(length(t)));
% 生成低通滤波器
h=fir1(60,[B_m*ts])；%设计低通滤波器频域相应，其中滤波器阶数60,截至频率
75Hz。
% 低通滤波
s_rec=filter(h,1,s_dem)；%理想低通滤波器filter滤除s_dem中的高频分量得恢复信号s_rec。
% 绘制恢复信号频谱
S_rec=fftshift(fft(s_rec,Lfft)/(length(t)))；
4 仿真结果曲线
4.1发送信号波形和频谱
-0.025-0.02-0.015-0.01-0.005
00.0050.010.0150.020.025
-2
-1.5
-1-0.500.511.5
2t(sec)
m (t )
图4发送信号时域波形
图4显示给出了用于调制的发送信号时域波形。图中横坐标和纵坐标分别对应表示时间和信号幅值。从图中可以明显看出发送信号为三角波信号。
-150
-100
-50
050
100
150
00.050.10.150.20.250.30.350.40.45
0.5f(Hz)
M (f )
图5 发送信号的频谱
图5显示给出了三角波发送信号对应的频谱。图中横坐标表示频率，纵坐标表示频谱幅值。从图中可以看出信号频谱主要集中在低频段，而且信号带宽较窄，大约在150Hz 。
4.2载波信号波形和频谱
-5
-4-3-2-1
012345
x 10
-3
-1.5-1
-0.5
0.51
1.5
t(sec)
s
c a r r
( t )
图6 载波信号时域波形
图6显示给出了发送端调制基带信号所使用的载波波形。图中横坐标表示载波持续时间，横坐标表示载波幅度。从图中可知本次调制使用的载波是幅度为1，周期为-3
2.010s ，频率为500Hz 的余玄波。
-600
-400
-200
0200
400
600
0.050.10.150.20.250.3
0.350.40.45
0.5f(Hz)
S
c a r r
( f )
图7 载波信号频谱
图7显示给出了系统发射信号所需载波信号频谱。图中横坐标表示频率，纵坐标表示信号频谱幅值。从图中可以明显看出载波频谱的中心频率分别为+500Hz ，-500Hz 。
4.3 AM 信号波形和频谱
-0.025-0.02-0.015-0.01-0.005
00.0050.010.0150.020.025
-2
-1.5
-1-0.500.5
11.5
2t(sec)
s
A M
( t )
图8 AM 信号波形(SNR=19.72dB)
图8显示给出了系统中经过AM 调制后信号的时域波形。图中横坐标和纵坐标分别对应表示时间和信号幅值。从图中可知，已调制信号的外包络仍然保持着跟调制信号相同的包络特性。
-600
-400
-200
0200
400
600
0.050.10.150.20.250.3
0.350.40.45
0.5f(Hz)
S
A M
( f )
图9 AM 信号频谱波形（SNR=19.72dB ）
图9显示给出了经历调制后信号所对应的频谱。频谱图中横坐标表示频率，纵坐标表示频谱幅值。从图中可知：原发射信号经过调制后，频谱明显从原来的低频部分搬移到载波频率对应的高频部分。但由于已调制信号中直流分量的作用，调制后信号频谱的幅值相比原调制信号频谱幅值并没有完全减半。
4.4相干解调波形和频谱
-0.025-0.02-0.015-0.01-0.005
00.0050.010.0150.020.025
-2
-1.5
-1-0.500.5
11.5
2t(sec)
s
r e c t
( t )
图10 经过相干解调后的时域信号波形
图10显示给出了相干解调中信号乘以跟发送载波同频的本地余玄波后的时域波形。图中横坐标为时间轴，纵坐标表示信号幅值。从上述图中可知：AM 信号完整的输入波形转化为单一极性输出，同时信号的频率增加。
-1000-500
05001000
0.1
0.2
0.3
0.4
0.5
f(Hz)
S
r e c t
( f )
图11 相干解调频谱
图11 显示给出了信号通过相干解调乘以本地载波后对应频谱图。图中横坐标表示频率变化范围，纵坐标表示频谱幅值。由上图对比图9可知，乘以本地载波后的信号频谱将AM 信号频谱又重新的搬移，而在高频段仍然保留频谱分量。
4.5恢复信号波形和频谱
-0.025-0.02-0.015-0.01-0.005
00.0050.010.0150.020.025
-0.5
0.5
1
t(sec)
s
r e c
( t )
图12 经过低通滤波器后恢复信号时域波形（SNR=19.72dB ）
图12显示给出了经过低通滤波器后恢复出的信号的时域波形。图中横坐标表示信号长度时间坐标轴，纵坐标表示恢复信号的幅值。从图中可看出，恢复的信号波形基本上跟发送端发送信号波形吻合。但由于受到噪声的影响，信号的包络发生了抖动。
-200
-150
-100
-50
050
100
150
200
00.050.10.150.20.250.3
0.350.40.45
0.5f(Hz)
S
r e c
( f )
图13 经历低通滤波器恢复信号频谱图
图13显示给出了相干解调信号经过低通滤波器后对应的频谱。上述图像中横坐标代表频率变化，纵坐标表示频谱幅度变化。对比原发射信号频谱可知，此恢复出的信号频谱基本跟发射信号频谱波形吻合并且带宽大约在150Hz ，说明能较好的恢复出原调制信号。
5总结
本报告使用MATLAB软件仿真了AM调制解调系统中各阶段信号的时域和频谱波形，更直观的研究了信号在通过系统中各阶段的波形和频谱搬移变化情况。
实验五BPSK调制及解调实验 一、实验目的 1、掌握BPSK调制和解调的基本原理； 2、掌握BPSK数据传输过程，熟悉典型电路； 3、了解数字基带波形时域形成的原理和方法，掌握滚降系数的概念； 4、熟悉BPSK调制载波包络的变化； 5、掌握BPSK载波恢复特点与位定时恢复的基本方法； 二、实验器材 1、主控&信号源、9号、13号模块各一块 2、双踪示波器一台 3、连接线若干 三、实验原理 1、BPSK调制解调（9号模块）实验原理框 PSK调制及解调实验原理框图 2、BPSK调制解调（9号模块）实验框图说明 基带信号的1电平和0电平信号分别与256KHz载波及256KHz反相载波相乘，叠加后得到BPSK调制输出；已调信号送入到13模块载波提取单元得到同步载波；已调信号与相干载波相乘后，经过低通滤波和门限判决后，解调输出原始基带信号。 四、实验步骤 实验项目一 BPSK调制信号观测（9号模块） 概述：BPSK调制实验中，信号是用相位相差180°的载波变换来表征被传递的信息。本项目通过对比观测基带信号波形与调制输出波形来验证BPSK调制原理。 1、关电，按表格所示进行连线。
2、开电，设置主控菜单，选择【主菜单】→【通信原理】→【BPSK/DBPSK数字调制解调】。将9号模块的S1拨为0000，调节信号源模块W3使256 KHz载波信号峰峰值为3V。 3、此时系统初始状态为：PN序列输出频率32KHz。 4、实验操作及波形观测。 （1）以9号模块“NRZ-I”为触发，观测“I”； （2）以9号模块“NRZ-Q”为触发，观测“Q”。 （3）以9号模块“基带信号”为触发，观测“调制输出”。 思考：分析以上观测的波形，分析与ASK有何关系？ 实验项目二 BPSK解调观测（9号模块） 概述：本项目通过对比观测基带信号波形与解调输出波形，观察是否有延时现象，并且验证BPSK解调原理。观测解调中间观测点TP8，深入理解BPSK解调原理。 1、保持实验项目一中的连线。将9号模块的S1拨为“0000”。 2、以9号模块测13号模块的“SIN”，调节13号模块的W1使“SIN”的波形稳定，即恢复出载波。 3、以9号模块的“基带信号”为触发观测“BPSK解调输出”，多次单击13号模块的“复位”按键。观测“BPSK解调输出”的变化。 4、以信号源的CLK为触发，测9号模块LPF-BPSK，观测眼图。 思考：“BPSK解调输出”是否存在相位模糊的情况？为什么会有相位模糊的情况？ 五、实验报告 1、分析实验电路的工作原理，简述其工作过程； 输入的基带信号由转换开关转接后分成两路，一路经过差分编码控制256KHz的载频，另一路经倒相去控制256KHz的载频。???解调采用锁相解调，只要在设计锁相环时，使它锁定在FSK的一个载频上此时对应的环路滤波器输出电压为零，而对另一载频失锁，则对应的环路滤波器输出电压不为零，那末在锁相环路滤波器输出端就可以获得原基带信号的信息。? 2、分析BPSK调制解调原理。 调制原理是：基带信号先经过差分编码得到相对码，再根据相对码进行绝对调相， 即将相对码的1电平和0电平信号分别与256K载波及256K反相载波相乘，叠加后得到DBPSK 调制输出。?
目录 一. FSK理论知识………………………………………………… 1.1FSK概念………………………………………………………………… 1.22FSK信号的波形及时间表示式………………………………………… 1.32FSK信号的产生方法…………………………………………………… 1.42FSK信号的功率谱密度………………………………………………… 1.52FSK信号的解调………………………………………………………… 1.6FSK的误码性能…………………………………………………………… 二.用MATLAB进行FSK原理及误码性能仿真……… 三、结论…………………………………………… 四、参考文献…………………………………………、 五、源程序……………………………………………
1、FSK理论知识 频率调制的最简单形式是二进制频率键控(FSK，frequency-shift keying)。FSK是调制解调器通过电话线路发送比特的方法。每个比特被转换为一个频率，0由较低的频率表示，1由较高的频率表示。 1.1、FSK概念 传“0”信号时，发送频率为f1的载波; 传“1”信号时，发送频率为f2的载波。可见，FSK是用不同频率的载波来传递数字消息的。 实现模型如下图： 1.2、2FSK信号的波形及时间表示式 根据上图模型的实现可以得到2FSK的信号波形如图：
2FSK信号的时间表达式为： 由以上表达式可见，2FSK信号由两个2ASK信号相加构成。 注意：2FSK有两种形式: （1）相位连续的2FSK； （2）相位不连续的2FSK。 在这里，我们只讨论相位不连续的频移键控信号，这样更具有普遍性。 1.3、2FSK信号的产生方法 2FSK信号的产生方法：2FSK信号可以两类方法来产生。 一是采用模拟调频的方法来产生（图1）；另一种方法是采用键控法（图2）； 图1.3-1 图1.3-2 1.4、2FSK信号的功率谱密度
数学实验报告 班级: 学号： 姓名： 实验序号：1 日期：年 月 日 实验名称：特殊函数与图形 ◆ 问题背景描述：绘图是数学中的一种重要手段，借助图形，可以使抽象的对象得到 明白直观的体现，如函数的性质等。同时，借助直观的图形，使初学者更容易接受新知识，激发学习兴趣。 ◆ 实验目的：本实验通过绘制一些特殊函数的图形，一方面展示这些函数的特点属性， 另一方面，就 Matlab 强大的作图功能作一个简单介绍。 实验原理与数学模型： 1、 球2222x y z R ++= ，x=Rsin φcos θ, y= Rsin φsin θ, z= cos φ, 0≤θ≤2π , 0≤φ≤π 环面 222222222()4(),(cos )cos ,x y z a r a x y x a r φθ+++-=+=- (cos )sin ,sin ,02,02y a r z r φθφφπθπ=-=≤≤≤≤ 2、 平面摆线：2 22 31150,(sin ),(1cos ),0233 x y x a t t y a t t π+-==-=-≤≤ 3、 空间螺线：（圆柱螺线）x=acost , y=asint , z=bt ;（圆锥螺线）22 cos ,sin ,x t t y t t z t === 4、 椭球面sin cos ,sin sin ,cos ,02,0x a y b z c φθφθφθπφπ===≤<≤≤ 双叶双曲面3 tan cos ,tan sin ,sec ,02,22 x a y b z c π φθφθφθπφπ===≤<- << 双曲抛物面2 sec ,tan 2 u x au y bu z θθ=== 实验所用软件及版本：mathematica(3.0) 主要内容（要点）： 1、 作出下列三维图形（球、环面） 2、 作出下列的墨西哥帽子 3、 作出球面、椭球面、双叶双曲面，单叶双曲面的图形 4、 试画出田螺上的一根螺线 5、 作出如图的马鞍面
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);
《通信原理》实验报告 实验一：抽样定理和PAM调制解调实验 系别：信息科学与工程学院 专业班级：通信工程1003班 学生姓名：陈威 同组学生：杨鑫 成绩： 指导教师：惠龙飞 （实验时间：2012 年 12 月 7 日——2012 年 12 月28日） 华中科技大学武昌分校
1、实验目的 1对电路的组成、波形和所测数据的分析，加深理解这种调制方法的优缺点。 2.通过脉冲幅度调制实验，使学生能加深理解脉冲幅度调制的原理。 2、实验器材 1、信号源模块 一块 2、①号模块 一块 3、60M 双踪示波器 一台 4、连接线 若干 3、实验原理 3.1基本原理 1、抽样定理 图3-1 抽样与恢复 2、脉冲振幅调制（PAM ） 所谓脉冲振幅调制，即是脉冲载波的幅度随输入信号变化的一种调制方式。如果脉冲载波是由冲激脉冲组成的，则前面所说的抽样定理，就是脉冲增幅调制的原理。 自然抽样 平顶抽样 ) (t m ) (t T
图3-3 自然抽样及平顶抽样波形 PAM方式有两种：自然抽样和平顶抽样。自然抽样又称为“曲顶”抽样，(t)的脉冲“顶部”是随m(t)变化的，即在顶部保持了m(t)变已抽样信号m s 化的规律（如图3-3所示）。平顶抽样所得的已抽样信号如图3-3所示，这里每一抽样脉冲的幅度正比于瞬时抽样值，但其形状都相同。在实际中，平顶抽样的PAM信号常常采用保持电路来实现，得到的脉冲为矩形脉冲。 四、实验步骤 1、将信号源模块、模块一固定到主机箱上面。双踪示波器，设置CH1通道为同步源。 2、观测PAM自然抽样波形。 （1）将信号源上S4设为“1010”，使“CLK1”输出32K时钟。 （2）将模块一上K1选到“自然”。 （3）关闭电源，连接 表3-1 抽样实验接线表 (5)用示波器观测信号源“2K同步正弦波”输出，调节W1改变输出信号幅度，使输出信号峰-峰值在1V左右。在PAMCLK处观察被抽样信号。CH1接PAMCLK（同步源），CH2接“自然抽样输出”（自然抽样PAM信号）。
摘要 在模拟通信系统中，由模拟信源产生的携带有信息的消息经过传感器转换成电信号。模拟基带信号在经过调制将低通频谱搬移到载波频率上适应信道，最终解调还原成电信号。本文应用了频率调制法产生调制解调信号。本论文中主要通过对SIMULINK工具箱的学习和使用，利用其丰富的模板以及本科对通信原理知识的掌握，完成了FM信号的调制与解调，以及用SIMULINK进行设计和仿真。首先利用简单的正玄波信号发生器作为信源，对模拟信号进行FM调制解调原理的仿真。 关键词：调制解调；FM ；MATLAB；SIMULINK仿真
Abstract In the simulation of communication systems, generated by the analog source carrying a message through the sensor into electrical signals. Analog baseband signal after the modul- -ation of the low pass spectrum to carrier frequency to adapt to the channel, the final reducti- -on into electrical signal demodulation. This paper applied the frequency modulation method to generate the signal modulation and demodulation. Mainly through the study and use of SIMULINK toolbox in this thesis, with its rich template and undergraduate course on comm--unication theory knowledge,the modulation and demodulation of FM signal, as well as the design and simulation with SIMULINK. Firstly, sine wave signal generator is simple as the source, simulation FM modulation anddemodulation principle of analogue signals. Then, using the song as the source. Keywords: modulation and demodulation；FM; MATLAB; SIMULINK simulation
实验一MATLAB操作基础 实验目的和要求： 1、熟悉MATLAB的操作环境及基本操作方法。 2、掌握MATLAB的搜索路径及设置方法。 3、熟悉MATLAB帮助信息的查阅方法 实验内容： 1、建立自己的工作目录，再设置自己的工作目录设置到MA TLAB搜索路径下，再试 验用help命令能否查询到自己的工作目录。 2、在MA TLAB的操作环境下验证课本；例1-1至例1-4，总结MATLAB的特点。 例1-1
例1-2 例1-3 例1-4
3、利用帮助功能查询inv、plot、max、round等函数的功能。 4、完成下列操作： （1）在matlab命令窗口输入以下命令： x=0:pi/10:2*pi; y=sin(x); （2）在工作空间窗口选择变量y，再在工作空间窗口选择回绘图菜单命令或在工具栏中单击绘图命令按钮，绘制变量y的图形，并分析图形的含义。
5、访问mathworks公司的主页，查询有关MATLAB的产品信息。 主要教学环节的组织： 教师讲授实验目的、开发环境界面、演示实验过程，然后同学上机练习。 思考题： 1、如何启动与退出MA TLAB集成环境？ 启动： （1）在windows桌面，单击任务栏上的开始按钮，选择‘所有程序’菜单项，然后选择MA TLAB程序组中的MA TLABR2008b程序选项，即可启动 MATLAB系统。 （2）在MA TLAB的安装路径中找到MA TLAB系统启动程序matlab.exe，然后运行它。 （3）在桌面上建立快捷方式后。双击快捷方式图标，启动MA TLAB。 退出： （1）在MA TLAB主窗口file菜单中选择exitMATLAB命令。 （2）在MA TLAB命令窗口中输入exit或quit命令。 （3）单击MATLAB主窗口的关闭按钮。 2、简述MATLAB的主要功能。 MATLAB是一种应用于科学计算领域的数学软件，它主要包括数值计算和符 号计算功能、绘图功能、编程语言功能以及应用工具箱的扩展功能。 3、如果一个MATLAB命令包含的字符很多，需要分成多行输入，该如何处理？
MATLAB 仿真实验报告 课题名称：MATLAB 仿真——图像处理 学院：机电与信息工程学院 专业：电子信息科学与技术 年级班级：2012级电子二班 一、实验目的 1、掌握MATLAB处理图像的相关操作，熟悉相关的函数以及基本的MATLAB语句。 2、掌握对多维图像处理的相关技能，理解多维图像的相关性质 3、熟悉Help 命令的使用，掌握对相关函数的查找，了解Demos下的MATLAB自带的原函数文件。 4、熟练掌握部分绘图函数的应用，能够处理多维图像。 二、实验条件
MATLAB调试环境以及相关图像处理的基本MATLAB语句，会使用Help命令进行相关函数查找 三、实验内容 1、nddemo.m函数文件的相关介绍 Manipulating Multidimensional Arrays MATLAB supports arrays with more than two dimensions. Multidimensional arrays can be numeric, character, cell, or structure arrays. Multidimensional arrays can be used to represent multivariate data. MATLAB provides a number of functions that directly support multidimensional arrays. Contents ： ●Creating multi-dimensional arrays 创建多维数组 ●Finding the dimensions寻找尺寸 ●Accessing elements 访问元素 ●Manipulating multi-dimensional arrays操纵多维数组 ●Selecting 2D matrices from multi-dimensional arrays从多维数组中选择二维矩 阵 (1)、Creating multi-dimensional arrays Multidimensional arrays in MATLAB are created the same way as two-dimensional arrays. For example, first define the 3 by 3 matrix, and then add a third dimension. The CAT function is a useful tool for building multidimensional arrays. B = cat(DIM,A1,A2,...) builds a multidimensional array by concatenating（联系起来）A1, A2 ... along the dimension DIM. Calls to CAT can be nested（嵌套）. (2)、Finding the dimensions SIZE and NDIMS return the size and number of dimensions of matrices. (3)、Accessing elements To access a single element of a multidimensional array, use integer subscripts（整数下标）. (4)、Manipulating multi-dimensional arrays
通信原理课程设计报告
一. 2DPSK基本原理 1.2DPSK信号原理 2DPSK方式即是利用前后相邻码元的相对相位值去表示数字信息的一种方式。现假设用Φ表示本码元初相与前一码元初相之差，并规定：Φ＝0表示0码，Φ＝π表示1码。则数字信息序列与2DPSK信号的码元相位关系可举例表示如2PSK信号是用载波的不同相位直接去表示相应的数字信号而得出的，在接收端只能采用相干解调，它的时域波形图如图2.1所示。 图1.1 2DPSK信号 在这种绝对移相方式中，发送端是采用某一个相位作为基准，所以在系统接收端也必须采用相同的基准相位。如果基准相位发生变化，则在接收端回复的信号将与发送的数字信息完全相反。所以在实际过程中一般不采用绝对移相方式，而采用相对移相方式。 定义?Φ为本码元初相与前一码元初相之差，假设： ?Φ=0→数字信息“0”； ?Φ=π→数字信息“1”。 则数字信息序列与2DPSK信号的码元相位关系可举例表示如下： 数字信息： 1 0 1 1 0 1 1 1 0 1
DPSK信号相位：0 π π 0 π π 0 π 0 0 π 或：π 0 0 π 0 0 π 0 π π 0 2. 2DPSK信号的调制原理 一般来说，2DPSK信号有两种调试方法，即模拟调制法和键控法。2DPSK 信号的的模拟调制法框图如图1.2.1所示，其中码变换的过程为将输入的单极性不归零码转换为双极性不归零码。 图1.2.1 模拟调制法 2DPSK信号的的键控调制法框图如图1.2.2所示，其中码变换的过程为将输入的基带信号差分，即变为它的相对码。选相开关作用为当输入为数字信息“0”时接相位0，当输入数字信息为“1”时接pi。 图1.2.2 键控法调制原理图 码变换相乘 载波 s(t)e o(t)
基于MATLAB的FSK调制解调 学生姓名：段斐指导老师：吴志敏 摘要本课程设计利用MATLAB集成环境下的M文件，编写程序来实现FSK 的调制解调，并绘制出解调前后的时域和频域波形及叠加噪声时解调前后的时频波形，并观察解调前后频谱有何变化以加深对F SK信号解调原理的理解。对信号叠加噪声，并迚行解调，绘制出解调前后信号的时频波形，改变噪声功率迚行解调，根据运行结果和波形来分析该解调过程的正确性及信道对信号传输的影响。完成整个FSK的调制解调过程。程序开发平台为MATLAB7.1，使用其自带的M文件实现。运行平台为Windows 2000。 关键词：程序设计；FSK ；调制解调；MATLAB7.1；M文件 1引言 本课程设计是利用MATLAB集成环境下的M文件，编写程序来实现FSK 的调制解调，并绘制出解调前后的时域和频域波形及叠加噪声时解调前后的时频波形，根据运行结果和波形来分析该解调过程的正确性及信道对信号传输的影响。 1.1课程设计目的 此次课程设计的目的是熟悉MATLAB中M文件的使用方法，编写M文件实现FSK的调制和解调，绘制出FSK信号解调前后在时域和频域中的波形，观察调解前后频谱的变化，再对信号迚行噪声叠加后解调同样绘制解调前后的
信号时频波形，最后改变噪声功率迚行调解，分析噪声对信号传输造成的影响，加深对FSK信号解调原理的理解。 1.2课程设计要求 熟悉MATLAB中M文件的使用方法，并在掌握FSK调制解调原理的基础上，编写出F SK调制解调程序。在M文件环境下运行程序绘制出F SK信号解调前后在时域和频域中的波形，观察波形在解调前后的变化，对其作出解释，同时对信号加入噪声后解调，得到解调后的时频波形，分析噪声对信号传输造成的影响。解释所得到的结果。 1.3课程设计步骤 本课程设计采用M文件编写的方法实现二迚制的FSK的调制与解调，然后在信号中叠加高斯白噪声。一，调用dmode函数实现FSK的解调，并绘制出F SK信号调制前后在时域和频域中的波形，两者比较。二，调用ddemod函数解调，绘制出F SK信号解调前后在时域和频域中的波形，两者比较。三，调用awgn函数在新海中叠加不同信噪比的噪声，绘制在各种噪声下的时域频域图。最后分析结果。 1.4设计平台简介 Matlab是美国MathWorks公司开发的用于概念设计,算法开发,建模仿真,实时实现的理想的集成环境。是目前最好的科学计算类软件。 作为和Mathematica、Maple并列的三大数学软件。其强项就是其强大的矩阵计算以及仿真能力。Matlab的由来就是Matrix + Laboratory = Matlab，这个软件在国内也被称作《矩阵实验室》。Matlab提供了自己的编译器：全面兼容C++以及Fortran两大语言。Matlab 7.1于2005.9最新发布-完整版,提供了
实验一 MATLAB 及仿真实验（控制系统的时域分析） 一、实验目的 学习利用MATLAB 进行控制系统时域分析，包括典型响应、判断系统稳定性和分析系统的动态特性； 二、预习要点 1、 系统的典型响应有哪些？ 2、 如何判断系统稳定性？ 3、 系统的动态性能指标有哪些？ 三、实验方法 （一） 四种典型响应 1、 阶跃响应： 阶跃响应常用格式： 1、)(sys step ；其中sys 可以为连续系统，也可为离散系统。 2、),(Tn sys step ；表示时间范围0---Tn 。 3、),(T sys step ；表示时间范围向量T 指定。 4、),(T sys step Y =；可详细了解某段时间的输入、输出情况。 2、 脉冲响应： 脉冲函数在数学上的精确定义：0 ,0)(1)(0 ?==?∞ t x f dx x f 其拉氏变换为：) ()()()(1)(s G s f s G s Y s f === 所以脉冲响应即为传函的反拉氏变换。 脉冲响应函数常用格式： ① )(sys impulse ； ② ); ,();,(T sys impulse Tn sys impulse ③ ),(T sys impulse Y = （二） 分析系统稳定性 有以下三种方法： 1、 利用pzmap 绘制连续系统的零极点图； 2、 利用tf2zp 求出系统零极点； 3、 利用roots 求分母多项式的根来确定系统的极点 （三） 系统的动态特性分析 Matlab 提供了求取连续系统的单位阶跃响应函数step 、单位脉冲响应函数impulse 、零输入响应函数initial 以及任意输入下的仿真函数lsim.
基于MATLAB的FSK的实验报告 姓 1.1
实现对FSK的MATLAB仿真. 重点研究问题: (1) 对FSK的概念、组成以及性能分析方法有深入的研究； (2) FSK调制与解调的原理及应用MATLAB软件实现仿真的方案. 1.2 FSK信号的调制方法 移频键控(FSK)：用数字调制信号的正负控制载波的频率。当数字信号的振幅为正时载波频率为f1，当数字信号的振幅为负时载波频率为 f2。有时也把代表两个以上符号的多进制频率调制称为移频键控。移频键控能区分通路，但抗干扰能力不如移相键控和差分移相键控。他的主要调制方法有以下两种: 方法一： 用一个矩形脉冲序列对一个载波进行调频。 图2-3 2FSK信号的产生(一) 方法二:键控法 图2-4 2FSK信号的产生(二) 键控法是利用矩形脉冲()t b来控制开关电路对两个不同的独立频率源进行选通。
1.3 FSK解调的方法 常见的FSK解调方法有两种:相干解调法与非相干解调法.现在我将对这两种解法。 1.4 设计总思路 如下图所示,我将FSK的调制与FSK的解调独立开作为两个子函数,其中FSK调制的输出即可作为FSK解调的输入信号.最后设计一主函数main将两个子函数同时调用完成整个仿真过程。 图3-1 设计总思路图 2.1 FSK调制的仿真设计 本文主要是对2FSK进行调制,而2FSK可看做是基带信号与载波频率的结合就可.FSK的产生思路参考的是键控法,如图4
图3-2 2FSK信号的产生(二) 2.2 FSK解调的仿真设计 如上图所示的FSK信号的相干检测原理图,FSK信号可以采用两个乘法检测器进行相干检测. 上图中输入信号为2FSK信号加上噪声组成 带通滤波器2的设计类似滤波器1,只是更改频率为fc2就可.
实验一小球做自由落体运动内容：一小球竖直方向做自由落体，并无损做往返运动。程序： theta=0:0.01:2*pi x=cos(theta) y=sin(theta) l=1 v=1 while l<10 for t=1:10 y=y+(-1)^l*v*t plot(x,y,[-1,1],[-56,2],'.') axis equal pause(0.1) end l=l+1 end 结果：
-50 -40 -30 -20 -10 收获：通过运用小球自由落体规律，及（-1）^n 来实现无损往 返运动！ 实验二 旋转五角星 内容：一个五角星在圆内匀速旋转 程序：x=[2 2 2 2 2 2] y=[0 4/5*pi 8/5*pi 2/5*pi 6/5*pi 0] y1=2*sin(y) x1=2*cos(y) theta=0:4/5*pi:4*pi
x2=2*cos(theta) y2=2*sin(theta) plot(x,y,x1,y1,x2,y2) axis equal theta1=theta+pi/10 x2=2*cos(theta1) y2=2*sin(theta1) plot(x2,y2) axis equal theta=0:4/5*pi:4*pi for rot=pi/10:pi/10:2*pi x=2*cos(theta+rot) y=2*sin(theta+rot) plot(x,y) pause(0.1) end 结果：
-2 -1.5-1-0.500.51 1.52 -2-1.5-1-0.500.511.5 2 收获：通过theta1=theta+pi/10，我们可以实现五角星在圆内匀速 旋转！ 实验三 转动的自行车 内容：一辆自行车在圆内匀速转动 程序：x=-4:0.08:4; y=sqrt(16-x.^2); theta1=-pi/2:0.01*pi:3*pi/2; x3=0.5*cos(theta1); y3=0.5*sin(theta1); theta=-pi/2+0.02*pi for k=1:100
MATLAB/Simulink 电力系统建模与仿真 实验报告 姓名：****** 专业：电气工程及其自动化 班级：******************* 学号：*******************
实验一无穷大功率电源供电系统三相短路仿真 1.1 无穷大功率电源供电系统仿真模型构建 运行MATLAB软件，点击Simulink模型构建，根据电路原理图，添加下列模块： （1）无穷大功率电源模块（Three-phase source） （2）三相并联RLC负荷模块（Three-Phase Parallel RLC Load） （3）三相串联RLC支路模块（Three-Phase Series RLC Branch） （4）三相双绕组变压器模块（Three-Phase Transformer (Two Windings)） （5）三相电压电流测量模块（Three-Phase V-I Measurement） （6）三相故障设置模块（Three-Phase Fault） （7）示波器模块（Scope） （8）电力系统图形用户界面（Powergui） 按电路原理图连接线路得到仿真图如下： 1.2 无穷大功率电源供电系统仿真参数设置 1.2.1 电源模块 设置三相电压110kV，相角0°，频率50Hz，接线方式为中性点接地的Y形接法，电源电阻0.00529Ω，电源电感0.000140H，参数设置如下图：
1.2.2 变压器模块 变压器模块参数采用标幺值设置，功率20MVA，频率50Hz，一次测采用Y型连接，一次测电压110kV，二次侧采用Y型连接，二次侧电压11kV，经过标幺值折算后的绕组电阻为0.0033，绕组漏感为0.052，励磁电阻为909.09，励磁电感为106.3，参数设置如下图： 1.2.3 输电线路模块 根据给定参数计算输电线路参数为：电阻8.5Ω，电感0.064L，参数设置如下图： 1.2.4 三相电压电流测量模块 此模块将在变压器低压侧测量得到的电压、电流信号转变成Simulink信号，相当于电压、电流互感器的作用，勾选“使用标签（Use a label）”以便于示波器观察波形，设置电压标签“Vabc”，电流标签“Iabc”，参数设置如下图：
基于MATLAB的ＡSK调制解调实现
————————————————————————————————作者：————————————————————————————————日期： ?
长沙理工大学 《通信原理》课程设计报告 学院专业 班级学号 学生姓名指导教师 课程成绩完成日期2016年1月8日
课程设计成绩评定 学院专业 班级学号 学生姓名指导教师 课程成绩完成日期20１6年1月8日 指导教师对学生在课程设计中的评价 评分项目优良中及格不及格课程设计中的创造性成果 学生掌握课程内容的程度 课程设计完成情况 课程设计动手能力 文字表达 学习态度 规范要求 课程设计论文的质量 指导教师对课程设计的评定意见 综合成绩指导教师签字２01６年1月8日
课程设计任务书 城南学院通信工程专业 课程名称通信原理课程设计时间２０15/2016学年第一学期１7~19 周 学生姓名指导老师 题目基于MＡTLAＢ的ＡＳK调制解调实现 主要内容: 利用MATＬAＢ集成环境下的M文件，编写程序来实现AＳＫ的调制解调， 要求采样频率为３６0ＨZ,并绘制出解调前后的时域和频域波形及叠加噪声时解 调前后的时频波形,根据运行结果和波形来分析该解调过程的正确性及信道对信 号传输的影响。 要求： 1)熟悉MATLAB中M文件的使用方法,并在掌握ASK调制解调原理 的基础上，编写出ASＫ调制解调程序。 ２)绘制出ASK信号解调前后在时域和频域中的波形,并观察解调前后频谱有何变化以加深对ＡSK信号解调原理的理解。 3）对信号叠加噪声，并进行解调，绘制出解调前后信号的时频波形,改变噪声功率进行解调,分析噪声对信号传输造成的影响。 4)在老师的指导下，要求独立完成课程设计的全部内容,并按要求编写课 程设计学年论文,能正确阐述和分析设计和实验结果。 应当提交的文件: （１）课程设计学年论文。 (2)课程设计附件。
实验一 Matlab基础知识 一、实验目的： 1.熟悉启动和退出Matlab的方法。 2.熟悉Matlab命令窗口的组成。 3.掌握建立矩阵的方法。 4.掌握Matlab各种表达式的书写规则以及常用函数的使 用。 二、实验内容： 1.求[100，999]之间能被21整除的数的个数。（rem） 2.建立一个字符串向量，删除其中的大写字母。（find） 3.输入矩阵，并找出其中大于或等于5的元素。（find） 4.不采用循环的形式求出和式 63 1 2i i= ∑ 的数值解。(sum) 三、实验步骤： ●求[100，199]之间能被21整除的数的个数。（rem） 1.开始→程序→Matlab 2.输入命令： ?m=100:999; ?p=rem(m,21); ?q=sum(p==0) ans=43 ●建立一个字符串向量，删除其中的大写字母。（find） 1.输入命令：
?k=input('’,’s’)； Eie48458DHUEI4778 ?f=find(k>=’A’&k<=’Z’)； f=9 10 11 12 13 ?k(f)=[ ] K=eie484584778 ●输入矩阵，并找出其中大于或等于5的元素。（find） 1.输入命令： ?h=[4 8 10;3 6 9; 5 7 3]； ?[i,j]=find(h>=5) i=3 j=1 1 2 2 2 3 2 1 3 2 3 ●不采用循环的形式求出和式的数值解。（sum） 1.输入命令： ?w=1:63； ?q=sum(2.^w) q=1.8447e+019
实验二 Matlab 基本程序 一、 实验目的： 1. 熟悉Matlab 的环境与工作空间。 2. 熟悉M 文件与M 函数的编写与应用。 3. 熟悉Matlab 的控制语句。 4. 掌握if,switch,for 等语句的使用。 二、 实验内容： 1. 根据y=1+1/3+1/5+……+1/（2n-1）,编程求：y<5时最大n 值以及对应的y 值。 2. 编程完成，对输入的函数的百分制成绩进行等绩转换，90~100为优，80~89为良，70~79为中，60~69为及格。 3. 编写M 函数文件表示函数 ，并分别求x=12和56时的函数值。 4. 编程求分段函数 2226;03 56;0532 1;x x x x y x x x x x x x +-<≠=-+≤<≠≠-+且且及其它，并求输入x=[-5.0,-3.0,1.0,2.0,2.5,3.0,3.5]时的输出y 。 三、 实验步骤： 根据y=1+1/3+1/5+……+1/（2n-1）,编程求：y<5时最大n 值以及对应的y 值。 1. 打开Matlab ，新建M 文件 2. 输入命令： 51022-+x
实验报告 课程名称：控制理论（乙） 指导老师：林峰 成绩：__________________ 实验名称：MATLAB 仿真实验 实验类型：________________同组学生姓名：__________ 一、实验目的和要求（必填） 二、实验内容和原理（必填） 三、主要仪器设备（必填） 四、操作方法和实验步骤 五、实验数据记录和处理 六、实验结果与分析（必填） 七、讨论、心得 实验九 控制系统的时域分析 一、 实验目的： 1．用计算机辅助分析的办法，掌握系统的时域分析方法。 2．熟悉Simulink 仿真环境。 二、实验原理及方法： 系统仿真实质上就是对系统模型的求解，对控制系统来说，一般模型可转化成某个微分方程或差分方程表示，因此在仿真过程中，一般以某种数值算法从初态出发，逐步计算系统的响应，最后绘制出系统的响应曲线，进而可分析系统的性能。控制系统最常用的时域分析方法是，当输入信号为单位阶跃和单位冲激函数时，求出系统的输出响应，分别称为单位阶跃响应和单位冲激响应。在MATLAB 中，提供了求取连续系统的单位阶跃响应函数step ，单位冲激响应函数impulse ，零输入响应函数initial 等等。 二、实验内容： 二阶系统，其状态方程模型为 ? 1x -0.5572 -0.7814 1x 1 = + u ? 2x 0.7814 0 2x 0 1x y = [1.9691 6.4493] +[0] u 2x 四、实验要求： 1．编制MATLAB 程序，画出单位阶跃响应曲线、冲击响应曲线、系统的零输入响应、斜坡输入响应； （1）画出系统的单位阶跃响应曲线； A=[-0.5572 -0.7814;0.7814 0 ]; B=[1;0];
通信原理课程设计报告 一. 2DPSK基本原理 1.2DPSK信号原理 2DPSK方式即是利用前后相邻码元的相对相位值去表示数字信息的一种方式。现假设用Φ表示本码元初相与前一码元初相之差，并规定：Φ＝0表示0码，
Φ＝π表示1码。则数字信息序列与2DPSK信号的码元相位关系可举例表示如2PSK信号是用载波的不同相位直接去表示相应的数字信号而得出的，在接收端只能采用相干解调，它的时域波形图如图2.1所示。 图1.1 2DPSK信号 在这种绝对移相方式中，发送端是采用某一个相位作为基准，所以在系统接收端也必须采用相同的基准相位。如果基准相位发生变化，则在接收端回复的信号将与发送的数字信息完全相反。所以在实际过程中一般不采用绝对移相方式，而采用相对移相方式。 定义?Φ为本码元初相与前一码元初相之差，假设： ?Φ=0→数字信息“0”； ?Φ=π→数字信息“1”。 则数字信息序列与2DPSK信号的码元相位关系可举例表示如下： 数字信息： 1 0 1 1 0 1 1 1 0 1 DPSK信号相位：0 π π 0 π π 0 π 0 0 π 或：π 0 0 π 0 0 π 0 π π 0 2. 2DPSK信号的调制原理 一般来说，2DPSK信号有两种调试方法，即模拟调制法和键控法。2DPSK 信号的的模拟调制法框图如图1.2.1所示，其中码变换的过程为将输入的单极性不归零码转换为双极性不归零码。
图1.2.1 模拟调制法 2DPSK信号的的键控调制法框图如图1.2.2所示，其中码变换的过程为将输入的基带信号差分，即变为它的相对码。选相开关作用为当输入为数字信息“0”时接相位0，当输入数字信息为“1”时接pi。 图1.2.2 键控法调制原理图 3. 2DPSK信号的解调原理 2DPSK信号最常用的解调方法有两种，一种是极性比较和码变换法，另一种是差分相干解调法。 (1) 2DPSK信号解调的极性比较法 它的原理是2DPSK信号先经过带通滤波器，去除调制信号频带以外的在信道中混入的噪声，再与本地载波相乘，去掉调制信号中的载波成分，再经过低通滤波器去除高频成分，得到包含基带信号的低频信号，将其送入抽样判决器中进行抽样判决的到基带信号的差分码，再经过逆差分器，就得到了基带信号。它的原理框图如图1.3.1所示。 码变换相乘 载波 s(t)e o(t) 相乘器低通滤波器抽样判决器2DPSK 带通滤波器 延迟T
（1）所用滤波器函数：巴特沃斯滤波器 % 注: wp（或Wp）为通带截止频率 ws(或Ws)为阻带截止频率 Rp为通带衰减 As为阻带衰减 %butterworth低通滤波器原型设计函数要求Ws>Wp>0 As>Rp>0 function [b,a]=afd_butt(Wp,Ws,Rp,As) N=ceil((log10((10^(Rp/10)-1)/(10^(As/10)-1)))/(2*log10(Wp/Ws))); %上条语句为求滤波器阶数 N为整数 %ceil 朝正无穷大方向取整 fprintf('\n Butterworth Filter Order=%2.0f\n',N) OmegaC=Wp/((10^(Rp/10)-1)^(1/(2*N))) %求对应于N的3db截止频率 [b,a]=u_buttap(N,OmegaC); （2）傅里叶变换函数 function [Xk]=dft(xn,N) n=[0:1:N-1]; k=[0:1:N-1]; WN=exp(-j*2*pi/N); nk=n'*k; WNnk=WN.^(nk); Xk=xn*WNnk; 设计部分： 1.普通AM调制与解调 %单音普通调幅波调制y=amod(x,t,fs,t0,fc,Vm0,ma)要求fs>2fc %x调制信号，t调制信号自变量,t0采样区间，fs采样频率， %fc载波频率，Vm0输出载波电压振幅，ma调幅度 t0=0.1;fs=12000; fc=1000;Vm0=2.5;ma=0.25; n=-t0/2:1/fs:t0/2; x=4*cos(150*pi*n); %调制信号 y2=Vm0*cos(2*pi*fc*n); %载波信号figure(1) subplot(2,1,1);plot(n,y2); axis([-0.01,0.01,-5,5]); title('载波信号'); N=length(x); Y2=fft(y2); subplot(2,1,2); plot(n,Y2); title('载波信号频谱'); %画出频谱波形y=Vm0*(1+ma*x/Vm0).*cos(2*pi*fc*n); figure(2) subplot(2,1,1);plot(n,x) title('调制信号'); subplot(2,1,2) plot(n,y) title('已调波信号'); X=fft(x);Y=fft(y);