matlab程序中功率谱分析的经典常用方法

一、直接法清除；clc全部关闭；%清除变量；清屏；关闭当前图形窗口Fs=1000 .t=0:1/Fs:1 .nfft=2048 .%改变nfft的值可对比不同采样值时的谱估计效果% * * * * * * * * * * * * * * * * * * * * *生成信号、噪声* * * * * * * * * * * x1=cos(2 * pi * 40 * t)3 * cos(2 * pi * 45 * t)；%信号x2=随机(大小(t)；%噪声x3=x1 x2%信号噪声Pxx，f=周期图(x3，窗口，nfft，Fs)；%直接法曲线图(f，10 * log 10(Pxx)；标题（直接法nfft=2048)；set(gca，xlim，1 120)；伊拉贝尔(调幅/分贝)；xlabel(频率/赫兹)；二、间接法Fs=1000 .%采样频率n=0:1/Fs:1 .%产生含有噪声的序列x1=cos(2 * pi * 40 * n)3 * cos(2 * pi * 45 * n)；%信号x2=随机(大小(n)；%噪声x3=x1 x2%信号噪声nfft=1024 .cxn=xcorr(x3)；%计算序列的自相关函数CXk=fft(cxn).Pxx=abs(CXk).指数=0:圆(nfft/2-1)；f=索引*Fs/nfft .plot_Pxx=10*log10(Pxx(索引1)；图(1)图(f，plot _ Pxx)；标题（间接法nfft=1024)；伊拉贝尔(调幅/分贝)；set(gca，xlim，1 120)；xlabel(频率/赫兹)；三、巴特利特法清除；clc全部关闭；%清除变量；清屏；关闭当前图形窗口Fs=1000 .t=0:1/Fs:1 .nfft=1024 .% * * * * * * * * * * * * * * * * * * * * *生成信号、噪声* * * * * * * * * * * * * * * * * * * * * * *x1=cos(2 * pi * 40 * t)3 * cos(2 * pi * 45 * t)；%信号x2=随机(大小(t)；%噪声x3=x1 x2%信号噪声窗口=汉明(512)；%海明窗no verlap=0；%数据无重叠p=0.9%置信概率Pxx，Pxxc=psd(x3，nfft，Fs，window，noverlap，p)；指数=0:圆(nfft/2-1)；k=索引*Fs/nfft .plot_Pxx=10*log10(Pxx(索引1)；plot_Pxxc=10*log10(Pxxc(索引1)；图(1)图(k，plot _ Pxx)；标题（巴特利特法海明窗)；set(gca，xlim，1120)；伊拉贝尔(调幅/分贝)；xlabel(频率/赫兹)；四韦尔奇法清除；clccloseall%清除变量；清屏；关闭当前图形窗口Fs=1000 .t=0:1/Fs:1 .nfft=1024 .% * * * * * * * * * * * * * * * * * * * * *生成信号、噪声* * * * * * * * * * * * * * * * * * * * * * *x1=cos(2 * pi * 40 * t)3 * cos(2 * pi * 45 * t)；%信号x2=随机(大小(t)；%噪声x3=x1 x2%信号噪声窗口=汉明(512)；%海明窗noverlap=128范围=单侧；%频率间隔为0Fs/2，只计算一半的频Pxx1，f=pwelch(x3，窗口、noverlap、nfft、Fs、范围)；plot _ Pxx 1=10 * log 10(Pxx 1)；图(1)；曲线图(f，曲线图_ Pxx 1)；标题（韦尔奇法海明窗)；伊拉贝尔(调幅/分贝)；set(gca，xlim，1120)；xlabel(频率/赫兹)；对所给的实测信号进行谱估计，本文采用了周期图法和韦尔奇法。其程序如下所示：一、周期图法清除；clc全部关闭；%清除变量；清屏；关闭当前图形窗口Fs=2000 .%采样频率加载8夏2数据。垫子；x1=香奈儿8夏2灰；x2=香奈儿8夏2 STIMlag=长度(x1).% * * * * * * * * * * * * * * * * * * * * *求功率谱密度* * * * * * * * * * * * * * * * * * * * * * *Pxx1，f=周期图(x1，窗口，长度(x1)，Fs)；%直接法Pxx2，f=周期图(x2，窗口，长度(x1)，Fs)；plot _ Pxx 1=10 * log 10(Pxx 1)；plot _ Pxx 2=10 * log 10(Pxx 2)；% * * * * * * * * * * * * * * * * * * * * *显示功率谱密度曲线* * * * * * * * * * * * * * * * * * * * * * *图(1)；曲线图(f，曲线图_Pxx1，b)；坐标轴(0，100，-50，40)；网格打开；标题（直接法)；set(gca，xlim，1 100)；伊拉贝尔(调幅/分贝)；坚持住。图(f，plot_Pxx2，r)；坐标轴(0，100，-50，40)；网格打开；标题（直接法)；xlabel(频率/赫兹)；伊拉贝尔(调幅/分贝)；二韦尔奇法清除；clc全部关闭；%清除变量；清屏；关闭当前图形窗口Fs=2000 .加载8夏2数据。垫子；x1=香奈儿8夏2灰；x2=香奈儿8夏2 STI窗口1=汉明(1024)；%海明窗noverlap=256%数据无重叠范围=单侧；%频率间隔为0 Fs/2，只计算一半频、% * * * * * * * * * * * * * * * * * * * * *求功率谱密度*% Pxx1，f=pwelch(x1，窗口1，noverlap，长度(x1)，Fs，范围)；Pxx2，f=pwelch(x2，窗口1，noverlap长度(x1)，Fs，范围)；plot _ Pxx 1=10 * log 10(Pxx 1)；plot _ Pxx 2=10 * log 10(Pxx 2)；% * * * * * * * * * * * * * * * * * * * * *显示功率谱密度曲线*%图(1)；曲线图(f，