clear;
h0=6.3;epsilon=0.3921;R=125;VR=1308.33;Sf=0.04558;h1=h0*(1-epsilon);
Delta_h=h0-h1;L=sqrt(R*Delta_h-(Delta_h^2)/4);
x=0:0.01:L;a=0.982;b=0.01;
z=Delta_h*x/(2*L)+h1/2;
delta_f=15+9.2*(abs(log(2*z/h0))).^0.11;
p1=10*6*b*L*VR/(h0-h1)*(((a*sqrt(h0*h1))/(4*(h1/2)^2)-(a*sqrt(h0*h1))./(4*z.^2))-(2/h1-1./z))+delta_f;
L1=plot(x,p1);
hold on
x=linspace(0,17.5,10);
z=Delta_h*x/(2*L)+h1/2;
delta_f=15+9.2*(abs(log(2*z/h0))).^0.11;
p1=10*6*b*L*VR/(h0-h1)*(((a*sqrt(h0*h1))/(4*(h1/2)^2)-(a*sqrt(h0*h1))./(4*z.^2))-(2/h1-1./z))+delta_f;
L2=plot(x,p1,'p');
hold on
delta_f=39;alpha=0.1407;
sita=0:0.0001:alpha;
hsita=h1+R*sita.^2;
p1=(delta_f*(pi/4*log(hsita/h1)+pi/4+sqrt(R/h1)*atan(sita*sqrt(R/h1)))).*(sita>=0&sita<0.0550) ...
+(delta_f*(pi/4*log(hsita/h0)+pi/4+sqrt(R/h1)*atan(alpha*sqrt(R/h1))-sqrt(R/h1)*atan(sita*sqrt(R/h1)))).*(sita>=0.0550&sita<=alpha);
L3=plot(R*sin(sita),p1);
hold on
sita=linspace(0,alpha,10);
hsita=h1+R*sita.^2;
p1=(delta_f*(pi/4*log(hsita/h1)+pi/4+sqrt(R/h1)*atan(sita*sqrt(R/h1)))).*(sita>=0&sita<0.0550) ...
+(delta_f*(pi/4*log(hsita/h0)+pi/4+sqrt(R/h1)*atan(alpha*sqrt(R/h1))-sqrt(R/h1)*atan(sita*sqrt(R/h1)))).*(sita>=0.0550&sita<=alpha);
L4=plot(R*sin(sita),p1,'d');
hold on
delta=0.27*2*L/Delta_h;
delta_f=40;x=0:0.01:L;
hx=h1/2+Delta_h./(2*L)*x;
p1=delta_f./delta.*((delta+1)*(2*hx/h1).^delta-1).*(x>=0&x<6.23) ...
+delta_f./delta.*((delta-1)*(h0./(2*hx)).^delta+1).*(x>=6.23&x<=L);
L5=plot(x,p1);
hold on
x=linspace(0,L,10);
hx=h1/2+Delta_h./(2*L)*x;
p1=delta_f./delta.*((delta+1)*(2*hx/h1).^delta-1).*(x>=0&x<6.23) ...
+delta_f./delta.*((delta-1)*(h0./(2*hx)).^delta+1).*(x>=6.23&x<=L);
L6=plot(x,p1,'p');
hold on
%L7=plot(p(:,1),p(:,2)); 这是另一个文件导入的数据先可以不要
set(gca,'position',[0.15,0.15,0.8,0.8]);
text(5,50,'n 0.982 \eta_p 0.01 K 40 \mu 0.27');
legend([L2,L4,L6,L7],'HM','R.B.Sims','А.И. Целиков','ex-p')
这应该是图形输出设置问题,如果方便给出具体绘图代码。
就图形绘制而言,MATLAB输出美观的图形不会有问题,只是设置比较麻烦。
建议期刊论文的图形绘制用Origin之类的软件来实现,这类软件特点是简单操作即可获得美观的图形。
不建议用Excel来绘制科技论文图,我还记得有次听领域内某顶级期刊主编做的关于科技论文写作图表规范的讲座,在讲座上该主编说 excel makes bad figures。