- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
可选性曲线程序.doc
x1=[3.46 7.33 9.47 10.48 11.19 11.79 12.78 20.50];
y1=[10.69 56.84 76.89 82.15 84.70 86.32 88.45 100];
x2=[3.46 8.23 15.50 25.5 34.28 42.94 52.91 79.64];
y2=[10.69 56.84 76.89 82.15 84.70 86.32 88.45 100];
x3=[20.5 22.54 37.85 57.40 66.64 72.04 75.48 79.64];
y3=[100 89.31 43.16 23.02 17.85 15.30 13.68 11.55];
x4=[1.3 1.4 1.5 1.6 1.7 1.8 1.9 ];
y4=[10.69 56.84 76.89 82.15 84.70 86.32 88.45 ];
x5=[ 1.4 1.5 1.6 1.7 1.8 1.9];
y5=[ 66.29 25.31 7.72 4.17 2.69 2.13];
plot(x1,y1,x1,y1,+,x2,y2,x2,y2,*);
xlabel(灰分Ad/%,fontsize,8);%标注并设置字体
ylabel(浮物产产率/%,fontsize,8);%标注并设置字体
text(16,90,\beta,fontsize,14);%在指定的坐标轴下注解曲线名称
text(12,57,\lambda,fontsize,14);%在指定的坐标轴下注解曲线名称
h1=gca;
h2=axes(position,get(h1,position));%保证坐标轴的统一性
plot(x3,y3,x3,y3,o);
set(h2,yaxislocation,right,color,none,xticklabel,[],xlim,get(h1,xlim),layer,top);
text(33,65,\theta,fontsize,14);
ylabel(沉物产率/%,fontsize,8);%标注并设置字体
h3=axes(position,get(h2,position));%保证坐标轴的统一性
plot(x4,y4,x4,y4,.,x5,y5,x5,y5,.);
axis(ij,[1.2 2.2 0 100]);
set(h3,xdir,reverse,xaxislocation,top,color,none);
title(密度/kg.cm-3,fontsize,8);%标注并设置字体
text(1.47,25,\epsilon,fontsize,14);%在指定的坐标轴下注解曲线名称
text(1.55,76,\delta,fontsize,14);%在指定的坐标轴下注解曲线名称
%坐标轴的统一,两套坐标体系有偏差,必须统一
xlimits=get(h1,xlim);
ylimits=get(h1,ylim);
xinc=(xlimits(2)-xlimits(1))/10;
yinc=(ylimits(2)-ylimits(1))/10;
set(h1,xtick,[xlimits(1):xinc:xlimits(2)],ytick,[ylimits(1):yinc:ylimits(2)])
grid on %打开栅格线 如果不想用栅格线,请把ON变成OFF
load x1;
load y1;
load x2;
load y2;
load x3;
load y3;
load x4;
load y4;
load x5;
load y5;
[m1,n1]=size(x1);
[m2,n2]=size(x2);
[m3,n3]=size(x3);
[m4,n4]=size(x4);
[m5,n5]=size(x5);
for i=1:1:n1-1
x11(1,i)=x1(1,i)+(x1(1,i+1)-x1(1,i))/2;
y11(1,i)=interpl(x1,y1,x11(1,i),spline)
end
for j=1:1:n2-1
x21(1,j)=x2(1,j)+(x2(1,j+1)-x2(1,j))/2;
y21(1,j)=interpl(x2,y2,x21(1,j),spline)
end
for k=1:1:n3-1
x31(1,k)=x3(1,k)+(x3(1,k+1)-x3(1,k))/2;
y31(1,k)=interpl(x3,y3,x31(1,k),spline)
end
for l=1:1:n4-1
x41(1,l)=x4(1,l)+(x4(1,l+1)-x4(1,l))/
文档评论(0)