浙江大学 信电系 信号与系统实验报告
《信号与系统》课程实验报告
《信号与系统》课程实验报告《信号与系统》课程实验报告一图1-1 向量表示法仿真图形2．符号运算表示法若一个连续时间信号可用一个符号表达式来表示，则可用ezplot命令来画出该信号的时域波形。
上例可用下面的命令来实现（在命令窗口中输入，每行结束按回车键）。
t=-10:0.5:10;f=sym('sin((pi/4)*t)');ezplot(f,[-16,16]);仿真图形如下：图1-2 符号运算表示法仿真图形三、实验内容利用MATLAB实现信号的时域表示。
三、实验步骤该仿真提供了7种典型连续时间信号。
用鼠标点击图0-3目录界面中的“仿真一”按钮，进入图1-3。
图1-3 “信号的时域表示”仿真界面图1-3所示的是“信号的时域表示”仿真界面。
界面的主体分为两部分：1) 两个轴组成的坐标平面（横轴是时间，纵轴是信号值）；2) 界面右侧的控制框。
控制框里主要有波形选择按钮和“返回目录”按钮，点击各波形选择按钮可选择波形，点击“返回目录”按钮可直接回到目录界面。
图1-4 峰值为8V，频率为0.5Hz，相位为180°的正弦信号图1-4所示的是正弦波的参数设置及显示界面。
在这个界面内提供了三个滑动条，改变滑块的位置，滑块上方实时显示滑块位置代表的数值，对应正弦波的三个参数：幅度、频率、相位；坐标平面内实时地显示随参数变化后的波形。
在七种信号中，除抽样函数信号外，对其它六种波形均提供了参数设置。
矩形波信号、指数函数信号、斜坡信号、阶跃信号、锯齿波信号和抽样函数信号的波形分别如图1-5～图1-10所示。
图1-5 峰值为8V，频率为1Hz，占空比为50%的矩形波信号图1-6 衰减指数为2的指数函数信号图1-7 斜率=1的斜坡信号图1-8 幅度为5V，滞后时间为5秒的阶跃信号图1-9 峰值为8V，频率为0.5Hz的锯齿波信号图1-10 抽样函数信号仿真途中，通过对滑动块的控制修改信号的幅度、频率、相位，观察波形的变化。
信号与系统测试实训报告
信号与系统测试实训报告信号与系统测试实训报告一、引言在信号与系统的学习过程中，理论知识的掌握是非常重要的，但是对于学生来说，实际操作和实验测试同样重要。
本报告旨在总结信号与系统测试实训的过程、方法和结果，并对实训中遇到的问题进行分析和解决。
二、实验目的本次实验主要目的是通过使用各种仪器设备和软件工具，对不同类型的信号进行测试和分析，验证信号与系统理论知识，并培养学生动手能力和问题解决能力。
三、实验设备和软件1. 实验设备：- 函数发生器：用于产生各种类型的基本信号。
- 示波器：用于观测信号波形和频谱。
- 模拟滤波器：用于对信号进行滤波处理。
- 数字示波器：用于观测高频数字信号。
- 电脑：用于运行软件工具。
2. 软件工具：- MATLAB：用于信号处理与分析。
- Simulink：用于模拟系统响应。
四、实验步骤及结果1. 实验一：基本周期信号测试首先使用函数发生器产生正弦、方波、三角波等基本周期信号，然后用示波器观测信号波形，并通过MATLAB进行频谱分析。
实验结果显示，正弦信号的频谱主要集中在一个频率上，而方波和三角波的频谱则包含了多个谐波分量。
2. 实验二：连续时间系统测试使用函数发生器产生输入信号，经过模拟滤波器后，用示波器观测输出信号，并通过MATLAB进行系统响应分析。
实验结果显示，滤波器对输入信号进行了平滑处理，并且系统响应与理论预期相符。
3. 实验三：离散时间系统测试使用函数发生器产生离散时间序列输入信号，经过数字滤波器后，用数字示波器观测输出信号，并通过MATLAB进行系统响应分析。
实验结果显示，数字滤波器对输入信号进行了去噪处理，并且系统响应与理论预期相符。
五、实验问题与解决1. 问题一：示波器观测到的信号与理论预期不符。
解决方法：首先检查示波器的设置是否正确，并确保连接线路无误。
检查函数发生器的输出是否稳定和准确。
检查示波器通道的控制参数是否正确设置。
2. 问题二：MATLAB频谱分析结果不准确。
信号与系统实验实验报告
信号与系统实验实验报告一、实验目的本次信号与系统实验的主要目的是通过实际操作和观察，深入理解信号与系统的基本概念、原理和分析方法。
具体而言，包括以下几个方面：1、掌握常见信号的产生和表示方法，如正弦信号、方波信号、脉冲信号等。
2、熟悉线性时不变系统的特性，如叠加性、时不变性等，并通过实验进行验证。
3、学会使用基本的信号处理工具和仪器，如示波器、信号发生器等，进行信号的观测和分析。
4、理解卷积运算在信号处理中的作用，并通过实验计算和观察卷积结果。
二、实验设备1、信号发生器：用于产生各种类型的信号，如正弦波、方波、脉冲等。
2、示波器：用于观测输入和输出信号的波形、幅度、频率等参数。
3、计算机及相关软件：用于进行数据处理和分析。
三、实验原理1、信号的分类信号可以分为连续时间信号和离散时间信号。
连续时间信号在时间上是连续的，其数学表示通常为函数形式；离散时间信号在时间上是离散的，通常用序列来表示。
常见的信号类型包括正弦信号、方波信号、脉冲信号等。
2、线性时不变系统线性时不变系统具有叠加性和时不变性。
叠加性意味着多个输入信号的线性组合产生的输出等于各个输入单独作用产生的输出的线性组合；时不变性表示系统的特性不随时间变化，即输入信号的时移对应输出信号的相同时移。
3、卷积运算卷积是信号处理中一种重要的运算，用于描述线性时不变系统对输入信号的作用。
对于两个信号 f(t) 和 g(t)，它们的卷积定义为：＼（f g)（t) ＝＼int_{＼infty}＾｛＼infty} f(＼tau) g(t ＼tau) d\tau ＼在离散时间情况下，卷积运算为：＼（f g)n ＝＼sum_{m ＝＼infty}＾｛＼infty} fm gn m ＼四、实验内容及步骤实验一：常见信号的产生与观测1、连接信号发生器和示波器。
2、设置信号发生器分别产生正弦波、方波和脉冲信号，调整频率、幅度和占空比等参数。
3、在示波器上观察并记录不同信号的波形、频率和幅度。
信号与系统实验报告
信号与系统实验报告一、实验目的(1) 理解周期信号的傅里叶分解，掌握傅里叶系数的计算方法；(2)深刻理解和掌握非周期信号的傅里叶变换及其计算方法；(3) 熟悉傅里叶变换的性质，并能应用其性质实现信号的幅度调制；(4) 理解连续时间系统的频域分析原理和方法，掌握连续系统的频率响应求解方法，并画出相应的幅频、相频响应曲线。
二、实验原理、原理图及电路图(1) 周期信号的傅里叶分解设有连续时间周期信号()f t ，它的周期为T ，角频率22fT，且满足狄里赫利条件，则该周期信号可以展开成傅里叶级数，即可表示为一系列不同频率的正弦或复指数信号之和。
傅里叶级数有三角形式和指数形式两种。
1）三角形式的傅里叶级数：01212011()cos()cos(2)sin()sin(2)2cos()sin()2n n n n a f t a t a t b t b t a a n t b n t 式中系数n a ，n b 称为傅里叶系数，可由下式求得：222222()cos(),()sin()T T T T nna f t n t dtb f t n t dtTT2）指数形式的傅里叶级数：()jn tn nf t F e式中系数n F 称为傅里叶复系数，可由下式求得：221()T jn tT nF f t edtT周期信号的傅里叶分解用Matlab进行计算时，本质上是对信号进行数值积分运算。
Matlab中进行数值积分运算的函数有quad函数和int函数。
其中int函数主要用于符号运算，而quad函数（包括quad8，quadl）可以直接对信号进行积分运算。
因此利用Matlab进行周期信号的傅里叶分解可以直接对信号进行运算，也可以采用符号运算方法。
quadl函数(quad系)的调用形式为：y＝quadl(‘func’,a,b)或y＝quadl(@myfun,a,b)。
其中func是一个字符串，表示被积函数的.m文件名（函数名）；a、b分别表示定积分的下限和上限。
浙江大学 信号与系统实验-软硬件结合操作实验
CopyrightAs one member of Information Science and Electronic Engineering Institute of Zhejiang University, I sincerely hope this will enable you to acquire more time to do whatever you like instead of struggling on useless homework. All the content you can use as you like. I wish you will have a meaningful journey on your college life.——W z h实验报告课程名称：信号与系统实验指导老师：周绮敏、史笑兴、李惠忠成绩：__ ___ 实验名称：MATLAB操作实验实验类型：电子电路实验名称：软硬件结合操作实验1——信号的卷积一、实验目的1、理解卷积的概念及物理意义。
2、通过实验的方法加深对卷积运算的图解方法及结果的理解。
二、实验仪器1、双踪示波器 1台2、信号源及频率计模块 1块3、数字信号处理模块 1块三、实验内容（一）脉冲信号的自卷积1）观察记录频率500Hz、幅度Vpp=1V、占空比为的矩形波自卷积2）减小矩形波的占空比，观察卷积结果；并记录占空比为时的矩形波自卷积3）分别观察记录频率500Hz、幅度Vpp=2V、占空比为、的矩形波自卷积（二）矩形信号与锯齿波信号的互卷积1）频率500Hz、幅度Vpp=1V、占空比为50%、25%的矩形波与幅度2V、频率500Hz锯齿波信号互卷积2）频率500Hz、幅度Vpp=2V、占空比为50%、25%的矩形波与幅度2V、频率500Hz锯齿波信号互卷积四、实验结果及分析研究（二）矩形信号与锯齿波信号的互卷积实验名称：软硬件结合操作实验2 __信号的分解及合成一、实验目的1、了解波形分解与合成原理，掌握用傅里叶级数进行谐波分析的方法。
信号与系统实验报告
信号与系统实验报告
实验名称：信号与系统实验
一、实验目的：
1.了解信号与系统的基本概念
2.掌握信号的时域和频域表示方法
3.熟悉常见信号的特性及其对系统的影响
二、实验内容：
1.利用函数发生器产生不同频率的正弦信号，并通过示波器观察其时域和频域表示。
2.通过软件工具绘制不同信号的时域和频域图像。
3.利用滤波器对正弦信号进行滤波操作，并通过示波器观察滤波前后信号的变化。
三、实验结果分析：
1.通过实验仪器观察正弦信号的时域表示，可以看出信号的振幅、频率和相位信息。
2.通过实验仪器观察正弦信号的频域表示，可以看出信号的频率成分和幅度。
3.利用软件工具绘制信号的时域和频域图像，可以更直观地分析信号的特性。
4.经过滤波器处理的信号，可以通过示波器观察到滤波前后的信号波形和频谱的差异。
四、实验总结：
通过本次实验，我对信号与系统的概念有了更深入的理解，掌
握了信号的时域和频域表示方法。
通过观察实验仪器和绘制图像，我能够分析信号的特性及其对系统的影响。
此外，通过滤波器的处理，我也了解了滤波对信号的影响。
通过实验，我对信号与系统的理论知识有了更加直观的了解和应用。
《信号与系统》课程实验报告材料
工业大学校区《信号与系统》课程实验报告专业班级学生《信号与系统》课程实验报告一实验名称一阶系统的阶跃响应姓名系院专业班级学号实验日期指导教师成绩一、实验目的1．熟悉一阶系统的无源和有源电路；2．研究一阶系统时间常数T的变化对系统性能的影响；3．研究一阶系统的零点对系统响应的影响。
二、实验原理1．无零点的一阶系统无零点一阶系统的有源和无源电路图如图2-1的(a)和(b)所示。
它们的传递函数均为：10.2s1G(s)＝+(a) 有源(b) 无源图2-1 无零点一阶系统有源、无源电路图2．有零点的一阶系统(|Z|<|P|)图2-2的(a)和(b)分别为有零点一阶系统的有源和无源电路图，它们的传递函数为：10.2s1)0.2(sG(s)++=，⎪⎪⎪⎪⎭⎫⎝⎛++=S611S161G(s）(a) 有源(b) 无源图2-2 有零点(|Z|<|P|)一阶系统有源、无源电路图3．有零点的一阶系统(|Z|>|P|)图2-3的(a)和(b)分别为有零点一阶系统的有源和无源电路图，它们的传递函数为：1s10.1sG(s)＝++(a) 有源(b) 无源图2-3 有零点(|Z|>|P|)一阶系统有源、无源电路图三、实验步骤1．打开THKSS-A/B/C/D/E型信号与系统实验箱，将实验模块SS02插入实验箱的固定孔中，利用该模块上的单元组成图2-1(a)(或(b))所示的一阶系统模拟电路。
2．实验线路检查无误后，打开实验箱右侧总电源开关。
3．将“阶跃信号发生器”的输出拨到“正输出”，按下“阶跃按键”按钮，调节电位器RP1，使之输出电压幅值为1V，并将“阶跃信号发生器”的“输出”端与电路的输入端“Ui”相连，电路的输出端“Uo”接到双踪示波器的输入端，然后用示波器观测系统的阶跃响应，并由曲线实测一阶系统的时间常数T。
4．再依次利用实验模块上相关的单元分别组成图2-2(a)(或(b))、2-3(a)(或(b))所示的一阶系统模拟电路，重复实验步骤3，观察并记录实验曲线。
浙江大学 信号与系统实验-拓展实验
CopyrightAs one member of Information Science and Electronic Engineering Institute of Zhejiang University, I sincerely hope this will enable you to acquire more time to do whatever you like instead of struggling on useless homework. All the content you can use as you like. I wish you will have a meaningful journey on your college life.——W z h实验报告课程名称：信号与系统实验指导老师：史笑兴、周绮敏、李惠忠成绩：__ ___实验名称：MATLAB拓展与应用实验实验类型：仿真设计（一、二为拓展实验，至少选做一题或都做；三、四为应用实验，至少选做一题或都做。
）一、周期信号的傅立叶级数分析与谐波合成1、Matlab代码与运行结果1.1 任务一产生周期为1，占空比为30%，高电平为1，低电平为0的周期方波，时间范围为[−2,2]，即能观察到4个完整的周期。
（时间变量的尺度变换）【代码】【结果】1.2 任务二产生一个周期方波，并计算其傅立叶系数，画出频谱分析图。
此任务可以分成三段程序来完成，分别是：①产生时域信号；②计算傅立叶系数；③绘制频谱分析图。
【代码】【结果】1.3 任务三周期方波的谐波合成。
（1）产生时域信号：产生一个方波信号，周期为1，高电平为1，低电平为0，占空比为50%，时间范围为[0,4]（此方波信号与【任务2】中周期方波相同，但是时间范围较长，包含4个周期）。
（2）根据【任务2】的计算结果，在[0,4] 的时间范围内，分别画出原始的方波信号以及以下各信号的波形进行对比。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
本科实验报告课程名称：信号与系统实验姓名：院系：专业：学号：指导教师：周绮敏、史笑兴、李惠忠年月日实验报告课程名称： 信号与系统实验 指导老师： 史笑兴、周绮敏、李惠忠 成绩： 实验名称： 实验一 MATLAB 基本实验 实验类型： 设计型 一、实验目的与要求1、掌握Matlab 的基本使用方法：操作界面、基本语法、基本函数等；2、掌握Matlab 在“信号与系统”课程中的应用： （1）信号的产生与表示；（2）产生基本信号的Matlab 函数； （3）卷积和的计算与表示；（4）离散时间信号的傅立叶变换的计算与表示，部分傅立叶变换的性质验证。
二、实验内容和原理1、P168，实验一，习题 1、2、3、4、5、7 2、P169，实验二，习题 1、1的补充题、2、3 3、P170，实验三，习题 1、2、3三、主要仪器设备带有matlab 的电脑一台四、操作方法和实验步骤4.1 新建文本文档并命名，在文本文档中写matlab 代码，另存为.m 的文件，此时注意名称开头不应为数字，标题中不应出现matlab 自带函数的名称； 4.2打开matlab 程序，运行已保存好的.m 文件，根据提示的错误信息修改代码直至无错误为止，运行代码，记录运行结果。
五、实验数据记录和处理 实验1.1 代码：clear;clf;clc; %Impulse sequence subplot(2,2,1); n=[-10:10]; z=(n==0); stem(n,z);axis([-10,10,-1,2]); title('Impulse sequence'); xlabel('n');ylabel('delta[n]');%Impulse function subplot(2,2,2); t=-10:0.1:10;专业： 姓名： 学号：日期：5.3—5.24 地点：东四420室实验时间： 周 第 节课 签到表编号：y=(t==0);plot(t,y,'r');axis([-10,10,-1,2]);title('Impulse function');xlabel('t');ylabel('\delta(t)');%Step Sequencesubplot(2,2,3);n=[-10:10];z=(n>=0);stem(n,z);axis([-10,10,-1,2]);title('Step sequence');xlabel('n');ylabel('u[n]');%Step functionsubplot(2,2,4);t=-10:0.1:10;y=(t>=0);plot(t,y,'b');axis([-10,10,-1,2]);title('Step function');xlabel('t');ylabel('u(t)');结果：实验1.2代码：%sine wavesubplot(2,2,1);t=-2*pi:pi/20:2*pi;plot(t,sin(t));axis([-2*pi,2*pi,-1.5,1.5]); title('Sine wave');xlabel('t');ylabel('sin(t)');%Triangular wavesubplot(2,2,2);t=-5*pi:pi/100:5*pi;x=sawtooth(t,0.5);plot(t,x);axis([-5*pi,5*pi,-1.5,1.5]); title('Triangular wave'); xlabel('t');ylabel('x');%Sawtooth wavesubplot(2,2,3);t=-5*pi:pi/100:5*pi;z=sawtooth(t,0);plot(t,z);axis([-5*pi,5*pi,-1.5,1.5]); title('Sawtooth wave t=0'); xlabel('t');ylabel('z');subplot(2,2,4);t=-5*pi:pi/100:5*pi;z=sawtooth(t,1);plot(t,z);axis([-5*pi,5*pi,-1.5,1.5]); title('Sawtooth wave t=1'); xlabel('t');ylabel('z');结果：实验1.3代码：clear;clf;clc;%Sinc functiont=-4*pi:pi/20:4*pi;subplot(2,1,1);plot(t,sinc(t));axis([-4*pi,4*pi,-0.5,1.5]);title('Sinc function');grid on;xlabel('t');ylabel('sinc(t)');%Sample functionsubplot(2,1,2);plot(t,sinc(t/pi));axis([-4*pi,4*pi,-0.5,1.5]);title('Sample function');grid on;xlabel('t');ylabel('sa(t)');结果：实验1.4代码：u1=stepseq(0,-5,20);u3=stepseq(5,-5,20);x=u1-u3;n=[-5:20];h=delta(n-5);subplot(3,1,1);stem(n,x);axis([-5,20,0,1.2]);title('x[n]'); xlabel('n');ylabel('x[n]');subplot(3,1,2);stem(n,h);axis([-5,20,0,1.2]);title('h[n]');xlabel('n');ylabel('h[n]');[y,ny]=conv_my(x,n,h,n)subplot(3,1,3);stem(ny,y);axis([-5,20,0,1.2]);title('y[n]');xlabel('n');ylabel('y[n]');结果：实验1.5代码：n=[-10:10];x1=delta(n);x2=delta(n+2);x3=delta(n-4);x4=2*delta(n+2)-delta(n-4);subplot(2,2,1);stem(n,x1);title('x_1[n]=\delta[n]');axis([-10,10,-1,2]);grid on;subplot(2,2,2);stem(n,x2);title('x_2[n]=\delta[n+2]');axis([-10,10,-1,2]);grid on;subplot(2,2,3);stem(n,x3);title('x_3[n]=\delta[n-4]');axis([-10,10,-1,2]);grid on;subplot(2,2,4);stem(n,x4);title('x_4[n]=2*\delta[n+2]-delta[n-4]'); axis([-10,10,-1,2]);grid on;结果：实验1.7代码：clear;clf;clc;n=[-15:15];x1=stepseq(0,-15,15);x2=stepseq(10,-15,15);x=x1-x2;[xe,xo]=evenodd(x);subplot(2,2,1);stem(n,x);title('Step sequence');xlabel('n');ylabel('x[n]');axis([-15.5,15.5,-1.2,2]);subplot(2,2,3);stem(n,xe);axis([-15,15,-1,2]);title('Even part'); xlabel('n');ylabel('xe[n]');axis([-15.5,15.5,-1.2,2]);subplot(2,2,4);stem(n,xo);title('Odd part');xlabel('n');ylabel('xe[n]');axis([-15.5,15,-1.2,1.2]);结果：实验2.1代码：clear;clf;clc;n=[-5:10];nx=[-5:10];nh=[-5:10];h=delta(n-5)*1.0;x=stepseq(0,-5,10)-stepseq(5,-5,10);[y,ny]=conv_m(x,nx,h,nh);subplot(3,1,1);stem(n,x);title('x[n]');xlabel('n');axis([-5,10,0,1.2] );grid on;subplot(3,1,2);stem(n,h);title('h[n]');xlabel('n');axis([-5,10,0,1.2] );grid on;subplot(3,1,3);stem(ny,y);title('y[n]');xlabel('n');axis([-10,20,0,1. 2]);grid on;结果：实验2.1补充题代码：u1=stepseq(5,-10,10);u3=stepseq(-5,-10,10);x=(abs(n)<=5)*1.0;n=[-10:10];nx=[-10:10];nh=[-10:10];h=delta(n+3)-delta(n-3);[y,ny]=conv_m(x,nx,h,nh);subplot(3,1,1);stem(n,x);axis([-10,10,-0.5,1.5]);title('x[n]'); xlabel('n');ylabel('x[n]');subplot(3,1,2);stem(n,h);axis([-10,10,-1.5,1.5]);title('h[n]');xlabel('n');ylabel('h[n]');subplot(3,1,3);stem(ny,y);axis([-20,20,-1.5,1.5]);title('y[n]'); xlabel('n');ylabel('y[n]');结果：实验2.2a 代码：n=-5:20;u1=stepseq(0,-5,20);u2=stepseq(6,-5,20);x=u1-u2;h=n.*x;subplot(3,1,1);stem(n,x);axis([-5,20,0,2]);title('Input Sequence');ylabel('x[n]');subplot(3,1,2);stem(n,h);axis([-5,20,0,6]);title('Impulse Sequence');ylabel('x[n]');[y1,ny]=conv_m(x,n,h,n);subplot(3,1,3);stem(ny,y1);title('Output Sequence');xlabel('n');ylabel('y_1[n]');结果：b 代码：n=-10:20;u1=stepseq(0,-10,20);u2=stepseq(6,-10,20);x=u1-u2;u3=stepseq(-5,-10,20);u4=stepseq(1,-10,20);x1=u3-u4;h=(n+5).*x1;subplot(3,1,1);stem(n,x);axis([-10,20,0,2]);title('Input Sequence');ylabel('x[n]');subplot(3,1,2);stem(n,h);axis([-10,20,0,6]);title('Impulse Response');ylabel('h[n+5]');[y2,ny]=conv_m(x,n,h,n);subplot(3,1,3);stem(ny,y2);title('OutputSequence');xlabel('n');ylabel('y_2[n]');结果：实验2.3代码：n=-5:50;u1=stepseq(0,-5,50);u2=stepseq(10,-5,50);x=u1-u2;h=(0.9).^n.*(n>=0);subplot(3,1,1);stem(n,x);axis([-5,20,-1,2]);title('x[n]'); ylabel('x[n]');subplot(3,1,2);stem(n,h);axis([-5,50,-1,2]);title('h[n]'); ylabel('h[n]');[y,ny]=conv_m(x,n,h,n);subplot(3,1,3);stem(ny,y);axis([-10,70,-2,8]);title('y[n]'); xlabel('n');ylabel('y[n]');结果：实验3.1代码：clear;clc;clf;n=[-20:20];x=stepseq(-5,-20,20)-stepseq(6,-20,20);w1=-2*pi;w2=2*pi;w_step=pi/400;w=w1:w_step:w2;ww=w/pi;xw=spec(x,n,w); %wsubplot(2,2,1);plot(ww,abs(xw));axis([-2,2,0,max(abs(xw))]);grid on; xlabel('w');title('Abs(X)');subplot(2,2,2);plot(ww,angle(xw));axis([-1.2,1.2,-1.2*pi,1.2*pi]); grid on;xlabel('w');title('Phase(X)');subplot(2,2,3);plot(ww,real(xw));axis([-2,2,min(real(xw)),max(real(xw))]);grid on;xlabel('w');title('Real(X)');subplot(2,2,4);plot(ww,imag(xw));axis([-2,2,min(imag(xw)),max(imag(xw))]);grid on;xlabel('w');title('Imag(X)');结果：实验3.2代码：n=0:100;x=cos(pi*n/2);py=exp(j*pi*n/4);y=py.*x;w = -pi : pi/800 : pi;ww=w/pi;xw=spec(x,n,w);yw=spec(y,n,w);subplot(2,2,1);plot(ww,abs(xw));title('Magnitude of X');grid; ylabel('|X|');xlabel('Normalized Freq.(¦Ð rad/s)');axis([ min(ww) , max(ww) , 0, max(abs(xw)) ]);subplot(2,2,2);plot(ww,angle(xw));title('Angle of X');grid; ylabel('Angle(X)');xlabel('Normalized Freq.(¦Ð rad/s)');axis([ min(ww) , max(ww) , -pi, pi ]);subplot(2,2,3);plot(ww,abs(yw));title('Magnitude of Y');grid; ylabel('|Y|');xlabel('Normalized Freq.(¦Ð rad/s)');axis([ min(ww) , max(ww) , 0, max(abs(yw)) ]);subplot(2,2,4);plot(ww,angle(yw));title('Angle of Y');grid; ylabel('Angle(Y)');xlabel('Normalized Freq.(¦Ð rad/s)');axis([ min(ww) , max(ww) , -pi, pi ]);结果：实验3.3代码：n=-5:10;x=sin(pi*n/2);w1=-pi;w2=pi;w_step=pi/800;w=w1:w_step:w2;X=spec(x,n,w);[xe,xo,m]=evenodd2(x,n);clf;XE=spec(xe,m,w);XO=spec(xo,m,w);XR=real(X);XI=imag(X);subplot(3,2,1);plot(w/pi,XR,'r');grid;axis([-1,1,min(XR),max(XR)]); xlabel('Normalized Freq.(\pi rad/sample)');ylabel('Re[X]');title('Real part of X');subplot(3,2,2);plot(w/pi,XI,'r');grid;axis([-1,1,min(XI),max(XI)]); xlabel('Normalized Freq.(\pirad/sample)');ylabel('Im(X)');title('Imaginary part of X');subplot(3,2,3);plot(w/pi,real(XE),'r');grid;axis([-1,1,min(real(XE)),max(real(XE))]);xlabel('Normalized Freq.(\pi rad/sample)');ylabel('XE or Re(XE)');title('DFT of even part-real');subplot(3,2,4);plot(w/pi,imag(XE),'r');grid;axis([-1,1,min(imag(XE)),max(imag(XE))]);xlabel('Normalized Freq.(\pi rad/sample)');ylabel('Im(XE)');title('DFT of even part-Imaginary');subplot(3,2,5);plot(w/pi,real(XO),'r');grid;axis([-1,1,min(real(XO)),max(real(XO))]);xlabel('Normalized Freq.(\pi rad/sample)');ylabel('Re(XO)');title('DFT of oddd part-Real');subplot(3,2,6);plot(w/pi,imag(XO),'r');grid;axis([-1,1,min(imag(XO)),max(imag(XO))]);xlabel('Normalized Freq.(\pi rad/sample)');ylabel('XO/j or Im(XO)');title('DFT of odd part-Imaginary');结果：子函数：1.stepseq.mfunction [x,n]=stepseq(n0,n1,n2)n=n1:n2;x=(n>=n0);函数的功能是产生阶跃函数，其中n1和n2是函数的横坐标边界，n0是跳变点。