实测数据轨道谱估计MATLAB程序

周期图法的功率谱图(直接法)
function [w]=tt()  %函数名称
l=1;    %波长下限
v=150/3.6;   %车速
fs=2*v/l;   %采样频率
f=fopen('0412.txt','r');   %打开04年12月轨道高低不平顺实测数据
y=fscanf(f,'%e');
%读入原始数据
nfft=2^15;    %FFT的点数
Xk=fft(y,nfft);    %快速傅立叶变换
%save fft.txt Xk -ascii;   %存储傅立叶变换值
p=abs(Xk).^2/length(Xk);   %计算功率谱
%save psd.txt p -ascii;   %存储功率谱值
index=0:round(nfft/2-1);   %时间坐标取一半
k=index.*fs/nfft;     %对应的频率
figure本文来自优.文,论-文·网原文请找腾讯752018766
w=loglog(k,p(index+1));    %画频率和功率谱密度的双对数坐标图
title('周期图法谱估计','fontsize',25)    %标题名,字号为25
xlabel('频率(Hz)','fontsize',20);   %x轴名称
grid onstm32检测不到RXNE标志
ylabel('功率谱密度(mm2)','fontsize',20);   %y轴名称
1.2.2  自相关法(间接法)
function [w]=indirect()   %函数名
l=1;   %波长下限
v=150/3.6;   %车速
fs=2*v/l;   %采样频率
f=fopen('0412.txt','r');   %打开04年12月轨道高低不平顺实测数据
y=fscanf(f,'%e');
%读取原始数据
nfft=2^15;    %FFT的点数
R=xcorr(y,'unbiased');     %自相关函数
F=fft(R,nfft);     %快速傅立叶变换
p=abs(F);     %求解功率谱密度
index=0:round(nfft/2-1);     %时间坐标取一半
k=index*fs/nfft;     %对应的频率
figure      %新窗口
w=loglog(k,p(index+1));     %画频率和功率谱密度的双对数坐标图
title('自相关法','fontsize',25)    %标题名,字号为25
xlabel('频率(Hz)','fontsize',20);     %x轴名称
grid on      %加网格2635 实测数据轨道谱估计MATLAB程序下载如图片无法显示或论文不完整，请联系qq752018766