当前位置：文档之家 › 信号处理与基础 实验报告

信号处理与基础 实验报告

  • 格式：doc
  • 大小：169.00 KB
  • 文档页数：16

河南城建学院实验报告课程名称：信号处理基础系：电气与电子工程系专业：自动化指导教师：梁成武姓名：学号：报告上交时间：2210 年12 月23 日教师评语成绩等级日期：2010年12月日实验一：MATLAB工具使用和信号的时域析一：实验目的1：对Matlab有初步认识2:掌握用Matlab软件产生基本信号的方法。

3：初步掌握用matlab对信号进行时域分析。

二：实验设备1：电脑一台2：Matlab软件三：实验原理与步骤MATLAB是矩阵实验室（Matrix Laboratory）的简称，是美国MathWorks公司出品的商业数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括MATLAB和Simulink两大部分。

MATLAB和Mathematica、Maple 并称为三大数学软件。

它在数学类科技应用软件中在数值计算方面首屈一指。

MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连matlab开发工作界面接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。

连续信号是指自变量的取值范围是连续的，且对于一切自变量的取值除了若干个不连续点之外，信号都有确定的取值与之对应。

严格来说，MATLAB并不能处理连续信号，而是用等时间间隔的样值来近似表示连续信号，当取样时间间隔足够小时，这些离散的样值就能较好的祭祀连续信号。

利用MA TLAB软件的信号处理工具箱中的函数可产生常见信号的函数如下表：1：产生正弦函数程序：t=（0：0.001:50）；y=sin(2*pi*50*y);plot(t（1：50)，y（1:50））2：产生sinc函数程序：x=linspace(-5,5);y=sinc（x）：plot（x，y）3：产生指数函数信号程序：x=(0:0.01:1);y=exp(-x);plot(x,y)在matlab中输入相应程序就可以得到相应的波形四：实验结果1：正弦函数图形2：sinc函数图形3：指数函数图形五：试验体会通过本次试验对MA TLAB有了初步的认识，了解了MA TLAB作图的基本方法。

并掌握了时域分析方法。

深刻体会到了作图和对信号分析的快速性，其强大的功能将对以后的学习有更大的帮助。

实验二：线性系统的时域特性一：实验目的1：掌握连续时间信号时域基本运算。

2：用MA TLAB实现信号的时移，反折，尺度变换以及平移，相乘等。

3：掌握线性是不变系统时域分析方法。

二：实验设备1：电脑一台2：Matlab软件三：实验原理与步骤连续时间信号在时域的一些基本运算——尺度变换，平移，翻转，叠加，相乘，微分，积分等不仅涉及信号的描述和分析，还进一步建立有关信号的基本概念和简化信号有着一定的意义。

1：尺度变换，平移和反转：尺度变换可分为幅度尺度变换和时间变换。

将信号一坐标轴为中心进行对称映射，就成了信号的翻转。

平移也称时移，对于信号x（t），考虑到大于零的常数t0 ，则平移信号x（t-t0）或x（t+t0）。

用matlab画出x（3t）和x（2-3t）的波形，其程序如下：%Program2_1_7t=-2:0.01:2;x1=tripuls(3*t,4,0.5);subplot(1,2,1);plot(t,x1);title('x(3t)');x2=tripuls((2-3*t),4,0.5);subplot(1,2,2);plot(t,x2);title('x(2-3t)')其图形如图1图12：叠加和相乘：信号的叠加和相乘运算在matlab中实现比较简单，直接使用算术运算符号“+”和“*”即可实现。

3微分和积分：在matlab中连续信号的微分用differ函数近似计算；定积分可用quad函数实现。

对采样函数进行微分和积分运算，其程序如下：%Program2_1_8h=0.001;t=(0:h:4);x=sinc(t);y1=diff(x)*1/h;plot(t(1:lenght(t)-1),y1);title('dSa(t)/dt');%Program2_1_9t=(0:0.01:4);y=sinc(x);plot(x,y)图形如图2图2四：实验结果与数据处理实验仿真图形与数据见图1 ，图2五：实验体会实验三：信号的频域分析一：实验目的1：掌握周期信号傅里叶级数分解与合成的计算公式。

2：掌握用matlab实现周期信号傅里叶级数分解与综合方法。

3：理解并掌握周期信号频谱特点。

二：实验设备1：电脑一台2：Matlab软件三：实验原理与步骤1：连续时间信号频域分析（傅里叶变换）傅里叶变换是对连续时间信号进行频谱分析的主要工具。

在matlab中可以利用函数fourier进行信号的傅里叶分析。

（1）：求信号x(t)=2te-的频谱。

程序如下：syms t w;x=exp(-t*t);Xw=fourier(x,t,w);ezplot(Xw)图形如下：（2）求单边指数信号x(t)=)e t-的频谱u(t程序如下：syms t w;ut=sym('Heaviside(t)'); x=exp(-t)*ut; Xw=fourier(x); XXP=abs(Xw); ezplot(XXP) 其图形如下（3）：求224)(a we w X -=的傅里叶反变换，其程序如下syms a w;X=exp((-w*w)/(4*a*a)); x=ifourier(X,w,t); x=simple(x)运行结果为：x =ifourier(exp(-1/4*w^2/a^2),w,t)2：连续时间系统频域分析一个信号除了可以分解为一系列冲击函数之和的形式外，还可以分解为一系列不同频率的正弦信号，为此，必须关心系统输出响应岁频率变化的规律.利用matlab 信号处理工具箱提供的函数freqs 可以计算系统的频率响应，其调用格式为：H=freqs （b ，a ，w ）例；已知描述某系统的微分方程为)()(2)()(3,,,t x t y t y t y =++， 求系统的幅频特性和相频特性解：对方程两边取傅里叶变换，得 ：)()(2)()(32w X w Y w jwY w Y jw =++）（得系统的频率特性函数为：2)(31)()()(2++==jw jw w X w Y w H其程序如下：w=linspace(0,6,300);b=[1]; a=[3 1 2]; H=freqs(b,a,w); subplot(1,2,1); plot(w,abs(H)); xlabel('\omega') ylabel('|H(omega)|') subplot(1,2,2); plot(w,angle(H)); xlabel('\omega') ylabel('\phi(\omega)')图形如下：：四：实验结果与数据处理实验仿真图形与数据见图1 ，图2五：实验体会通过本次试验掌握了周期信号傅里叶级数分解与合成的计算公式和利用matlab 实现周期信号傅里叶级数分解与综合方法，理解了信号频谱特点。

实验四：信号的复频域分析一、实验目的1.学会用MATLAB 进行部分分式展开；2.学会用MATLAB 分析LTI 系统的特性；3.学会用MATLAB 进行Laplace 正、反变换。

二：实验设备1：电脑一台 2：Matlab 软件三、实验原理与步骤1．用MATLAB 进行部分分式展开用MATLAB 函数residue 可以得到复杂有理分式F(s)的部分分式展开式，其调用格式为 [],,(,)r p k residue num den =其中，num,den 分别为F(s)的分子和分母多项式的系数向量，r 为部分分式的系数，p 为极点，k 为F(s)中整式部分的系数，若F(s)为有理真分式，则k 为零。

例：用部分分式展开法求F(s)的反变换 322()43s F s s s s+=++其MATLAB 程序为 format rat; num=[1,2]; den=[1,4,3,0];[r,p]=residue(num,den)程序中format rat 是将结果数据以分数形式显示F(s)可展开为 210.536()13F s s s s --=++++ 所以，F(s)的反变换为 3211()()326t t f t e e u t --⎡⎤=--⎢⎥⎣⎦2．用MATLAB 分析LTI 系统的特性系统函数H （s ）通常是一个有理分式，其分子和分母均为多项式。

计算H （s ）的零极点可以应用MATLAB 中的roots 函数，求出分子和分母多项式的根，然后用plot 命令画图。

在MATLAB 中还有一种更简便的方法画系统函数H （s ）的零极点分布图，即用pzmap 函数画图。

其调用格式为pzmap(sys)sys 表示LTI 系统的模型，要借助tf 函数获得，其调用格式为sys=tf(b,a)式中，b 和a 分别为系统函数H （s ）的分子和分母多项式的系数向量。

如果已知系统函数H （s ），求系统的单位冲激响应h(t)和频率响应H ω（j ）可以用以前介绍过的impulse 和freqs 函数。

例： 已知系统函数为 321221s s s +++H(s)=试画出其零极点分布图，求系统的单位冲激响应h(t)和频率响应H ω（j ），并判断系统是否稳定。

解：其MATLAB 程序如下： num=[1]; den=[1,2,2,1]; sys=tf(num,den); figure(1);pzmap(sys); t=0:0.02:10;h=impulse(num,den,t); figure(2);plot(t,h) title('Impulse Response') [H,w]=freqs(num,den); figure(3);plot(w,abs(H)) xlabel('\omega')title('Magnitude Response') 结果如下图：3．用MATLAB 进行Laplace 正、反变换MATLAB 的符号数学工具箱提供了计算Laplace 正、反变换的函数Laplace 和ilaplace,其调用格式为()()F laplace f f ilaplace F ==上述两式右端的f 和F 分别为时域表示式和s 域表示式的符号表示，可以应用函数sym实现，其调用格式为S=sym(A)式中，A 为待分析表示式的字符串，S 为符号数字或变量。

合集下载

信号处理实验报告

信号处理实验报告

信号处理实验报告实验目的：通过实验了解信号处理的基本原理和方法，并掌握使用MATLAB进行信号处理的基本操作。

实验原理：信号处理是指对模拟信号或数字信号进行分析、处理、提取有用信息的过程。

信号处理包括信号的采集、滤波、降噪、特征提取等核心内容。

MATLAB 是一种功能强大的数学软件，也是信号处理的常用工具。

通过使用MATLAB，可以对信号进行快速、准确的处理和分析。

实验过程：1. 使用MATLAB生成一个正弦信号，频率为100Hz，幅值为1，时长为1s。

matlabt = 0:0.001:1;f = 100;x = sin(2*pi*f*t);2. 绘制该信号的时域图像。

matlabfigure;plot(t, x);xlabel('时间（s）');ylabel('幅值');title('正弦信号的时域图像');3. 使用MATLAB进行频谱分析。

matlabN = length(x);f = (0:N-1)*(1/N);X = fft(x);P = abs(X).^2/N;figure;plot(f,P);xlabel('频率（Hz）');ylabel('功率谱密度');title('信号的频谱图像');4. 对信号进行滤波，去除高频成分。

matlabfs = 1000;Wp = 200/(fs/2);Ws = 300/(fs/2);Rp = 3;Rs = 60;[n,Ws] = cheb2ord(Wp,Ws,Rp,Rs); [b,a] = cheby2(n,Rs,Ws);y = filter(b,a,x);figure;plot(t,y);xlabel('时间（s）');ylabel('幅值');title('去除高频成分后的信号');5. 对滤波后的信号进行降噪处理。

信号分析与处理实验报告

信号分析与处理实验报告

《信号分析与处理》实验报告华北电力大学前言1．实验总体目标通过实验，巩固掌握课程的讲授内容，使学生对信号分析与线性系统分析的基本理论及分析方法有一个感性认识和更好地理解，使学生在分析问题与解决问题的能力及实践技能方面有所提高。

2．适用专业自动化专业本科生3．先修课程信号分析与处理4．实验课时分配5需要配置微机及MATLAB工具软件。

6．实验总体要求1、掌握信号分解的基本思想及信号在时域、频域和变换域进行分解的基本理论及描述方法，用MATLAB编程语言实现基本信号的表示及可视化，计算和分析信号的频谱；2、掌握在时域、频域和变换域分析LTI系统的方法，及系统在时域、频域和变换域的描述方法，用MATLAB编程语言实现LTI系统的时域分析及频率分析。

3、掌握信号的调制与解调，用MATLAB编程语言仿真分析信号的调制与解调。

⒎ 本实验的重点、难点及教学方法建议实验通过MATLAB编程语言来实现基本信号的表示及可视化，计算分析信号的频谱，实现LTI系统的时域分析及频率分析，并仿真分析信号的调制与解调，使学生对信号分析与线性系统分析的基本理论及分析方法有一个感性认识和更好地理解。

实验的重点及难点是：掌握基本信号的数学表示，信号的频谱特点，计算LTI系统的典型响应，掌握信号的调制与解调。

在这样的理论基础上，学会用MATLAB编程语言来实现对信号与系统响应的可视化及对数字滤波器进行设计。

教学建议：打好理论基础，熟练编程语言。

目录实验一信号的时域与频域分析 3实验二信号的时域与频域处理 4实验三数字滤波器的设计 5实验一一、实验目的1、熟悉MATLAB 平台，高效的数值计算及符号计算功能；2、实现基本信号的表示及可视化计算；3、分析信号的频谱。

二、 实验类型验证型 三、 实验仪器微机，MATLAB 工具软件。

四、 实验原理MATLAB 是功能强大的数学软件，它提供了计算周期连续函数和周期离散序列的频谱的一系列函数。

信号处理实验报告

信号处理实验报告

一、实验目的本次实验旨在通过MATLAB软件平台，对数字信号处理的基本概念、原理和方法进行学习和实践。

通过实验，加深对以下内容的理解：1. 离散时间信号的基本概念和性质；2. 离散时间系统及其特性；3. 离散傅里叶变换（DFT）及其性质；4. 离散傅里叶逆变换（IDFT）及其应用；5. 窗函数及其在信号处理中的应用。

二、实验内容1. 离散时间信号的产生与性质（1）实验步骤：1.1 利用MATLAB生成以下离散时间信号：- 单位脉冲序列：δ[n]；- 单位阶跃序列：u[n]；- 矩形序列：R[n]；- 实指数序列：a^n；- 复指数序列：e^(jωn)。

1.2 分析并比较这些信号的性质，如自相关函数、功率谱密度等。

（2）实验结果：实验结果显示，不同类型的离散时间信号具有不同的性质。

例如，单位脉冲序列的自相关函数为δ[n]，功率谱密度为无穷大；单位阶跃序列的自相关函数为R[n]，功率谱密度为有限值；矩形序列的自相关函数为R[n]，功率谱密度为无穷大；实指数序列和复指数序列的自相关函数和功率谱密度均为有限值。

2. 离散时间系统及其特性（1）实验步骤：2.1 利用MATLAB构建以下离散时间系统：- 线性时不变系统：y[n] = x[n] a^n；- 非线性时不变系统：y[n] = x[n]^2；- 线性时变系统：y[n] = x[n] (1 + n)。

2.2 分析并比较这些系统的特性，如稳定性、因果性、线性时不变性等。

（2）实验结果：实验结果显示，不同类型的离散时间系统具有不同的特性。

例如，线性时不变系统的输出与输入之间存在线性关系，且满足时不变性；非线性时不变系统的输出与输入之间存在非线性关系，但满足时不变性；线性时变系统的输出与输入之间存在线性关系，但满足时变性。

3. 离散傅里叶变换（DFT）及其性质（1）实验步骤：3.1 利用MATLAB对以下离散时间信号进行DFT变换：- 单位脉冲序列：δ[n]；- 单位阶跃序列：u[n]；- 矩形序列：R[n]。

数字信号处理基础实验报告一概要

数字信号处理基础实验报告一概要

Laboratory Exercise 1 (2 class hours) DISCRETE-TIME SIGNALS: TIME-DOMAIN REPRESENTATION Project 1.1 Unit impulse and unit step sequencesA copy of Program P1_1 is given below.% Program P1_1% Generation of a Unit impulse Sequenceclf;% Generate a vector from -10 to 20n = -10:20;% Generate the unit impulse sequencedelta = [zeros(1,10) 1 zeros(1,20)];% Plot the unit impulse sequencestem(n,delta);xlabel('Time index n');ylabel('Amplitude');title('Unit Impulse Sequence');axis([-10 20 0 1.2]);Answers:Q1.1The unit impulse sequence δ[n] generated by running Program P1_1 is shown below:Q1.2The modified Program P1_1 to generate a delayed unit sample sequence δd[n]with a delay of 11 samples is given below along with the sequence generated by running this program.% Program Q1.2% Generation of a Unit impulse Sequence with a delay of 11 samplesclf;n = -10:20;delta = [zeros(1,21) 1 zeros(1,9)];stem(n,delta);xlabel('Time index n');ylabel('Amplitude');title('Unit Impulse Sequence with a delay of 11 samples');axis([-10 20 0 1.2]);Q1.3The modified Program P1_1 to generate a unit step sequence u[n]is given below along with the sequence generated by running this program.% Program Q1.3% Generation of a unit step sequenceclf;n = -10:20;delta = [zeros(1,10),ones(1,21)];stem(n,delta);xlabel('Time index n');ylabel('Amplitude');title('Unit Step Sequence');axis([-10 20 0 1.2]);Q1.4The modified Program P1_1 to generate a unit step sequence ud[n] with an advance of 7 samples is given below along with the sequence generated by running this program.% Program Q1.4% Generation of a unit step sequence with an advance of 7 samplesclf;n = -10:20;delta = [zeros(1,17),ones(1,14)];stem(n,delta);xlabel('Time index n');ylabel('Amplitude');title('Unit Step Sequence with an advance of 7 samples');axis([-10 20 0 1.2]);Project 1.2 Exponential signalsA copy of Programs P1_2 and P1_3 are given below.% Program P1_2% Generation of a complex exponential sequence clf;c = -(1/12)+(pi/6)*i;K = 2;n = 0:40;x = K*exp(c*n);subplot(2,1,1);stem(n,real(x));xlabel('Time index n');ylabel('Amplitude'); title('Real part');subplot(2,1,2);stem(n,imag(x));xlabel('Time index n');ylabel('Amplitude'); title('Imaginary part');% Program P1_3% Generation of a real exponential sequenceclf;n = 0:35; a = 1.2; K = 0.2;x = K*a.^n;stem(n,x);xlabel('Time index n');ylabel('Amplitude');Answers:Q1.5The complex-valued exponential sequence generated by running Program P1_2 is shown below:Q1.6The purpose of the operator real is: Complex real part.确定一个复数或矩阵的实部；The purpose of the operator imag is: Complex image part. 确定一个复数或矩阵的虚部.Q1.7The purpose of the command subplot is:Create axes in tiled positions.（将图形窗口分成多个矩形窗格来显示多个图形）Q1.8The real-valued exponential sequence generated by running Program P1_3 is shown below:Q1.9The difference between the arithmetic operators ^ and .^ is : ^数值乘方，而.^矩阵乘方。