matlab划分网格
《matlab划分网格》由会员分享,可在线阅读,更多相关《matlab划分网格(2页珍藏版)》请在装配图网上搜索。
1、clearclc %划分矩形网格的时候,底边的一定是偶数份,侧边无所谓 x_low=0,4;y_low=0,0;% 下边界 x_r=4,4;y_r=0,4;% 右边界 x_up=4,0;y_up=4,4;% 上边界x_l=0,0;y_l=4,0;% 左边界 x=x_low,x_r,x_up,x_l;%边界坐标数组 xy=y_low,y_r,y_up,y_l;%边界坐标数组 yplot(x,y,m);% 画出图形 hold onn1=input( 下边界等分段数 n1=);x_lowcp=linspace(x_low(1),x_low(2),n1+1);%n 等分下边界y_lowcp=linsp
2、ace(y_low(1),y_low(2),n1+1); plot(x_lowcp,y_lowcp,o)% 画出等分点 hold onn2=input( 右边界等分段数 n2=);x_rcp=linspace(x_r(1),x_r(2),n2+1);%n1 等分右边界 y_rcp=linspace(y_r(1),y_r(2),n2+1);n3=input( 上边界等分段数 n3=);x_upcp=linspace(x_up(1),x_up(2),n3+1);%n1 等分上边界 y_upcp=linspace(y_up(1),y_up(2),n3+1);n4=input( 左边界等分段数 n4=
3、); x_lcp=linspace(x_l(1),x_l(2),n4+1);%n1 等分左边界 y_lcp=linspace(y_l(1),y_l(2),n4+1);x_cp=x_lowcp,x_rcp,x_upcp,x_lcp;% 等分点坐标数组y_cp=y_lowcp,y_rcp,y_upcp,y_lcp;plot(x_cp,y_cp,o)% 画出等分点hold on for i=2:n2+1连接底边中点plot(x_lowcp(n1/2+1),x_rcp(i),y_lowcp(n1/2+1),y_rcp(i),m);% 和右边界等分点end连接底边中for i=2:n3+1 plot(x
4、_lowcp(n1/2+1),x_upcp(i),y_lowcp(n1/2+1),y_upcp(i),m);% 点和上边界等分点end连接底边中点for i=2:n4+1 plot(x_lowcp(n1/2+1),x_lcp(i),y_lowcp(n1/2+1),y_lcp(i),m);% 和左边界等分点endhold on for i=2:n1/2 y_vertical(i-1)=(x_lowcp(n1/2+1)-x_lowcp(i)/(x_lowcp(n1/2+1)-x_lowcp(1)*(y_l(1)-y_l(2);%左bu竖线的纵坐标画出左边竖线画出plot(x_lowcp(i),x_lowcp(i),y_lowcp(i),y_vertical(i-1);%plot(x_lowcp(i),x_lowcp(n1-i+2),y_vertical(i-1),y_vertical(i-1);%横线endfor i=n1:-1:n1/2+2y_vertical(i)=(x_lowcp(n1/2+1)-x_lowcp(i)/(x_lowcp(n1+1)-x_lowcp(n1/2+1)*(画出右边竖线y_r(1)-y_r(2);% 右边竖线的纵坐标plot(x_lowcp(i),x_lowcp(i),y_lowcp(i),y_vertical(i);%end
- 温馨提示:
1: 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
2: 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
3.本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
5. 装配图网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。