《matlab划分网格.doc》由会员分享，可在线阅读，更多相关《matlab划分网格.doc（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