huanyuan.m

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

%a=xlsread('d:\12345','sheet3','a1:b95');
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];
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	87	94	83
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	81	23	42];
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)<23
            b(i,j)=a(i,j)+1;
        elseif a(i,j)>23&&a(i,j)<36
            b(i,j)=a(i,j)+2;
        elseif a(i,j)>36&&a(i,j)<99
            b(i,j)=a(i,j)+3;
        else
            b(i,j)=a(i,j)+4;
        end
    end
end
%xlswrite('d:\12345',b,'sheet3','f1');
for i=1:95
    X(1)=x(b(1,i));
    X(2)=x(b(2,i));
    Y(1)=y(b(1,i));
    Y(2)=x(b(2,i));
    plot(X,Y);
    hold on;
end