基于Matlab的数字信号处理基本操作
- 格式:doc
- 大小:262.92 KB
- 文档页数:25
温州大学物理与电子信息工程学院Matlab 仿真及其应用 实验报告实验一Matlab 基本功能和基础知识操作 [实验目的和要求]1、 熟练掌握Matlab 的启动与退出2、 熟悉Matlab 的命令窗口、常用命令、帮助系统3、 熟悉Matlab 的数据类型、基本矩阵操作、运算符和字符串处理[实验内容]1、 用逻辑表达式球下列分段函数的值 22201112,=0:0.5:2.52123t t y t t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩其中2、 求[100,999]之间能被32整除的数的个数3、 建立一个字符串向量,删除其中的小写字母。
4、 输入矩阵1234514789A ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦,求出此矩阵的行列式,逆和特征根,并找出A 中大于5和小于9的元素,用行列式表示。
5、 不采用循环的形式求出和式63230034ii i i S ===+∑∑6、 给定矩阵E=rand (4,4),计算C+E ,C*E ,C\E实验结果及分析:经过Matlab 软件的程序编辑和测试分析,得出以下实验结果: 详见程序代码、注释及屏幕截图:【题1】程序代码:t=0:0.5:2.5y=t.^2.*((t>=0)&(t<1))+(t.^2-1).*((t>=1)&(t<2))+(t.^2-2*t+1).*((t>=2)&(t<3)) 效果截图:【题2】程序代码:p=rem([100:999],32)==0;sum(p)效果截图:【题3】程序代码:ch='dfghjGUIJKVC',k=find(ch>'a'&ch<='z'),ch(k)=[]效果截图:【题4】程序代码:A=[1 2 3;4 5 14;7 8 9];[i,j]=find(A>5&A<9) %定位for n=1:length(i)m(n)=A(i(n),j(n))endDA=det(A) %行列式IA=inv(A) %逆矩阵EA=eig(A) %特征根效果截图:【题5】程序代码:E=rand(4,4); %产生随机数0~1 C=rand(4,4);B1=C+E;B2=C*E;B3=C/E;B1B2B3效果截图:【题6】程序代码:E=rand(4,4); %产生随机数0~1 C=rand(4,4);B1=C+E;B2=C*E;B3=C/E;B1B2B3效果截图:4、心得:通过本次Matlab课程实验,我已熟练Matlab的命令窗口、常用命令、帮助系统,并掌握Matlab的数据类型、基本矩阵操作、运算符和字符串处理。
Matlab中的数字信号处理方法与实例数字信号处理是一门研究数字信号在数字域中分析、处理和改变的学科。
Matlab是一种强大的数值计算工具,被广泛应用于信号处理领域。
本文将介绍一些在Matlab中常用的数字信号处理方法与实例,并通过实例来展示它们的应用。
1. 信号的采样与重构信号采样是指将连续时间信号转化为离散时间信号的过程。
在Matlab中,我们可以使用“sample”函数对信号进行采样,并使用“hold”函数对采样后的信号进行重构。
下面是一个示例:```matlabfs = 100; % 采样频率t = 0:1/fs:1; % 时间序列x = sin(2*pi*5*t); % 原始信号subplot(2,1,1);plot(t,x);title('原始信号');xlabel('时间');ylabel('幅值');subplot(2,1,2);stem(t,x);title('采样和重构后的信号');xlabel('时间');ylabel('幅值');```在这个例子中,我们生成了一个频率为5Hz的正弦信号,然后对该信号进行采样和重构。
从结果可以看出,原始信号和重构后的信号基本上是一致的。
2. 信号的频谱分析频谱分析是指将信号从时域转换到频域的过程,可以用来分析信号的频率成分。
在Matlab中,我们可以使用“fft”函数对信号进行傅里叶变换,并使用“abs”函数获取信号的幅度谱。
下面是一个示例,演示如何对信号进行频谱分析:```matlabfs = 100; % 采样频率t = 0:1/fs:1; % 时间序列x = sin(2*pi*5*t); % 原始信号N = length(x); % 信号长度X = fft(x); % 傅里叶变换f = (0:N-1)*(fs/N); % 频率坐标plot(f,abs(X));title('信号的频谱');xlabel('频率');ylabel('幅度');```在这个示例中,我们同样生成了一个频率为5Hz的正弦信号,然后对该信号进行傅里叶变换,并绘制出信号的频谱图。
Matlab数字信号处理实验报告本次实验使用MATLAB进行数字信号处理操作,目的是熟悉MATLAB中数字信号处理的相关工具箱,并进一步理解数字信号处理的基本概念和算法。
一、实验内容1.信号的生成与显示2.时域分析和频域分析3.滤波器设计4.数字滤波器性能分析二、实验步骤在MATLAB中,使用sawtooth函数生成一个锯齿波信号,并使用plot函数进行时域波形的显示。
代码如下:f = 1000;fs = 40000;t = 0:1/fs:0.01;y = sawtooth(2*pi*f*t);plot(t,y);xlabel('Time (s)');ylabel('Amplitude');title('Sawtooth Wave');时域分析包括波形的观察和参数分析,如幅值、均值、方差等。
频域分析则是对信号进行傅里叶变换,得到其频谱图,包括频率分布和强度分布。
%时域分析amp = max(y)-min(y);mean_y = mean(y);var_y = var(y);设计一个低通滤波器,将高于1kHz的频率成分滤掉。
对滤波后的信号进行时域分析和频域分析,比较滤波前后信号的特征参数和频谱特征,并绘制原始信号、滤波后信号及其频谱图。
subplot(2,2,1);plot(t,y);xlabel('Time (s)');ylabel('Amplitude');title('Sawtooth Wave');subplot(2,2,2);plot(t,y_filt);xlabel('Time (s)');ylabel('Amplitude');title('Sawtooth Wave After Filter');subplot(2,2,3:4);plot(f2,fft_y_filt,'r',f,fft_y,'g');xlabel('Frequency (Hz)');ylabel('Amplitude');title('Sawtooth Wave Spectrum Comparison');legend('After Filter','Before Filter');三、实验结果与分析通过生成并显示一段锯齿波信号,并对其进行时域和频域分析,可以得到该信号的关键信息,如幅值、均值、方差和频率分布特性。
使用Matlab进行数字信号处理的技巧与注意事项1. 引言数字信号处理(Digital Signal Processing,简称DSP)是指通过对离散时间信号进行算法处理以达到某种目的的一种技术。
在现代科技和工程领域中,DSP被广泛应用于通信、图像处理、音频处理、医学诊断、雷达和控制系统等众多领域。
而Matlab作为一种强大的数学计算工具,具备优秀的信号处理和算法库,成为众多工程师和研究人员进行数字信号处理的首选之一。
本文将介绍一些使用Matlab 进行数字信号处理时的技巧与注意事项。
2. 信号处理基础知识在使用Matlab进行数字信号处理之前,有一些基础知识是必备的。
首先是对信号的了解,信号可以分为连续时间信号和离散时间信号。
连续时间信号指的是信号在所有时间上都有定义,而离散时间信号则只在某些时间点上有定义。
数字信号处理主要针对离散时间信号进行。
此外,还需要了解采样定理、变换、滤波器以及噪声等基本概念。
3. 信号与信号处理在进行数字信号处理时,首先需要得到待处理的信号。
Matlab提供了多种方法来生成信号,比如使用波形发生器函数、加载文件以及使用模型等。
根据具体情况选择合适的方法生成待处理信号。
4. 信号的可视化与分析在开始处理信号之前,可以使用Matlab中的图形工具对信号进行可视化和分析。
例如,使用plot函数可以绘制信号的时域波形图,使用spectrogram函数可以绘制信号的频谱图,利用histogram函数可以绘制信号的直方图等。
这些图像可以帮助我们更好地理解信号的特征和性质。
5. 信号的滤波处理滤波是数字信号处理中常用的操作之一,用于去除信号中的噪声、增强信号的频率特征等。
Matlab提供了丰富的滤波器设计函数,包括有限冲激响应(FIR)滤波器和无限冲击响应(IIR)滤波器等。
通过选择合适的滤波器类型、阶数和截止频率等参数,可以实现对信号的滤波处理。
6. 时频分析时频分析用于分析信号在时间和频率上的变化情况,帮助我们更全面地认识信号的特性。
如何使用Matlab进行信号处理导言信号处理是指从源信号中提取或改变信息的过程,而Matlab是一种功能强大的数学计算和编程软件,广泛应用于各个领域的数据分析和处理。
本文将介绍如何使用Matlab进行信号处理,包括信号读入、滤波、频谱分析和信号合成等几个方面。
一、信号读入在进行信号处理之前,首先需要将信号读入到Matlab中。
Matlab提供了多种方式来读入信号,常用的有以下几种:1. 读取音频信号使用Matlab的`audioread`函数可以读取多种音频格式的信号文件,例如:```[x, fs] = audioread('audio.wav');```其中,`x`是读入的音频信号,`fs`是采样率。
2. 读取图像信号使用Matlab的`imread`函数可以读取多种图像格式的信号文件,例如:```x = imread('image.jpg');```其中,`x`是读入的图像信号。
3. 生成模拟信号如果需要生成模拟信号进行处理,可以使用Matlab的信号生成函数,例如生成正弦信号:```fs = 1000; % 采样率t = 0:1/fs:1; % 时间向量f = 10; % 信号频率x = sin(2*pi*f*t); % 正弦信号```二、滤波滤波是信号处理中常用的技术,可以去除信号中的噪声或者提取感兴趣的频率成分。
Matlab提供了丰富的滤波函数,常用的有以下几种:1. 低通滤波低通滤波可以去除高频噪声,保留低频信号。
使用Matlab的`lowpass`函数可以设计低通滤波器并滤波信号,例如将采样率为1000Hz的信号x通过一个截止频率为100Hz的低通滤波器:```Fc = 100; % 截止频率Wn = Fc / (fs/2); % 归一化截止频率b = fir1(50, Wn, 'low'); % 设计低通滤波器y = filter(b, 1, x); % 低通滤波```2. 高通滤波高通滤波可以去除低频噪声,保留高频信号。
基于MATLAB的数字信号处理1MATLAB的综合实验一、实验目的及要求培养学生利用Matlab解决专业问题的能力。
二、实验设备(环境)及要求1.计算机2.Matlab软件编程实验平台三、实验内容1、编程实现一个数字信号处理的仿真系统。
要求具有界面并实现以下功能:1)能产生(得到)并选择各种数字信号(sin、方波、三角波、语音、噪声及其叠加);2)具有DFT、DCT和DWT变换功能,并对各种信号进行变换;3)设计滤波器实现低通、高通、带通滤波,得到输出信号的频域特性和时间序列;4)输入一段叠加了噪声的语音信号,显示其频谱特性,通过变换或滤波对其降噪,得到输出信号的频域特性和时间序列。
四、设计思想本系统包含有三个主要部分:信号产生与变换模块,滤波器模块和语音噪声处理。
信号产生与变换通过输入信号频率和采样频率实现正弦、方波、三角波、语音信号的产生以及噪声的叠加,系统设定信号持续时间为0.05s,语音信号为截取了一段2s的声音信号。
同时对各个信号进行DFT,DCT和DWT变换,且变换点数N=256,同时设定DWT变换时的小波类型为db1。
滤波器模块设计了四个IIR 滤波器(巴特沃斯、切比雪夫Ⅰ型,切比雪夫Ⅱ型和椭圆滤波器),并分别实现低通,高通和带通。
界面设计了各种滤波器所需参数的输入模块。
系统设定待滤波信号持续时间为0.05s,包含有3个频率成分,S=sin(2*pi*f*t)+ sin(2*pi*5*f*t)+sin(2*pi*8*f*t),其中f为输入信号频率,S通过低通、带通、高通滤波器之后,分别得到频率为f,5f和8f的正弦信号,实现信号滤波。
语音噪声处理部分是一个复选框按钮,通过巴特沃斯低通滤波器对其进行降噪,设计中通过观察噪声语音信号的频谱得到低通滤波器的截止频率和阻带起始频率,并合理输入通带衰减与阻带衰减,最终得到理想的降噪结果。
数字滤波器设计过程中用到了如下的一些matlab设计函数:buttord、butter,cheb1ord、cheby1,cheb2ord、cheby2,ellipord、ellip。
Matlab中的模拟和数字信号处理方法引言:Matlab是一种强大的计算软件工具,广泛应用于科学、工程和数学等领域。
在信号处理领域,Matlab提供了丰富的模拟和数字信号处理方法,极大地方便了信号处理的研究和应用。
本文将介绍一些主要的模拟和数字信号处理方法,以及它们在Matlab中的实现。
一、模拟信号处理方法:1. Fourier变换Fourier变换是一种重要的信号分析方法,可以将信号从时间域转换到频率域,从而揭示信号的频谱特性。
在Matlab中,可以使用fft函数进行傅里叶变换,ifft 函数进行逆傅里叶变换。
通过傅里叶变换,我们可以分析信号的频谱,包括频率成分、功率谱密度等。
2. 滤波滤波是信号处理中常用的方法,可以消除信号中的噪声或者选择感兴趣的频率成分。
在Matlab中,提供了丰富的滤波函数,包括低通滤波器、高通滤波器、带通滤波器等。
通过设计滤波器,我们可以选择不同的滤波方式,如巴特沃斯滤波、切比雪夫滤波等。
3. 时域分析时域分析是对信号在时间域上的特性进行研究,包括信号的振幅、频率、相位等。
在Matlab中,我们可以使用时域分析函数来计算信号的均值、方差、自相关函数等。
通过时域分析,可以更好地了解信号的时间特性,比如周期性、正弦信号等。
二、数字信号处理方法:1. 数字滤波器数字滤波器是将连续时间的信号转换为离散时间的信号,并对其进行滤波处理的一种方法。
在Matlab中,我们可以使用fir1、fir2等函数设计数字滤波器,以满足不同的滤波需求。
数字滤波器可以消除离散信号中的噪声,提取感兴趣的频率成分。
2. 频谱分析频谱分析是对离散信号的频谱进行研究,可以了解信号在频域上的特性。
在Matlab中,可以使用fft函数进行快速傅里叶变换,得到离散信号的频谱。
通过频谱分析,我们可以掌握信号的频率成分、频率幅度等信息。
3. 信号编码信号编码是将模拟信号转换为数字信号的过程,以进行数字信号处理和传输。
实验一Matlab的基本操作一、实验目的:熟悉MATLAB软件布局,掌握MA TLAB基本命令和操作。
二、实验仪器:电脑一台,MATLAB6.5或更高级版本软件一套。
三、实验内容:1.MATLAB简介:MATLAB的名称源自Matrix Laboratory,1984年由美国Mathworks 公司推向市场。
它是一种科学计算软件,专门以矩阵的形式处理数据。
MA TLAB将高性能的数值计算和可视化集成在一起,并提供了大量的内置函数,从而被广泛地应用于科学计算、控制系统、信息处理等领域的分析、仿真和设计工作。
MATLAB软件包括五大通用功能:数值计算功能(Nemeric);符号运算功能(Symbolic);数据可视化功能(Graphic);数据图形文字统一处理功能(Notebook)和建模仿真可视化功能(Simulink)。
该软件有三大特点:一是功能强大;二是界面友善、语言自然;三是开放性强。
目前,Mathworks公司已推出30多个应用工具箱。
MATLAB在线性代数、矩阵分析、数值及优化、数理统计和随机信号分析、电路与系统、系统动力学、信号和图像处理、控制理论分析和系统设计、过程控制、建模和仿真、通信系统、以及财政金融等众多领域的理论研究和工程设计中得到了广泛应用。
MATLAB在信号与系统中的应用主要包括符号运算和数值计算仿真分析。
由于信号与系统课程的许多内容都是基于公式演算,而MATLAB借助符号数学工具箱提供的符号运算功能能基本满足信号与系统课程的需求。
例如,解微分方程、傅里叶正反变换、拉普拉斯正反变换、z正反变换等。
MATLAB在信号与系统中的另一主要应用是数值计算与仿真分析,主要包括函数波形绘制、函数运算、冲激响应与阶跃响应仿真分析、信号的时域分析、信号的频谱分析、系统的S域分析、零极点图绘制等内容。
数值计算仿真分析可以帮助学生更深入理解信号与系统的理论知识,并为将来使用MA TLAB进行信号处理领域的各种分析和实际应用打下基础。
使用MATLAB进行信号处理和音频分析的基本教程第一章信号处理基础信号处理是指对信号进行获取、加工和分析的过程。
MATLAB作为一种强大的计算工具,提供了丰富的信号处理函数和工具箱。
在本章中,将介绍信号的概念、信号的表示和MATLAB中常用的信号处理函数。
1.1 信号的概念信号是指随着时间、空间或者其他变量而变化的物理量。
常见的信号类型包括连续时间信号和离散时间信号,以及模拟信号和数字信号。
1.2 信号的表示MATLAB使用向量或矩阵来表示信号。
向量表示一维信号,矩阵表示多维信号。
可以使用MATLAB中的数组操作函数来创建和操作信号。
1.3 信号处理函数MATLAB提供了丰富的信号处理函数,可以用于信号滤波、谱分析、频域变换等。
常用的信号处理函数包括filter、fft、ifft等。
第二章音频处理基础音频处理是指对声音信号进行分析、过滤和增强的过程。
MATLAB提供了强大的音频处理工具箱和函数库。
本章将介绍音频信号的特点、音频处理的基本原理和MATLAB中的音频处理函数。
2.1 音频信号的特点音频信号是由声音振动引起的连续变化的电信号。
它的特点包括频率、幅度、相位等。
2.2 音频处理的基本原理音频处理的基本原理包括滤波、均衡、静音检测、音量控制等。
MATLAB提供了相关函数和工具箱,可以方便地实现这些音频处理功能。
2.3 音频处理函数MATLAB提供了丰富的音频处理函数,包括从音频文件中读取数据、音频信号的滤波、语音识别等。
常用的音频处理函数包括audioread、audiowrite、speechrecognition等。
第三章信号处理实例本章将通过实例演示如何使用MATLAB进行信号处理和音频分析。
具体包括信号滤波、谱分析和音频处理等。
3.1 信号滤波以滤波为例,介绍如何使用MATLAB对信号进行滤波处理。
首先,使用filter函数设计滤波器,然后将信号输入滤波器,最后绘制滤波后的信号波形图。
基于Matlab的数字信号处理基本操作 一、实验目的 1、熟悉Matlab软件的的运行环境, 2、了解Matlab基本命令与操作,熟悉数据输入、输出与编辑方法; 3、用Matlab函数产生离散信号; 4、用Matlab进行离散信号的表示与运算; 5、运用基本函数绘制常用函数的图形。 二、实验环境: 1、Matlab 6.5 2、PC机 三、实验方法: Matlab操作界面的默认窗口如图1.1所示,Matlab的通用操作界面包括9个常用的窗口,分别是:命令窗口、历史命令窗口、当前目录浏览器窗口、工作空间浏览器窗口、数组编辑窗口、交互界面分类目录窗口、M文件编辑/调试器窗口、帮助导航/浏览器窗口和程序性能剖析窗口。 图1-1 进入Matlab 6.5命令窗口(Command Window),命令窗口会默认出现在Matlab界面的右侧,在命令窗口可键入各种Matlab的命令、函数和表达式,并显示除图形以外的所有运算结果。 命令串口也可单独显示,选择菜单“VIEW”→“Undock Command Window”命令,可使单独的命令窗口返回Matlab界面。 Matlab 6.5版运行时,命令窗口中的每个命令行前会出现提示符“〉〉”。与以前的版本不同,命令窗口内显示的字符和数值采用不同的颜色,在默认的情况下,输入的命令、表达式以及计算结果等采用黑色字体,字符串采用红色,“if”、“for”等关键词采用蓝色。 命令窗口中命令行的编辑,由于Matlab把命令窗口中输入的所有命令都记录在专门的“历史命令”(Command History)空间中,因此Matlab命令窗口不仅可以对输入的命令进行编辑和运行,而且可以对已输入的命令进行回调、编辑和重运行。 命令窗编辑功能:键入和修改程序的方法与通常的文字处理相仿。特殊的功能键为: ESC:恢复命令输入的空白状态 ↓:调出下一行命令 ↑:调出上一行(历史)命令 命令窗编辑功能在程序调试时十分有用。对于已执行过的命令,如要做些修改后重新执行,就可不必重新键入,用↑键调出原命令做修改即可。 主菜单中的编辑项功能:用它可以把屏幕上加深了的文字裁减(Cut)或复制(Copy)下来,然后粘贴(Paste)到任一其他视窗的任何位置上去。 主菜单中的视图项功能:用它可以改变屏幕上显示的视窗布局。例如,我们希望只显示命令窗,使它占整个屏幕,依次引出View的下拉菜单,即【View】→【Desktop Layout】→【Command Window Only】 键入“help”,屏幕上将显示系统中已装入的函数库(即子目录)的名称。 §键入help子目录名,如help elfun,即得出elfun库中各函数名。 §键入help函数名,如help tan2,即得到tan2函数的意义及用法。 退出MATLAB有两种方法。一种是键入exit或quit,另一种是用鼠标双击左上角的小方块或单击右上角的×号,后者是非正常退出,该次的输入命令将不记录在“历史命令”中,所以应当尽量避免使用。 在命令窗口中,默认情况下数值计算结果的显示格式为:当数值为整数,以整数显示;当数值为实数,以小数后4位的精度情况显示,即以“短(Short)”格式显示;如果数值的有效数字超出了这一范围,则以科学计数法显示结果。 注意:数值的显示精度并不代表数值的存储精度。 在命令窗口可以直接输入操作命令,对命令窗口进行控制。 例如: (1)用于清空命令窗口中的显示内容 语法:clc (2)在命令窗口中控制其后每页的显示内容行数。 语法:more 通常,只要执行了一种绘图命令,就会自动产生图形窗,以后的绘图都在这一个图形窗中进行。如想再建一个或几个图形窗,则可键入figure,MATLAB会新建一个图形窗,并自动给它依次排序。如果要认为规定新图为图3,则可键入figure(3)。如果调看已经存在的图形窗3,也应键入figure(3)。 MATLAB程序编制有两种方式。一种称为行命令方式,这就是在命令窗中一行一行地输入程序,计算机每次对一行命令作出反应,像计算器那样。这只能编辑简单的程序,在入门时可以用这种方式。程序稍复杂一些,就应把程序写成一个有多行语句组成的文件,让MATLAB来执行这个文件。编写和修改这种文件程序就要用到文本编辑器。 命令窗上方最左边的按钮是用来打开文本编辑器空白页的,左边第二个按钮是用来打开原有程序文件的。 在命令窗中键入demo,将出现MATLAB的演示窗,演示窗的左侧是库目录。双击目录中的名称或选中该项目后点击右下角的【Run „】方框,即出现该项目的演示界面。通常,演示画面的右侧是一些功能按钮,左上半部是图形,而左下半部则是相应的MATLAB程序语句。还可以在界面上直接修改这些语句并重新执行。因此演示程序也是一个很好的学习过程。 四、实验用MATLAB语言工具函数简介 本次试验在熟悉软件的运行环境的基础上,要求完成信号的产生以及图形的输出。 MATLAB可以根据给出的数据,用绘图命令在屏幕上画出其图形,通过图形对科学计算进行描述。 Plot命令用来绘制x-y坐标中的曲线。它是一个功能很强的命令。输入变量不同,可以产生很多不同的结果。 1、 plot(y):输入一个数组的情况 如果y是一个数组,函数plot(y)给出线性直角坐标的二维图,以y中元素的下标作为X坐标,y中元素的值作为Y坐标,一一画出在X-Y坐标平面图上,而且将各点以直线相连。例如,要画出是个随机数的曲线。可列出: y=5*(rand(1,10)-.5) 由Rand函数产生的随机数的最大值为1,最小数为0,平均值为0.5。所以y的最大值为2.5,最小值为-2.5,平均值为0。键入plot(y),MATLAB会产生一个图形窗,自动规定最合适的坐标比例绘图。X方向是下标,从1~10,Y方向范围则是-4~4,并自动标出刻度。可以用title命令给图加上标题,用xlabel,ylabel命令给坐标轴加上说明,用text或gext命令可在图上任何位置加标注,也可用grid命令在图上打上坐标网格线。 MATLAB实现程序和形成的图如下: y=5*(rand(1,10)-.5) y = Columns 1 through 9 2.2506 -1.3443 0.5342 -0.0701 1.9565 1.3105 -0.2177 -2.4075 1.6070 Column 10 -0.2765 >> plot(y) >> title('my first plot') >> xlabel('x'),ylabel('y') >> grid
2、 plot(x,y):输入两个数组的情况 如果数组x和y具有相同长度,命令plot(x,y)将绘出以x元素为横坐标,y元素为纵坐标的曲线。 3、 用plot(t,[y1,y2,„])命令 该语句中t是向量,y=[y1,y2,„]是矩阵,若t是列(行)向量,则y的列(行)长与t长度相同。Y的行(列)数就是曲线的根树。例如, >> t=0:0.5:4*pi; >> y=exp(-0.1*t).*sin(t); >> y1=exp(-0.1*t).*sin(t+1); >> plot(t,[y;y1])
4、 用hold命令 在画完前一张图后用hold命令保持住,再画下一张曲线。如键入 Plot(t,y),hold on,plot(t,y1,‘g’) 执行此命令时,图形窗产生第一幅图形,同时,命令屏幕显示Current plot held,图形处于保持状态。再执行plot(t,y1,‘g’),就把第二幅图以绿色的曲线迭合在同一张图上。 5、 图形控制 (1)figure:打开图形窗口。MATLAB中的第一幅图随plot命令自动打开,以后的plot命令都画在同一张图上。如要画在另一张新图上,就要用figure命令打开新的图形窗口。有了顺序为1,2,3,„的几个图形窗后,再用plot语句,即键入figure(i),表示打开第i幅图。否则,所有的图都会画在最后显示的那张图上。 (2)clf:清除当前图形窗的内容。 (3)hold:保持当前图形窗的内容,再键入hold,就解除冻结。这种拉线开关式的控制有时会造成混乱,可以用hold on和hold off命令以得到确定的状态。 (4)close:关闭当前图形窗。close all:关闭所有图形窗。 (5)subplot(n,m,p)命令:将图形窗口分为n╳m个子图,在第p个子图处绘制图形。 五、实验内容: 1、利用Matlab产生单位脉冲序列)(n,并画图。 2、利用Matlab产生单位阶跃序列)(n,并画图。 3、利用Matlab产生实指数序列:nX)9.0(,并画图。 4、产生被噪声污染的正弦信号。 5、实现一个由50 HZ和120HZ正弦信号构成的和信号。
6、周期为N的正弦序列)2sin(nN,且10Nn。 7、矩形序列)(nRN 选出4—5个题目作为实验题目。 六、实验报告 1、列出实验程序清单,并附上必要的程序说明。 2、记录实验结果,实验后,对结果进行分析。 实验二 用Matlab进行信号的一维快速傅立叶变换(FFT) 一、试验目的 1、了解傅立叶变换的物理意义:序列的离散傅立叶变换就是序列的离散频谱; 2、了解FFT,掌握FFT函数的调用格式; 3、应用FFT计算信号的频谱。 二、实验原理: 1、离散傅里叶变换(DFT)及其主要性质 DFT表示离散信号的离散频谱,DFT的主要性质中有奇偶对称特性,虚实特性等。通过实验可以加深理解。 例如:实序列的DFT具有偶对称的实部和奇对称的虚部。这可以证明如下: 由定义
10)()(Nnkn
NWnxkX
1010)2sin()()2cos()(NnNnknNnxjknNnx
10)()()(NnnkN
NWnxkNX
10)(NnknNNnWWnx
10)(Nnkn
NWnx
1010)2sin()()2cos()(NnNnknNnxjknNnx
)(*)(kNXkX 实序列DFT的这个特性,在本实验中可以通过实指数序列及三角序列看出来。 对于单一频率的三角序列来说它的DFT谱线也是单一的,这个物理意义我们可以从实验中得到验证,在理论上可以推导如下: 设: