BT法求功率谱程序BT法求功率谱程序
%自相关函数
clc; %清空
clear all;%清除所有变量
close all;%关闭所有窗口
load('C:\Users\caolili\Desktop\FileDat01_2.mat')
fs=4096;%设定采样频率
N=4096; %采样点数
n=0:N-1;
t=n/fs;
Lag=300;%延迟样点数
[c,lags]=xcorr(s,Lag,'unbiased');%对原始信号进行无偏自相关估计
subplot(1,2,1); plot(n,s);...
BT法求功率谱程序
%自相关函数
clc; %清空
clear all;%清除所有变量
close all;%关闭所有窗口
load('C:\Users\caolili\Desktop\FileDat01_2.mat')
fs=4096;%设定采样频率
N=4096; %采样点数
n=0:N-1;
t=n/fs;
Lag=300;%延迟样点数
[c,lags]=xcorr(s,Lag,'unbiased');%对原始信号进行无偏自相关估计
subplot(1,2,1); plot(n,s);%绘制原信号的时域波形 xlabel('时间t(1/4096s)'); ylabel('信号s');
title('带噪声的信号波形');
grid on;
subplot(1,2,2); plot(lags/fs,c);%绘制自相关函数图象 xlabel('时间t');
ylabel('Rx(t)'); title('带噪声的信号自相关函数'); grid on;
%功率谱
clc; %清空
clear all; %清除所有变量
close all; %关闭所有窗口
load('C:\Users\caolili\Desktop\FileDat01_2.mat')
fs=4096;%设定采样频率
N=4096; %采样点数
n=0:N-1;
t=n/fs;
Lag=300;%延迟样点数
[c,lags]=xcorr(s,Lag,'unbiased');%求信号的自相关函数
fy=fft(c,N); %对自相关函数做FFT变换 t1=0:round(N/2-1); a=t1*fs/N;
P=10*log10(fy(t1+1));%纵坐标为相对功率谱密度，单位dB/Hz
figure(gcf);
plot(a,P);
ylabel('功率谱密度dBw/Hz'); title('信号的功率谱');
grid;
本文档为【BT法求功率谱程序】，请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑，
图片更改请在作品中右键图片并更换，文字修改请直接点击文字进行修改，也可以新增和删除文档中的内容。
该文档来自用户分享，如有侵权行为请发邮件ishare@vip.sina.com联系网站客服，我们会及时删除。
[版权声明] 本站所有资料为用户分享产生，若发现您的权利被侵害，请联系客服邮件isharekefu@iask.cn，我们尽快处理。
本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权，请谨慎使用。
网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传，仅限个人学习分享使用，禁止用于任何广告和商用目的。
浙公网安备 33021202002483