随机信号实验.docVIP

实验一 信号的产生及傅立叶分析（设计性） 一 实验目的 1 学会利用计算机仿真信号。 2 理解信号采样思想。DFT）: 正变换 反变换 和都是点数为N的有限长序列。实质上有限长序列都是作为周期序列的一个周期来表示的，都隐含有周期性意义。 三 实验内容 1 几种常用序列（如正弦、矩形、指数序列等）的产生。 1．用stem函数来画出序列的波形，通过改变N和s得值来改变时间长度、抽样频率 N=500;k=0:N; s=0.03; X=5*sin(s*pi*k); plot(k,X,k,zeros(1,N+1)); Xlabel(k); Ylabel(X[k]); title(余弦序列); 2．指数序列 clear；clc %c：指数序列的幅度 %a：指数序列的底数 %k1：绘制序列的起始序号 %k2：绘制序列的终止序号 c=1;a=0.75;k1=0;k2=20;k=k1:k2; x=c*(a.^k);stem(k,x);%filled Xlabel(k); Ylabel(x);title(??êyDòáD); 3各种序列 t=0:0.01:1;k=1:200; x1=0.1*exp(-2*t); %指数序列 x2=2*cos(2*pi*4*t); %余弦序列 x3=[ones(1,10) zeros(1,90) ones(1,10) zeros(1,90)]; subplot(3,1,1); plot(t,x1); title(指数序列); subplot(3,1,2); plot(t,x2); title(余弦序列); subplot(3,1,3); plot(k,x3); title(矩形序列); 4 编程实现序列的离散傅里叶变换(DFT)，输入x(n)，输出X(k)并且对于不同序列（如矩形序列等）做DFT. clear; clc k=0:31; x1=2*((0.75).^k); subplot 321; stem(k,x1); title(指数序列); y1=fft(x1,32); subplot 322; stem(k,y1); title(指数序列DFT); k=0:31; x2=sin(k); subplot 323; stem(k,x2); title(正弦序列Sin(k)); y2=fft(x2,32); subplot 324; stem(k,y2); title(正弦序列DFT); x3=[ones(1,8) zeros(1,8) ones(1,8) zeros(1,8)]; subplot 325; stem(k,x3); title(矩形序列); y3=fft(x3,32); subplot 326; stem(k,y3); title(矩形序列DFT); 实验二 快速傅立叶变换FFT及频谱分析（设计性） 一 实验目的 1 进一步加深对DFT算法原理和基本性质的理解 2 熟悉FFT算法原理和FFT的子程序应用 3 学习用FFT对连续时间信号进行频谱分析的方法，了解可能出现的分析误差及原因 二 实验原理（参考P187，P189） FFT只是DFT的一种快速算法，利用FFT可减少运算量，提高速度。 按时间抽选的基2FFT算法原理： 将输入序列X(n)按照奇偶进行分组，基本运算就是蝶形运算。 对于输入序列点数为N=2L，共有L级蝶形，每级都由N/2个蝶形运算组成。 以N=8时，基2的按时间抽选（输入倒序，输出顺序）的FFT流图为例： 三 实验内容 1 编程实现基2的按时间抽选（输入倒序，输出顺序）的FFT。 2 利用FFT对给定正弦信号进行频谱分析，信号截断采用矩形窗。 给定正弦序列，分别用DFT和FFT求x(n)的频谱, 比较二者的效率，可用运行时间或运算量（程序中）表示。 3 再分别用凯泽窗和哈明窗对给定信号进行频谱分析。 四 实验报告要求 分析DFT和FFT求频谱的区别。 分析两者求频谱的效率。 对于不同类型序列，用DFT FFT求频谱，分析结果。 验收程序。 实验程序： clear; M=30; F=5; N=70; T=1000; n=0:N-1; x=sin(F*2*pi*n/N); window=[boxcar(M),zeros(1,N-M)]; %window=[kaiser(M),zeros(1,N-M)]; %凯泽窗 %window=[hamming(M),zeros(1,N-M)]; %哈明窗 t=-T/2:(T-1)/2; subplot(2,2,1); stem(x); title(正弦序列); subplot(2,2,2); y=ff