基于MATLAB实现的对数据进行傅立叶分析，从时域转换到频域，观察信号的频谱结构，从分析信号的特点+使用说明文档.rar

%用fft对求信号的幅频和相频曲线 function fftfenxi clear;clc;clf; %N=2048; %fft默认计算的信号是从0开始的 %t=linspace(1,2,N);deta=t(2)-t(1);fs=1/deta; %x=5*sin(2*pi*10*t)+5*sin(2*pi*35*t); % N1=256;N2=512;w1=0.2*2*pi;w2=0.3*2*pi;w3=0.4*2*pi; % x=(t>=-200&t<=-200+N1*deta).*sin(w1*t)+(t>-200+N1*deta&t<=-200+N2*deta).*sin(w2*t)+(t>-200+N2*deta&t<=200).*sin(w3*t); fs=12000; % y=load('g:\毕业设计\故障数据及代码\正常\ZB4409.TXT'); % y=load('g:\毕业设计\故障数据及代码\内圈剥落\NB3009.TXT'); y=load('g:\毕业设计\故障数据及代码\外圈剥落\WB2709.TXT'); % y=load('g:\毕业设计\故障数据及代码\滚动体剥落\GB3509.TXT'); N=length(y); deta=1/fs; m=0:N-1; % f=(1:N-2)/N*(fs/2); % f=1./(N*deta)*m; f=m*fs/N %可以查看课本就是这样定义横坐标频率范围的 %下面计算的Y就是x(t)的傅里叶变换数值 %Y=exp(i*4*pi*f).*fft(y)%将计算出来的频谱乘以exp(i*4*pi*f)得到频移后[-2,2]之间的频谱值 Y=fft(y,N); % z=2*abs(Y); z=sqrt(Y.*conj(Y))/fs;%conj(Y)求Y的共轭 figure(7) % plot(f(1:100),z(1:100)); plot(f,z); title('幅频曲线') % xiangwei=angle(Y); % figure(8) % plot(f,xiangwei) % title('相频曲线') % figure(9) % plot(1:N,y,'r') % title('原始信号')