MATLAB课件第六章201404汇编.ppt

第六章 MATLAB在信号分析与处理中的应用 信号及其表示 连续时间信号 离散时间信号 信号分析及处理 序列的操作 信号频谱分析 图像线性变换 例1. 产生一个方波，正信号40% t 0:0.01:6*pi; x square t,40 plot t,x axis [0,20,-1.1,1.1] 例2. 产生一个1.5s的50Hz的锯齿波，采样率10kHz. fs 10000; t 0:1/fs:1.5; x sawtooth 2*pi*50*t ; plot t,x axis [0,0.2,-1,1] 例3. 产生一个sinc信号. t linspace -5,5 ; x sinc t ; plot t,x 例4. 产生一个脉冲调制信号 t 0:1/50e3:10e-3; %抽样频率50kHz,抽样时间10ms d [0:1/1e3:10e-3;0.8.^ 0:10 ]’; %第一列说明每一脉冲的延迟时间，第二列说明每一次幅值衰减 x pulstran t,d,’gauspuls’,10e3,0.5 ; %调制波为高斯噪声，频率1kHz,带宽50% plot t,x 6.1.2 离散时间信号的表示 例. 绘制离散时间信号的棒状图。其中x -1 -1, x 0 1, x 1 2, x 2 1, x 3 0, x 4 -1。 n -3:5; %定位时间变量 x [0,0,-1,1,2,1,0,-1,0]; stem n,x ; grid; line [-3,5],[0,0] ; xlabel n ; ylabel x[n] 6.1.3 几种常用离散时间信号的表示 6.2 信号分析与处理 6.2.1 序列的操作 1. 信号加 x n x1 n +x2 n MATLAB实现： x x1+x2; 2. 信号乘 x n x1 n x2 n MATLAB实现： x x1.*x2; 3. 改变比例 y n k*x n MATLAB实现： y k*x; 4. 折叠y n x -n MATLAB实现： y fliplr x ; 5. 抽样和 MATLAB实现： y sum x n1:n2 ; 6. 抽样积 MATLAB实现： y prod x n1:n2 ; 7. 信号能量 MATLAB实现： y sum abs x .^2 ; 8. 信号功率 MATLAB实现： y sum abs x .^2 /N ; 例5 产生一个叠加信号，并计算它的抽样和、抽样积、信号能量功率等。 t 0:0.001:1 ’; x sin 2*pi*50*t +2*sin 2*pi*120*t ; y1 sum x 1:50 ; %抽样和 y2 prod x 1:50 ; %抽样积 Ex sum abs x .^2 ; %信号能量 Px sum abs x .^2 /length x ; %信号功率 6.2.2 信号频谱分析 将时域信号变换至频域加以分析的方法称为频谱分析。频谱分析的目的是把复杂的时间历程波形，经过傅里叶变换分解为若干单一的谐波分量来研究，以获得信号的频率结构以及各谐波和相位信息。 MATLAB频谱分析函数 y fft x ，返回向量x的不连续fourier变换 信号频谱分析举例 例6 利用函数生成一组数据用以模拟电力现场测量数据，并利用MATLAB编程实现其频谱分析。其中，x 2 sin 2pf1t +?sin 2pf2t ，f1? ?100，f2? ?300。在信号x中加入随机信号，用以模拟信号噪声。 fs 1000; t 0:1/fs:0.6 %设置采样频率为1000 Hz f1 100;f2 300;x sin 2*pi*f1*t +sin 2*pi*f2*t ; subplot 411 ;plot x ; title f1 100Hz \f2 300Hz xlabel n ; grid on; number 512 %设置用于FFT计算的数据点数为512个； y fft x,number ; %对x信号进行512点的傅里叶变换； n 0:length y -1; f fs*n/length y %设置频率轴坐标 subplot 412 plot f,abs y title FFT xlabel Hz grid on x x+randn size x ; %在信号中加入噪声 subplot 413 ,plot x title 含随机噪声的信号 xlabel 序列N ,grid on y fft x,number ; %对含噪声信号进行512点的傅立叶变换 n 0:length y