毕业论文论文范文课程设计实践报告法律论文英语论文教学论文医学论文农学论文艺术论文行政论文管理论文计算机安全
您现在的位置: 毕业论文 >> 报告 >> 正文

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

更新时间:2014-5-16:  来源:毕业论文
周期图法的功率谱图(直接法)
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
设为首页 | 联系站长 | 友情链接 | 网站地图 |

copyright©lwfree.cn 六维论文网 严禁转载
如果本毕业论文网损害了您的利益或者侵犯了您的权利,请及时联系,我们一定会及时改正。