基本信息
源码名称：频域积分的matlab实现
源码大小：6.08KB
文件格式：.m
开发语言：MATLAB
更新时间：2021-11-22
友情提示：（无需注册或充值，赞助后即可获取资源下载链接）
嘿，亲！知识可是无价之宝呢，但咱这精心整理的资料也耗费了不少心血呀。小小地破费一下，绝对物超所值哦！如有下载和支付问题，请联系我们QQ(微信同号)：813200300
本次赞助数额为： 2 元×
微信扫码支付：2 元
×
请留下您的邮箱，我们将在2小时内将文件发到您的邮箱
源码介绍
频域积分的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; % 建立分析频率序列
频域积分的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; % 建立分析频率序列