实验一 离散时间信号分析
实验一 离散时间信号的频域分析
实验一离散时间信号的频域分析一实验目的：信号的变换域分析是信号处理中一种有效的工具。
在离散信号的时域分析中，我们通常将信号表示成单位采样序列δ[n]的线性组合，而在频域中，我们将信号表示成复变量e-jwn或e-j(2*pi/N)n的线性组合，通过这样的表示，可以将时域的离散序列映射到频域以便于进一步的处理。
在本实验中，将学习利用MATLAB计算离散时间信号的DTFT和DFT，并加深对其相互关系的理解。
二实验原理：（1） DTFT和DFT的定义及其相互关系：序列x[n]的DTFT定义：X（e jw）=∑x[n]e-jnw（n从负无穷到正无穷）它是关于自变量w的复函数，且是以2*pi为周期的连续函数。
X（e jw）可以表示为：X（e jw）=X re（e jw）+ jX im（e jw）其中，X re（e jw）和X im（e jw）分别是X （e jw）的实部和虚部；还可以表示为：X（e jw）= |X（e jw）|e jØ(w)其中。
|X（e jw）|和Ø(w)=arg{ X（e jw）}分别是X（e jw）的幅度函数和相位函数；它们都是w的实函数，也是以2*pi为周期的周期函数。
序列x[n]的N点DFT定义：X[k]=X(e j(2*pi/N)k)=∑x[n]e-j(2*pi/N)kn (0<=n<=N-1)X[k]是周期为N的序列。
X（e jw）与X[k]的关系：X[k]是对X（e jw）在一个周期中的谱的等间隔N点采样，即：X[k]= X（e jw）|w=(2*pi/N)k而X（e jw）可以通过对X[k]内插获得。
（2）使用MATLAB命令：A.基于DTFT离散时间信号分析函数:freqz,real,imag,abs,angle,unware.函数freqz可以用来计算一个以e jw的有理分式形式给出的序列的DTFT值。
Freqz的形式多样，常见的有H=freqz(num,den,w),其中num表示序列有理分式DTFT的分子多项式系数，den表示分母多项式系数（均按z的降幂排列），矢量w表示在0到2*pi中给定的一系列频率集合点。
实验一 离散时间信号的时域分析
实验一 离散时间信号的时域分析实验1 序列的产生1． 目的：熟悉C 语言产生和绘制，熟悉MATLAB 中产生信号和绘制信号的基本命令。
2． 具体实验：2.1 单位样本和单位阶跃序列。
Q1.1 运行程序P1.1 ，以产生单位样本序列u[n]并显示它。
答：如图1-1所示。
Q1.2 命令clf , axis , title , xlabel 和ylabel 的作用是什么？ 答：clf ：擦除当前图形窗口中的图形。
Axis ：调整坐标轴X 轴Y 轴的范围。
Title:给绘制的图形加上标题。
Xlabel:给X 轴加上标注。
Ylabel: 给Y 轴加上标注。
Q1.3 修改程序P1.1 以产生带有延时11个样本的延迟单位样本序列ud[n]。
运行修改的程序并显示产生的序列。
答：如图1-2所示。
Q1.4修改程序P1.1 以产生单位步长序列s[n]。
运行修改的程序并显示产生的序列。
答：如图1-3所示。
Q1.5修改程序P1.1 以产生带有超前7个样本的延迟单位样本序列sd[n]。
运行修改的程序并显示产生的序列。
答：如图1-4所示。
Figure 1-2 The unit sample sequence ud[n]Figure 1-1 The shifted unit sample sequence u[n]单位样本序列公式如下所示：Time index nA m p l i t u d eUnit Sample Sequence u[n]Time index nA m p l i t u d eShifted Unit Sample Sequence ud[n]1 , n=0 1 , n=k δ[n]= δ[n-k]=0 , 0≠0 0 , 0≠kFigure 1-3 The unit step sequence s[n] Figure 1-4 The shifted unit step sequence sd[n]单位阶跃序列公式如下所示：1 , n ≥0 1 , n ≥k μ[n]= μ[n-k]=0 , n ＜0 0 , n ＜k2.2 指数信号Q1.6 运行程序P1.2 ，以产生复数值的指数序列。
离散时间信号的时域分析实验报告
离散时间信号的时域分析实验报告实验报告：离散时间信号的时域分析一、实验目的本实验旨在通过MATLAB软件，对离散时间信号进行时域分析，包括信号的显示、基本运算（如加法、减法、乘法、反转等）、以及频域变换（如傅里叶变换）等，以加深对离散时间信号处理的基本概念和原理的理解。
二、实验原理离散时间信号是在时间轴上离散分布的信号，其数学表示为离散时间函数。
与连续时间信号不同，离散时间信号只能在特定的时间点取值。
离散时间信号的时域分析是研究信号的基本属性，包括幅度、时间、频率等。
通过时域分析，我们可以对信号进行各种基本运算和变换，以提取有用的信息。
三、实验步骤1.信号生成：首先，我们使用MATLAB生成两组简单的离散时间信号，一组为正弦波，另一组为方波。
我们将这些信号存储在数组中，以便后续分析和显示。
2.信号显示：利用MATLAB的绘图功能，将生成的信号在时域中显示出来。
这样，我们可以直观地观察信号的基本属性，包括幅度和时间关系。
3.基本运算：对生成的信号进行基本运算，包括加法、减法、乘法、反转等。
将这些运算的结果存储在新的数组中，并绘制出运算后的信号波形。
4.傅里叶变换：使用MATLAB的FFT（快速傅里叶变换）函数，将信号从时域变换到频域。
我们可以得到信号的频谱，进而分析信号的频率属性。
5.结果分析：对上述步骤得到的结果进行分析，包括比较基本运算前后的信号波形变化，以及傅里叶变换前后的频谱差异等。
四、实验结果1.信号显示：通过绘制图形，我们观察到正弦波和方波在时域中的波形特点。
正弦波呈现周期性的波形，方波则呈现明显的阶跃特性。
2.基本运算：通过对比基本运算前后的信号波形图，我们可以观察到信号经过加法、减法、乘法、反转等运算后，其波形发生相应的变化。
例如，两个信号相加后，其幅度和时间与原信号不同。
反转信号则使得波形在时间轴上反向。
3.傅里叶变换：通过FFT变换，我们将时域中的正弦波和方波转换到频域。
正弦波的频谱显示其频率为单一的直流分量，方波的频谱则显示其主要频率分量是直流分量和若干奇数倍的谐波分量。
实验一离散时间信号的时域分析
实验一离散时间信号的时域分析离散时间信号是一种离散的信号形式，其具有离散的时间间隔。
这种信号在数字信号处理中得到了广泛的应用。
时域分析是分析信号的一种方法，它通常包括分析信号的幅度、相位、频率等参数，并从中获得信号的特征。
在本实验中，我们将探讨离散时间信号的时域分析方法。
1.实验目的• 了解离散时间信号的基本概念和性质。
• 熟悉MATLAB软件的使用，理解信号处理工具箱的使用方法。
2.实验原理离散时间信号是一种在离散时间点上定义的数列。
它通常用序列来表示，序列的元素是按照一定的时间间隔离散采样得到的。
离散时间信号的采样频率通常表示为Fs，单位是赫兹。
离散时间信号可以写成如下的形式：x(n) = [x(0),x(1),x(2),...,x(N-1)]其中，n表示离散时间点的下标，N表示离散时间信号的长度。
• 幅度分析：指分析离散时间信号的振幅大小。
离散时间信号的幅度、相位、频率的分析通常使用傅里叶变换、离散傅里叶变换等变换方法来实现。
3.实验步骤3.1 生成离散时间信号使用MATLAB编写程序，生成一个离散时间信号。
例如，我们可以生成一个正弦信号：t = 0:0.01:1;x = sin(2*pi*100*t);其中，t表示时间向量，x表示正弦信号。
将信号进行离散化，得到离散时间信号：其中，fs表示采样频率，n表示采样时间点，xn表示采样后的信号。
使用MATLAB的plot函数，绘制离散时间信号的时域图像。
figure(1);plot(n, xn);xlabel('Time');ylabel('Amplitude');其中，figure(1)表示创建一个新的窗口，用于显示图像。
xlabel和ylabel用于设置图像的横轴和纵轴标签。
3.3 使用FFT进行幅度分析X = fft(xn);n = length(X);f = (0:n-1)*(fs/n);power = abs(X).^2/n;其中，X表示离散时间信号的傅里叶变换结果，n表示离散时间信号的长度，f表示频率向量，power表示幅度谱。
《数字信号处理》上机实验指导书
《数字信号处理》上机实验指导书实验1 离散时间信号的产生1．实验目的数字信号处理系统中的信号都是以离散时间形态存在，所以对离散时间信号的研究是数字信号处理的基本所在。
而要研究离散时间信号，首先需要产生出各种离散时间信号。
MATLAB 是一套功能强大的工程计算及数据处理软件，广泛应用于工业，电子，医疗和建筑等众多领域。
使用MATLAB软件可以很方便地产生各种常见的离散时间信号，而且它还具有强大的绘图功能，便于用户直观地输出处理结果。
通过本实验，学生将学习如何用MATLAB产生一些常见的离散时间信号，并通过MATLAB中的绘图工具对产生的信号进行观察，加深对常用离散信号的理解。
2．实验要求本实验要求学生运用MATLAB编程产生一些基本的离散时间信号，并通过MATLAB的几种绘图指令画出这些图形，以加深对相关教学内容的理解，同时也通过这些简单的函数练习了MATLAB的使用。
3．实验原理（1）常见的离散时间信号1）单位抽样序列，或称为离散时间冲激，单位冲激：?(n)???1?0n?0 n?0如果?(n)在时间轴上延迟了k个单位，得到?(n?k)即：?1n?k ?(n?k)??0n?0?2）单位阶跃序列n?0?1 u(n)?n?0?0在MATLAB中可以利用ones( )函数实现。
x?ones(1,N);3）正弦序列x(n)?Acos(?0n??)这里，A,?0,和?都是实数，它们分别称为本正弦信号x(n)的振幅，角频率和初始相位。
f0??02?为频率。
x(n)?ej?n4）复正弦序列5）实指数序列x(n)?A?n（2）MATLAB编程介绍MATLAB是一套功能强大，但使用方便的工程计算及数据处理软件。
其编程风格很简洁，没有太多的语法限制，所以使用起来非常方便，尤其对初学者来说，可以避免去阅读大量的指令系统，以便很快上手编程。
值得注意得就是，MATLAB中把所有参与处理的数据都视为矩阵，并且其函数众多，希望同学注意查看帮助，经过一段时间的训练就会慢慢熟练使用本软件了。
数字信号处理实验报告一二
数字信号处理课程实验报告实验一 离散时间信号和系统响应一. 实验目的1. 熟悉连续信号经理想采样前后的频谱变化关系，加深对时域采样定理的理解2. 掌握时域离散系统的时域特性3. 利用卷积方法观察分析系统的时域特性4. 掌握序列傅里叶变换的计算机实现方法，利用序列的傅里叶变换对离散信号及系统响应进行频域分析二、实验原理1. 采样是连续信号数字化处理的第一个关键环节。
对采样过程的研究不仅可以了解采样前后信号时域和频域特性的变化以及信号信息不丢失的条件，而且可以加深对离散傅里叶变换、Z 变换和序列傅里叶变换之间关系式的理解。
对连续信号()a x t 以T 为采样间隔进行时域等间隔理想采样，形成采样信号： 式中()p t 为周期冲激脉冲，()a x t 为()a x t 的理想采样。
()a x t 的傅里叶变换为()a X j Ω：上式表明将连续信号()a x t 采样后其频谱将变为周期的，周期为Ωs=2π/T 。
也即采样信号的频谱()a X j Ω是原连续信号xa(t)的频谱Xa(jΩ)在频率轴上以Ωs 为周期，周期延拓而成的。
因此，若对连续信号()a x t 进行采样，要保证采样频率fs ≥2fm ，fm 为信号的最高频率，才可能由采样信号无失真地恢复出原模拟信号ˆ()()()a a xt x t p t =1()()*()21()n a a a s X j X j P j X j jn T π∞=-∞Ω=ΩΩ=Ω-Ω∑()()n P t t nT δ∞=-∞=-∑计算机实现时，利用计算机计算上式并不方便，因此我们利用采样序列的傅里叶变换来实现，即而()()j j n n X e x n e ωω∞-=-∞=∑为采样序列的傅里叶变换2. 时域中，描述系统特性的方法是差分方程和单位脉冲响应，频域中可用系统函数描述系统特性。
已知输入信号，可以由差分方程、单位脉冲响应或系统函数求出系统对于该输入信号的响应。
实验一离散时间信号的分析报告
工程大学信号分析与处理实验一专业：通信02班学生：瑶华学号：**********完成时间：2022年4月27日实验一: 离散时间信号的分析一、实验目的1.认识常用的各种信号，理解其数学表达式和波形表示。
2.掌握在计算机中生成及绘制数字信号波形的方法。
3.掌握序列的简单运算及计算机实现与作用。
4.理解离散时间傅立叶变换、Z 变换及它们的性质和信号的频域特性。
二、实验设备计算机，MATLAB 语言环境。
三、实验基础理论1.序列的相关概念2.常见序列● 单位取样序列⎩⎨⎧≠==0n 0,0n 1n ，）（δ ● 单位阶跃序列⎩⎨⎧<≥=0,00,1)(n n n u ● 单位矩形序列⎩⎨⎧-≤≤=其他,010,1)(N n n R N ● 实指数序列)()(n u a n x n =● 复指数序列n jw e n x )(0)(+=σ● 正弦型序列)n sin()(0ϕ+=w A n x3.序列的基本运算● 移位 y(n)=x(n-m)● 反褶 y(n)=x(-n)● 和 )()()(21n x n x n y +=● 积 )()()(21n x n x n y •=● 标乘 y(n)=mx(n)● 累加∑-∞==nm m x n y )()( ● 差分运算 ⎩⎨⎧--=∇-+=∆)1()()()()1()(x n x n x n x n x n x n 后相差分前向差分 4.离散傅里叶变换的相关概念● 定义 ∑+∞-∞=-=n jwn jwe n x e X )(）（● 两个性质1) [])2()2()2()()(,2)(ππππ++∞-∞=+-+--===∑w j n nw j jw n w j jwn jw e X e n x e X e ew e X 故有。
由于的周期函数，周期为是 2) 当x （n ）为实序列时，)(jw e X 的幅值)(jw e X 在π20≤≤w 区间是偶对称函数，相位)(arg jw e X 是奇对称函数。
实验一离散时间信号与系统时域分析
实验一离散时间信号与系统时域分析实验目的1学习MATLAB语言编程和调试技巧2学会简单的矩阵输入和图形表示法3掌握简单的绘图命令一实验目的1学习MATLAB语言编程和调试技巧2学会简单的矩阵输入和图形表示法3掌握简单的绘图命令二、实验原理本实验主要为了熟悉MATLAB环境，重点掌握简单的矩阵（信号）输入和绘图命令，特别是绘图命令tem()和plot()。
实验内容中涉及到信号的无失真采样、离散卷积运算和差分方程求解这三个主要的问题。
其基本原理分别如下：对一个模拟信号某(t)进行采样离散化某(n)，为了不失真地从采样信号某(n)中恢复原始信号某(t)，采样时必须满足采样定理，即采样频率必须大于等于模拟信号中最高频率分量的2倍。
一个离散时间系统，输入信号为某(n)，输出信号为y(n)，运算关系用T[﹒]表示，则输入与输出的关系可表示为y(n)＝T[某(n)]。
（1）线性时不变(LTI)系统的输入输出关系可通过h(n)表示：y(n)=某(n)某h(n)=式中某表示卷积运算。
（2）LTI系统的实现可物理实现的线性时不变系统是稳定的、因果的。
这种系统的单位脉冲响应是因果的（单边）且绝对可和的，即：h(n)0,n0;nh(n)0在MATLAB语言中采用conv实现卷积运算，即：Y=conv(某,h),它默认从n=0开始。
常系数差分方程可以描述一个LTI系统，通过它可以获得系统的结构，也可以求信号的瞬态解。
利用MATLAB 自带的filter（），可以代替手工迭代运算求解系统的差分方程，求解的过程类似于对输入信号进行滤波处理。
三、实验内容1、试画出如下序列的波形（1）某(n)3(n3)(n2)2(n1)4(n1)2(n2)3(n3)（2）某(n)0.5R10(n)解：用MATLAB描述波形1(1)某=[3120-42-3];%矩阵输入某n=-3:1:3;%输入自变量n,以间隔为1从-3到3变化n实验目的1学习MATLAB语言编程和调试技巧2学会简单的矩阵输入和图形表示法3掌握简单的绘图命令tem(n,某);%tem()函数绘制火柴杆图，注意n,某元素个数必须相等某label('n');%横坐标显示nylabal('某(n)');%纵坐标显示某(n)grid;%绘制网格1(2)n=0:9;某=0.5.^n;tem(n,某);某label('n');ylabel('某(n)');gri实验目的1学习MATLAB语言编程和调试技巧2学会简单的矩阵输入和图形表示法3掌握简单的绘图命令2、用MATLAB计算序列{-201–13}和序列{120-1}的离散卷积，即计算某(n)2(n)(n2)(n3)3(n4)与h(n)(n)2(n1)(n3)解：用MATLAB描述波形。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
n0=0;n1=-3;n2=3; impseq(n0,n1,n2); 1.1.2.单位阶跃序列 u(n)
源程序：
function [x,n] =stepseq(n0,n1,n2) %产生x(n)=u(n-n0);n1<n0<n2 if((n0<n1)|(n0>n2)|(n1>n2)) error('参数不满足n1<n0<n2') end n=n1:n2
锯齿波 源程序： function [y,t] = saw( t ) t=0:0.001:10; y=t-fix(t);%取 t 的小数部分 plot(t,y)
2. 总结 matlab 编程和 C 编程的相同点和不同点？ （1）相同点：都是可编程软件；很多编程语句和规则都是相同的（如赋值、变量命名、循 环语句、条件语句、函数嵌套调用等）；都是实现结构化程序块编程。 （2）不同点：MATLAB 是一种数值计算环境和编程语言，它可以进行矩阵运算、绘制函数和 数据、实现算法、创建用户界面、连接其他编程语言的程序等，主要应用于工程计算、控制 设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。Matlab 中很 多函数可以直接使用而非定义声明之后调用。C 语言，是一种通用的、程序式的编程语言， 广泛用于系统与应用软件的开发。具有高效、灵活、功能丰富、表达力强和较高的移植性等 特点，也就是说，C 语言的目的是用与软件开发，它有着丰富的应用支持。C 语言中调用函 数必须事先自定义和声明。 3. 总结叙述 matlab 编程的一个新应用领域，工具箱或功能 目前，MATLAB 已经把工具箱延伸到了科学研究和工程应用的诸多领域，如数据采集、数据 库接口、概率统计、样条拟合、优化算法、偏微分方程求解、神经网络、小波分析、信号处 理、图像处理、系统辨识、控制系统设计、线性矩阵不等式控制、鲁棒控制、模型预测、模 糊逻辑、金融分析、地图工具、非线性控制设计、实时快速原型及半物理仿真、嵌入式系统 开发、定点仿真、DSP 与通讯、电力系统仿真等领域。 模型预测控制工具箱提供了一系列用于模型预测控制的分析、 设计和仿真的函数。 这些函数 的类型主要有： ① 系统模型辨识函数： 主要功能包括通过多变量线性回归方法计算 MISO 脉冲响应模型和 阶跃响应模型及测量数据的归一化等。 ② 模型建立和转换函数：主要功能包括建立模型预测控制工具箱使用的 MPC 状态空间模 型 及状态空间模型与 MPC 状态空间模型、 阶跃响应模型、 脉冲响应模型之间的转换。 ③ 模型预测控制器设计和仿真工具主要功能包括面向阶跃响应模型的预测控制器设计 与仿真函数和面向 MPC 状态空间模型的设计和仿真函数两类。 ④ 系统分析工具主要功能包括计算模型预测控制系统频率响应、极点和奇异值的有关函 数。 ⑤ 其他功能函数主要功能包括绘图和矩阵计算函数等。
源程序： function [x,t] = X2( A,a,f,t1,t2 ) %产生函数 x(t)=e^(-at) sin(2pift),t1<t2 if(t1>t2) error('参数不满足 t1<t2') end t=t1:0.01:t2; x=exp(-a.*t).* sin(2.*pi.*f.*t); plot(t,x) A=2;a=0.5;f=2;t1=-2;t2=2; x=X2(A,a,f,t1,t2) 1.4 谐波信号
（4）求和：将以上所有对应点的乘积累加起来，即得y(n)。
三、主要实验仪器及材料
PC 机、Matlab7.0。
四、实验内容
1．知识准备 认真复习以上基础理论，理解本实验所用到的实验原理。 2．离散时间信号（序列）的产生 利用 MATLAB 产生和绘制下列有限长序列： （1） 单位脉冲序列δ (n)，单位阶跃序列 u(n)，矩形序列 R8(n) （2） 在一幅图上绘出曲线 x1(t)=Asin(5 t + 3 );x2(t)=Acos(t-3 );x3(t)=2sin(3t)cos(2t-3 ). （3） 绘出正弦衰减信号：x2(t)=Ae−α t sin 2πft , A = 2, α = 0.5, f = 2Hz （ 4 ） 绘 出 谐 波 信 号 ： x3(t)= A 1sin(2 π ft)+A2sin(2 π f2t)+A3sin(2 π f3t);A1=1,A2=0.5,A3=0.2,f=2Hz. 3．序列的运算 生成下列序列： x(n)=δ (n)+2δ (n-1)+3δ (n-2)+4δ (n-3)+5δ (n-4) h(n)= δ (n)+2δ (n-1)+δ (n-2)+2δ (n-3) (1) 利用 MATLAB 编程完成上述两序列的卷积，并绘制运算后序列的波形。 （2）编写一个完成卷积的函数，输入是两个序列，输出是结果。 4．采样定理的研究 分别令采样周期 Ts 为不同值，绘出不同采样周期下 x(t) =sin(t)的频谱（直接用 FFT 函数） ， 观察频谱混叠现象，和防止混叠的采样周期。 5．上机实验并完成实验报告。请将实验报告保存为压缩文件（.rar， 或.zip 格 式） ，包括源程序（.M 文件）和文档（.doc） ，统一名为“xxx——实验 x.rar”或 “xxx——实验 x.zip”
N=8;ns=-2;nf=10; x=R(N,ns,nf) 1.2 在同一图中绘制 x1(t),x2(t),x3(t)
源程序： function [x1,x2,x3,t ] = X1(A,a,b )
t=a:0.001:b;a<b if(a>=b) error('参数不满足 a<b') end x1=A.*sin(pi./5.*t+ pi./3); x2=A.*cos(t-pi./3); x3=2.*sin(3.*t).*cos(2.*t-pi./3); plot(t,x1,x2,x3) A=2;a=-pi;b=pi; [x1,x2,x3]=X1(A,a,b)； 1.3 正弦衰减信号
(2)fs=π /2
(3)fs=2
(4)fs=10
Байду номын сангаас
由图可以看出：当采样周期 T<=π /2 时正弦函数 x=sin(t)不产生混叠。
源程序： function [ ] =sample(fs,N) T=1/fs; %采样周期 n=0:N-1; t=n*T; %采样时间点 x=sin(t); %生成正弦信号 subplot(1,2,1) plot(t,x);%画正弦信号的时域波形 y=fft(x,N); %fft 变换 mag=sqrt(y.*conj(y)); %求幅值 f=(0:N-1)*fs/N; %进行对应的频率转换 subplot(1,2,2); plot(f,mag); %画频谱图 axis([-10,50,0,80]);
π π π π
五、实验结果
1.1.1 单位脉冲序列δ （n）
源程序： function [x,n] =impseq(n0,n1,n2) %产生 x(n)=δ (n-n0);n1<=n0<=n2 if((n0<n1)|(n0>n2)|(n1>n2)) error('参数不满足 n1<=n0<=n2') end n=n1:n2; x=[(n-n0)==0]; %x=[zeros(1,n0-n1),1,zeros(1,(n2-n0))] stem(n,x);
源程序： function [ x,t ] =X3( A1,A2,A3,f,t1,t2 ) %产生 x3(t)=A1.*sin(2.*pi.*f.*t)+A2.*sin(2.*pi.*f.*2.*t)+A3.*sin(2.*pi.*f.*3.* %t),t1<t2 if(t1>t2) error('参数不满足 t1<t2') end t=t1:0.0001:t2; x=A1.*sin(2.*pi.*f.*t)+A2.*sin(2.*pi.*f.*2.*t)+A3.*sin(2.*pi.*f.*3.*t); plot(t,x) A1=1;A2=0.5; A3=0.2;f=2;t1=-1;t2=1; x=X3(A1,A2,A3,f,t1,t2)； 2 矩阵的运算
x=[(n-n0)>=0] %x=[zeros(1,(n0-n1)),ones(1,(n2-n1))] stem(n,x);
n0=0;n1=-3;n2=3; stepseq(n0,n1,n2); 1.1.3.矩形序列
源程序： function [ x,n ] = R( N,ns,nf ) %产生序列x(n)=u(n-n0)-u(n-n0-N)n=ns:nf;np=0; x=stepseq(np,ns,nf)-stepseq(N,ns,nf); %x=[0<=n<7] stem(n,x);
源程序： function [y,ny] = convwthn(x,nx,h,nh) ny1=nx(1)+nh(1); ny2=nx(end)+nh(end)
%x*h 在 nx(1)+nh(1)到 nx(end)+nh(end)之间取值； %nx(1)、nh(1)为最小下标，end 为最大下标； y=conv(x,h); ny=ny1:ny2; stem(ny,y); end x=[1,2,3,4,5]; nx=0:4; h=[1,2,1,2]; nh=0:3; [y,ny]=convwthn(x,nx,h,nh); 3 采样定理的研究 （1） fs=0.5
五、思考题
1. 如何产生方波信号序列和锯齿波信号序列？ 方波 源程序： function [y,t ] = squ( T ) %方波 T=0:0.001:2*pi;%方波的时间向量 y=square(6*T);%方波 w=6,周期为 T=2*pi/6 plot(T,y);%画方波； axis([T(1)-1,T(end)+1,-2,2]);%坐标轴区域