gonglp1.m

内燃机转子仿真

%function mpyy1=gonglp1()
clear all;
clc;
close all;
format long;
load 'b2neiranji';
load 'b2neiranji1';
load 'b2neiranji2';
load 'b2neiranji3';
load 'b2neiranji4';
load 'b2neiranji5';

ip=7;
x=(xyz(ip,:)-mean(xyz(ip,:)))*180/pi;
x1=(xyz1(ip,:)-mean(xyz1(ip,:)))*180/pi;
x2=(xyz2(ip,:)-mean(xyz2(ip,:)))*180/pi;
x3=(xyz3(ip,:)-mean(xyz3(ip,:)))*180/pi;
x4=(xyz4(ip,:)-mean(xyz4(ip,:)))*180/pi;
x5=(xyz5(ip,:)-mean(xyz5(ip,:)))*180/pi;
%x=aaii(ip,2000:80000);
%x1=aaii1(ip,2000:80000);
%x2=aaii2(ip,2000:80000);
%x3=aaii3(ip,2000:80000);
%x4=aaii4(ip,2000:80000);
%x5=aaii5(ip,2000:80000);



step=h;

sf=1/step;
%t=0:step:80;
%y1=sin(2*pi*10*t);
%y2=sin(2*pi*5*t);
ng=length(x);
%nfft=2^nextpow2(ng);
nfft=2^15;
%w1=xcorr(y1,z1,200);
%w2=xcorr(y2,z2,200);
%x=-nn*step:step:nn*step;
yy=fft(x,nfft);
pyy=yy.*conj(yy)/nfft;
yy1=fft(x1,nfft);
pyy1=yy1.*conj(yy1)/nfft;
yy2=fft(x2,nfft);
pyy2=yy2.*conj(yy2)/nfft;
yy3=fft(x3,nfft);
pyy3=yy3.*conj(yy3)/nfft;
yy4=fft(x4,nfft);
pyy4=yy4.*conj(yy4)/nfft;
yy5=fft(x5,nfft);
pyy5=yy5.*conj(yy5)/nfft;

ff=(0:sf/nfft:sf/2-sf/nfft);%%%%%%%%%%%%%%%%%%%%%%%%
nn=1:length(ff);
upx=0;%0-120 120-220 220-400 400-600 600-800
lox=150;
upx1=150;%0-120 120-220 220-400 400-600 600-800
lox1=300;
upx2=300;%0-120 120-220 220-400 400-600 600-800
lox2=450;
upx3=450;%0-120 120-220 220-400 400-600 600-800
lox3=600;
upx4=600;%0-120 120-220 220-400 400-600 600-800
lox4=800;

ulimy=3;
llimy=0;
figure(1);
plot(ff(nn),pyy(nn),'k',ff(nn),pyy1(nn),'-.k',ff(nn),pyy2(nn),'.-k',ff(nn)...
    ,pyy3(nn),'--k',ff(nn),pyy4(nn),':k',ff(nn),pyy5(nn),'-kx');
hl=legend('e=0','e=0.1','e=0.2','e=0.3','e=0.4','e=0.45',1);
xlabel('频率(Hz)');
ylabel('幅值');
xlim([upx,lox]);
%ylim([llimy ulimy]);
set(hl,'Interpreter','none');
%xlim('auto');
ylim('auto');

title('额定转速第九轴段频率在[0 150]Hz扭角位移自谱图');

figure(2);
plot(ff(nn),pyy(nn),'k',ff(nn),pyy1(nn),'-.k',ff(nn),pyy2(nn),'.-k',ff(nn)...
    ,pyy3(nn),'--k',ff(nn),pyy4(nn),':k',ff(nn),pyy5(nn),'-kx');
hl=legend('e=0','e=0.1','e=0.2','e=0.3','e=0.4','e=0.45',1);
xlabel('频率(Hz)');
ylabel('幅值');
xlim([upx1,lox1]);
%ylim([llimy ulimy]);
set(hl,'Interpreter','none');
%xlim('auto');
ylim('auto');


title('额定转速第九轴段频率在[150 300]Hz扭角位移自谱图');

figure(3);
plot(ff(nn),pyy(nn),'k',ff(nn),pyy1(nn),'-.k',ff(nn),pyy2(nn),'.-k',ff(nn)...
    ,pyy3(nn),'--k',ff(nn),pyy4(nn),':k',ff(nn),pyy5(nn),'-kx');
hl=legend('e=0','e=0.1','e=0.2','e=0.3','e=0.4','e=0.45',1);
xlabel('频率(Hz)');
ylabel('幅值');
xlim([upx2,lox2]);
%ylim([llimy ulimy]);
set(hl,'Interpreter','none');
%xlim('auto');
ylim('auto');


title('额定转速第九轴段频率在[300 450]Hz扭角位移自谱图');

figure(4);
plot(ff(nn),pyy(nn),'k',ff(nn),pyy1(nn),'-.k',ff(nn),pyy2(nn),'.-k',ff(nn)...
    ,pyy3(nn),'--k',ff(nn),pyy4(nn),':k',ff(nn),pyy5(nn),'-kx');
hl=legend('e=0','e=0.1','e=0.2','e=0.3','e=0.4','e=0.45',1);
xlabel('频率(Hz)');
ylabel('幅值');
xlim([upx3,lox3]);
%ylim([llimy ulimy]);
set(hl,'Interpreter','none');
%xlim('auto');
ylim('auto');


title('额定转速第九轴段频率在[450 600]Hz扭角位移自谱图');

figure(5);
plot(ff(nn),pyy(nn),'k',ff(nn),pyy1(nn),'-.k',ff(nn),pyy2(nn),'.-k',ff(nn)...
    ,pyy3(nn),'--k',ff(nn),pyy4(nn),':k',ff(nn),pyy5(nn),'-kx');
hl=legend('e=0','e=0.1','e=0.2','e=0.3','e=0.4','e=0.45',1);
xlabel('频率(Hz)');
ylabel('幅值');
xlim([upx4,lox4]);
%ylim([llimy ulimy]);
set(hl,'Interpreter','none');
%xlim('auto');
ylim('auto');


title('额定转速第九轴段频率在[600 800]Hz扭角位移自谱图');