inputdata1.m

内燃机转子仿真

clear all;
clc;
format short;
%a=[1 2 3;4 5 6;7 8 9];
load 'sjjl1';%numak
load 'zhongxinchafen';%x
load 'xierbotedata';%xierbote
load 'sjjl';%v
load 'duibi';%wilson xy
load 'jqj';%jingquejie
%numak=numak';
%x=(x(2:3,:))';
%xierbote=(xierbote(2:3,:))';
%v=(v(2:3,:))';
%xy=(xy(2:3,:))';
a=[qx(1:2,1:1000);v(2,:);x(2,:);xierbote(2,:);xy(2,:);numak(2,:)];
fid=fopen('ls1.txt','wt+');
for i=1:1:1000
    for j=1:7
        if j==7 
      fprintf(fid,'  %.7f\n',a(j,i));
        else
       fprintf(fid,'  %.7f',a(j,i));  
        end
    end
end
fclose(fid);
plot(qx(1,1:200),qx(2,1:200),'k',qx(1,1:200),v(2,1:200),'-.*k',qx(1,1:200),x(2,1:200),'-.k',qx(1,1:200),xierbote(2,1:200),':k',qx(1,1:200),xy(2,1:200),'--k',qx(1,1:200),numak(2,1:200),'-.ko')
h=legend('精确解','PIM','中心差分法','希尔伯特法','威尔逊法','纽马克法',1);
xlabel('时间(s)');
ylabel('位移x1');
set(h,'box','on');