MATLAB程序集.docVIP

1．数据预处理主程序（m文件名：wz_main_090531.m）: %wz_main %建立日期: 2009.05.01 %修改日期: 2009.05.31 %清除窗口及变量空间 zz %===========参数设置区======================================== sta_K=3; %数据开始位置 end_K=7; %数据结束位置 %============================================================= w_BT=请选择你要装载的数据文件; [DataLoad_1 WJ_Name] = wzfun_LoadFile(w_BT); disp([装载的文件名为： WJ_Name]) JiPin=input(请输入基频频率(Hz)：) disp(幅值归一化后结果) K=1024; [Data_hang Data_lie]=size(DataLoad_1); DataLoad(:,1)=DataLoad_1(:,1); DataLoad(:,2:Data_lie)=DataLoad_1(:,2:Data_lie)/1000; for w1_i=2:Data_lie Data_Yuan=DataLoad(sta_K*K+1:end_K*K,[1 w1_i]); Data_Lvbo=wzfun_LvBo(Data_Yuan(:,2),1024,10,400); w_td_str=[第 num2str(w1_i-1) 通道信号]; wzfun_Graph(Data_Yuan,5*(w1_i-2)+1,1,原始信号时域波形图,时间(s),振动加速度(m/s^2),b,w_td_str) Data_Lvbo_t=[Data_Yuan(:,1) Data_Lvbo]; wzfun_Graph(Data_Lvbo_t,5*(w1_i-2)+2,1,滤波后信号时域波形图,时间(s),振动加速度(m/s^2),r,w_td_str) [f_Yuan,Pin_A_Yuan]= wzfun_PinPu(Data_Yuan); [f_Lvbo,Pin_A_Lvbo]= wzfun_PinPu(Data_Lvbo_t); Data_Yuan_Pin_A=[f_Yuan,Pin_A_Yuan]; Data_Lvbo_Pin_A=[f_Lvbo,Pin_A_Lvbo]; wzfun_Graph(Data_Yuan_Pin_A,5*(w1_i-2)+3,1,原始信号频谱（幅值谱）图,频率(Hz),振动加速度(m/s^2),b,w_td_str) wzfun_Graph(Data_Lvbo_Pin_A,5*(w1_i-2)+4,1,滤波后信号频谱（幅值谱）图,频率(Hz),振动加速度(m/s^2),r,w_td_str) %求基频及前四次谐波所对应的幅值 %=============================================================== %[x1_zb y1_zz]=wzfun_ZuiZhi(f_Lvbo,Pin_A_Lvbo); f_max=max(f_Lvbo); munm=size(f_Lvbo); for i=1:4 fw_min=(JiPin)*i*0.8; fw_max=(JiPin)*i*1.2; ni=round(fw_min*munm/f_max+1); if ni=1; ni=1; end na=round(fw_max*munm/f_max+1); if na=munm na=munm; end [x_zb(i) y_zz(i)]=wzfun_ZuiZhi(f_Lvbo(ni:na),Pin_A_Lvbo(ni:na)); end %[pn,minp,maxp,tn,mint,maxt] = premnmx(p,t) %对p,t进行归一化 y_zz_guiyi=norm(y_zz); shujubiao=[x_zb y_zz/y_zz_guiyi]; h=figure(5*(w1_i-2)+5); set(h,Color,White); bar(shujubiao(:,2),0.5); ti