细化谱
data=load('D:\私人\毕业设计\实验数据\7月份数据\滚动体点蚀\滚动体点蚀6k\07519102360001144321.txt');long=length(data);
n=1:8192;
fs=6000*2.56;
t=n/fs;
%时域图
figure(1);
plot(t,data);
title('时域图与频谱图');
xlabel('时间(s)');
ylabel('幅值(m/s^2)');
%幅值谱(频谱)
nfft=8192;
fuzhi=abs(fft(data,nfft));
figure(2);
plot((0:nfft/2-1)/nfft*fs,fuzhi(1:nfft/2));
xlabel('频率(HZ)');
ylabel('幅值(m/s^2)');
title('频谱');
%移频
yipin=data.*exp(-2*j*pi*0*t);%移频后时域
pin=abs(fft(yipin,nfft));%移频后频域
plot((0:nfft/2-1)/nfft*fs,pin(1:nfft/2));
title('信号移频');
%滤波器
wp=0.065*pi;%通带截止频率
%过度频带意义和作用?
ws=0.09*pi;%阻带截止频率
rs=40;%阻带最小衰减
DB=ws-wp;
beta=0.5842*(rs-21)^0.4+0.07886*(rs-21);
M=ceil((rs-8)/2.285/DB);
wc=(wp+ws)/2/pi;
hn=fir1(M,wc,kaiser(M+1,beta));
y=fftfilt(hn,yipin);%滤波后时域波形
zhihou=abs(fft(y,nfft));%滤波后频域波形
figure(3);
subplot(211);
plot((0:nfft/2-1)/nfft*fs,zhihou(1:nfft/2));
title('信号滤波');
%重采样
ccy=resample(y,1,10);%重采样时域放大倍数
ccyp=abs(fft(ccy,nfft));%重采样频域
subplot(212);
plot((0:nfft/2-1)/nfft*fs,ccyp(1:nfft/2));
title('信号重采样');
这是我自己写的一段细化谱的程序,但是标红的地方我觉得有问题,因为细化之后的谱图和细化之前的横坐标是相同的。哪为高手可以帮忙看看 支持原创
页:
[1]