- 1、本文档共6页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
matlab常语句整理
设置横坐标间隔
使用绘图函数后(例如使用了plot(x,y)),可以使用下面语句进行横纵坐标间距的设置。
set(gca,xtick,-4:0.5:2.5)
以下是设置坐标的详细代码:
axis([0 24 0 900]);set(gca,xtick,[0 1 2 ...24]);set(gca,ytick,[0 100 200 300 400 500 600 700 800 900]);
如何使三维图的颜色更漂亮
比如使用了三维绘图函数plot3,(plot3绘制三维图的例子,详见三维图绘图教程)
shading interp????????????%灯光处理a=load(final.txt);
5. 佛山南海区项目用到的程序:
yxieliang=interp1(shuiwei,xieliuliang, 17,spline)
%由水位查泄流量,其中,shuiwei是水位的数组,xieliuliang是泄流量数组
yshuiwei=interp1(kurong,shuiweiku, 794.0969,spline)
%由库容查水位
Q=1*0.32*5*(2*9.8)^0.5*H0.^1.5
6. 生成类似abc1、abc2这样的顺序变量名
for i=1:5
a=[i,i+3];
str=([Var,num2str(i) = a, ;]);
eval(str)
end
这是matlab运行的结果
7. 求泄量过程线使用的程序(2011-4下旬使用的)
下面程序运行的条件:
在matlab中创建两个变量:shuiweiku,kurong其中shuiweiku表示水位-库容的水位数组,kurong表示水位-库容的库容数组,创建后,复制黏贴下面代码,就可以运行。
e=1/3; %e是表示水闸的开启程度,e=1/3表示开三分之一
ashuiwei(1)=16.48; %ashuiwei(1)=16.48是指从16.48开始泄水
%ashuiwei(i)表示时段i的实时水位
akurong(1)= interp1(shuiweiku,kurong, ashuiwei(1),spline);
e2=(17-15.5)*e; %17-15.5表示水闸的高度
for i=1:327
if ashuiwei(i)16.47 %16.47是高溢洪道的高度
Qh(i)=1*0.32*10*(2*9.8)^0.5*(ashuiwei(i)-16.47)^1.5; %Qh(i)表示高溢洪道的流量
else
Qh(i)=0; %如果水位低于16.47,则高溢洪道流量为0
end
Qy(i)=e*0.32*5*(2*9.8)^0.5*(ashuiwei(i)-15.5)^1.5; %宽顶堰泄流公式
H=ashuiwei(i)-15.5; %H表示图1的H
H1=H-e2; %H1表示图1 的H1
if H1=0 %如果H1《=0,则不用乘以那个系数
Ql(i)=Qy(i);
else
Ql(i)=Qy(i)*((H^1.5-H1^1.5)/H^1.5); %Ql(i)表示低溢洪道的流量,乘以那个系数
end
Q(i)=Qh(i)+Ql(i); %总流量=高、低溢洪道流量之和
axieliang(i)=Q(i); %axieliang表示实时泄流量
akurong(i+1)=akurong(i)-axieliang(i)*60*60/10000; %akurong表示实时库容
ashuiwei(i+1)=interp1(kurong,shuiweiku, akurong(i+1),spline);
if ashuiwei(i+1)=15.5
break
end
end
%运行下面句子,再查看total1和total2的值,
total1(1,:)=ashuiwei(2:255);
%total(1,:)表示total矩阵的第一行的值是ashuiwei第2到第255的值
total1(2,:)=axieliang(1:254);
total1(3,:)=akurong(2:255);
total2(1,:)=ashuiwei(256:328);
total2(2,:)=axieliang(255:327);
total2(3,:)=akurong(256:328);
5.1.2检验17.31是否合理所使用的程序
xzshuiwei2(1)=16.47;
xzxieliang2(1)= interp1(shuiwei,xieliuliang,xzshuiwei2(1), linear);
xz
文档评论(0)