abc1=abc(abc(:,1)>0& abc(:,1) <18,:);
li1=li(li(:,1)>0& li(:,1) <18,:);
fd=[abc1(:,2),abc1(:,3),abc1(:,4),li1(:,2)-0.927,li1(:,3)+0.24,li1(:,4)+0.31,li1(:,5)+0.24,li1(:,6)-0.344,li1(:,7)+1.375];
% fj=ones(size(fd,1),9);
% fj(:,1)=fd(:,1);fj(:,2)=fd(:,2);fj(:,3)=fd(:,3);
% fj(:,4)=0.927;fj(:,5)=-0.24;fj(:,6)=-0.31;fj(:,7)=-0.24;fj(:,8)=0.344;fj(:,9)=-1.375;
csvwrite('f31819d.dat',fd);
% csvwrite('f31819j.dat',fj);
% [m,n]=size(fd);
% a=zeros(m,1);
% b=zeros(m,1);
% for i=1:1:m
% a(i,1)=(cos(fd(i,2))*sin(fd(i,1))*cos(fd(i,3))-sin(fd(i,2))*cos(fd(i,3)))/(cos(fd(i,1))*cos(fd(i,3)));
% a(i,1)=atan(a(i,1))*180/pi;
% b(i,1)=sin(fd(i,1))*sin(fd(i,2))*cos(fd(i,3))+sin(fd(i,3))*cos(fd(i,2));
% b(i,1)=asin(b(i,1))*180/pi;
% end
% plot(a);hold on;plot(b);
filname=input('请输入文件名(实验车次):','s'); %数据文件名
[m,n]=size(fd);
Cy=zeros(m,6);
Ct=zeros(m,6);
Cy(:,1)=fd(:,4);
Cy(:,2)=fd(:,5);
Cy(:,3)=fd(:,6);
Cy(:,4)=fd(:,7);
Cy(:,5)=fd(:,8);
Cy(:,6)=fd(:,9);
for i=1:1:m
Ct(i,1)=6.206057*Cy(i,1);
Ct(i,2)=29.4335*Cy(i,2);
Ct(i,3)=11.09471*Cy(i,3);
Ct(i,4)=1.222321*Cy(i,4);
Ct(i,5)=0.974101*Cy(i,5);
Ct(i,6)=2.362205*Cy(i,6);
end
modfil=fopen('scharect.inf','r');
str1=fscanf(modfil,'%s',1);
s=fscanf(modfil,'%f',1);%模型参考面积(m2)
str2=fscanf(modfil,'%s',1);
l=fscanf(modfil,'%f',1);%模型展长(m)
str3=fscanf(modfil,'%s',1);
ba=fscanf(modfil,'%f',1);%模型平均气动弦长(m)
str4=fscanf(modfil,'%s',1);
dx=fscanf(modfil,'%f',1);%模型重心距天平中心位置(m)
str5=fscanf(modfil,'%s',1);
dy=fscanf(modfil,'%f',1);%模型重心距天平中心位置(m)
str6=fscanf(modfil,'%s',1);
dz=fscanf(modfil,'%f',1);%模型重心距天平中心位置(m)
str7=fscanf(modfil,'%s',1);
v=fscanf(modfil,'%f',1);%风速(m/s)
qfilname='Q.dat';qfilname=strcat(filname,qfilname);
tfilname='T.dat';tfilname=strcat(filname,tfilname);
Ctw=zeros(m,6);
for i=1:1:m
Ctw(i,1)=Ct(i,1);
Ctw(i,2)=Ct(i,2);
Ctw(i,3)=Ct(i,3);
Ctw(i,4)=Ct(i,4)+Ct(i,3)*dy-Ct(i,2)*dz;
Ctw(i,5)=Ct(i,5)-Ct(i,1)*dz-Ct(i,3)*dx;
Ctw(i,6)=Ct(i,6)+Ct(i,1)*dy+Ct(i,2)*dx;
end
Cq=zeros(m,6);
for i=1:1:m
Cq(i,1)=cosd(fd(i,1))*cosd(fd(i,2))*Ctw(i,1)+sind(fd(i,1))*cosd(fd(i,2))*Ctw(i,2)-sind(fd(i,2))*Ctw(i,3);
Cq(i,2)=-sind(fd(i,1))*Ctw(i,1)+cosd(fd(i,1))*Ctw(i,2);
Cq(i,3)=cosd(fd(i,1))*sind(fd(i,2))*Ctw(i,1)+sind(fd(i,1))*sind(fd(i,2))*Ctw(i,2)+cosd(fd(i,2))*Ctw(i,3);
Cq(i,4)=cosd(fd(i,1))*cosd(fd(i,2))*Ctw(i,4)-sind(fd(i,1))*cosd(fd(i,1))*Ctw(i,5)+sind(fd(i,2))*Ctw(i,6);
Cq(i,5)=sind(fd(i,1))*Ctw(i,4)+cosd(fd(i,1))*Ctw(i,5);
Cq(i,6)=-cosd(fd(i,1))*sind(fd(i,2))*Ctw(i,4)+sind(fd(i,1))*sind(fd(i,2))*Ctw(i,5)+cosd(fd(i,2))*Ctw(i,6);
end
%气动力转换为气动力系数
for i=1:1:m
Ctw(i,1)=Ctw(i,1)*9.8/(0.5*1.225*v^2*s);
Ctw(i,2)=Ctw(i,2)*9.8/(0.5*1.225*v^2*s);
Ctw(i,3)=Ctw(i,3)*9.8/(0.5*1.225*v^2*s);
Ctw(i,4)=Ctw(i,4)*9.8/(0.5*1.225*v^2*s*l);
Ctw(i,5)=Ctw(i,5)*9.8/(0.5*1.225*v^2*s*l);
Ctw(i,6)=Ctw(i,6)*9.8/(0.5*1.225*v^2*s*ba);
Cq(i,1)=Cq(i,1)*9.8/(0.5*1.225*v^2*s);
Cq(i,2)=Cq(i,2)*9.8/(0.5*1.225*v^2*s);
Cq(i,3)=Cq(i,3)*9.8/(0.5*1.225*v^2*s);
Cq(i,4)=Cq(i,4)*9.8/(0.5*1.225*v^2*s*l);
Cq(i,5)=Cq(i,5)*9.8/(0.5*1.225*v^2*s*l);
Cq(i,6)=Cq(i,6)*9.8/(0.5*1.225*v^2*s*ba);
end
% T=[fd(:,1),fd(:,2),fd(:,3),Ctw]; Q=[fd(:,1),fd(:,2),fd(:,3),Cq];
% csvwrite('x3097T.dat',T);
% csvwrite('x3097Q.dat',Q);
tfilname=fopen(tfilname,'wt');
fprintf(tfilname,' %d %d\n',[m,n]);
fprintf(tfilname,' α\t\t\t\t\tβ\t\t\t\t\tγ\t\t\t\t\tCx\t\t\t\t\tCy\t\t\t\t\tCz\t\t\t\t\tCmx\t\t\t\t\tCmy\t\t\t\t\tCmz\n');
fprintf(tfilname,' % -8.4f % -8.4f % -8.4f % -8.4f % -8.4f % -8.4f % -8.4f % -8.4f % -8.4f\n',[fd(:,1),fd(:,2),fd(:,3),Ctw]');
fclose(tfilname);
qfilname=fopen(qfilname,'wt');
fprintf(qfilname,' %d %d\n',[m,n]);
fprintf(qfilname,' α\t\t\t\t\tβ\t\t\t\t\tγ\t\t\t\t\tCx\t\t\t\t\tCy\t\t\t\t\tCz\t\t\t\t\tCmx\t\t\t\t\tCmy\t\t\t\t\tCmz\n');
fprintf(qfilname,' % -8.4f % -8.4f % -8.4f % -8.4f % -8.4f % -8.4f % -8.4f % -8.4f % -8.4f\n',[fd(:,1),fd(:,2),fd(:,3),Cq]');
fclose(qfilname);
plot(abc1(:,2),Cq(:,2));
xlabel('α(°)');ylabel('升力系数Cy');
grid on;