huang.m

内燃机转子仿真

clear all;
clc;
fomart long;
fname=input('the name of the data=','s');
fid=fopen(fname,'r');
%%%%%%%%%%%%
de=fscanf(fid,'%f',1);
qa=fscanf(fid,'%f',1);
ra=fscanf(fid,'%f',1);
e=fscanf(fid,'%f',1);
%%%%%%%%%%%%%%%%%%%%

dd2=fscanf(fid,'%f',1);
v=fscanf(fid,'%f',1);
rf=fscanf(fid,'%f',1);
%%%%%%%%%%%%%%%%%%%%%%%%%%
r2=fscanf(fid,'%f',1);
r1=fscanf(fid,'%f',1);
pf=fscanf(fid,'%f',1);
qi=fscanf(fid,'%f',1);
dd=fscanf(fid,'%f',1);
k=fscanf(fid,'%f',1);
fclose(fid);

dp1=de*e*(1-qa^2)/4/ra/qa^2;%3
p2=e*dd2/2/((1+qa^2)/(1-qa^2)+v);%4
p11=(e*dd2/2*(r2^2-r1^2)/yf2+((1-v)*r2^2+(1+v)*r1^2)*p2)/2/r1^2;%5
p1=dp1+p11;%6
ddt=2*pf*rf*((1+qa^2)/(1-qa^2)+(1+qi^2)/(1-qi^2))/e%7
dd1t=2*r2*(2*r1^2*p1+((v-1)*r2^2-(v+1)*r1^2)*p2)+pf*rf*((1+qi^2)/(1-qi^2)-v)/e;%8
dd1=ddt+dd1t-dd;%9
s=(dd1+dd)/k;%10
fid=fopen(output,'wt');
fprintf(fid,'%f','','');
status=fclose(fid);