仅供参考
KPa10=[0 71 4.9 56;4.9 56 10 45;10 45 15.1 37;15.1 37 20.1 29;20.1 29 25 21;25 21 30.1 14];
c1=KPa10(:,1);
c2=KPa10(:,2);
c3=KPa10(:,3);
c4=KPa10(:,4);
a=[]; b=[]; syms x y ;
for i=1:6
[solx,soly]=solve([((0.1-y*pi*c1(i)/2000).*(0.08-y*pi*c1(i)/2000)).*x*10000==c2(i), ((0.1-y*pi*c3(i)/2000).*(0.08-y*pi*c3(i)/2000)).*x*10000==c4(i)], [x y]);
solx=vpa(solx,4);
soly=vpa(soly,4);
a(i,:)=double(solx);
b(i,:)=double(soly);
fid1=fopen('3.txt','a');
fprintf(fid1,'a(%d) = %d \r\nb(%d) = %d \r\n',i,a(i),i,b(i));
fclose(fid1);
end