频域积分的matlab实现

【实例简介】频域积分的matlab实现
【实例截图】
【核心代码】f_01=5.98;f_02=30.64;f_03=61.26;f_04=91.12;f_05=121.15; % 频率
f=[f_02,f_03,f_04,f_05,f_01];
fai_01=43.15;fai_02=60.26;fai_03=30.37;fai_04=80.54;fai_05=100.27;% 相位
fai=[fai_02,fai_03,fai_04,fai_05,fai_01];
fai=fai-(fai>180)*360;
a_01=0.037;a_02=0.260;a_03=0.180;a_04=0.070;a_05=0.029;% 幅值
a=[a_02,a_03,a_04,a_05,a_01];
%加速度 、速度、位移
x=a_01*cos(2*pi*f_01*n1/fs fai_01*pi/180) a_02*cos(2*pi*f_02*n1/fs fai_02*pi/180) a_03*cos(2*pi*f_03*n1/fs fai_03*pi/180) a_04*cos(2*pi*f_04*n1/fs fai_04*pi/180) a_05*cos(2*pi*f_05*n1/fs fai_05*pi/180) 0.15;
v_theory=a_01/(2*pi*f_01)*1000*cos(2*pi*f_01*n1/fs fai_01*pi/180-0.5*pi) a_02/(2*pi*f_02)*1000*cos(2*pi*f_02*n1/fs fai_02*pi/180-0.5*pi) a_03/(2*pi*f_03)*1000*cos(2*pi*f_03*n1/fs fai_03*pi/180-0.5*pi) a_04/(2*pi*f_04)*1000*cos(2*pi*f_04*n1/fs fai_04*pi/180-0.5*pi) a_05/(2*pi*f_05)*1000*cos(2*pi*f_05*n1/fs fai_05*pi/180-0.5*pi);
disp_theory=a_01/((2*pi*f_01)^2)*1000000*cos(2*pi*f_01*n1/fs fai_01*pi/180-pi) a_02/((2*pi*f_02)^2)*1000000*cos(2*pi*f_02*n1/fs fai_02*pi/180-pi) a_03/((2*pi*f_03)^2)*1000000*cos(2*pi*f_03*n1/fs fai_03*pi/180-pi) a_04/((2*pi*f_04)^2)*1000000*cos(2*pi*f_04*n1/fs fai_04*pi/180-pi) a_05/((2*pi*f_05)^2)*1000000*cos(2*pi*f_05*n1/fs fai_05*pi/180-pi);


y2=awgn(x,10,'measured');%高斯白噪声


x_accr=y2;  %加入噪声之后的信号,作为加速度信号
x_accr=x_accr-mean(x_accr);
D=length(x_accr); % 数据长度
time_series=t; %分析时间序列
% 加汉宁窗 
w=hann(N,'periodic');
% fft 变换
fft_complex=fft(x_accr.*w');
fft_complex=fft_complex(1:N/2)/N*4;
freq_series=sf*(0:N/2-1)/N; %  建立分析频率序列