产生离散衰减正弦序列()π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通信原理仿真 学号： 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);
数字基带传输实验 实验报告
一、实验目的 1、提高独立学习的能力； 2、培养发现问题、解决问题和分析问题的能力； 3、学习Matlab 的使用； 4、掌握基带数字传输系统的仿真方法； 5、熟悉基带传输系统的基本结构； 6、掌握带限信道的仿真以及性能分析； 7、通过观测眼图和星座图判断信号的传输质量。 二、系统框图及编程原理 1.带限信道的基带系统模型（连续域分析） ?输入符号序列―― ?发送信号―― ――比特周期，二进制码元周期 ?发送滤波器―― 或或 ?发送滤波器输出――
?信道输出信号或接收滤波器输入信号 （信道特性为1） ?接收滤波器―― 或或 ?接收滤波器的输出信号 其中 （画出眼图） ?如果位同步理想，则抽样时刻为 ?抽样点数值为（画出星座图） ?判决为 2.升余弦滚降滤波器 式中称为滚降系数，取值为, 是常数。时，带宽为Hz；时，带宽为Hz。此频率特性在内可以叠加成一条直线，故系统无码间干扰传输的最小符号间隔为s，或无码间干扰传输的最大符号速率为Baud。
相应的时域波形为 此信号满足 在理想信道中，，上述信号波形在抽样时刻上无码间干扰。 如果传输码元速率满足，则通过此基带系统后无码间干扰。 3.最佳基带系统 将发送滤波器和接收滤波器联合设计为无码间干扰的基带系统，而且具有最佳的抗加性高斯白噪声的性能。 要求接收滤波器的频率特性与发送信号频谱共轭匹配。由于最佳基带系统的总特性是确定的，故最佳基带系统的设计归结为发送滤波器和接收滤波器特性的选择。 设信道特性理想，则有
（延时为0） 有 可选择滤波器长度使其具有线性相位。 如果基带系统为升余弦特性，则发送和接收滤波器为平方根升余弦特性。 由模拟滤波器设计数字滤波器的时域冲激响应 升余弦滤波器（或平方根升余弦滤波器）的带宽为，故其时域抽样速率至少为，取，其中为时域抽样间隔，归一化为1。 抽样后，系统的频率特性是以为周期的，折叠频率为。故在一个周期内 以间隔抽样，N为抽样个数。频率抽样为，。 相应的离散系统的冲激响应为 将上述信号移位，可得因果系统的冲激响应。 5．基带传输系统（离散域分析） ?输入符号序列―― ?发送信号―― ――比特周期，二进制码元周期 ?发送滤波器――
实验1：上采样与内插 一、实验目的 1、了解上采样与内插的基本原理和方法。 2、掌握上采样与内插的matlab程序的设计方法。 二、实验原理 上采样提高采样频率。上采样使得周期降低M倍，即新采样周期Tu和原有采样周期Ts的关系是T u=T s/M,根据对应的连续信号x(t),上采样过程从原有采样值x(kT s)生成新采样值x(kT u)=x(kT s/M)。操作的结果是在每两个采样值之间放入M-1个零值样点。 更实用的内插器是线性内插器，线性内插器的脉冲响应定义如下： 上采样值x(kT u)=x(kT s/M)通过与线性内插器的脉冲响应的卷积来完成内插。 三、实验内容 仿真正弦波采样和内插，通过基本采样x（k）,用M=6产生上采样x u(k)，由M=6线性内插得到样点序列x i（k）。 四、实验程序 % File: c3_upsampex.m M = 6; % upsample factor h = c3_lininterp(M); % imp response of linear interpolator t = 0:10; % time vector tu = 0:60; % upsampled time vector x = sin(2*pi*t/10); % original samples xu = c3_upsamp(x,M); % upsampled sequence subplot(3,1,1) stem(t,x,'k.') ylabel('x') subplot(3,1,2) stem(tu,xu,'k.') ylabel('xu') xi = conv(h,xu); subplot(3,1,3) stem(xi,'k.') ylabel('xi') % End of script file. % File: c3_upsample.m function out=c3_upsamp(in,M)
通信原理实验报告 实验一 数字基带传输实验 一、实验目的 1、提高独立学习的能力； 2、培养发现问题、解决问题和分析问题的能力； 3、学习Matlab 的使用； 4、掌握基带数字传输系统的仿真方法； 5、熟悉基带传输系统的基本结构； 6、掌握带限信道的仿真以及性能分析； 7、通过观测眼图和星座图判断信号的传输质量。 二、实验原理 1. 带限信道的基带系统模型（连续域分析） 输入符号序列 ————｛al ｝ 发送信号 ————1 0()()L l d t al t lTb δ-==-∑ Tb 是比特周期，二进制码元周期 发送滤波器 ————GT(w)或GT （t ）
发送滤波器输出 ———— 11 00()()*()()*()()L L l b T l T b T l l x t d t t a t lT g t a g t lT g δ--====-=-∑∑ 信道输出信号或接收滤波器输入信号()()()y t x t n t =+ 接收滤波器 ()R G ω或()R G f 接收滤波器输出信号 1 0()()*()()*()*()()*()()()L R T R R l b R l r t y t g t d t g t g t n t g t a g t lT n t -===+=-+∑ 其中2()()()j ft T R g t G f G f e df π∞ -∞=? 如果位同步理想，则抽样时刻为b l T ? 0 1l L =- 判决为 '{}l a 2. 升余弦滚降滤波器 1()||2s s H f T f T α-=≤； ()H f =111[1cos (||)]||2222s s s s s T T f f T T T παααα--++-<≤ ()H f = 10||2s f T α+> 式中α 称为滚降系数，取值为0 <α ≤1, T s 是常数。α = 0时，带宽为1/ 2T s Hz ；α =1时， 带宽为1/T s Hz 。此频率特性在(?1/(2T s ),1/(2T s ))内可以叠加成一条直线，故系统无码间干 扰传输的最小符号间隔为T s s ，或无码间干扰传输的最大符号速率为1/T s Baud 。相应的时 域波形h (t )为 222sin /cos /()/14/s s s s t T t T h t t T t T παππα=?- 此信号满足
《信号与系统》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验证低通抽样定理 二、实验目的 1、掌握抽样定理的工作原理。 2、通过MATLAB编程实现对抽样定理的验证，加深抽样定理的理解。同时训练应用计算机分析问题的能力。 3、了解MATLAB软件，学习应用MATLAB软件的仿真技术。它主要侧重于某些理论知识的灵活运用，以及一些关键命令的掌握，理解，分析等。 4、计算在临界采样、过采样、欠采样三种不同条件下恢复信号的误差，并由此总结采样频率对信号恢复产生误差的影响，从而验证时域采样定理。 三、实验步骤及原理 1、对连续信号进行等间隔采样形成采样信号，采样信号的频谱是原连续信号的频谱以采样频率为周期进行周期性的延拓形成的。 2、设连续信号的的最高频率为Fmax，如果采样频率Fs>2Fmax，那么采样信号可以唯一的恢复出原连续信号，否则Fs<=2Fmax会造成采样信号中的频谱混叠现象，不可能无失真地恢复原连续信号。 四、实验内容 1、画出连续时间信号的时域波形及其幅频特性曲线，信号为 x=cos(4*pi*t)+1.5*sin(6*pi*t)+0.5*cos(20*pi*t) 2、对信号进行采样，得到采样序列，画出采样频率分别为10Hz，20 Hz，50 Hz时的采样序列波形； 3、对不同采样频率下的采样序列进行频谱分析，绘制其幅频曲线，对比各频率下采样序列和的幅频曲线有无差别。 4、对信号进行谱分析，观察与3中结果有无差别。 5、由采样序列恢复出连续时间信号，画出其时域波形，对比与原连续时间信号的时域波形。 五、实验仿真图 (1) x=cos(4*pi*t)+1.5*sin(6*pi*t)+0.5*cos(20*pi*t)的时域波形及幅频特性曲线。clear; close all; dt=0.05; t=-2:dt:2 x=cos(4*pi*t)+1.5*sin(6*pi*t)+0.5*cos(20*pi*t); N=length(t); Y=fft(x)/N*2; fs=1/dt; df=fs/(N-1); f=(0:N-1)*df;
习题三 绘制典型信号及其频谱图 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 )');
通信原理实验报告实验名称：采样定理 实验时间： 201211日年12月 指导老师：应娜 学院：计算机学院 级：班 学号： 姓名：
通信原理实验报告 一、实验名称 MATLAB验证低通抽样定理 二、实验目的 1、掌握抽样定理的工作原理。 2、通过MATLAB编程实现对抽样定理的验证，加深抽样定理的理解。同时训练应用计算机分析问题的能力。 3、了解MATLAB软件，学习应用MATLAB软件的仿真技术。它主要侧重于某些理论知识的灵活运用，以及一些关键命令的掌握，理解，分析等。 4、计算在临界采样、过采样、欠采样三种不同条件下恢复信号的误差，并由此总结采样频率对信号恢复产生误差的影响，从而验证时域采样定理。 三、实验步骤 1、画出连续时间信号的时域波形及其幅频特性曲线，信号为 f(x)=sin(2*pi*80*t)+ cos(2*pi*30*t)； 2、对信号进行采样，得到采样序列，画出采样频率分别为80Hz，110 Hz，140 Hz时的采样序列波形； 3、对不同采样频率下的采样序列进行频谱分析，绘制其幅频曲线，对比各频率下采样序列和的幅频曲线有无差别。 4、对信号进行谱分析，观察与3中结果有无差别。 5、由采样序列恢复出连续时间信号，画出其时域波形，对比与原连续时间信号的时域波形。 四、数据分析 (1)部分程序分析： f=[fs0*k2/m2,fs0*k1/m1]; %设置原信号的频率数组 axis([min(t),max(t),min(fx1),max(fx1)]) %画原信号幅度频谱 f1=[fs*k2/m2,fs*k1/m1]; %设置采样信号的频率数组 fz=eval(fy); %获取采样序列 FZ=fz*exp(-j*[1:length(fz)]'*w); %采样信号的离散时间傅里叶变换 TMN=ones(length(n),1)*t-n'*T*ones(1,length(t)); 由采样信号恢复原信号fh=fz*sinc(fs*TMN); %． (2)原信号的波形与幅度频谱：
实验一 设计任务： 用MatLib仿真一个BFSK通信系统，基本参数： 1）fc=1000Hz； 2）Rb=100bps； 3）信息序列:“Hello world”的ASCII 实验与报告基本要求： 1）Matlab程序，要点旁注（可打印后手写）； 2）绘出信号波形，绘出信号PSD； 3）给出解调后的信息序列； 4）将信息重复3遍以上,FSK信号保存为WAV文件格式，使用音频播放，聆听；M文件： wave.m function[t,mt]=wave(m,dt,fs) l=length(m); mt=[]; ddt=1/fs; n=floor(dt*fs); m_add=ones(1,n); for i=1:l if(m(i)) mt=[mt,m(i),m_add]; else mt=[mt,m(i),m_add*0]; end t=(1:((n+1)*l))*ddt; end my_filter.m function[num,den]=my_filter(wp,ws,ap,as) if nargin<4 as=15; end if nargin<4 ap=3; end [n,wn]=buttord(wp,ws,ap,as); [num,den]=butter(n,wn); end 代码：
f0=800;%‘0’码载波频率 f1=1200;%‘1’码载波频率 fs=4000;%采样频率 Rb=100;%比特率 dt=1/Rb;%一个比特发送时间 A0=2;%调制幅度 A1=2;%相干解调幅度 miu=0;sigma=0.3;%miu:高斯白噪声均值，sigma：高斯白噪声均方差 str='Hello world';%信号字符串 m_dec=abs(str);%将信号字符串转换成ASCII码(十进制) m_bin=dec2bin(m_dec,8); m_bin=abs(m_bin)-48;%将十进制转换成8比特二进制矩阵 m=[]; for i=1:size(m_bin,1) m=[m,m_bin(i,:)]; end%将二进制转换成行向量 [t,m]=wave(m,dt,fs);%对信号采样 mt_f1=m.*cos(2*pi*f1*t)*A0;%频率f1调制 mt_f0=(~m).*cos(2*pi*f0*t)*A0;%频率f0调制 mt=mt_f1+mt_f0;%发送信号 l=length(mt); subplot(2,1,1);plot(t,mt); grid on;xlabel('t/s');title('m(t)');%发送信号波形subplot(2,1,2);periodogram(mt,[],l,fs);grid on;%发送信号PSD
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光滑得多。
《Matlab与通信仿真》实验指导书(上) 刘毓杨辉徐健和煦黄庆东吉利萍编著 通信与信息工程学院 2011-1
目录 第一章 MALTAB基础知识 (1) 1.1MATLAB基础知识 (1) 1.2MATLAB基本运算 (2) 1.3MATLAB程序设计 (7) 第二章 MATLAB计算结果可视化和确知信号分析 (13) 2.1计算结果可视化 (13) 2.2确知信号分析 (17) 第三章随机信号与数字基带仿真 (23) 3.1基本原理 (23) 3.2蒙特卡罗算法 (30) 第四章模拟调制MATLAB实现 (34) 4.1模拟调制 (34) 4.2信道加性高斯白噪声 (35) 4.3AM调制解调的MATLAB实现 (36) 第五章模拟信号的数字传输 (45) 5.1脉冲编码调制 (45) 5.2低通抽样定理 (45) 5.3均匀量化原理 (46) 5.4非均匀量化 (48) 第六章数字频带传输系统 (52) 6.1数字频带传输原理 (52) 6.2数字频带传输系统的MATLAB实现 (53) 第七章通信系统仿真综合实验 (67) 7.1基本原理 (67) 7.2实验内容 (67)
第一章 MALTAB基础知识 本章目标 ●了解MATLAB 程序设计语言的基本特点，熟悉MATLAB软件运行环境 ●掌握创建、保存、打开m文件及函数的方法 ●掌握变量等有关概念，具备初步的将一般数学问题转化为对应的计算机模型并进行处理的能力 1.1 MATLAB基础知识 1.1.1 MATLAB程序设计语言简介 MATLAB，Matrix Laboratory的缩写，是由MathWorks公司开发的一套用于科学工程计算的可视化高性能语言，具有强大的矩阵运算能力。与大家常用的Fortran和C等高级语言相比，MATLAB的语法规则更简单，更贴近人的思维方方式，被称为“草稿纸式的语言”。MATLAB软件主要由主包、仿真系统（simulink）和工具箱（toolbox）三大部分组成。 1.1.2 MATLAB界面及帮助 MATLAB基本界面如图1-1所示，命令窗口包含标题栏、菜单栏、工具栏、命令行区、状态栏、垂直和水平波动条等区域。 图1-1 MATLAB基本界面 （1）菜单栏
通信原理软件实验报告 学院：信息与通信工程学院 班级： 一、通信原理Matlab仿真实验 实验八 一、实验内容 假设基带信号为m(t)=sin(2000*pi*t)+2cos(1000*pi*t),载波频率为20kHz，请仿真出AM、DSB-SC、SSB信号，观察已调信号的波形和频谱。 二、实验原理 1、具有离散大载波的双边带幅度调制信号AM 该幅度调制是由DSB-SC AM信号加上离散的大载波分量得到，其表达式及时间波形图为： 应当注意的是，m(t)的绝对值必须小于等于1，否则会出现下图的过调制： AM信号的频谱特性如下图所示： 由图可以发现，AM信号的频谱是双边带抑制载波调幅信号的频谱加上离散的大载波分量。 2、双边带抑制载波调幅（DSB—SC AM）信号的产生 双边带抑制载波调幅信号s(t)是利用均值为0的模拟基带信号m(t)和正弦载波 c(t)相乘得到，如图所示： m(t)和正弦载波s(t)的信号波形如图所示：
若调制信号m(t)是确定的，其相应的傅立叶频谱为M(f)，载波信号c(t)的傅立叶频谱是C(f),调制信号s(t)的傅立叶频谱S(f)由M(f)和C(f)相卷积得到，因此经过调制之后，基带信号的频谱被搬移到了载频fc处，若模拟基带信号带宽为W，则调制信号带宽为2W，并且频谱中不含有离散的载频分量，只是由于模拟基带信号的频谱成分中不含离散的直流分量。 3、单边带条幅SSB信号 双边带抑制载波调幅信号要求信道带宽B=2W, 其中W是模拟基带信号带宽。从信息论关点开看，此双边带是有剩余度的，因而只要利用双边带中的任一边带来传输，仍能在接收机解调出原基带信号，这样可减少传送已调信号的信道带宽。 单边带条幅SSB AM信号的其表达式： 或 其频谱图为： 三、仿真设计 1、流程图：
MATLA仿真实验报告 学院：计算机与信息学院 课程：—随机信号分析 姓名： 学号： 班级： 指导老师： 实验一
题目：编写一个产生均值为1,方差为4的高斯随机分布函数程序, 求最大值，最小值，均值和方差，并于理论值比较。 解：具体的文件如下，相应的绘图结果如下图所示 G仁random( 'Normal' ,0,4,1,1024); y=max(G1) x=mi n(G1) m=mea n(G1) d=var(G1) plot(G1);
实验二 题目：编写一个产生协方差函数为CC)=4e":的平稳高斯过程的程序，产生样本函数。估计所产生样本的时间自相关函数和功率谱密度，并求统计自相关函数和功率谱密度，最后将结果与理论值比较。 解：具体的文件如下，相应的绘图结果如下图所示。 N=10000; Ts=0.001; sigma=2; beta=2; a=exp(-beta*Ts); b=sigma*sqrt(1-a*a); w=normrnd(0,1,[1,N]); x=zeros(1,N); x(1)=sigma*w(1); for i=2:N x(i)=a*x(i-1)+b*w(i); end %polt(x); Rxx=xcorr(x0)/N; m=[-N+1:N-1]; Rxx0=(sigma A2)*exp(-beta*abs(m*Ts)); y=filter(b,a,x) plot(m*Ts,RxxO, 'b.' ,m*Ts,Rxx, 'r');
periodogram(y,[],N,1/Ts); 文件旧硯化）插入（1〕 ZMCD 克闻〔D ］窗口曲） Frequency (Hz) 50 100 150 200 250 300 350 400 450 500 NH---.HP)&UO 二 balj/ 」- □歹
实验一、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]);
实验教程
目录 实验一：连续时间信号与系统的时域分析-------------------------------------------------6 一、实验目的及要求---------------------------------------------------------------------------6 二、实验原理-----------------------------------------------------------------------------------6 1、信号的时域表示方法------------------------------------------------------------------6 2、用MATLAB仿真连续时间信号和离散时间信号----------------------------------7 3、LTI系统的时域描述-----------------------------------------------------------------11 三、实验步骤及内容--------------------------------------------------------------------------15 四、实验报告要求-----------------------------------------------------------------------------26 实验二：连续时间信号的频域分析---------------------------------------------------------27 一、实验目的及要求--------------------------------------------------------------------------27 二、实验原理----------------------------------------------------------------------------------27 1、连续时间周期信号的傅里叶级数CTFS---------------------------------------------27 2、连续时间信号的傅里叶变换CTFT--------------------------------------------------28 3、离散时间信号的傅里叶变换DTFT -------------------------------------------------28 4、连续时间周期信号的傅里叶级数CTFS的MATLAB实现------------------------29 5、用MATLAB实现CTFT及其逆变换的计算---------------------------------------33 三、实验步骤及内容----------------------------------------------------------------------34 四、实验报告要求-------------------------------------------------------------------------48 实验三：连续时间LTI系统的频域分析---------------------------------------------------49 一、实验目的及要求--------------------------------------------------------------------------49 二、实验原理----------------------------------------------------------------------------------49 1、连续时间LTI系统的频率响应-------------------------------------------------------49 2、LTI系统的群延时---------------------------------------------------------------------50 3、用MATLAB计算系统的频率响应--------------------------------------------------50 三、实验步骤及内容----------------------------------------------------------------------51 四、实验报告要求-------------------------------------------------------------------------58 实验四：调制与解调以及抽样与重建------------------------------------------------------59 一、实验目的及要求--------------------------------------------------------------------------59 二、实验原理----------------------------------------------------------------------------------59
题目: 数字信号处理MATLAB仿真实验 姓名 学院 专业 班级 学号 班内序号
实验一：数字信号的 FFT 分析 1、实验内容及要求 (1) 离散信号的频谱分析： 设信号 此信号的0.3pi 和 0.302pi 两根谱线相距很近，谱线 0.45pi 的幅度很小，请选择合适的序列长度 N 和窗函数，用 DFT 分析其频谱，要求得到清楚的三根谱线。 (2) DTMF 信号频谱分析 用计算机声卡采用一段通信系统中电话双音多频（DTMF ）拨号数字 0～9的数据，采用快速傅立叶变换（FFT ）分析这10个号码DTMF 拨号时的频谱。 2、实验目的 通过本次实验，应该掌握： (a) 用傅立叶变换进行信号分析时基本参数的选择。 (b) 经过离散时间傅立叶变换（DTFT ）和有限长度离散傅立叶变换（DFT ） 后信号频谱上的区别，前者 DTFT 时间域是离散信号，频率域还是连续的，而 DFT 在两个域中都是离散的。 (c) 离散傅立叶变换的基本原理、特性，以及经典的快速算法（基2时间抽选法），体会快速算法的效率。 (d) 获得一个高密度频谱和高分辨率频谱的概念和方法，建立频率分辨率和时间分辨率的概念，为将来进一步进行时频分析（例如小波）的学习和研究打下基础。 (e) 建立 DFT 从整体上可看成是由窄带相邻滤波器组成的滤波器组的概念，此概念的一个典型应用是数字音频压缩中的分析滤波器，例如 DVD AC3 和MPEG Audio 。 3、程序代码 (1) N=5000; n=1:1:N; x=0.001*cos(0.45*pi*n)+sin(0.3*pi*n)-cos(0.302*pi*n-pi/4); y=fft(x,N); magy=abs(y(1:1:N/2+1)); k=0:1:N/2; w=2*pi/N*k; stem(w/pi,magy) axis([0.25,0.5,0,50]) （2） column=[1209,1336,1477,1633]; line=[697,770,852,941]; fs=10000; N=1024; 00010450303024().*cos(.)sin(.)cos(.)x n n n n ππππ=+--
实验一 基本信号的产生与运算 一、 实验目的 学习使用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('杨婕婕朱艺星');
补充内容：模拟调制系统的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; %%根据抽样时间间隔进行抽样