基带信号眼图实验matlab仿真
数字基带信号的眼图实验——m a t l a b
仿真
一、实验目的
1、掌握无码间干扰传输的基本条件和原理，掌握基带升余弦滚降系统的实现方法；
2、通过观察眼图来分析码间干扰对系统性能的影响，并观察在输入相同码率的NRZ 基带信号下，不同滤波器带宽对输出信号码间干扰大小的影响程度；
3、熟悉MATLAB 语言编程。
二、实验预习要求
1、复习《数字通信原理》第七章节——奈奎斯特第一准则内容；
2、复习《数字通信原理》第七章节——数字基带信号码型内容；
3、认真阅读本实验内容，熟悉实验步骤。
三、实验原理和电路说明
1、基带传输特性
基带系统的分析模型如图3-1所示，要获得良好的基带传输系统，就应该
图3-1 基带系统的分析模型
抑制码间干扰。设输入的基带信号为
()n s n a t nT δ-∑，s T 为基带信号的码元周期，则经过基带传输系统后的输出码元为()n s
n
a h t nT -∑。其中 1()()2j t h t H e d ωωωπ+∞-∞=? （3-1）
理论上要达到无码间干扰，依照奈奎斯特第一准则，基带传输系统在时域应满足：
10()0,s k h kT k =?=??
，为其他整数 (3-2)
频域应满足： ()0,s s T T H πωωω?≤?=???
，其他 (3-3)
图3-2 理想基带传输特性 此时频带利用率为2/Baud Hz ,这是在抽样值无失真条件下，所能达到的最高频率利用率。
由于理想的低通滤波器不容易实现，而且时域波形的拖尾衰减太慢，因此在得不到严格定时时，码间干扰就可能较大。在一般情况下，只要满足：
222(),s i s s s s i H H H H T T T T T ππππ
ωωωωω??????+=-+++=≤ ? ? ???????∑ (3-4)
基带信号就可实现无码间干扰传输。这种滤波器克服了拖尾太慢的问题。
从实际的滤波器的实现来考虑，采用具有升余弦频谱特性()H ω时是适宜的。
(1)(1)1sin (),2(1)()1,0(1)0,s s s s s s T T T T H T T ππαπαωωαπαωωπαω???-+--≤≤???????-?=≤≤???+>???
(3-5)
这里α称为滚降系数，01α≤≤。
所对应的其冲激响应为： ()222sin
cos()()14s s s s
t T t T h t t t T T παππα=- (3-6) 此时频带利用率降为2/(1)Baud/Hz α+，这同样是在抽样值无失真条件下，所能达到的最高频率利
用率。换言之，若输入码元速率'1/s s R T >，则该基带传输系统输出码元会产生码间干扰。
2、眼图
所谓眼图就是将接收滤波器输出的，未经再生的信号，用位定时以及倍数作为同步信号在示波器上重复扫描所显示的波形（因传输二进制信号时，类似人的眼睛）。干扰和失真所产生的畸变可以很清楚的从眼图中看出。眼图反映了系统的最佳抽样时间，定时的灵敏度，噪音容限，信号幅度的畸变范围以及判决门限电平，因此通常用眼图来观察基带传输系统的好坏。
图3-3 眼图示意图
四、仿真环境
Windows NT/2000/XP/Windows 7/VISTA ；
MATLAB 以上。
五、仿真程序设计
1、程序框架
图3-4 程序框架
首先，产生M进制双极性NRZ码元序列，并根据系统设置的抽样频率对该NRZ码元序列进行
抽样，再将抽样序列送到升余弦滚降系统，最后画出输出码元序列眼图。
2、参数设置
该仿真程序应具备一定的通用性，即要求能调整相应参数以仿真不同的基带传输系统，并观察输出眼图情况。因此，对于NRZ码元进制M、码元序列长度Num、码元速率Rs，采样频率Fs、升余弦滚降滤波器参考码元周期Ts、滚降系数alpha、在同一个图像窗口内希望观测到的眼图个数Eye_num
等均应可以进行合理设置。
3、实验内容
根据现场实验题目内容，设置仿真程序参数，编写仿真程序，仿真波形，并进行分析给出结论。
4、仿真结果参考
参考例程参数设置如下：
无码间干扰时：
Ts=1e-2; %升余弦滚降滤波器的理想参考码元周期，单位s
Fs=1e3; %采样频率，单位Hz。注意：该数值过大将
%严重增加程序运行时间
Rs=50; %输入码元速率，单位Baud
M=2; %输入码元进制
Num=100; %输入码元序列长度。注意：该数值过大将
%严重增加程序运行时间
Eye_num=2; %在一个窗口内可观测到的眼图个数。
图3-5(a) 仿真参考结果图（1）
图3-5(b) 仿真参考结果图（2）
图3-5(c) 仿真参考结果图（3）
从眼图张开程度可以得出没有发生码间干扰，这是因为基带信号的码元速率Rs为50Baud，而升余弦滚降滤波器和FIR滤波器的等效带宽B=60Hz（Ts=10ms），Rs<2B，满足了奈奎斯特第一准则的条件。
有码间干扰时：
Ts=5*(1e-2); %升余弦滚降滤波器的参考码元周期，单位s
Fs=1e3; %采样频率，单位Hz。注意：该数值过大将
%严重增加程序运行时间
Rs=50; %输入码元速率，单位Baud
M=2; %输入码元进制
Num=100; %输入码元序列长度。注意：该数值过大将
%严重增加程序运行时间
Eye_num=2; %在一个窗口内可观测到的眼图个数。
图3-5(d) 仿真参考结果图（4）
眼图基本闭合，存在较为严重的码间干扰，这是因为码元速率Rs虽然仍为50Baud，但滤波器等效带宽已经变为12Hz（Ts=50ms），Rs>2B不再满足奈奎斯特第一准则。
多进制码元情况：
图3-6 四进制NRZ码元眼图
六、实验报告要求
1、整理实验数据，画出相应的波形。
2、结合奈奎斯特第一准则，分析波形，表述出码间干扰ISI现象与滤波器的等效带宽设定值之间的关系，给出原因。
3、结合奈奎斯特第一准则，分析波形，表述出码间干扰ISI现象与滤波器的滚降系数设定值之间的关系，给出原因。
七、思考题
1、自行编写升余弦滚降滤波器冲激响应函数，特别注意当公式中分子分母均为0时的特殊情况。
2、采用MATLAB自带眼图函数eyediagram来观察眼图。
八、参考程序
close all;
alpha=; %设置滚降系数，取值范围在[0,1]
Ts=1e-2; %升余弦滚降滤波器的参考码元周
%期, Ts=10ms,无ISI。
% Ts=2*(1e-2); %Ts=20ms,已经出现ISI（临界点）
% Ts=5*(1e-2); %Ts=50ms,出现严重ISI
Fs=1e3; %采样频率，单位Hz。注意：该数
%值过大将严重增加程序运行时间
Rs=50; %输入码元速率，单位Baud
% M=2;
M=4; %输入码元进制
Num=100; %输入码元序列长度。注意：该数值
%过大将严重增加程序运行时间。Samp_rate=Fs/Rs %采样率，应为大于1的正整数，即
%要求Fs,Rs之间呈整数倍关系
% Eye_num=2; %在一个窗口内可观测到的眼图个数。Eye_num=4; %在一个窗口内可观测到的眼图个数。%产生双极性NRZ码元序列
NRZ=2*randint(1,Num,M)-M+1;
figure(1);
stem(NRZ);
xlabel('时间');
ylabel('幅度');
hold on;
grid on;
title('双极性NRZ码元序列');
%对双极性NRZ码元序列进行抽样
k=1;
for ii=1:Num
for jj=1:Samp_rate
Samp_data(k)=NRZ(ii);
k=k+1;
end
end
%基带升余弦滚降系统冲激响应
[ht,a] = rcosine(1/Ts,Fs,'fir',alpha);
%画出基带升余弦滚降系统冲激响应波形
figure(2);
subplot(2,1,1);
plot(ht);
xlabel('时间');
ylabel('冲激响应');
hold on;
grid on;
title('升余弦滚降系统冲激响应,滚降因子\alpha=');
%将信号送入基带升余弦滚降系统，即做卷积操作
st = conv(Samp_data,ht)/(Fs*Ts);
subplot(2,1,2);
plot(st);
xlabel('时间');
ylabel('信号幅度');
hold on;
grid on;
title('经过升弦滚降系统后的码元')
%画眼图，在同一个图形窗口重复画出一个或若干个码元figure(3);
for k = 10:floor(length(st)/Samp_rate)-10
%不考虑过渡阶段信号，只观测稳定阶段
ss = st(k*Samp_rate+1:(k+Eye_num)*Samp_rate);
plot(ss);
hold on;
end
xlabel('时间');
ylabel('信号幅度');
hold on;
grid on;
title('基带信号眼图');
% eyediagram(st,Samp_rate);
% xlabel('时间');
% ylabel('信号幅度');
% hold on;
% grid on;
% title('基带信号眼图');
信号与系统matlab实验及答案
产生离散衰减正弦序列()π0.8sin 4n x n n ?? = ??? , 010n ≤≤，并画出其波形图。 n=0:10; x=sin(pi/4*n).*0.8.^n; stem(n,x);xlabel( 'n' );ylabel( 'x(n)' ); 用MATLAB 生成信号()0sinc at t -, a 和0t 都是实数,410t -<<，画波形图。观察并分析a 和0t 的变化对波形的影响。 t=linspace(-4,7); a=1;
t0=2; y=sinc(a*t-t0); plot(t,y); t=linspace(-4,7); a=2; t0=2; y=sinc(a*t-t0); plot(t,y);
t=linspace(-4,7); a=1; t0=2; y=sinc(a*t-t0); plot(t,y);
三组对比可得a 越大最大值越小，t0越大图像对称轴越往右移 某频率为f 的正弦波可表示为()()cos 2πa x t ft =，对其进行等间隔抽样，得到的离散样值序列可表示为()()a t nT x n x t ==，其中T 称为抽样间隔，代表相邻样值间的时间间隔，1 s f T = 表示抽样频率，即单位时间内抽取样值的个数。抽样频率取40 Hz s f =，信号频率f 分别取5Hz, 10Hz, 20Hz 和30Hz 。请在同一张图中同时画出连续信号()a x t t 和序列()x n nT 的波形图，并观察和对比分析样值序列的变化。可能用到的函数为plot, stem, hold on 。 fs = 40; t = 0 : 1/fs : 1 ; % ?μ?ê·?±e?a5Hz,10Hz,20Hz,30Hz f1=5; xa = cos(2*pi*f1*t) ; subplot(1, 2, 1) ;
基于MATLAB的GMSK调制与解调课设报告
基于Matlab的GMSK调制与解调 1.课程设计目的 （1）加深对GMSK基本理论知识的理解。 （2）培养独立开展科研的能力和编程能力。 （3）通过SIMULINK对BT=0.3的GMSK调制系统进行仿真。 2.课程设计要求 （1）观察基带信号和解调信号波形。 （2）观察已调信号频谱图。 （3）分析调制性能和BT参数的关系。 3.相关知识 3.1GMSK调制 调制原理图如图2.2，图中滤波器是高斯低通滤波器，它的输出直接对VCO 进行调制，以保持已调包络恒定和相位连续。 非归零数字序 高斯低通滤 波器频率调制器 （VCO） GMSK已 调信号 图3.1GMSK调制原理图 为了使输出频谱密集，前段滤波器必须具有以下待性： 1.窄带和尖锐的截止特性，以抑制FM调制器输入信号中的高频分量； 2.脉冲响应过冲量小，以防止FM调制器瞬时频偏过大； 3.保持滤波器输出脉冲响应曲线下的面积对应丁pi／2的相移。以使调制指数为1／2。前置滤波器以高斯型最能满足上述条件，这也是高斯滤波器最小移频键控(GMSK)的由来。
GMSK 信号数据 3.2GMSK 解调 GMSK 本是MSK 的一种，而MSK 又是是FSK 的一种，因此，GMSK 检波也可以采用FSK 检波器，即包络检波及同步检波。而GMSK 还可以采用时延检波，但每种检波器的误码率不同。 GMSK 非相干解调原理图如图2.3，图中是采用FM 鉴频器（斜率鉴频器或相位鉴频器）再加判别电路，实现GMSK 数据的解调输出。 图3.2GMSK 解调原理图 4.课程设计分析 4.1信号发生模块 因为GMSK 信号只需满足非归零数字信号即可，本设计中选用（Bernoulli Binary Generator）来产生一个二进制序列作为输入信号。 图4.1GMSK 信号产生器 该模块的参数设计这只主要包括以下几个。其中probability of a zero 设置为0.5表示产生的二进制序列中0出现的概率为0.5；Initial seed 为61表示随机数种子为61；sample time 为1/1000表示抽样时间即每个符号的持续时为0.001s。当仿真时间固定时，可以通过改变sample time 参数来改变码元个数。例如仿真时间为10s，若sample time 为1/1000，则码元个数为10000。 带通滤 波器限幅器判决器鉴频器GMSK 信号 输出
matlab信号与系统实验报告
实验一 基本信号的产生与运算 一、 实验目的 学习使用MATLAB 产生基本信号、绘制信号波形、实现信号的基本运算。 二、 实验原理 MATLAB 提供了许多函数用于产生常用的基本信号：如阶跃信号、脉冲信号、指数信号、正弦信号和周期方波等等。这些信号是信号处理的基础。 1、 利用MATLAB 产生下列连续信号并作图。 （1）51),1(2)(<<---=t t u t x （2）300),3 2sin()(3.0<<=-t t e t x t （3）1.01.0,3000cos 100cos )(<<-+=t t t t x （4）2000),8.0cos()1.0cos()(<<=t t t t x ππ 答：（1）、 >> t=-1:0.02:5; >> x=(t>1); >> plot(t,-2*x); >> axis([-1,5,-3,1]); >> title('杨婕婕 朱艺星'); >> xlabel('x(t)=-2u(t-1)'); （2）、 >> t=0:0.02:30; >> x=exp(-0.3*t).*sin(2/3*t); >> plot(t,x);
>> title('杨婕婕朱艺星'); >> xlabel('x(t)=exp(-0.3*t).*sin(2/3*t)'); 因为原函数在t=15后x(t)取值接近于零，所以将横坐标改成0到15，看得更清晰 axis([0,15,-0.2,0.6]); （3）>> t=-0.1:0.01:0.1;x=cos(100*t)+cos(3000*t);plot(t,x); >> title('杨婕婕朱艺星'); >>xlabel('x=cos(100*t)+cos(3000*t)'); 因为t的间隔取太大，以至于函数不够准确，缩小t的间隔： t=-0.1:0.002:0.2;x=cos(100*t)+cos(3000*t); plot(t,x);title('杨婕婕') >> t=-0.1:0.0001:0.1; x=cos(100*t)+cos(3000*t); >> plot(t,x);title('杨婕婕朱艺星'); >> xlabel('x=cos(100*t)+cos(3000*t)');
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)实验指导书
《信号与系统》实验指导书 张建奇骆崇编写 浙江工业大学之江学院信息工程分院 2012年2月
目录 实验一MATLAB的基本使用 (1) 一、实验目的 (1) 二、实验原理 (1) 三、实验内容与要求 (8) 四、实验报告 (9) 实验二时域波形的MATLAB实现 (10) 一、实验目的 (10) 二、预习要求 (10) 三、实验原理 (10) 四、实验内容与要求 (18) 五、实验报告 (19) 实验三用MATLAB对系统时域分析 (20) 一、实验目的 (20) 二、预习要求 (20) 三、实验原理 (20) 四、实验内容与要求 (29)
实验一MATLAB的基本使用 一、实验目的 1、了解和掌握MATLAB的基本操作 2、了解MATLAB的库函数 3、会用MATLAB进行简单的操作。 二、实验原理 1、界面操作 MATLAB是“MATrix LABoratory”的缩写（矩阵实验室），它是由美国Mathworks公司于1984年正式推出的一种科学计算软件，由于其强大的功能，在欧美的一些大学里MATLAB已经成为许多诸如数字信号处理、自动控制理论等高级教程的主要工具软件，同时也成为理工科学生，必须掌握的一项基本技能。 当需要运行程序时，只需选择桌面上（或开始）中的MATLAB6.5应用程序图标即可 通常情况下，MATLAB的工作环境主要由一下几个窗口组成： 命令窗口（Command Window）
工作区间浏览器（Workspace） 历史命令窗口（Command History） 图形窗口(Figure) 文本编辑窗口(Editor) 当前路径窗口（Current Directory） MATLAB的命令窗与命令操作 当用户使用命令窗口进行工作时，在命令窗口中可以直接输入相应的命令，系统将自动显示信息。 例如在命令输入提示符“>>”后输入指令： >>t=[1,2,3;4,5,6;7,8,9]; 按回车键(Enter)后，系统即可完成对变量t的赋值。 MATALB提供了非常方便的在线帮助命令(help)，它可提供各个函数的用法指南，包括格式、参数说明、注意事项及相关函数等内容。 2、图形窗 MATLAB图形窗(Figure)主要用于显示用户所绘制的图形。 通常，只要执行了任意一种绘图命令，图形窗就会自动产生。
信号与系统MATLAB实验报告
《信号与系统》MATLAB实验报告 院系：专业： 年级：班号： 姓名：学号： 实验时间： 实验地点：
实验一 连续时间信号的表示及可视化 实验题目： )()(t t f δ=；)()(t t f ε=；at e t f =)(（分别取00<>a a 及）； )()(t R t f =；)()(t Sa t f ω=；)2()(ft Sin t f π=（分别画出不同周期个数 的波形）。 解题分析： 以上各类连续函数，先运用t = t1: p:t2的命令定义时间范围向量，然后调用对应的函数，建立f 与t 的关系，最后调用plot （）函数绘制图像，并用axis （）函数限制其坐标范围。 实验程序： （1） )()(t t f δ= t=-1:0.01:3 %设定时间变量t 的范围及步长 f=dirac(t) %调用冲激函数dirac （） plot(t,f) %用plot 函数绘制连续函数 axis([-1,3,-0.5,1.5]) %用axis 函数规定横纵坐标的范围 （2） )()(t t f ε= t=-1:0.01:3 %设定时间变量t 的范围及步长 f=heaviside(t) %调用阶跃函数heaviside （） plot(t,f) %用plot 函数绘制连续函数 title('f(t)=heaviside(t)') %用title 函数设置图形的名称 axis([-1,3,-0.5,1.5]) %用axis 函数规定横纵坐标的范围 （3） at e t f =)( a=1时： t=-5:0.01:5 %设定时间变量t 的范围及步长 f=exp(t) %调用指数函数exp （）
基于matlab的通信信道及眼图的仿真 通信原理课程设计
通信原理课程设计 基于matlab的通信信道及眼图的仿真 作者: 摘要 由于多径效应和移动台运动等影响因素，使得移动信道对传输信号在时间、频率和角度上造成了色散，即时间色散、频率色散、角度色散等等，因此多径信道的特性对通信质量有着重要的影响，而多径信道的包络统计特性则是我们研究的焦点。根据不同无线环境，接收信号包络一般服从几种典型分布，如瑞利分布、莱斯分布等。因此我们对瑞利信道、莱斯信道进行了仿真并针对服从瑞利分布的多径信道进行模拟仿真。由于眼图是实验室中常用的一种评价基带传输系统的一种定性而方便的方法，“眼睛”的张开程度可以作为基带传输系统性能的一种度量，它不但反映串扰的大小，而且也可以反映信道噪声的影响。为此，我们在matlab上进行了仿真,加深对眼图的理解。 关键词：瑞利信道莱斯信道多径效应眼图 一、瑞利信道 在移动通信系统中，发射端和接收端都可能处于不停的运动状态之中，这种相对运动将产生多普勒频移。在多径信道中，发射端发出的信号通过多条路径到达接收端，这些路径具有不同的延迟和接收强度，它们之间的相互作用就形成了衰落。MATLAB中的多径瑞利衰落信道模块可以用于上述条件下的信道仿真。 多径瑞利衰落信道模块用于多径瑞利衰落信道的基带仿真，该模块的输入信号为复信号，可以为离散信号或基于帧结构的列向量信号。无线系统中接收机与发射机之间的相对运动将引起信号频率的多普勒频移，多普勒频移值由下式决定： 其中v是发射端与接收端的相对速度，θ是相对速度与二者连线的夹角，λ是信号的波长。
Fd的值可以在该模块的多普勒平移项中设置。由于多径信道反映了信号在多条路径中的传输，传输的信号经过不同的路径到达接收端，因此产生了不同的时间延迟。当信号沿着不同路径传输并相互干扰时，就会产生多径衰落现象。在模块的参数设置表中，Delay vector（延迟向量）项中，可以为每条传输路径设置不同的延迟。如果激活模块中的Normalize gain vector to 0 dB overall gain,则表示将所有路径接收信号之和定为0分贝。信号通过的路径的数量和Delay vector（延迟向量）或Gain vector(增益向量)的长度对应。Sample time(采样时间)项为采样周期。离散的Initial seed(初始化种子)参数用于设置随机数的产生。 1.1、Multipath Rayleigh Fading Channel(多径瑞利衰落信道)模块的主要参数 参数名称参数值 Doppler frequency(Hz) 40/60/80 Sample time 1e-6 Delay vector(s) [0 1e-6] Gain vector(dB) [0 -6] Initial seed 12345 使能 Normalize gain vector to 0 dB overall gain Bernoulli Random Binary Generator(伯努利二进制随机数产生器)的主要参数 参数名称参数值 Probability of a zero0.5 Initial seed54321
实验1 用MATLAB进行信号频谱分析(推荐文档)
实验1 用MATLAB 进行信号频谱分析 一、实验目的 ㈠ 初步掌握MATLAB 产生常用离散时间信号的编程方法。 ㈡ 学习编写简单的FFT 算法程序，对离散信号进行幅频谱分析。 ㈢ 观察离散时间信号频谱的特点。 二、实验原理 ㈠ 常用的离散时间信号 在 MATLAB 语言主要是研究离散信号的。常用的离散信号有： 1．单位取样序列 ???≠==000 1)(n n n δ 2．单位阶跃序列 ?? ?<≥=0 01 )(n n n u 3．实指数序列 R a n a n x n ∈?=;)( 4．复指数序列 n e n x n j ?=+)(0)(ωσ 5．正(余)弦序列 )c o s ()(0θω+=n n x n ? 6．周期序列 n N n x n x ?+=)()( ㈡ 离散信号的产生 离散信号的图形显示使用stem 指令。 在 MATLAB 中的信号处理工具箱中,主要提供的信号是离散信号。
由于MATLAB 对下标的约定为从1开始递增，例如x=[5,4,3,2,1,0],表示x(1)=5,x(2)=4,X(3)=3… 因此要表示一个下标不由1开始的数组x(n)，一般应采用两个矢量，如 n=[-3，-2，-1，0，l ，2，3，4，5]； x=[1，-l ，3，2，0，4，5，2，1]； 这表示了一个含9个采样点的矢量：X(n)={x(-3)，x(-2)，x(-1)，x(0)，x(1)，x(2)，x(3)，x(4)，x(5)}。 1．单位取样序列 ?? ?≠==δ0 001)(n n n 这一函数实现的方法有二： 方法一：可利用MATLAB 的zeros 函数。 x=zeros(1，N)； %建立一个一行N 列的全零数组 x(1)=1； %对X （1）赋1 方法二：可借助于关系操作符实现 n=1:N; x=[n==1]; %n 等于1时逻辑关系式结果为真，x=1；n 不等于1时为假，x=0 如要产生 ?????≤<<=≤≤=-δ2 0210 100)(10)(n n n n n n n n n n n n 则可采用MATLAB 实现： n=n1:n2； x=[(n-n0)==0]；%n=n0时逻辑关系式结果为真，x=1；n ≠n0时为
(完整word版)信号与系统matlab实验
习题三 绘制典型信号及其频谱图 1.更改参数，调试程序，绘制单边指数信号的波形图和频谱图。观察参数a对信号波形 及其频谱的影响。 程序代码： close all; E=1;a=1; t=0:0.01:4; w=-30:0.01:30; f=E*exp(-a*t); F=1./(a+j*w); plot(t,f);xlabel('t');ylabel('f(t)'); figure; plot(w,abs(F));xlabel('\omega');ylabel('|F(\omega)|';
E=1，a=1,波形图频谱图更改参数E=2，a=1;
更改参数a，对信号波形及其频谱的影响。（保持E=2）上图为a=1图像 a=2时
a=4时 随着a的增大，f(t)曲线变得越来越陡，更快的逼近0，而对于频谱图，随着a增大，图像渐渐向两边张开，峰值减小，陡度减小，图像整体变得更加平缓。 2.矩形脉冲信号 程序代码： close all; E=1;tao=1; t=-4:0.1:4; w=-30:0.1:30;
f=E*(t>-tao/2&tao/2)+0*(t<=-tao/2&t>=tao/2); F=(2*E./w).*sin(w*tao/2); plot(t,f);xlabel('t');ylabel('f(t)'); figure; plot(w,abs(F));xlabel('\omega');ylabel('|F(\omega)|') ; figure; plot(w,20*log10(abs(F))); xlabel('\omega');ylabel('|F(\omega)| in dB'); figure; plot(w,angle(F));xlabel('\omega');ylabel('\phi(\omega )');
MATLAB通信系统仿真实验报告1
MATLAB通信系统仿真实验报告
实验一、MATLAB的基本使用与数学运算 目的：学习MATLAB的基本操作，实现简单的数学运算程序。 内容： 1-1要求在闭区间[0，2π]上产生具有10个等间距采样点的一维数组。试用两种不同的指令实现。 运行代码：x=[0:2*pi/9:2*pi] 运行结果： 1-2用M文件建立大矩阵x x=[0.10.20.30.40.50.60.70.80.9 1.11.21.31.41.51.61.71.81.9 2.12.22.32.42.52.62.72.82.9 3.13.23.33.43.53.63.73.83.9] 代码：x=[0.10.20.30.40.50.60.70.80.9 1.11.21.31.41.51.61.71.81.9 2.12.22.32.42.52.62.72.82.9 3.13.23.33.43.53.63.73.83.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=[56;78]B=[910;1112]x1=A+B X2=A-B X3=A*B X4=A.*B X5=A^3 X6=A.^3X7=A/B X8=A\B
运行结果： 1-4任意建立矩阵A，然后找出在[10,20]区间的元素位置。 程序代码及运行结果： 代码：A=[1252221417;111024030;552315865]c=A>=10&A<=20运行结果： 1-5总结：实验过程中，因为对软件太过生疏遇到了些许困难，不过最后通过查书与同学交流都解决了。例如第二题中，将文件保存在了D盘，而导致频频出错，最后发现必须保存在MATLAB文件之下才可以。第四题中，逻辑语言运用到了ij，也出现问题，虽然自己纠正了问题，却也不明白错在哪了，在老师的讲解下知道位置定位上不能用ij而应该用具体的整数。总之第一节实验收获颇多。
matlab信号与系统实验报告
实验一 基本信号的产生与运算 一、 实验目的 学习使用MATLAB 产生基本信号、绘制信号波形、实现信号的基本运算。 二、 实验原理 MATLAB 提供了许多函数用于产生常用的基本信号：如阶跃信号、脉冲信号、指数信号、正弦信号和周期方波等等。这些信号是信号处理的基础。 1、 利用MATLAB 产生下列连续信号并作图。 （1）51),1(2)(<<---=t t u t x （2）300),3 2 sin()(3.0<<=-t t e t x t （3）1.01.0,3000cos 100cos )(<<-+=t t t t x （4）2000),8.0cos()1.0cos()(<<=t t t t x ππ 答：（1）、 >> t=-1:0.02:5; >> x=(t>1); >> plot(t,-2*x); >> axis([-1,5,-3,1]); >> title('杨婕婕 朱艺星'); >> xlabel('x(t)=-2u(t-1)');
（2）、 >> t=0:0.02:30; >> x=exp(-0.3*t).*sin(2/3*t); >> plot(t,x); >> title('杨婕婕朱艺星'); >> xlabel('x(t)=exp(-0.3*t).*sin(2/3*t)');
因为原函数在t=15后x(t)取值接近于零，所以将横坐标改成0到15，看得更清晰 axis([0,15,-0.2,0.6]);
（3）>> t=-0.1:0.01:0.1;x=cos(100*t)+cos(3000*t);plot(t,x); >> title('杨婕婕朱艺星'); >>xlabel('x=cos(100*t)+cos(3000*t)'); 因为t的间隔取太大，以至于函数不够准确，缩小t的间隔：t=-0.1:0.002:0.2;x=cos(100*t)+cos(3000*t); plot(t,x);title('杨婕婕')
信号与系统MATLAB实验
2016-2017学年第一学期 信号与系统实验报告 班级： 姓名： 学号： 成绩： 指导教师：
实验一常见信号的MATLAB 表示及运算 一．实验目的 1．熟悉常见信号的意义、特性及波形 2．学会使用MATLAB 表示信号的方法并绘制信号波形 3. 掌握使用MATLAB 进行信号基本运算的指令 4. 熟悉用MATLAB 实现卷积积分的方法 二．实验原理 信号一般是随时间而变化的某些物理量。按照自变量的取值是否连续，信号分为连续时间信号和离散时间信号，一般用()f t 和()f k 来表示。若对信号进行时域分析，就需要绘制其波形，如果信号比较复杂，则手工绘制波形就变得很困难，且难以精确。MATLAB 强大的图形处理功能及符号运算功能，为实现信号的可视化及其时域分析提供了强有力的工具。 根据MATLAB 的数值计算功能和符号运算功能，在MATLAB 中，信号有两种表示方法，一种是用向量来表示，另一种则是用符号运算的方法。在采用适当的MATLAB 语句表示出信号后，就可以利用MATLAB 中的绘图命令绘制出直观的信号波形了。下面分别介绍连续时间信号和离散时间信号的MATLAB 表示及其波形绘制方法。 1.连续时间信号 所谓连续时间信号，是指其自变量的取值是连续的，并且除了若干不连续的点外，对于一切自变量的取值，信号都有确定的值与之对应。从严格意义上讲，MATLAB 并不能处理连续信号。在MATLAB 中，是用连续信号在等时间间隔点上的样值来近似表示的，当取样时间间隔足够小时，这些离散的样值就能较好地近似出连续信号。在MATLAB 中连续信号可用向量或符号运算功能来表示。 ⑴ 向量表示法 对于连续时间信号()f t ，可以用两个行向量f 和t 来表示，其中向量t 是用形如12::t t p t 的命令定义的时间范围向量，其中，1t 为信号起始时间，2t 为终止时间，p 为时间间隔。向量f 为连续信号()f t 在向量t 所定义的时间点上的样值。 说明：plot 是常用的绘制连续信号波形的函数。 严格说来，MATLAB 不能表示连续信号，所以，在用plot()命令绘制波形时，要对自变量t 进行取值，MATLAB 会分别计算对应点上的函数值，然后将各个数据点通过折线连接起来绘制图形，从而形成连续的曲线。因此，绘制的只是近似波形，而且，其精度取决于t 的取样间隔。t 的取样间隔越小，即点与点之间的距离越小，则近似程度越好，曲线越光滑。例如：图1-1是在取样间隔为p=0.5时绘制的波形，而图1-2是在取样间隔p=0.1时绘制的波形，两相对照，可以看出图1-2要比图1-1光滑得多。
通信原理数字基带传输系统习题及其答案
第四章（数字基带传输系统）习题及其答案 【题4－1】设二进制符号序列为，试以矩形脉冲为例，分别画出相应的单极性码型，双极性码波形，单极性归零码波形，双极性归零码波形，二进制差分码波形。 【答案4－1】 【题4－2】设随机二机制序列中的0和1分别由()g t 和()g t -组成，其出现概率分别为p 和(1)p -： 1）求其功率谱密度及功率； 2）若()g t 为图（a ）所示的波形，s T 为码元宽度，问该序列存在离散分量 1 s f T =否 3）若()g t 改为图（b ）所示的波形，问该序列存在离散分量1 s f T =否 【答案4－2】 1）随机二进制序列的双边功率谱密度为 由于 可得： 式中：()G f 是()g t 的频谱函数。在功率谱密度()s P ω中，第一部分是其连续谱成分，第二部分是其离散谱成分。 随机二进制序列的功率为 2）当基带脉冲波形()g t 为 ()g t 的付式变换()G f 为 因此 式中： 1 s s f T = 。 所以，该二进制序列不存在离散分量。 3）当基带脉冲波形()g t 为
()g t 的付式变换()G f 为 因此 式中： 1s s f T = 。 所以，该二进制序列存在离散分量。 【题4－3】设二进制数字基带信号的基本脉冲序列为三角形脉冲，如下图所示。图中s T 为码元宽度，数字信息1和0分别用()g t 的有无表示，且1和0出现的概率相等： 1）求数字基带信号的功率谱密度； 2）能否重该数字基带信号中提取同步所需的频率1 s s f T =的分量若能，计 算该分量的功率。 【答案4－3】 1）由图得 ()g t 的频谱函数()G ω为 由题设可知 所以 代入二进制数字基带信号的双边功率谱密度函数表达式，可得 2）二进制数字基带信号的离散谱分量()v P ω为 当1m =±时，s f f =±，代入上式可得 因为该二进制数字基带信号中存在1s s f T =的离散分量，所以能从该数字基带信号中提取码元同步所需的频率1s s f T =的分量。 该频率分量的功率为 【题4－5】已知信息代码为，求相应的AMI 码、HDB3码、PST 码及双相码。 【答案4－5】 AMI 码： +10000 0000 –1＋1 HDB3码； +1000+V-B00-V0+1-1
信号与系统MATLAB实验总汇
实验一、MATLAB 编程基础及典型实例 一、实验目的 （1）熟悉MATLAB 软件平台的使用； （2）熟悉MATLAB 编程方法及常用语句； （3）掌握MATLAB 的可视化绘图技术； （4）结合《信号与系统》的特点，编程实现常用信号及其运算。 示例一：在两个信号进行加、减、相乘运算时，参于运算的两个向量要有相同的维数，并且它们的时间变量范围要相同，即要对齐。编制一个函数型m 文件，实现这个功能。function [f1_new,f2_new,n]=duiqi(f1,n1,f2,n2) a=min(min(n1),min(n2)); b=max(max(n1),max(n2)); n=a:b; f1_new=zeros(1,length(n)); f2_new=zeros(1,length(n)); tem1=find((n>=min(n1))&(n<=max(n1))==1); f1_new(tem1)=f1; tem2=find((n>=min(n2))&(n<=max(n2))==1); f2_new(tem2)=f2; 四、实验内容与步骤 （2）绘制信号x(t)=)3 2sin(2t e t ?的曲线，t 的范围在0~30s ，取样时间间隔为0.1s 。t=0:0.1:30; y=exp(-sqrt(2)*t).*sin(2*t/3); plot(t,y);
（3）在n=[-10:10]范围产生离散序列：?? ?≤≤?=Other n n n x ,033,2)(,并绘图。n=-10:1:10; z1=((n+3)>=0); z2=((n-3)>=0); x=2*n.*(z1-z2); stem(n,x);（4）编程实现如下图所示的波形。 t=-2:0.001:3; f1=((t>=-1)&(t<=1)); f2=((t>=-1)&(t<=2)); f=f1+f2; plot(t,f); axis([-2,3,0,3]);
(完整word版)使用matlab绘制眼图.docx
使用 matlab 绘制数字基带信号的眼图实验 一、实验目的 1、掌握无码间干扰传输的基本条件和原理，掌握基带升余弦滚降系统的实现方法； 2、通过观察眼图来分析码间干扰对系统性能的影响，并观察在输入相同码率的NRZ 基带信号下，不同滤波器带宽对输出信号码间干扰大小的影响程度； 3、熟悉 MATLAB语言编程。 二、实验原理和电路说明 1、基带传输特性 基带系统的分析模型如图3-1 所示，要获得良好的基带传输系统，就应该 a n t nT s 基带传输a n h t nT s n n抽样判决 H ( ) 图 3-1基带系统的分析模型 抑制码间干扰。设输入的基带信号为a n t nT s， T s为基带信号的码元周期，则经过 n 基带传输系统后的输出码元为a n h t nT s。其中 n h(t )1H ()e j t d（3-1 ） 2 理论上要达到无码间干扰，依照奈奎斯特第一准则，基带传输系统在时域应满足： ，k 0 h( kT s)(3-2) 0,k为其他整数 频域应满足： T s， T s(3-3) H ( ) 0,其他
H ( ) T s T s T s 图 3-2 理想基带传输特性 此时频带利用率为 2Baud / Hz , 这是在抽样值无失真条件下，所能达到的最高频率利用率。 由于理想的低通滤波器不容易实现， 而且时域波形的拖尾衰减太慢， 因此在得不到严格 定时时，码间干扰就可能较大。在一般情况下，只要满足： 2 i H 2 2 , (3-4) H H ( ) H T s i T s T s T s T s 基带信号就可实现无码间干扰传输。这种滤波器克服了拖尾太慢的问题。 从实际的滤波器的实现来考虑，采用具有升余弦频谱特性 H ( ) 时是适宜的。 1 sin T s ( ) , (1 ) (1 ) 2 T s T s T s H ( ) 1, (1 ) 0 (3-5) T s 0, (1 ) T s 这里 称为滚降系数， 1。 所对应的其冲激响应为： sin t cos( t T s ) h(t ) T s (3-6) t 1 4 2t 2 T s 2 T s 此时频带利用率降为 2 / (1 ) Baud/ Hz ，这同样是在抽样值无失真条件下， 所能达到的最 高频率利用率。换言之，若输入码元速率 R s ' 1/ T s ，则该基带传输系统输出码元会产生码
基带信号眼图实验——matlab仿真
基带信号眼图实验——matlａb 仿真
————————————————————————————————作者：————————————————————————————————日期： ?
数字基带信号的眼图实验——matla ｂ仿真 一、实验目的 1、掌握无码间干扰传输的基本条件和原理，掌握基带升余弦滚降系统的实现方法； 2、通过观察眼图来分析码间干扰对系统性能的影响，并观察在输入相同码率的ＮRZ 基带信号下，不同滤波器带宽对输出信号码间干扰大小的影响程度； 3、熟悉MATL ＡB 语言编程。 二、实验预习要求 1、复习《数字通信原理》第七章７.1节——奈奎斯特第一准则内容； 2、复习《数字通信原理》第七章7．2节——数字基带信号码型内容； 3、认真阅读本实验内容,熟悉实验步骤。 三、实验原理和电路说明 １、基带传输特性 基带系统的分析模型如图3-1所示,要获得良好的基带传输系统，就应该 () n s n a t nT δ-∑() H ω() n s n a h t nT -∑基带传输抽样判决 图3-１?基带系统的分析模型 抑制码间干扰。设输入的基带信号为()n s n a t nT δ-∑，s T 为基带信号的码元周期,则经过基 带传输系统后的输出码元为 ()n s n a h t nT -∑。其中 1 ()()2j t h t H e d ωωωπ +∞ -∞ = ? ?（3-1） 理论上要达到无码间干扰,依照奈奎斯特第一准则,基带传输系统在时域应满足： 10()0,s k h kT k =?=? ? ， 为其他整数 ?? ?(３－２) 频域应满足：
数字基带信号的传输码型
基带数字信号的表示和传输 图1-1：基带传输模型图 1）信号形成器：产生适合于信道传输的基带信号波形。 2）信道：允许基带信号通过的媒介。 3）接收滤波器：用来接收信号，尽可能滤除信道噪声和其他干扰，对信道特性进行均衡，使输出的基带信号有利于判决。 4）抽样判决器：在传输特性不理想及噪声背景下，在规定时刻（由位定时脉冲控制）对将接收滤波器的输出波形进行抽样判决，以恢复或再生基带信号。 5）定时脉冲和同步提取：用来抽样的位定时脉冲依靠同步提取电路从接收信号中提取，位定时的准确性将直接影响判决效果。 2.常见的数字基带传输码型 （1）AMI AMI(Alternative Mark Inversion)码的全称是传号交替反转码，其编码规则是三元码，“1”交替地变换为“+1”和“-1”，“0”保持不变采用归零码，脉冲宽度为码元宽度之半“0”，“1”不等概时也无直流；零频附近的低频分量小；频率集中在1/2码速处；编解码电路简单，且可以利用传号极性交替这一规律观察五码情况；整流成归零码之后，从中可以提取定时分量。 连0码多时，AMI 整流后的RZ 码连0也多，不利于提取高质量的位同步信号。AMI 码的波形图如图1-6所示： 1 011100000000111 +1-1000000000+1+1+1-1-1二进制码 二进制波形AMI 波形 AMI 码 图1-4 数字基带传输系统模型
图1-6 AMI 码波形 （2）HDB 3码 HDB 3(3nd Order High Density Bipolar)码的全称是三阶高密度双极性码， 是AMI 码的一种改进，保持了AMI 码的优点，使“0”连续不超过3个。其编码规则为：“1”交替地变换为＋1与－1的半占空归零码，但连“0”数小于或者等于3。当连“0”数等于4时，用取代节“000V ”或者“B00V ”代替，“V ”的极性与前一个非零符号的极性相同(这破坏了极性交替的规则，所以V 又称为破坏脉冲)；并要求相邻的“V ”也满足极性必须交替。V 的取值为+1或-1.B 的取值可以是0、+1、-1，以使V 同时满足(3)中的要求。 HDB3码波形如下： 1 011100000000111 +1-10000+V -B 00-V +1-1 +1-1+1二进制码 二进制波形HDB3码波形 HDB3码图1-7 HDB 3码波形 （3）双相码 双向码又称为曼彻斯特(Manchester)码，用一个周期的正负对称方波表示“0”，而用其反相波形表示“1”，其编码规则：“1”用“10”表示，“0”用“01”表示，是一种双极性不归零波形，只有极性相反的两个电平；每个码元中心都有电平跳变，含有丰富的定时信息，且没有直流分量，编码过程也简单；缺点是占用带宽加宽，使频率利用率降低。双相码波形如下： 1 01110000011 二进制码 二进制波形双相码波形 双相码1001100110101010 01010101 图1-8 双相码波形
根据Matlab的信号与系统实验指导材料
基于Matlab 的信号与系统实验指导 实验一 连续时间信号在Matlab 中的表示 一、实验目的 1、学会运用Matlab 表示常用连续时间信号的方法 2、观察并熟悉这些信号的波形和特性 二、实验原理及实例分析 1、信号的定义与分类 2、如何表示连续信号？ 连续信号的表示方法有两种；符号推理法和数值法。 从严格意义上讲，Matlab 数值计算的方法不能处理连续时间信号。然而，可利用连续信号在等时间间隔点的取样值来近似表示连续信号，即当取样时间间隔足够小时，这些离散样值能被Matlab 处理，并且能较好地近似表示连续信号。 3、Matlab 提供了大量生成基本信号的函数。如： （1）指数信号：K*exp(a*t) （2）正弦信号：K*sin(w*t+phi)和K*cos(w*t+phi) （3）复指数信号：K*exp((a+i*b)*t) （4）抽样信号：sin(t*pi) 注意：在Matlab 中用与Sa(t)类似的sinc(t)函数表示，定义为：)t /()t (sin )t (sinc ππ= （5）矩形脉冲信号：rectpuls(t,width) （6）周期矩形脉冲信号：square(t,DUTY)，其中DUTY 参数表示信号的占空比
DUTY%，即在一个周期脉冲宽度（正值部分）与脉冲周期的比值。占空比默认为0.5。 （7）三角波脉冲信号：tripuls(t, width, skew)，其中skew 取值范围在-1~+1之间。 （8）周期三角波信号：sawtooth(t, width) （9）单位阶跃信号：y=(t>=0) 三、实验内容 1、验证实验内容 直流及上述9个信号 2、程序设计实验内容 （1）利用Matlab 命令画出下列连续信号的波形图。 （a ）)4/3t (2cos π+ （b ） )t (u )e 2(t -- （c ）)]2()(u )][t (cos 1[--+t u t π （2）利用Matlab 命令画出复信号) 4/t (j 2e )t (f π+=的实部、虚部、模和辐角。 四、实验报告要求 1、格式：实验名称、实验目的、实验原理、实验环境、实验内容、实验思考等 2、实验内容：程序设计实验部分源代码及运行结果图示。
基于MATLAB的QAM 眼图和星座图
南昌大学信息工程学院 《随机信号分析》课程作业 题目：QAM调制信号的眼图及星座图仿真指导老师：虞贵财 作者：毕圣昭 日期：2011-12-05
QAM调制信号的眼图及星座图仿真 1. 眼图 眼图是在数字通信的工程实践中测试数字传输信道质量的一种应用广泛、简单易行的方法。实际上它的一个扫描周期是数据码元宽度1~2倍并且与之同步的示波器。对于二进制码元，显然1和0的差别越大，接受判别时错判的可能性就越小。由于传输过程中受到频带限制，噪声的叠加使得1和0的差别变小。在接收机的判决点，将“1”和“0”的差别用眼图上“眼睛”张开的大小来表示，十分形象、直观和实用。MATLAB工具箱中有显示眼图和星座图的仪器，下面通过具体的例子说明它们的应用。 图1-1所示是MATLAB Toolbox\Commblks中的部分内容，展示了四进制随机数据通过基带QPSK调制、升余弦滤波（插补）及加性高斯白噪声传输环境后信号的眼图。 图1-1 通过QPSK基带调制升余弦滤波及噪声环境后观察眼图的仿真实验系统 图1-2所示是仿真运行后的两幅眼图，上图是I（同相）信号，下图是Q（正交）信号。 图1-2 通过QPSK基带调制及噪声传输环境后观察到的眼图
2. 星座图 星座图是多元调制技术应用中的一种重要的测量方法。它可以在信号空间展示信号所在的位置，为系统的传输特性分析提供直观的、具体的显示结果。 为了是系统的功率利用率、频带利用率得到充分的利用，在特定的调制方式下，在信号空间中如何排列与分布信号？在传输过程中叠加上噪声以后，信号之间的最小距离是否能保证既定的误码率的要求这些问题的研究用星座图仪十分直观方便。多元调制都可以分解为In-phase（同相）分量及Quadrature（正交）分量。将同相分量用我们习惯的二维空间的Ｘ轴表示，正交分量用Ｙ轴表示。信号在Ｘ－Ｙ平面（同相－正交平面）的位置就是星座图。MATLAB通信系统的工具箱里有着使用方便、界面美观的星座图仪。 图1-3所示是随机数据通过基带QAM调制及噪声环境传输后，观察星座图的仿真系统。 图1-3 通过基带QAM调制及噪声环境传输后观察星座图的仿真系统图1-4所示是运行仿真后的星座图 图1-4 通过基带QAM调制及噪声环境传输后观察到的星座图