chuantongfenxi.m

内燃机转子仿真

%xinhao de yiban fenxi
clear all;
clc;
format long;
close all;
load 'bneiranji';
load 'b1neiranji1';
load 'b1neiranji2';
load 'b1neiranji3';
load 'b1neiranji4';
load 'b1neiranji5';
%ni=20000;
for i=2:9
%bx=xyz(i,:)-mean(xyz(i,:));
bx=xyz(i,:);
%bx=vv(i,:);
ni=length(bx);
fbx=max(bx);%feng zhi
jdpjzbx=sum(abs(bx))/ni;%jue dui ping jun zhi
fgfzbx=(sum(sqrt(abs(bx)))/ni)^2;%fang gen fu zhi
jfgbx=sqrt(sum(bx.*bx)/ni);;%jun fang gen zhi
fzzbbx=fbx/jdpjzbx;%feng zhi zhi biao
bxzbbx=jfgbx/jdpjzbx;%bo xing zhi bao
mczbbx=fbx/jdpjzbx;%mai chong zhi biao
ydzbbx=fbx/fgfzbx;%yu du zhi biao

%cx1=xyz1(i,:)-mean(xyz1(i,:));
cx1=xyz1(i,:);
%cx1=vv1(i,:);
fcx1=max(cx1);
jdpjzcx1=sum(abs(cx1))/ni;
fgfzcx1=(sum(sqrt(abs(cx1)))/ni)^2;
jfgcx1=sqrt(sum(cx1.*cx1)/ni);
fzzbcx1=fcx1/jdpjzcx1;
bxzbcx1=jfgcx1/jdpjzcx1;
mczbcx1=fcx1/jdpjzcx1;
ydzbcx1=fcx1/fgfzcx1;

%cx2=xyz2(i,:)-mean(xyz2(i,:));
cx2=xyz2(i,:);
%cx2=vv2(i,:);
fcx2=max(cx2);
jdpjzcx2=sum(abs(cx2))/ni;
fgfzcx2=(sum(sqrt(abs(cx2)))/ni)^2;
jfgcx2=sqrt(sum(cx2.*cx2)/ni);
fzzbcx2=fcx2/jdpjzcx2;
bxzbcx2=jfgcx2/jdpjzcx2;
mczbcx2=fcx2/jdpjzcx2;
ydzbcx2=fcx2/fgfzcx2;

%cx3=xyz3(i,:)-mean(xyz3(i,:));
cx3=xyz3(i,:);
%cx3=vv3(i,:);
fcx3=max(cx3);
jdpjzcx3=sum(abs(cx3))/ni;
fgfzcx3=(sum(sqrt(abs(cx3)))/ni)^2;
jfgcx3=sqrt(sum(cx3.*cx3)/ni);
fzzbcx3=fcx3/jdpjzcx3;
bxzbcx3=jfgcx3/jdpjzcx3;
mczbcx3=fcx3/jdpjzcx3;
ydzbcx3=fcx3/fgfzcx3;

%cx4=xyz4(i,:)-mean(xyz4(i,:));
cx4=xyz4(i,:);
%cx4=vv4(i,:);
fcx4=max(cx4);
jdpjzcx4=sum(abs(cx4))/ni;
fgfzcx4=(sum(sqrt(abs(cx4)))/ni)^2;
jfgcx4=sqrt(sum(cx4.*cx4)/ni);
fzzbcx4=fcx4/jdpjzcx4;
bxzbcx4=jfgcx4/jdpjzcx4;
mczbcx4=fcx4/jdpjzcx4;
ydzbcx4=fcx4/fgfzcx4;


cx5=xyz5(i,:);
%cx4=vv4(i,:);
fcx5=max(cx5);
jdpjzcx5=sum(abs(cx5))/ni;
fgfzcx5=(sum(sqrt(abs(cx5)))/ni)^2;
jfgcx5=sqrt(sum(cx5.*cx5)/ni);
fzzbcx5=fcx5/jdpjzcx5;
bxzbcx5=jfgcx5/jdpjzcx5;
mczbcx5=fcx5/jdpjzcx5;
ydzbcx5=fcx5/fgfzcx5;



fengz(:,i)=[fbx;fcx1;fcx2;fcx3;fcx4;fcx5];
fengz(:,i)=fengz(:,i)/max(fengz(:,i));
fzzbi(:,i)=[fzzbbx;fzzbcx1;fzzbcx2;fzzbcx3;fzzbcx4;fzzbcx5];
fzzbi(:,i)=fzzbi(:,i)/max(fzzbi(:,i));
bxzbi(:,i)=[bxzbbx;bxzbcx1;bxzbcx2;bxzbcx3;bxzbcx4;bxzbcx5];
bxzbi(:,i)=bxzbi(:,i)/max(bxzbi(:,i));
mczbi(:,i)=[mczbbx;mczbcx1;mczbcx2;mczbcx3;mczbcx4;mczbcx5];
mczbi(:,i)=mczbi(:,i)/max(mczbi(:,i));
ydzbi(:,i)=[ydzbbx;ydzbcx1;ydzbcx2;ydzbcx3;ydzbcx4;ydzbcx5];
ydzbi(:,i)=ydzbi(:,i)/max(ydzbi(:,i));




pbx=skewness(bx);%pian du
pdzbbx=pbx/jfgbx;%pian du zhi biao

pcx1=skewness(cx1);
pdzbcx1=pcx1/jfgcx1;

pcx2=skewness(cx2);
pdzbcx2=pcx2/jfgcx2;

pcx3=skewness(cx3);
pdzbcx3=pcx3/jfgcx3;

pcx4=skewness(cx4);
pdzbcx4=pcx4/jfgcx4;

pcx5=skewness(cx5);
pdzbcx5=pcx5/jfgcx5;

pd=[pbx;pcx1;pcx2;pcx3;pcx4;pcx5];
xdpd=pd/max(pd);
pdc=[pbx-pbx;pcx1-pbx;pcx2-pbx;pcx3-pbx;pcx4-pbx;pcx5-pbx];
pdzbxyz(:,i)=[pdzbbx;pdzbcx1;pdzbcx2;pdzbcx3;pdzbcx4;pdzbcx5];
pdzbxyz(:,i)=pdzbxyz(:,i)/max(pdzbxyz(:,i));

pdi(:,i)=xdpd;
pdci(:,i)=pdc;

qbx=kurtosis(bx);%qiao du
qdzbbx=qbx/jfgbx;%qiao du zhi biao

qcx1=kurtosis(cx1);
qdzbcx1=qcx1/jfgcx1;

qcx2=kurtosis(cx2);
qdzbcx2=qcx2/jfgcx2;

qcx3=kurtosis(cx3);
qdzbcx3=qcx3/jfgcx3;

qcx4=kurtosis(cx4);
qdzbcx4=qcx4/jfgcx4;

qcx5=kurtosis(cx5);
qdzbcx5=qcx5/jfgcx5;

qd=[qbx;qcx1;qcx2;qcx3;qcx4;qcx5];
xdqd=qd/max(qd);
qdc=[qbx-qbx;qcx1-qbx;qcx2-qbx;qcx3-qbx;qcx4-qbx;qcx5-qbx];
qdzbxyz(:,i)=[qdzbbx;qdzbcx1;qdzbcx2;qdzbcx3;qdzbcx4;qdzbcx5];
qdzbxyz(:,i)=qdzbxyz(:,i)/max(qdzbxyz(:,i));
qdi(:,i)=xdqd;
pdci(:,i)=qdc;
end
ii=9;
ifengz=fengz';
ifzzbi=fzzbi';
ibxzbi=bxzbi';
imczbi=mczbi';
iydzbi=ydzbi';
ipdi=pdi';
ipdzbxyz=pdzbxyz';
iqdi=qdi';
iqdzbxyz=qdzbxyz';
bardata=[ifengz(ii,:);ifzzbi(ii,:);imczbi(ii,:);ipdi(ii,:);ipdzbxyz(ii,:);iqdi(ii,:);iqdzbxyz(ii,:)];
bar(bardata);
ylim([0.5 1.02]);
%legl=legend('e=0','e=0.05','e=0.1','e=0.2','e=0.35');
%set(legl,'Interpreter','none');
grid on;
if ii==2
    title('轴段2时域参数随裂纹深度的变化');
    set(gca,'XTickLabel',{'峰值';'峰值指标';'脉冲指标';'偏度';'偏度指标';'峭度';'峭度指标'});
elseif ii==3
    title('轴段3时域参数随裂纹深度的变化');
    set(gca,'XTickLabel',{'峰值';'峰值指标';'脉冲指标';'偏度';'偏度指标';'峭度';'峭度指标'});
elseif ii==4
    title('轴段4时域参数随裂纹深度的变化');
    set(gca,'XTickLabel',{'峰值';'峰值指标';'脉冲指标';'偏度';'偏度指标';'峭度';'峭度指标'});
elseif ii==5
    title('轴段5时域参数随裂纹深度的变化');
    set(gca,'XTickLabel',{'峰值';'峰值指标';'脉冲指标';'偏度';'偏度指标';'峭度';'峭度指标'});
elseif ii==6
    title('轴段6时域参数随裂纹深度的变化');
    set(gca,'XTickLabel',{'峰值';'峰值指标';'脉冲指标';'偏度';'偏度指标';'峭度';'峭度指标'});
elseif ii==7
    title('轴段7时域参数随裂纹深度的变化');
    set(gca,'XTickLabel',{'峰值';'峰值指标';'脉冲指标';'偏度';'偏度指标';'峭度';'峭度指标'});
elseif ii==8
    title('轴段8时域参数随裂纹深度的变化');
    set(gca,'XTickLabel',{'峰值';'峰值指标';'脉冲指标';'偏度';'偏度指标';'峭度';'峭度指标'});
elseif ii==9
    title('轴段9时域参数随裂纹深度的变化');
    set(gca,'XTickLabel',{'峰值';'峰值指标';'脉冲指标';'偏度';'偏度指标';'峭度';'峭度指标'});
end