matlab划分网格.docx
- 文档编号:1639415
- 上传时间:2022-10-23
- 格式:DOCX
- 页数:2
- 大小:12.91KB
matlab划分网格.docx
《matlab划分网格.docx》由会员分享,可在线阅读,更多相关《matlab划分网格.docx(2页珍藏版)》请在冰豆网上搜索。
clear
clc
%划分矩形网格的时候,底边的一定是偶数份,侧边无所谓
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];%边界坐标数组x
y=[y_low,y_r,y_up,y_l];%边界坐标数组y
plot(x,y,'m');%画出图形
holdon
n1=input('下边界等分段数n1=');
x_lowcp=linspace(x_low
(1),x_low
(2),n1+1);%n等分下边界
y_lowcp=linspace(y_low
(1),y_low
(2),n1+1);
plot(x_lowcp,y_lowcp,'o')%画出等分点
holdon
n2=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=');
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')%画出等分点
holdon
fori=2:
n2+1
plot([x_lowcp(n1/2+1),x_rcp(i)],[y_lowcp(n1/2+1),y_rcp(i)],'m');%连接底边中点和右边界等分点
end
fori=2:
n3+1
plot([x_lowcp(n1/2+1),x_upcp(i)],[y_lowcp(n1/2+1),y_upcp(i)],'m');%连接底边中点和上边界等分点
end
fori=2:
n4+1
plot([x_lowcp(n1/2+1),x_lcp(i)],[y_lowcp(n1/2+1),y_lcp(i)],'m');%连接底边中点和左边界等分点
end
holdon
fori=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)]);%画出横线
end
fori=n1:
-1:
n1/2+2
y_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
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- matlab 划分 网格