2 下载《标准实验报告》 1 3 下载 实验题目 4 提交 实验报告 切记：请按时上传作业！到时将自动关机！ 13:26:23.

Presentation on theme: "2 下载《标准实验报告》 1 3 下载 实验题目 4 提交 实验报告 切记：请按时上传作业！到时将自动关机！ 13:26:23."— Presentation transcript:

1 2 下载《标准实验报告》 1 3 下载 实验题目 4 提交 实验报告 切记：请按时上传作业！到时将自动关机！ 13:26:23

2 程序、图、结果 13:26:23

3 实验一 FT及互相关的初步应用 教育网 http://222.197.183.243/wlxt
2019/8/4 实验一 FT及互相关的初步应用 教育网 公网 13:26:23

4 一、实验目的 1. 熟悉频谱（幅度谱）的特性 2. 熟悉互相关的初步应用 13:26:23

5 三、实验内容 1.利用FFT/IFFT实现简单的频域滤波 要求与思路： 1）构造一低频离散信号f；
2）在此低频信号上加随机噪声得到混合信号g； 3）对信号g进行FFT变换得到其幅度谱和相位谱； 4）对高于某个频率的幅度谱全部置为0； 5）利用IFFT对频谱做反变换得到滤波后的信号f’; 6）画图比较 f 和 f’。 部分提示见下页…. 13:26:24

6 y1 = sin(x1*0.3); %自己构造有用的信号/序列 noise = rand(1,nx) - 0.5;
x1 = [1:nx]; y1 = sin(x1*0.3); %自己构造有用的信号/序列 noise = rand(1,nx) - 0.5; y2 = y1 + noise; figure(1); plot(x1,y1,'r', x1,y2,'b:'); XXX = fft(y2); figure(2); plot(abs(XXX)); title('信号的幅度谱'); %改造频谱，实现滤波 for i = 10 : nx-10 %自己选择合适滤波范围 XXX(i) = 0; end f1 = ifft(XXX); figure(3); plot(x1,y1,'r',x1,f1,'b:'); title('滤波的效果'); 13:26:24

7 幅度 13:26:24 样点序号

8 2019/8/4 幅度 设置为0 π (π) (2π) 13:26:24

9 13:26:24 样点序号

10 目的：利用互相关运算寻找被噪声淹没的有用信号的位置。
2. 互相关的一个简单应用 目的：利用互相关运算寻找被噪声淹没的有用信号的位置。 13:26:24

11 y1[n] y1[n]出现在y2[n]的某个位置 y2[n] 在y2[n]上加随机噪声 y3[n] 13:26:24 样点序号

12 y1[n] y3[n] r[n] = xcorr(y3,y1); 13:26:24 样点序号

13 %寻找序列r中的最大值maxR及其所在的位置(序号)id： [maxR, id] = max(r);
… %自己构造序列y3[n]和y1[n] r = xcorr(y3,y1); nx = length(r); nn = (nx-1)/2; plot([-nn:nn], r); …… 注：设y3[n], y1[n]两个序列的最大长度为M，则xcorr(y3, y1)返回的互相关序列r的长度为2*M-1，其样点号范围是[-(M-1):(M-1)]. %寻找序列r中的最大值maxR及其所在的位置(序号)id： [maxR, id] = max(r); （用在此题目中，需再根据上面的注明进一步 调整id，以得到信号所在的正确位置） 13:26:24

14 XCORR Cross-correlation function estimates.
>> help xcorr XCORR Cross-correlation function estimates. C = XCORR(A,B), where A and B are length M vectors (M>1), returns the length 2*M-1 cross-correlation sequence C. If A and B are of different length, the shortest one is zero-padded. …… The zeroth lag of the output correlation is in the middle of the sequence, at element or row M. 13:26:24

15 y1[n] y3[n] M 2*M-1 13:26:24