tianranqi4.m

天然气管道连接问题原程序

x=[95.0129,23.1139,60.6843,48.5982,89.1299,76.2097,45.6468,1.8504,82.1407,44.4703,61.5432,79.1937,92.1813,73.8207,17.6266,40.5706,93.547,91.6904,41.027,89.365,5.7891,35.2868,81.3166,0.9861,13.8891,20.2765,19.8722,60.3792,27.2188,19.8814,1.5274,74.6786,44.5096,93.1815,46.5994,41.8649,84.6221,52.5152,20.2647,67.2137,83.8118,1.964,68.1277,37.9481,83.1796,50.2813,70.9471,42.8892,30.4617,18.9654,19.3431,68.2223,30.2764,54.1674,15.0873,69.7898,37.8373,86.0012,85.3655,59.3563,49.6552,89.9769,82.1629,64.491,81.7974,66.0228,34.1971,28.9726,34.1194,53.4079,72.7113,30.929,83.8496,56.8072,37.0414,70.274,54.6571,44.488,69.4567,62.131,79.4821,95.6843,52.259,88.0142,17.2956,97.9747,27.1447,25.2329,87.5742,73.7306,13.6519,1.1757,89.3898,19.9138,29.8723,66.1443,28.4409,46.9224,6.4781,98.8335];
y=[58.2792,42.3496,51.5512,33.3951,43.2907,22.595,57.9807,76.0365,52.9823,64.0526,20.9069,37.9818,78.3329,68.0846,46.1095,56.7829,79.4211,5.9183,60.2869,5.0269,41.5375,30.4999,87.4367,1.5009,76.795,97.0845,99.0083,78.8862,43.8659,49.8311,21.3963,64.3492,32.0036,96.0099,72.6632,41.1953,74.4566,26.7947,43.9924,93.338,68.3332,21.256,83.9238,62.8785,13.3773,20.7133,60.7199,62.9888,37.0477,57.5148,45.1425,4.3895,2.7185,31.2685,1.2863,38.3967,68.3116,9.2842,3.5338,61.2395,60.854,1.576,1.6355,19.0075,58.6918,5.7581,36.7568,63.1451,71.7634,69.2669,8.4079,45.4355,44.1828,35.325,15.3606,67.5645,69.9213,72.7509,47.8384,55.4842,12.1047,45.0754,71.5883,89.2842,27.3102,25.4769,86.5603,23.235,80.4872,90.8398,23.1894,23.9313,4.9754,7.8384,64.0815,19.0887,84.3869,17.39,17.0793,99.4295];
x1=[3.206	17.4571	4.7576 3.206];
y1=[12.9166	19.3377	20 12.9166];
x2=[50	53.7465	46.9222	43.1123 33.3207	50];
y2=[30	48.449	57.1195	 56.3187 39.805 30];
x3=[54.6982	53.7465	46.9222 54.6982];
y3=[70	90	80 70];
x4=[90	80	70 90];
y4=[75	95	80 75];
for i=1:3
    k1(i)=(y1(i+1)-y1(i))/(x1(i+1)-x1(i));
end
for i=1:5
    k2(i)=(y2(i+1)-y2(i))/(x2(i+1)-x2(i));
end
for i=1:3
    k3(i)=(y3(i+1)-y3(i))/(x3(i+1)-x3(i));
end
for i=1:3
    k4(i)=(y4(i+1)-y4(i))/(x4(i+1)-x4(i));
end
for i=1:100
    if k1(1)*(x(i)-x1(1))+y1(1)-y(i)<0 & k1(2)*(x(i)-x1(2))+y1(2)-y(i)>0 & k1(3)*(x(i)-x1(3))+y1(3)-y(i)>0
        fprintf('第%d个点在第一障碍区\t\n',i);
    elseif k2(1)*(x(i)-x2(1))+y2(1)-y(i)<0 & k2(2)*(x(i)-x2(2))+y2(2)-y(i)>0 & k2(3)*(x(i)-x2(3))+y2(3)-y(i)>0 & k2(4)*(x(i)-x2(4))+y2(4)-y(i)>0 & k2(5)*(x(i)-x2(5))+y2(5)-y(i)<0
        fprintf('第%d个点在第二障碍区\t\n',i);
    elseif k3(1)*(x(i)-x3(1))+y3(1)-y(i)>0 & k3(2)*(x(i)-x3(2))+y3(2)-y(i)>0 & k3(3)*(x(i)-x3(3))+y3(3)-y(i)<0
        fprintf('第%d个点在第三障碍区\t\n',i);
    elseif k4(1)*(x(i)-x4(1))+y4(1)-y(i)>0 & k4(2)*(x(i)-x4(2))+y4(2)-y(i)>0 & k4(3)*(x(i)-x4(3))+y4(3)-y(i)<0
        fprintf('第%d个点在第四障碍区\t\n',i);
    end
end