噪声中正弦信号的经典法频谱分析
声音信号的频谱分析与频率测量方法
声音信号的频谱分析与频率测量方法声音是我们日常生活中不可或缺的一部分，我们通过声音来交流、表达情感，甚至通过声音来判断事物的性质。
然而，声音是如何产生的？我们如何对声音进行分析和测量呢？本文将介绍声音信号的频谱分析与频率测量方法。
声音信号是由空气中的振动引起的，当物体振动时，会产生压力波，通过空气传播出去，我们就能听到声音。
声音信号可以通过振动的频率和振幅来描述，其中频率是指振动的周期性，而振幅则是指振动的强度。
频谱分析是一种将声音信号分解成不同频率成分的方法。
它可以帮助我们了解声音信号的频率分布情况，从而更好地理解声音的特性。
频谱分析的基本原理是将声音信号转换为频域表示，即将信号从时域转换为频域。
这可以通过傅里叶变换来实现。
傅里叶变换是一种将时域信号转换为频域信号的数学方法。
它将信号分解成一系列正弦波的叠加，每个正弦波都有不同的频率和振幅。
通过傅里叶变换，我们可以得到声音信号的频谱图，从而了解声音信号中不同频率成分的贡献程度。
频谱图通常以频率为横轴，振幅或能量为纵轴，通过不同的颜色或灰度表示不同频率成分的强度。
频谱图可以直观地展示声音信号的频率分布情况，帮助我们分析声音的特性。
例如，在音乐领域，频谱分析可以用来研究音乐的音色特点，判断乐器的类型等。
除了频谱分析，频率测量是对声音信号进行定量分析的重要方法。
频率是声音信号中最基本的特征之一，它决定了声音的音调高低。
频率测量可以通过多种方法实现，其中一种常用的方法是自相关法。
自相关法是一种基于信号自身的周期性特点进行频率测量的方法。
它通过计算信号与自身的延迟版本之间的相似程度来确定信号的周期性。
具体而言，自相关法将信号与其自身进行延迟，然后计算它们之间的相关性。
通过寻找最大相关性的延迟值，我们可以得到信号的主要频率成分。
除了自相关法，还有一些其他的频率测量方法，如峰值检测法、零交叉法等。
这些方法在不同的应用场景下有着各自的优势和适用性。
例如，峰值检测法适用于测量周期性信号的频率，而零交叉法适用于测量非周期性信号的频率。
一种提取噪声中正弦信号的总体最小二乘法_喻胜
∑i =1 ∑j =1 |X ij -Y ij |
( E)
m
n
/2 2 1
( 15)
文献[ 4] 的定理二指出 : 在以上( 15) 式定义的 Frobenius 量最小的意义下 , 与( 14) 式实际观测 矩阵 H 最接近的秩为 2K 的增强数据矩阵 H 可按如下的迭代方法求出( 算法 2) : ① 将 Hankel 数据矩阵 H 作奇异值分解展开并取前 2K 个最大奇异值而将其余奇异值 置0 : H(E) = ② 将( 16) 式所得 H
2k
-k )=0 n =2K +1 , …, N
( 4) ( 5) ( 6) ( 7)
式中的线性预测系数{ ak } 具有如下的中心元素对称性 : a 2K -k =a k k = 1 , … , K 且 a 0 =a 2K = 1 若能由观测序列求出{ ak } 的估计值 , 则以下 z 的多项式 Χ ( z )=
2 参数 。 这里假定 e ( n) 为具有零均值 、 方差为 σ 的宽带噪声 。 将纯谐波过程
本文于 1999 年 2 月收到 。 喻 胜 : 博士研究生 ; 闫 波 : 讲师 ; 陈光
: 教授 , 博导 。
第 2 期 一种提取噪声中正弦信号的总体最小二乘法 · 7 · y( n)= 写作 y( n) =
k
ω kn + k ) ∑k =1 Ak cos(
n
k
( 2) ( 3)
= 1 ( bkzk ∑ k
+bk zk ) n = 1, …, N
＊
＊
式中 bk = ( 1/2) Ak exp( j k) ; zk =exp( jω k) 。 很明显 , 当被测过程的各频率量 ω k 求出后 , 以上 N 个方程组成的方程组( 3) 是复振幅 bk 的线性方程 , 不难由传统的线性最小二乘法求出诸 { bk } 。 所以 , 本文的讨论仅限于由观测序列{ x( n) } 估计各正弦信号频率 。 我们知道 , 对任意频率的 K 个纯净实正弦信号( 2) 可以用一个 AR( 2K ) 线性预测模型来 描述 : n ∑k =0 aky(
用FFT对信号作频谱分析
用FFT对信号作频谱分析快速傅立叶变换（FFT）是一种在信号处理中常用于频谱分析的方法。
它是傅立叶变换的一种快速算法，通过将信号从时间域转换到频域，可以提取信号的频率信息。
FFT算法的原理是将信号分解为不同频率的正弦波成分，并计算每个频率成分的幅度和相位。
具体而言，FFT将信号划分为一系列时间窗口，每个窗口内的信号被认为是一个周期性信号，然后对每个窗口内的信号进行傅立叶变换。
使用FFT进行频谱分析可以得到信号的频率分布情况。
频谱可以显示信号中各个频率成分的强度。
通过分析频谱可以识别信号中的主要频率成分，判断信号中是否存在特定频率的干扰或噪声。
常见的应用包括音频信号处理、图像处理、通信系统中的滤波和解调等。
使用FFT进行频谱分析的步骤如下：1.首先，获取待分析的信号，并确保信号是离散的，即采样频率与信号中的最高频率成分满足奈奎斯特采样定理。
2.对信号进行预处理，包括去除直流分量和任何不需要的干扰信号。
3.对信号进行分段，分段后的每个窗口长度在FFT算法中通常为2的幂次方。
常见的窗口函数包括矩形窗、汉明窗等。
4.对每个窗口内的信号应用FFT算法，将信号从时间域转换到频域，并计算每个频率成分的幅度和相位。
5.对所有窗口得到的频谱进行平均处理，以得到最终的频谱分布。
在使用FFT进行频谱分析时需要注意的问题有：1.噪声的影响：FFT对噪声敏感，噪声会引入幅度偏差和频率漂移。
可以通过加窗等方法来减小噪声的影响。
2.分辨率的选择：分辨率是指在频谱中能够分辨的最小频率间隔。
分辨率与信号长度和采样频率有关，需要根据需求进行选择。
3.漏泄效应：当信号中的频率不是FFT长度的整数倍时，会出现漏泄效应。
可以通过零填充等方法来减小漏泄效应。
4.能量泄露：FFT将信号限定在一个周期内进行计算，如果信号过长，则可能导致部分频率成分的能量泄露到其他频率上。
总之，FFT作为信号处理中常用的频谱分析方法，能够提取信号中的频率信息，广泛应用于多个领域。
利用Matlab绘制正弦信号的频谱图并做相关分析
利用Matlab绘制正弦信号的频谱图并做相关分析一、作业要求：1、信号可变（信号的赋值、相位、频率可变）；2、采样频率fs可变；3、加各种不同的窗函数并分析其影响；4、频谱校正；5、频谱细化。
二、采用matlab编写如下程序：clear;clf;fs=100;N=1024; %采样频率和数据点数A=20;B=30;C=0.38;n=0:N-1;t=n/fs; %时间序列x=A*sin(2*pi*B*t+C); %信号y=fft(x,N); %对信号进行傅里叶变换yy=abs(y); %求得傅里叶变换后的振幅yy=yy*2/N; %幅值处理f=n*fs/N; %频率序列subplot(3,3,1),plot(f,yy); %绘出随频率变化的振幅xlabel('频率/\itHz');ylabel('振幅');title('图1：fs=100，N=1024');grid on;%两种信号叠加，x=A*sin(2*pi*B*t+C)+2*A*sin(2*pi*1.5*B*t+2.5*C); %信号y=fft(x,N); %对信号进行傅里叶变换yy=abs(y); %求得傅里叶变换后的振幅yy=yy*2/N; %幅值处理f=n*fs/N; %频率序列subplot(3,3,2),plot(f,yy); %绘出随频率变化的振幅xlabel('频率/\itHz');ylabel('振幅');title('图2：fs=100,N=1024，两种信号叠加');grid on;%加噪声之后的图像x=A*sin(2*pi*B*t+C)+28*randn(size(t));y=fft(x,N);yy=abs(y);yy=yy*2/N; %幅值处理subplot(3,3,3),plot(f(1:N/2.56),yy(1:N/2.56));xlabel('频率/\itHz');ylabel('振幅');title('图3：fs=100,N=1024混入噪声');grid on;%改变采样点数N=128N=128;n=0:N-1;t=n/fs; %时间序列x=A*sin(2*pi*B*t+C); %信号y=fft(x,N); %对信号进行傅里叶变换yy=abs(y); %求得傅里叶变换后的振幅yy=yy*2/N; %幅值处理f=n*fs/N; %频率序列subplot(3,3,4),plot(f(1:N/2.56),yy(1:N/2.56)); %绘出随频率变化的振幅xlabel('频率/\itHz');ylabel('振幅');title('图4：fs=100，N=128');grid on;%改变采样频率为200Hz时的频谱fs=400;N=1024;n=0:N-1;t=n/fs;x=A*sin(2*pi*B*t+C); %信号y=fft(x,N); %对信号进行快速傅里叶变换yy=abs(y); %求取傅里叶变换的振幅yy=yy*2/N; %幅值处理f=n*fs/N;subplot(3,3,5),plot(f(1:N/2.56),yy(1:N/2.56)); %绘出随频率变化的振幅xlabel('频率/\itHz');ylabel('振幅');title('图5：fs=400，N=1024');grid on;%加三角窗函数fs=100;N=1024; %采样频率和数据点数n=0:N-1;t=n/fs; %时间序列x=A*sin(2*pi*B*t+C); %信号window=triang(N);%生成三角窗函数x=x.*window';%加窗函数y=fft(x,N); %对信号进行傅里叶变换yy=abs(y); %求得傅里叶变换后的振幅yy=yy*2/N; %幅值处理f=n*fs/N; %频率序列subplot(3,3,6),plot(f(1:N/2.56),2*yy(1:N/2.56)); %绘出随频率变化的振幅xlabel('频率/\itHz');ylabel('振幅');title('图6：fs=100,N=1024,加三角窗函数');grid on;%加海明窗函数后的频谱fs=100;N=1024;n=0:N-1;t=n/fs;x=A*sin(2*pi*B*t+C); %信号window=hamming(N);%生成海明窗函数x=x.*window';%加窗函数y=fft(x,N); %对信号进行快速傅里叶变换yy=abs(y); %求取傅里叶变换的振幅yy=yy*2/N; %幅值处理f=n*fs/N;subplot(3,3,7),plot(f(1:N/2.56),1.852*yy(1:N/2.56)); %绘出随频率变化的振幅xlabel('频率/\itHz');ylabel('振幅');title('图7：fs=100，N=1024,加海明窗函数');grid on;%加汉宁窗函数后的频谱fs=100;N=1024;n=0:N-1;t=n/fs;x=A*sin(2*pi*B*t+C); %信号window=hanning(N);%生成汉宁窗函数x=x.*window';%加窗函数y=fft(x,N); %对信号进行快速傅里叶变换yy=abs(y); %求取傅里叶变换的振幅yy=yy*2/N; %幅值处理f=n*fs/N;subplot(3,3,8),plot(f(1:N/2.56),2*yy(1:N/2.56)); %绘出随频率变化的振幅xlabel('频率/\itHz');ylabel('振幅');title('图8：fs=100，N=1024,加汉宁窗函数');grid on;三、运行结果如下：四、分析与结论：1）从所做图像可以看出，信号的幅值均小于真实值，说明在截断信号时存在泄露。
几种基于FFT的频率估计方法精度分析_齐国清
的是 ΡCRB 与 FFT 的频率分辨率 ∃f 的比值。当N =
2 048, 信噪比为 40 dB 时, 标准差下限约为 10- 4∃f ;
当 N = 256, 信 噪 比 为 10 dB 时 标 准 差 下 限 约 为
10- 2∃f 。
∆δ =
±
1
Α +
Α
(4)
式 (4) 中的符号根据次大值是在最大值的左侧还是
摘要: 介绍了噪声背景中正弦信号频率估计的方差下限, 对利用 FFT 主瓣内两条幅度最大谱线进行插值的频率估 计方法 (R ife2J ane 方法和 Q u inn 方法) 以及利用 FFT 相位进行频率插值的方法 (分段 FFT 相位差法和重叠 FFT 相位差法) 的方差进行了理论分析, 推导了 Q u inn 方法的频率估计方差计算公式, 提出了通过滤波进一步提高分段 FFT 相位差法的频率估计精度的方法。通过计算机M on te Carlo 模拟实验对上述各种方法的频率估计精度以及加 窗函数的影响进行了分析并与理论下限进行了比较, 指出了每种方法所能达到的估计精度。
介绍了噪声背景中正弦信号频率估计的方差下限对利用fft主瓣内两条幅度最大谱线进行插值的频率估计方法rife2jane方法和quinn方法以及利用fft相位进行频率插值的方法分段fft相位差法和重叠fft相位差法的方差进行了理论分析推导了quinn方法的频率估计方差计算公式提出了通过滤波进一步提高分段fft相位差法的频率估计精度的方法
引 起较大的频率估计误差。 图 2 (a) 所示为 N =
图 1 频率估计标准差下限随信噪比和 FFT 长度变化情况
1. 2 R ife-Jane 频率估计方法
s (n) 的 N 点 FFT 记为 S (k ) , 鉴于实序列 F F T 的对称性, 只考虑离散频谱的前 N 2 点, 在不加窗 (即加矩形窗) 情况下有
频谱分析(完整版)
Matlab 信号处理工具箱 帮助文档 谱估计专题翻译：无名网友 & Lyra频谱分析Spectral estimation （谱估计）的目标是基于一个有限的数据集合描述一个信号的功率（在频率上的）分布。
功率谱估计在很多场合下都是有用的，包括对宽带噪声湮没下的信号的检测。
从数学上看，一个平稳随机过程n x 的power spectrum （功率谱）和correlation sequence （相关序列）通过discrete-time Fourier transform （离散时间傅立叶变换）构成联系。
从normalized frequency （归一化角频率）角度看，有下式()()j mxx xxm S R m eωω∞-=-∞=∑注：()()2xx S X ωω=，其中()/2/2lim N j n n N N X x e ωω=-=∑πωπ-<≤。
其matlab 近似为X=fft(x,N)/sqrt(N)，在下文中()L X f 就是指matlab fft 函数的计算结果了使用关系2/s f f ωπ=可以写成物理频率f 的函数，其中s f 是采样频率()()2/sjfm f xx xxm S f R m eπ∞-=-∞=∑相关序列可以从功率谱用IDFT 变换求得：()()()/22//22sss f jfm f j m xx xx xx sf S e S f e R m d df f πωππωωπ--==⎰⎰序列n x 在整个Nyquist 间隔上的平均功率可以表示为()()()/2/202ss f xx xx xx sf S S f R d df f ππωωπ--==⎰⎰ 上式中的()()2xx xx S P ωωπ=以及()()xx xx sS f P f f =被定义为平稳随机信号n x 的power spectral density (PSD)（功率谱密度） 一个信号在频带[]1212,,0ωωωωπ≤<≤上的平均功率可以通过对PSD 在频带上积分求出从上式中可以看出()xx P ω是一个信号在一个无穷小频带上的功率浓度，这也是为什么它叫做功率谱密度。
基于matlab的正弦信号叠加噪声信号的频谱分析(实验报告)
基于MATLAB的正弦信号叠加噪声信号的频谱分析——《随机信号分析基础》课题设计报告学院：弘深学院班级：电子信息实验班学号： 20136927姓名：文政指导老师：欧翔2015年4月3日基于MATLAB的正弦信号叠加噪声信号的频谱分析目录一、课题目的 (2)二、课题要求 (2)三、设计原理 (2)1.瑞丽分布 (2)2.工具软件使用 (3)四、实验过程 (4)1.产生瑞丽分布的声音信号 (4)2.绘制直方图、均值图和方差图 (4)3.绘制噪声信号的密度函数图 (5)4.绘制自相关函数图 (5)5.绘制幅频特性图、相频特性图和频谱图 (5)6.产生频率为10Hz的正弦信号，并绘制合成图像 (6)7.改变sigma大小，分析噪声方差对正弦函数的影响 (7)五、实验结果及分析 (7)六、结论 (13)附录 (14)(MATLAB 源程序代码) (14)一、课题目的1.熟悉MATLAB语言的基本语法；2.掌握MATLAB语言中求信号的期望、方差及自相关函数的方法；3.掌握信号和噪声产生及叠加的方法；4.掌握MATLAB语言中信号自相关函数图和频谱图的绘制方法二、课题要求1.产生一个符合瑞利分布的噪声信号并绘制出其直方图、均值图、方差图；2.绘制出所产生的噪声信号的概率密度函数图及自相关函数图；3.绘制出所产生的噪声信号的幅频特性、相频特性图以及频谱图；4.产生一个频率为10Hz的正弦信号并绘制出图形；5.将噪声信号叠加到所产生的正弦信号上并绘制出图形；6.分析改变方差时，噪声对正弦信号产生的影响，以图举例说明。
三、设计原理1.瑞丽分布瑞利分布是一种最常见的用于描述平坦衰落信号接收包络或独立多径分量接受包络统计时变特性的一种分布类型，它是在高斯正态分布的基础上，优化组合得到的。
设相互独立的随机变量ξ、η服从一维正态分布ξ~N(0,σ2), η~N(0,σ2);其密度函数为f(ξ)=−ξ22σ2;f(η)=2πσ−η22σ2;则其联合分布密度函数为f(ξ,η)=f(ξ)∙f(η)=12πσ2e−ξ2+η22σ2;其分布密度函数图像为：图1 二维联合正态分布概率密度图若将ξ−η平面坐标系转为r−θ极坐标系，则r=√ξ2+η2, tanθ=ξηr与θ相互独立，则有：f(r)=rσe−r22σ2∙u(r)即为随机变量r的瑞利分布概率密度函数，其图像为图2 瑞利分布概率密度函数图2.工具软件使用本文使用MATLAB（使用版本MATLAB R2014b）软件对服从瑞利分布的噪声进行相关的分析处理。
[重点]对正弦信号的采样频谱分析
![[重点]对正弦信号的采样频谱分析](https://img.taocdn.com/s3/m/4f2f80ccac51f01dc281e53a580216fc700a530f.png)
一、题目要求：给定采样频率fs，两个正弦信号相加，两信号幅度不同、频率不同。
要求给定正弦信号频率的选择与采样频率成整数关系和非整数关系两种情况，信号持续时间选择多种情况分别进行频谱分析。
二、题目原理与分析：本题目要对正弦信号进行抽样，并使用fft对采样信号进行频谱分析。
因此首先对连续正弦信号进行离散处理。
实际操作中通过对连续信号间隔相同的抽样周期取值来达到离散化的目的。
根据抽样定理，如果信号带宽小于奈奎斯特频率（即采样频率的二分之一），那么此时这些离散的采样点能够完全表示原信号。
高于或处于奈奎斯特频率的频率分量会导致混叠现象。
设抽样周期为TS(抽样角频率为ωS)，则可见抽样后的频谱是原信号频谱的周期性重复，当信号带宽小于奈奎斯特频率的二分之一时不会产生频谱混叠现象。
因此，我们对采样频率的选择采取fs>2fo,fs=2fo,fs<2fo三种情况进行分析。
对信号采样后，使用fft函数对其进行频谱分析。
为了使频谱图像更加清楚，更能准确反映实际情况并接近理想情况，我们采用512点fft。
取512点fft不仅可以加快计算速度，而且可以使频谱图更加精确。
若取的点数较少，则会造成频谱较大的失真。
三、实验程序：本实验采用matlab编写程序，实验中取原信号为ft=sin(2πfXt)+2sin(10πfXt)，取频率f=1kHz，实验程序如下：f=1000;fs=20000;Um=1;N=512;T=1/fs;t=0:1/fs:0.01;ft=Um*sin(2*pi*f*t)+2*Um*sin(10*pi*f*t);subplot(3,1,1);plot(t,ft);grid on;axis([0 0.01 1.1*min(ft) 1.1*max(ft)]);xlabel('t'),ylabel('ft');title('抽样信号的连续形式');subplot(3,1,2);stem(t,ft);grid on;axis([0 0.01 1.1*min(ft) 1.1*max(ft)]);xlabel('t'),ylabel('ft');title('实际抽样信号');k=0:N-1;Fw=fft(ft,N);subplot(3,1,3);plot(k,abs(Fw));grid on;axis([0 550 -0.2 65*pi]);title('抽样信号幅度谱')在实际操作过程中，对于信号频率与采样频率所成整数倍与非整数倍关系时，信号持续时间不同时，只需改变程序中的相关语句即可。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实验报告
一、实验名称
噪声中正弦信号的经典法频谱分析
二、实验目的
通过对噪声中正弦信号的经典法频谱分析，来理解和掌握经典谱估计的知识，以及学会应用经典谱估计的方法。
三、基本原理
1．周期图法：又称直接法。
把随机信号)(n x 的N 点观察数据)(n x N 视为一能量有限信号，直接取)(n x N 的傅里叶变换，得)(jw N e X ，然后再取其幅值的平方，并除以N ，作为对)(n x 真
实的功率谱)(jw e P 的估计，以)(ˆjw PER
e P 表示用周期图法估计出的功率谱，则2)(1)(ˆw X N
w P n
PER =。
2．自相关法：又称为间接法功BT 法。
先由)(n x N 估计出自相关函数)(ˆm r
，然后对)(ˆm r 求傅里叶变换得到)(n x N 的功率谱，记之为)(ˆw P BT
，并以此作为对)(w P 的估计，即1,)(ˆ)(ˆ-≤=--=∑N M e
m r w P jwm
M
M
m BT。
3．Bartlett 法：对L 个具有相同的均值μ和方差2σ的独立随机变量1X ，2X ，…，L X ，新随机变量L X X X X L /)(21+++= 的均值也是μ，但方差是L /2σ，减小了L 倍。
由此得
到改善)(ˆw P PER
方差特性的一个有效方法。
它将采样数据)(n x N 分成L 段，每段的长度都是M ，即N=LM ，第i 段数据加矩形窗后，变为L i e n x
M
w x
M n jwn i N
I
PER ≤≤=∑-=-1,)(1)(ˆ2
10。
把)(ˆw P PER
对应相加，再取平均，得到平均周期图2
1110
)(1)(ˆ1)(∑∑∑==-=-==L i L i M n jwn i
N i PER PER e n x ML w P L w P 。
4．Welch 法：它是对Bartlett 法的改进。
改进之一是，在对)(n x N 分段时，可允许每一段的数据有部分的交叠。
改进之二是，每一段的数据窗口可以不是矩形窗口，例如使用汉宁窗或汉明窗，记之为)(2n d 。
这样可以改善由于矩形窗边瓣较大所产生的谱失真。
然后按Bartlett
法求每一段的功率谱，记之为)(ˆw P P E R
，即210
2)()(1)(ˆ∑-=-=M n j w n
i N
i PER e
n d n x
MU
w P ，式中
∑-==
10
22
)(1M n n d
M
U 是归一化因子。
四、主要编程步骤
（一）构造一个频率为100Hz 的正弦信号，再构造一个方差为0.01的高斯白噪声，将正弦信号与噪声相加，得到信号xn=sin(2*pi*100*n)+0.1*randn(size(n))； （二）进行功率谱估计 1．周期图法
① 对xn 加汉明窗,调用函数[Pxx,f]=periodogram(xn,window,nfft,Fs)，其中采样频率Fs=1000Hz ，分别取不同的信号长度128，256，512，1024时，分别进行功率谱估计，画功率谱图。
② 对xn 分别加矩形窗，汉宁窗，海明窗，blackman 窗。
固定信号长度为256，采样频率为1000Hz, 调用函数[Pxx,f]=periodogram(xn,window,nfft,Fs), 分别进行功率谱估计，画功率谱图。
2．间接法
用cxn=xcorr(xn,'unbiased')来计算xn 的自相关函数，然后对其进行傅里叶变换，便得到它的功率谱，画图。
3．Bartlett 法 对xn 加矩形窗，调用函数[Pxx,Pxxc]=psd(xn,nfft,Fs,window,noverlap,p)，进行功率谱估计，画功率谱图。
4．Welch 法
对xn 加海明窗，调用函数[Pxx,f]=pwelch(xn,window,noverlap,nfft,Fs,range)，进行功率谱估计，画功率谱图。
五、实验结果及分析
1．周期图法，
采样频率Fs=1000Hz ，信号频率f=100Hz ，高斯白噪声方差为0.01。
① 取不同的信号长度128，256，512，1024时 程序运行结果如下图显示：
分析：当N逐渐增大时，曲线的起伏也逐渐加剧，这是因为N增大，使互不相关的点增多，才导致了曲线起伏的加剧。
②对xn分别加矩形窗，汉宁窗，海明窗，blackman窗。
固定信号长度N=256，调用函数[Pxx,f]=periodogram(xn,window,nfft,Fs), 分别进行功率谱估计，画功率谱图。
程序运行结果如下图显示：
③改变信噪比
固定信号长度N=512,信噪比为-20dB,0dB,20dB时
分析：信噪比越大，效果越好。
④改变窗函数的长度
2．间接法
采样频率Fs=1000Hz,采样点数N=1024，信号频率=100Hz，高斯白噪声方差为0.01。
3.Bartlett法
采样频率Fs=1000Hz,采样点数N=1024，信号频率=100Hz，高斯白噪声方差为0.01。
4．Welch法
采样频率Fs=1000Hz,采样点数N=1024，信号频率=100Hz，高斯白噪声方差为0.01。
六、结论
根据分析结果，给出明确的，与分析相一致的结论。