tianranqi.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];
nodenumber=[1	2	3	4	5	6	7	8	9	10	11	12	13	14	15	16	17	18	19	20	21	22	23	24	25	26	27	28	29	30	31	32	33	34	35	36	37	38	39	40	41	42	43	44	45	46	47	48	49	50	51	52	53	54	55	56	57	58	59	60	61	62	63	64	65	66	67	68	69	70	71	72	73	74	75	76	77	78	79	80	81	82	83	84	85	86	87	88	89	90	91	92	93	94	95	96	97	98	99	100];
plot(x,y,'*');
for i=1:100
text(x(i)+0.2,y(i)+0.2,num2str(nodenumber(i)));
end
hold on
x1=[3.206	17.4571	4.7576 3.206];
y1=[12.9166	19.3377	20 12.9166];
nodenumber1=[1	2	3];
plot(x1,y1,'r-');
hold on
x2=[50	53.7465	46.9222	43.1123 33.3207	50];
y2=[30	48.449	57.1195	 56.3187 39.805 30];
nodenumber2=[1 2 3 4 5];
plot(x2,y2,'r-');
hold on
x3=[54.6982	53.7465	46.9222 54.6982];
y3=[70	90	80 70];
nodenumber3=[1 2 3];
plot(x3,y3,'r-');
hold on
x4=[90	80	70 90];
y4=[75	95	80 75];
nodenumber4=[1 2 3];
plot(x4,y4,'r-');
hold on
grid on
grid minor
a=[90	39	92	74	48	93	16	45	25	48	75	90	94	89	80	63	45	36	90	61	18	73	60	64	30	78	41	69	77	56	28	27	95	51	71	58	55	66	86	18	44	70	55	54	88	62	80	68	38	64	77	43	96	79	87	46	78	47	70	75	91	92	81	86	87	85	70	62	40	31	53	76	76	35	58	62	95	74	93	35	78	65	91	88	23	85	84	20	79	94	52	72	94	83	12
19	29	65	67	36	61	12	9	24	14	33	17	84	29	67	49	18	2	59	10	15	13	56	46	13	42	18	27	3	19	14	2	43	35	51	6	42	54	12	6	30	4	19	41	82	8	33	49	34	21	57	35	32	4	37	27	68	28	11	54	52	41	32	81	40	82	8	30	26	21	11	53	3	31	57	38	72	26	5	10	5	47	50	39	7	21	24	14	1	66	22	50	23	42	1];
for i=1:2
    for j=1:95
        if a(i,j)<4
            b(i,j)=a(i,j);
        elseif a(i,j)>=4&&a(i,j)<22
            b(i,j)=a(i,j)+1;
        elseif a(i,j)>=22&&a(i,j)<34
            b(i,j)=a(i,j)+2;
        elseif a(i,j)>=34&&a(i,j)<96
            b(i,j)=a(i,j)+3;
        elseif a(i,j)>=96
            b(i,j)=a(i,j)+4;
        end
    end
end
%xlswrite('d:\12345',b,'sheet3','f1');
for i=1:95
    X(1,i)=x(b(1,i));
    X(2,i)=x(b(2,i));
    Y(1,i)=y(b(1,i));
    Y(2,i)=y(b(2,i));
    plot(X,Y,'b-');
    hold on;
end