clc;clear;close all;
sm1=[0 180.73 456.15 1147.75];
sa1=[241.3 147.87 152.05 0];
l1=241.3;
l2=233.5141747;
l3=480.8243182;
l4=1147.75;
arfa1=90/180*pi;
arfa2=39.3/180*pi;
arfa3=18.4/180*pi;
arfa4=0/180*pi;
ll1=21.2856;
ll2=19.4835;
ll3=41.3846;
ll4=61.2515;
csm21=l1-ll1;
csm22=l2-ll2;
csm23=l3-ll3;
csm24=l4-ll4;
sm21=csm21*cos(arfa1);
sa21=csm21*sin(arfa1);
sm22=csm22*cos(arfa2);
sa22=csm22*sin(arfa2);
sm23=csm23*cos(arfa3);
sa23=csm23*sin(arfa3);
sm24=csm24*cos(arfa4);
sa24=csm24*sin(arfa4);
csm31=l1+ll1;
csm32=l2+ll2;
csm33=l3+ll3;
csm34=l4+ll4;
sm31=csm31*cos(arfa1);
sa31=csm31*sin(arfa1);
sm32=csm32*cos(arfa2);
sa32=csm32*sin(arfa2);
sm33=csm33*cos(arfa3);
sa33=csm33*sin(arfa3);
sm34=csm34*cos(arfa4);
sa34=csm34*sin(arfa4);
sm2=[sm21 sm22 sm23 sm24];
sa2=[sa21 sa22 sa23 sa24];
sm3=[sm31 sm32 sm33 sm34];
sa3=[sa31 sa32 sa33 sa34];
[x1Data, y1Data] = prepareCurveData( sm1, sa1 );
[x2Data, y2Data] = prepareCurveData( sm2, sa2 );
[x3Data, y3Data] = prepareCurveData( sm3, sa3 );
ft1 = 'pchipinterp';
pict1=fit( x1Data, y1Data, ft1);
pict2=fit( x2Data, y2Data, ft1);
pict3=fit( x3Data, y3Data, ft1);
figure( 'Name', '王海航第二题作业' );
h1 = plot( pict1,'k',sm1, sa1,'ko')
set(h1,'LineWidth',1.5,'MarkerSize',4,{'MarkerFaceColor'},get(h1,'color'));
hold on
h2 = plot( pict2,'b--',sm2, sa2,'b^')
set(h2,'LineWidth',1.5,'MarkerSize',4,{'MarkerFaceColor'},get(h1,'color'));
hold on
h3 = plot( pict3,'b--',sm3, sa3,'b^')
set(h3,'LineWidth',1.5,'MarkerSize',4,{'MarkerFaceColor'},get(h1,'color'));
hh1=legend('\sigma_m_1-\sigma_a_1','R0.5','\sigma_m_2-\sigma_a_2','R0.99','\sigma_m_3-\sigma_a_3','R0.99');%,'Unlubricated','Location','NorthWest');
set(hh1,'Fontname','Times New Roman','Fontsize',10,'Box','off');
%以下为美化图形程序
fig_FontSize = 10;
gcf1_position=[600 500 480 360];
xlim([0,1200]);ylim([0, 300]);
xlabel('\sigma_m(MPa)','Fontname','Times New Roman');
ylabel('\sigma_a(MPa)','Fontname','Times New Roman');
set(gcf,'Position',gcf1_position,'color','white');
set(gca,'Position',[0.13 0.16 0.80 0.77],'fontsize',fig_FontSize,'Fontname', 'Times New Roman','FontWeight','Normal','Xtick',0:200:1200);
set(get(gca,'XLabel'),'FontSize',fig_FontSize,'Vertical','top');
set(get(gca,'YLabel'),'FontSize',fig_FontSize,'Vertical','bottom');