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(100