function?[xy_Lxy_Rxy_Uxy_D]?=?linear_boundary（x_maxx_miny_maxy_minnum_node1num_node2n1n2）
%?構造線性矩形邊界
%?主要用到線性插值函數linspace
%?INPUT:
%??????x_max?|?x_min?|?y_max?|?y_min?分別是矩形外邊界的范圍，注意要包含翼型
%??????num_node1?矩形邊界x方向節點數，num_node2為y向節點數
%??????n1?為矩形下邊界節點數，n2為矩形左邊界節點數
%?OUTPUT：
%??????依次為??【左，右，上，下】邊界坐標，每個均為2列（包含xy?方向）的列向量
%%?導入翼型坐標并處理
xy0?=?importdata（‘CH10?（smoothed）.txt‘）;?%?該文件由profili軟件得到，去除里面的文字，只保留點坐標
%?plot（xy0（:1）xy0（:2））;
num_xy0?=?size（xy01）;?%?應該為偶數上下面點數相同

a?=?xy0（?2?:?num_xy0/2??:?）;?%?處理數據格式
b?=?xy0（?num_xy0/2+1?:?end?:）;
b?=?flipud（b）;
xy0?=?[b;a];???????????????????%?得到逆時針排序的翼型點坐標
%?plot（xy0（:1）xy0（:2））;
num_xy0?=?size（xy01）;
clear?a?b;
%%?構造邊界
%?構造下邊界xy_D
%?num_node1?=?89;
n_temp?=?（num_node1?-?num_xy0）/2?;
a_temp?=?zeros（n_temp+12）;
b_temp?=?a_temp;
a_temp（:1）?=?linspace（x_max100n_temp+1）;
a_temp（:2）?=?0;
b_temp（:1）?=?linspace（100x_maxn_temp+1）;
b_temp（:2）?=?0;
xy_D?=?[a_temp;xy0（2:end-1:）;b_temp];?
clear?a_temp?b_temp?xy0;???????????????????%?得到xy_D
%?構造上邊界xy_U
xy_U?=?xy_D;
%?n1?=?39;
%?n2?=?10;
n3?=?num_node1?-?n1?-?n2;
xy_U（1:n11）?=?linspace（x_maxx_minn1）;
xy_U（1:n12）?=?y_min;

xy_U（n1:n1+n22）?=?linspace（y_miny_maxn2+1）;
xy_U（n1:n1+n21）?=?x_min;

xy_U（n1+n2:end1）?=?linspace（x_minx_maxn3+1）;
xy_U（n1+n2:end2）?=?y_max;
%?構造左邊界xy_L
%?num_node2?=?6;
xy_L?=?zeros（num_node22）;
xy_L（:2）?=?linspace（y_min0num_node2）;
xy_L（:1）?=?x_max;
%?構造右邊界
xy_R?=?xy_L;
xy_R（:2）?=?linspace（y_max0num_node2）;
xy_R（:1）?=?x_max;
clear?num_node1?num_node2?;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1533??2010-05-08?15:33??劃分網格\CH10?（smoothed）.txt

?????文件???????1771??2010-05-09?00:07??劃分網格\linear_boundary.m

?????文件???????1301??2010-05-09?10:43??劃分網格\Main_proc.m

?????文件????????127??2010-05-09?09:00??劃分網格\plot_boundary.m

?????文件???????2081??2010-05-09?00:12??劃分網格\solve_4_3_13.m

?????目錄??????????0??2014-04-11?17:10??劃分網格

-----------?---------??----------?-----??----

?????????????????6813????????????????????6