2012年全国大学生数学建模竞赛A题葡萄酒的评价（含word论文和源代码资源）

亦是远方

已于 2024-08-04 18:44:00 修改

阅读量1.7k

点赞数 1

分类专栏：数学建模文章标签：数学建模

于 2023-06-01 13:22:24 首次发布

本文链接：https://blog.csdn.net/qq_35500719/article/details/130986159

版权

数学建模专栏收录该内容

36 篇文章 17 订阅

订阅专栏

文章目录

零、相关内容
一、部分题目
二、部分论文
三、部分源代码
四、完整word版论文和源代码

零、相关内容

2012年全国大学生数学建模竞赛A题“葡萄酒的评价”旨在通过建立数学模型，分析葡萄酒评价结果的显著性差异，并探讨利用理化指标评估葡萄酒质量的可能性。

在解答这一问题时，需要对提供的评分数据进行统计分析，特别是使用假设检验来确定两组评酒员的评分是否存在显著性差异。根据相关博客和论文的分析，可以采用T检验或Wilcoxon符号秩检验等方法来判断两组评酒数据的一致性[¹][³]。这些检验可以帮助确定哪一组评酒员的评价更加可信，从而为后续的质量评价提供可靠的数据基础。

另外，还需要对酿酒葡萄和葡萄酒的理化指标进行分析，以探究它们与葡萄酒质量之间的关系。这可以通过回归分析、主成分分析(PCA)等方法来实现。例如，可以使用逐步回归分析来筛选出对葡萄酒质量影响最大的理化指标，并通过多级嵌套的方法改进回归模型，以提高预测的准确性[³]。

在确定了理化指标与葡萄酒质量之间的关系后，还可以尝试建立一个综合评价模型，该模型可以基于理化指标来预测葡萄酒的质量等级。这类模型的建立可以参考聚类分析、K-means算法等无监督学习方法，通过比较不同聚类方法的效果，选择最佳的分级策略[⁵]。

综上所述，2012年全国大学生数学建模竞赛A题不仅要求参赛者具备扎实的数学建模能力，还要求能够灵活运用多种数据分析方法来解决实际问题。通过对评酒员评分的统计分析、理化指标与葡萄酒质量之间关系的探索以及构建综合评价模型，可以实现对葡萄酒质量的科学评价。

一、部分题目

2012高教社杯全国大学生数学建模竞赛题目

A题葡萄酒的评价

确定葡萄酒质量时一般是通过聘请一批有资质的评酒员进行品评。每个评酒员在对葡萄酒进行品尝后对其分类指标打分，然后求和得到其总分，从而确定葡萄酒的质量。酿酒葡萄的好坏与所酿葡萄酒的质量有直接的关系，葡萄酒和酿酒葡萄检测的理化指标会在一定程度上反映葡萄酒和葡萄的质量。附件1给出了某一年份一些葡萄酒的评价结果，附件2和附件3分别给出了该年份这些葡萄酒的和酿酒葡萄的成分数据。请尝试建立数学模型讨论下列问题：

1. 分析附件1中两组评酒员的评价结果有无显著性差异，哪一组结果更可信？
2. 根据酿酒葡萄的理化指标和葡萄酒的质量对这些酿酒葡萄进行分级。
3. 分析酿酒葡萄与葡萄酒的理化指标之间的联系。
4．分析酿酒葡萄和葡萄酒的理化指标对葡萄酒质量的影响，并论证能否用葡萄和葡萄酒的理化指标来评价葡萄酒的质量？
附件1：葡萄酒品尝评分表（含4个表格）
附件2：葡萄和葡萄酒的理化指标（含2个表格）
附件3：葡萄和葡萄酒的芳香物质（含4个表格）

二、部分论文

在这里插入图片描述

三、部分源代码

clc;
clear;

% Φ -> A   纬度
% δ -> F   太阳赤道纬度夹角
% ω -> C   太阳时角
%  h -> Oh  太阳高度角
% λ -> D   经度
%

% 10月22日北京时间9:00-15:00

% 10月22日是一年的第 295 天
n = 295.0;

%地理位置 北纬39度54分26秒,东经116度23分29秒
% 纬度
A = 39 + 54/60 + 26/3600;
% 经度
D = 116 + 23/60 + 29/3600;
% 时区经度
Dm = 120;

% 太阳赤纬夹角（度）
F = 23.45*sin(2*pi*(284+n)/365);

% 太阳时
B = 2*pi*(n -81)/364;
E = 9.87*sin(2*B) - 7.53*cos(B) - 1.5*sin(B); 

ii = 0;

LsArr = [];

X = [];

Larr = [3 4 5 6 7 8 9 10]

% T0： m 时 n 分  
% 9:00-15:00
% m = [9, 15]  n = [0, 59]
for m = 9:1:14
    for nn = 0:10:59
        ii = ii +1;
        X(ii) = m+nn/60;
        
        T0 = m + nn/60;
        Ts = T0 + E/60 + (D - Dm)/15;
        % 太阳时角（度）
        C = 15*(Ts - 12);

        % 太阳高度角
        Oh = asin(sin(A*pi/180)*sin(F*pi/180) + cos(A*pi/180)*cos(F*pi/180)*cos(C*pi/180));

        % 杆长 L = 3 m
        L = 3;

        % 影长 Ls
        Ls = L/tan(Oh);
        LsArr(ii) = Ls;
    end
end




plot(X, LsArr, 'red');
%axis([9 15 22 40]);
xlabel('时间');
ylabel('影子长度');
title('北京3米杆长9-15点影长随时间变化曲线');

问题2

（0）第二问主函数

clc;
clear;

% Φ -> A   纬度
% δ -> F   太阳赤道纬度夹角
% ω -> C   太阳时角
%  h -> Oh  太阳高度角
% λ -> D   经度

dB = [0.4555 	0.4409 	0.4247 	0.4136 	0.3986 	0.3919 	0.3777 	0.3656 	0.3582 	0.3481 	0.3438 	0.3305 	0.3264 	0.3169 	0.3120 	0.3069 	0.2987 	0.2928 	0.2876 	0.2853 	];

dLs = [1.1496 	1.1822 	1.2153 	1.2491 	1.2832 	1.3180 	1.3534 	1.3894 	1.4262 	1.4634 	1.5015 	1.5402 	1.5799 	1.6201 	1.6613 	1.7033 	1.7462 	1.7901 	1.8350 	1.8809 	1.9279 ];

% 4月18日是一年的第 108 天
n = 108;

% 太阳赤纬夹角（度）
C = 23.45*sin(2*pi*(284+n)/365);

% 太阳时
B = 2*pi*(n -81)/364;
E = 9.87*sin(2*B) - 7.53*cos(B) - 1.5*sin(B); 

Jm = 120;

hour = [14	14	14	14	14	14	15	15	15	15	15	15	15	15	15	15	15	15	15	15	15];
minutes = [42	45	48	51	54	57	0	3	6	9	12	15	18	21	24	27	30	33	36	39	42];

value1 = 0;
value2 = 0;

min = 100000000;



j = 1;

X = [];
Y = [];

MinArr = ones(12,5);
dFsArr = [];

LsArr = [];

DD = ones(12, 20);

LsArrr = ones(12, 21);
% 杆长
for L = 0:0.1:3
    % 纬度
    for W = 15:0.1:25
        % 经度
        for J = 105:0.1:115
            % 时间
            for i = 1:1:21
                
                T0 = hour(i) + minutes(i)/60;
                Ts = T0 + E/60 + (J - Jm)/15;
            
                % 太阳时角（度）
                S = 15*(Ts - 12);

                % 太阳高度角
                Oh = asin(sin(W*pi/180)*sin(C.*pi/180) + cos(W*pi/180)*cos(C*pi/180)*cos(S*pi/180));

                % 太阳方位角
                if(S <0)
                    Fs = acos(  (sin(C*pi/180) - sin(Oh)*sin(W*pi/180)) / (cos(Oh)*cos(W*pi/180)));
                    FsArr(i) = Fs;
                else
                    Fs = 2*pi - acos(  (sin(C*pi/180) - sin(Oh)*sin(W*pi/180)) / (cos(Oh)*cos(W*pi/180)));
                    FsArr(i) = Fs;
                end
                
                if(i >= 2)
                    value1 = value1 + (  FsArr(i-1) - FsArr(i) - dB(i - 1) )^2;
                    dFsArr(i-1) = FsArr(i-1) - FsArr(i);
                end
                
                % 影长 Ls
                Ls = L / tan(Oh);
                
                LsArr(i) = Ls;
                
                value2 = value2 + (Ls - dLs(i))^2;


                i = i +1;
                
            end
            
            value = value1/20*value2/21;
            %if(value < min)
                %min = value;
                %MinArr = [L W J]
            %end
            
            if (value < 0.0000003)
                X(j) = j;
                j
                Y(j) = value;
                LsArrr(j, :) = LsArr;
                DD(j, :) = dFsArr*180/pi;
                MinArr(j , 1:5) = [j L W J value];
                j = j+1;
            end
            
            
            
            value1 = 0;
            value2 = 0;
            value = 0;
            
               
        end
    end
end




plot(X, Y,'*');

（1）经度灵敏度分析

clc;
clear;

% Φ -> A   纬度
% δ -> F   太阳赤道纬度夹角
% ω -> C   太阳时角
%  h -> Oh  太阳高度角
% λ -> D   经度
%

% 10月22日北京时间9:00-15:00

% 10月22日是一年的第 295 天
n = 295.0;




%地理位置 北纬39度54分26秒,东经116度23分29秒
% 纬度
A = 40;
% 经度

D2 = 116;
D1 = D2 - 2;
D3 = D2 + 2;
% 时区经度
Dm = 120;

% 太阳赤纬夹角（度）
F = 23.45*sin(2*pi*(284+n)/365);

% 太阳时
B = 2*pi*(n -81)/364;
E = 9.87*sin(2*B) - 7.53*cos(B) - 1.5*sin(B); 

i = 0;

LsArr1 = [];
LsArr2 = [];
LsArr3 = [];



% T0： m 时 n 分  
% 9:00-15:00
% m = [9, 15]  n = [0, 59]
for m = 9:1:14
    for nn = 0:10:59
        i = i +1;
        I(i) = m+nn/60;
        
        T0 = m + nn/60;
        
        Ts1 = T0 + E/60 + (D1 - Dm)/15;
        Ts2 = T0 + E/60 + (D2 - Dm)/15;
        Ts3 = T0 + E/60 + (D3 - Dm)/15;
        % 太阳时角（度）
        C1 = 15*(Ts1 - 12);
        C2 = 15*(Ts2 - 12);
        C3 = 15*(Ts3 - 12);


        % 太阳高度角
        Oh1 = asin(sin(A*pi/180)*sin(F*pi/180) + cos(A*pi/180)*cos(F*pi/180)*cos(C1*pi/180));
        Oh2 = asin(sin(A*pi/180)*sin(F*pi/180) + cos(A*pi/180)*cos(F*pi/180)*cos(C2*pi/180));
        Oh3 = asin(sin(A*pi/180)*sin(F*pi/180) + cos(A*pi/180)*cos(F*pi/180)*cos(C3*pi/180));



        % 杆长 L = 3 m
        L = 3;
        
        
        
        L3 = 3.1;

        % 影长 Ls
        Ls1 = L/tan(Oh1);
        Ls2 = L/tan(Oh2);
        Ls3 = L/tan(Oh3);
        LsArr1(i) = Ls1;
        LsArr2(i) = Ls2;
        LsArr3(i) = Ls3;
    end
end


plot(I, LsArr1);
hold on;
plot(I, LsArr2);
hold on;
plot(I, LsArr3);
%axis([9 15 22 40]);
xlabel('时间');
ylabel('影子长度');
title('经度灵敏度分析');
legend('114°E' ,'116°E','118°E');

（2）维度灵敏度分析

clc;
clear;

% Φ -> A   纬度
% δ -> F   太阳赤道纬度夹角
% ω -> C   太阳时角
%  h -> Oh  太阳高度角
% λ -> D   经度
%

% 10月22日北京时间9:00-15:00

% 10月22日是一年的第 295 天
n = 295.0;




%地理位置 北纬39度54分26秒,东经116度23分29秒
% 纬度
A1 = 39;
A2 = 40;
A3 = 41;


% 经度
D = 116;

% 时区经度
Dm = 120;

% 太阳赤纬夹角（度）
F = 23.45*sin(2*pi*(284+n)/365);

% 太阳时
B = 2*pi*(n -81)/364;
E = 9.87*sin(2*B) - 7.53*cos(B) - 1.5*sin(B); 

i = 0;

LsArr1 = [];
LsArr2 = [];
LsArr3 = [];



% T0： m 时 n 分  
% 9:00-15:00
% m = [9, 15]  n = [0, 59]
for m = 9:1:14
    for nn = 0:10:59
        i = i +1;
        I(i) = m+nn/60;
        
        T0 = m + nn/60;
        
        Ts = T0 + E/60 + (D - Dm)/15;
        % 太阳时角（度）
        C = 15*(Ts - 12);


        % 太阳高度角
        Oh1 = asin(sin(A1*pi/180)*sin(F*pi/180) + cos(A1*pi/180)*cos(F*pi/180)*cos(C*pi/180));
        Oh2 = asin(sin(A2*pi/180)*sin(F*pi/180) + cos(A2*pi/180)*cos(F*pi/180)*cos(C*pi/180));
        Oh3 = asin(sin(A3*pi/180)*sin(F*pi/180) + cos(A3*pi/180)*cos(F*pi/180)*cos(C*pi/180));



        % 杆长 L = 3 m
        L = 3;
        
        
        
        L3 = 3.1;

        % 影长 Ls
        Ls1 = L/tan(Oh1);
        Ls2 = L/tan(Oh2);
        Ls3 = L/tan(Oh3);
        LsArr1(i) = Ls1;
        LsArr2(i) = Ls2;
        LsArr3(i) = Ls3;
    end
end


plot(I, LsArr1, '*');
%axis([9 15 22 40]);
xlabel('时间');
ylabel('影子长度');
title('纬度灵敏度分析');
legend('38°N' ,'40°N','42°N');

问题3

clc;
clear;

% Φ -> A   纬度
% δ -> F   太阳赤道纬度夹角
% ω -> C   太阳时角
%  h -> Oh  太阳高度角
% λ -> D   经度
%

% 时间
hour = [13	13	13	13	13	13	13	13	13	13	13	13	13	13	13	13	13	14	14	14	14
];
minutes = [9	12	15	18	21	24	27	30	33	36	39	42	45	48	51	54	57	0	3	6	9
];


Fs = [0.909541043	0.906049014	0.900872887	0.897969489	0.893046127	0.887167053	0.884149304	0.876983995	0.872571801	0.866083489	0.860234734	0.853663332	0.847169536	0.841329236	0.833261206	0.827336861	0.819957272	0.813116998	0.805282613	0.798637273
];

Ls = [3.533142184	3.546768029	3.561797643	3.578100715	3.595750783	3.61493428	3.635425983	3.657218272	3.680541115	3.705167836	3.731278025	3.758917911	3.788087888	3.818701015	3.850809619	3.88458522	3.919911828	3.956875992	3.99553479	4.035750835	4.077863059
];

FsArr = [0.880834641	0.877278202	0.873429185	0.869296509	0.864889603	0.86021836	0.855293081	0.850124425	0.844723357	0.839101093	0.833269051	0.827238796	0.821021996	0.81463037	0.808075644	0.801369507	0.794523573	0.787549336	0.78045814	0.773261144
0.942777771	0.938816179	0.93450787	0.929864071	0.924896687	0.919618231	0.914041753	0.908180764	0.902049163	0.895661161	0.889031212	0.882173936	0.875104057	0.867836335	0.860385502	0.852766202	0.844992939	0.837080025	0.829041531	0.820891243
0.857881568	0.854766987	0.851416033	0.847835574	0.844032851	0.840015435	0.835791197	0.831368273	0.826755025	0.821960009	0.81699194	0.811859655	0.806572083	0.801138208	0.795567044	0.789867595	0.784048838	0.778119686	0.772088968	0.765965404
0.811583635	0.808971923	0.806162644	0.803160921	0.799972154	0.796602007	0.793056378	0.789341381	0.785463315	0.781428647	0.777243984	0.772916051	0.768451664	0.763857715	0.759141139	0.754308902	0.749367974	0.744325313	0.739187842	0.733962434
0.918172984	0.914354463	0.910206395	0.905739192	0.900963889	0.895892085	0.890535871	0.884907771	0.879020667	0.87288774	0.8665224	0.859938224	0.853148893	0.846168134	0.83900966	0.831687121	0.824214048	0.81660381	0.808869568	0.801024239
0.829406698	0.826764921	0.823932874	0.820915692	0.817718772	0.814347759	0.810808516	0.807107108	0.803249773	0.799242903	0.795093017	0.790806742	0.786390788	0.781851926	0.777196968	0.772432747	0.767566095	0.762603822	0.757552707	0.752419469
1.16684882	1.161542901	1.155756961	1.149510206	1.142822972	1.135716575	1.128213157	1.12033553	1.112107026	1.103551347	1.094692424	1.085554272	1.076160866	1.066536017	1.056703255	1.04668573	1.036506114	1.026186518	1.015748419	1.005212594
0.776438352	0.774972254	0.773400039	0.771723843	0.769945926	0.768068664	0.766094543	0.764026149	0.761866166	0.759617365	0.757282599	0.754864792	0.752366936	0.749792082	0.74714333	0.744423826	0.741636753	0.738785324	0.735872772	0.732902351
0.986400514	0.982380236	0.978059427	0.973448721	0.968559276	0.963402715	0.95799106	0.952336676	0.946452202	0.940350494	0.934044566	0.92754753	0.920872543	0.914032748	0.907041234	0.899910979	0.892654811	0.885285366	0.877815052	0.870256013
0.908336067	0.905061021	0.901561466	0.897844594	0.893917931	0.889789305	0.885466806	0.880958755	0.876273669	0.871420224	0.866407222	0.86124356	0.855938194	0.850500108	0.844938288	0.83926169	0.833479212	0.827599673	0.821631784	0.815584131
0.851786397	0.848670393	0.845315849	0.841729649	0.837919054	0.833891659	0.829655364	0.825218336	0.820588975	0.815775875	0.810787793	0.80563361	0.800322302	0.794862902	0.789264472	0.783536068	0.777686718	0.771725388	0.765660957	0.759502198
1.15080961	1.145488865	1.139714088	1.133503829	1.126877646	1.119855966	1.112459944	1.104711324	1.096632294	1.088245354	1.079573184	1.070638516	1.061464021	1.052072191	1.042485241	1.032725012	1.022812889	1.012769719	1.002615754	0.992370583
0.766445573	0.765410907	0.764302819	0.763122633	0.761871748	0.760551639	0.759163853	0.757710001	0.756191761	0.754610869	0.752969117	0.75126835	0.749510462	0.747697392	0.745831119	0.743913659	0.74194706	0.739933402	0.737874789	0.735773346
0.949103285	0.945024226	0.940573036	0.935761792	0.930603338	0.925111202	0.919299514	0.91318292	0.906776503	0.900095693	0.89315619	0.885973881	0.878564761	0.870944862	0.86313018	0.85513661	0.84697988	0.838675499	0.830238699	0.821684391
0.872482366	0.869707792	0.86675363	0.863625236	0.860328206	0.856868357	0.8532517	0.849484418	0.845572847	0.841523448	0.837342789	0.83303752	0.828614353	0.82408004	0.819441355	0.814705071	0.809877945	0.804966698	0.799977998	0.794918447
1.165164421	1.159898002	1.154165688	1.147986183	1.141379263	1.134365631	1.126966769	1.119204796	1.111102318	1.10268229	1.093967877	1.084982325	1.075748834	1.066290444	1.056629926	1.046789685	1.036791665	1.026657277	1.01640732	1.006061927
0.77193515	0.77043976	0.768835757	0.767125332	0.7653108	0.763394597	0.761379272	0.759267475	0.757061957	0.754765559	0.752381202	0.749911884	0.747360669	0.744730682	0.742025097	0.739247136	0.736400056	0.733487143	0.730511706	0.727477072
0.905602462	0.901903458	0.897896226	0.89359042	0.888996252	0.884124439	0.878986138	0.873592896	0.867956582	0.862089335	0.856003504	0.849711589	0.84322619	0.836559953	0.829725517	0.822735471	0.815602302	0.808338359	0.800955811	0.793466614
0.898373814	0.89494397	0.891252409	0.887307273	0.883117137	0.878690961	0.874038046	0.869167994	0.864090661	0.85881611	0.853354573	0.847716406	0.841912046	0.835951973	0.829846672	0.823606596	0.817242129	0.810763559	0.804181042	0.797504579
0.835269284	0.832968563	0.830521393	0.827931731	0.825203718	0.822341661	0.819350021	0.816233394	0.8129965	0.809644165	0.80618131	0.802612929	0.798944084	0.79517988	0.791325461	0.787385989	0.783366636	0.779272569	0.775108935	0.770880857
0.756511095	0.755043914	0.753460159	0.751762019	0.749951824	0.748032031	0.746005227	0.743874111	0.741641495	0.739310289	0.736883498	0.734364211	0.731755593	0.729060877	0.726283358	0.723426379	0.720493331	0.717487637	0.714412751	0.711272144
0.923487372	0.919571891	0.915303424	0.91069311	0.905752792	0.900494943	0.89493259	0.889079246	0.882948824	0.876555573	0.869913998	0.86303879	0.855944755	0.848646747	0.841159605	0.833498093	0.825676839	0.817710288	0.80961265	0.801397857
0.933325196	0.929582055	0.925546146	0.921227055	0.916634881	0.911780179	0.906673905	0.901327362	0.895752142	0.889960068	0.883963149	0.877773514	0.871403374	0.864864961	0.858170489	0.851332105	0.84436185	0.837271617	0.83007312	0.822777855
0.796767395	0.794780956	0.792656098	0.790396082	0.788004342	0.785484471	0.78284021	0.780075436	0.777194148	0.774200459	0.771098578	0.7678928	0.764587496	0.761187096	0.757696082	0.754118971	0.750460307	0.746724651	0.742916565	0.739040608
1.243502726	1.237845026	1.231594153	1.224773755	1.217409131	1.209527011	1.20115534	1.192323051	1.183059844	1.173395968	1.163362015	1.152988713	1.142306738	1.131346539	1.120138171	1.108711145	1.0970943	1.085315677	1.073402422	1.061380695
1.175450272	1.169844136	1.16378715	1.157298795	1.150399515	1.14311057	1.13545389	1.127451928	1.119127519	1.110503745	1.101603797	1.092450854	1.083067965	1.073477935	1.063703229	1.053765874	1.043687383	1.03348868	1.023190034	1.012811008
0.775855532	0.77455824	0.773168066	0.771686813	0.770116388	0.768458793	0.766716127	0.764890575	0.762984405	0.760999963	0.758939664	0.756805988	0.754601477	0.752328723	0.749990369	0.747589097	0.745127626	0.742608706	0.740035112	0.737409637
0.992982651	0.988966022	0.984669573	0.980103525	0.975278558	0.970205754	0.964896541	0.959362637	0.953615993	0.947668739	0.94153313	0.935221494	0.928746183	0.922119525	0.915353778	0.908461092	0.901453464	0.89434271	0.887140422	0.879857946
0.872535306	0.86976433	0.866813662	0.863688655	0.860394899	0.856938208	0.853324587	0.849560219	0.845651434	0.841604693	0.837426562	0.833123688	0.828702783	0.8241706	0.81953391	0.814799488	0.80997409	0.805064439	0.800077205	0.795018989
0.920623503	0.917130591	0.913362934	0.909329138	0.905038279	0.900499854	0.895723737	0.890720126	0.885499497	0.880072556	0.874450187	0.868643412	0.86266334	0.856521126	0.850227927	0.843794865	0.837232987	0.830553228	0.823766383	0.816883073
0.999820957	0.995561315	0.990929205	0.98593744	0.980599596	0.974929927	0.968943278	0.962654994	0.956080835	0.949236887	0.942139478	0.934805094	0.927250298	0.919491657	0.911545667	0.903428684	0.895156866	0.886746113	0.878212012	0.869569793
0.974520982	0.9702612	0.965578247	0.960485875	0.954998796	0.949132578	0.942903541	0.936328652	0.929425422	0.922211795	0.914706048	0.906926689	0.898892358	0.890621736	0.882133453	0.873446007	0.864577684	0.855546491	0.846370085	0.83706572
0.851079291	0.848245613	0.845214512	0.841991661	0.838583007	0.834994755	0.831233335	0.827305381	0.823217701	0.818977256	0.814591131	0.81006651	0.805410652	0.800630869	0.795734499	0.790728886	0.785621359	0.780419209	0.775129674	0.769759915
0.863351547	0.860115825	0.856624763	0.852885826	0.848906894	0.844696225	0.840262415	0.835614356	0.830761194	0.825712295	0.820477198	0.815065578	0.809487211	0.80375193	0.797869596	0.791850057	0.785703121	0.779438521	0.773065891	0.766594732
0.912638993	0.909284623	0.90567833	0.901827969	0.897741805	0.893428469	0.888896919	0.884156396	0.879216385	0.874086572	0.868776803	0.863297045	0.857657343	0.851867789	0.845938479	0.839879482	0.833700809	0.827412376	0.821023983	0.814545282
1.102038464	1.097081067	1.091756867	1.086080936	1.080069034	1.073737509	1.067103203	1.06018335	1.052995478	1.045557323	1.037886728	1.030001567	1.021919654	1.013658672	1.0052361	0.996669145	0.987974686	0.979169214	0.970268793	0.961289011
0.79139496	0.789566326	0.787611916	0.785534605	0.783337414	0.781023507	0.77859618	0.776058849	0.773415042	0.770668385	0.767822597	0.764881475	0.761848884	0.758728751	0.755525047	0.752241787	0.748883011	0.74545278	0.741955164	0.738394237
0.992724323	0.988455597	0.983770915	0.978683987	0.973209457	0.967362804	0.961160235	0.954618585	0.947755212	0.940587894	0.933134727	0.925414025	0.917444225	0.909243797	0.900831153	0.892224571	0.883442115	0.874501567	0.865420364	0.85621554
0.773376087	0.772333578	0.771223491	0.770047126	0.76880585	0.767501098	0.766134369	0.764707217	0.763221257	0.761678155	0.760079626	0.758427431	0.756723375	0.7549693	0.753167084	0.751318638	0.7494259	0.747490834	0.745515426	0.743501679
0.923412893	0.919490277	0.915214889	0.910597884	0.905651117	0.900387073	0.894818792	0.888959791	0.882823994	0.876425654	0.869779282	0.86289957	0.855801326	0.848499406	0.841008648	0.833343811	0.825519522	0.817550223	0.809450116	0.801233129
1.034975228	1.030607771	1.025822418	1.020633403	1.015055898	1.009105917	1.002800206	0.996156133	0.989191588	0.981924866	0.974374572	0.966559511	0.958498596	0.950210748	0.941714815	0.933029481	0.924173196	0.915164102	0.90601997	0.896758141
0.785562609	0.783395484	0.781069396	0.778588103	0.775955568	0.773175946	0.770253571	0.767192936	0.763998683	0.760675584	0.757228523	0.753662489	0.749982551	0.746193848	0.742301571	0.738310953	0.734227249	0.730055728	0.725801653	0.721470275
0.875951652	0.873183945	0.870240906	0.867127823	0.863850215	0.860413813	0.856824533	0.853088456	0.849211809	0.84520094	0.841062294	0.836802399	0.832427838	0.827945231	0.823361217	0.818682434	0.813915503	0.809067005	0.804143473	0.799151371
0.968341856	0.964191426	0.959645226	0.954716096	0.949417734	0.943764604	0.93777184	0.931455161	0.924830766	0.91791525	0.910725507	0.903278639	0.895591874	0.887682478	0.879567677	0.871264585	0.86279013	0.854160996	0.84539356	0.83650384
0.823061176	0.820214297	0.817157136	0.813895483	0.810435439	0.806783377	0.802945928	0.798929942	0.794742469	0.790390725	0.785882069	0.781223972	0.776423993	0.771489753	0.766428909	0.761249127	0.755958065	0.750563345	0.745072535	0.739493129
0.795622377	0.794309436	0.79292292	0.791464512	0.789935972	0.788339128	0.786675874	0.784948164	0.78315801	0.781307476	0.779398672	0.777433756	0.77541492	0.773344395	0.771224439	0.769057338	0.766845401	0.764590951	0.762296328	0.75996388
1.04039069	1.035881605	1.03104631	1.025897463	1.020448299	1.014712556	1.008704395	1.002438322	0.995929118	0.989191758	0.982241345	0.975093037	0.967761984	0.960263263	0.95261182	0.944822419	0.936909589	0.928887579	0.920770317	0.912571375
0.95499944	0.95088855	0.94637825	0.941481282	0.936211257	0.930582567	0.924610292	0.918310106	0.91169818	0.904791093	0.897605735	0.890159218	0.88246879	0.874551747	0.866425355	0.858106777	0.849612999	0.840960768	0.832166528	0.823246374
1.097463232	1.092605612	1.087282617	1.08151135	1.075310019	1.068697805	1.061694729	1.054321514	1.046599443	1.038550234	1.030195898	1.02155862	1.012660632	1.003524104	0.994171032	0.984623138	0.974901785	0.965027888	0.955021843	0.944903463
1.134447758	1.12948065	1.12401669	1.118074226	1.111672849	1.104833251	1.097577065	1.089926712	1.081905247	1.073536207	1.06484346	1.055851065	1.046583133	1.037063696	1.027316592	1.017365346	1.007233072	0.996942382	0.986515302	0.9759732
];

LsArr =[3.533527414	3.547020355	3.561858157	3.578051826	3.595613547	3.614556735	3.634896098	3.656647695	3.679829014	3.704459044	3.730558354	3.75814919	3.787255564	3.817903362	3.850120453	3.883936811	3.919384646	3.956498542	3.99531561	4.035875651	4.07822133
3.534475715	3.547686811	3.562292509	3.578301742	3.595724463	3.614571704	3.634855628	3.656589587	3.679788194	3.704467389	3.730644518	3.758338416	3.787569491	3.818359826	3.850733273	3.884715566	3.920334438	3.957619743	3.996603594	4.037320503	4.079807537
3.532933249	3.546860501	3.562080563	3.578606619	3.596453168	3.615636081	3.63617267	3.65808175	3.681383723	3.706100656	3.732256372	3.759876549	3.788988824	3.819622911	3.851810722	3.885586503	3.92098698	3.958051516	3.996822287	4.037344462	4.079666408
3.534510095	3.548323175	3.56340881	3.579782284	3.597460383	3.616461455	3.636805486	3.658514181	3.681611044	3.706121481	3.7320729	3.75949482	3.788419001	3.818879574	3.850913185	3.884559157	3.91985966	3.956859902	3.995608329	4.036156854	4.078561096
3.533936913	3.547178293	3.561798398	3.577806811	3.595214191	3.614032326	3.634274188	3.655953996	3.679087283	3.703690966	3.729783428	3.757384599	3.786516048	3.817201082	3.849464847	3.883334442	3.918839041	3.956010022	3.994881105	4.035488501	4.077871079
3.530398238	3.544576367	3.560010361	3.576715602	3.594708935	3.614008738	3.634634991	3.656609354	3.679955254	3.704697976	3.730864768	3.758484949	3.787590027	3.818213836	3.850392675	3.884165462	3.919573903	3.956662675	3.99547963	4.036076005	4.078506666
3.53345517	3.546621747	3.561232455	3.577288964	3.594793279	3.613747786	3.634155299	3.656019114	3.679343061	3.704131565	3.730389699	3.758123253	3.787338792	3.818043725	3.850246373	3.883956039	3.91918308	3.955938984	3.994236443	4.034089435	4.075513306
3.53452323	3.548435418	3.563574086	3.57995946	3.597613672	3.616560851	3.636827221	3.658441209	3.681433564	3.705837492	3.731688793	3.759026025	3.787890677	3.818327353	3.850383988	3.884112074	3.91956691	3.956807879	3.995898754	4.036908026	4.079909281
3.53026871	3.544390283	3.559816918	3.57655728	3.594620945	3.614018436	3.634761282	3.656862067	3.680334488	3.705193416	3.731454968	3.759136569	3.788257036	3.818836656	3.850897272	3.884462376	3.919557209	3.956208862	3.994446397	4.034300958	4.075805908
3.530584706	3.545093692	3.56084406	3.577847973	3.596118745	3.615670887	3.636520169	3.658683678	3.682179885	3.707028719	3.733251642	3.760871735	3.789913787	3.820404394	3.852372064	3.885847332	3.92086288	3.957453671	3.995657091	4.035513105	4.077064423
3.533564657	3.547443652	3.562620091	3.579107178	3.596919442	3.616072795	3.636584595	3.658473718	3.681760636	3.706467494	3.732618211	3.760238569	3.789356326	3.820001329	3.85220564	3.88600367	3.921432331	3.958531192	3.99734265	4.037912125	4.080288255
3.531214699	3.544830691	3.559867207	3.576326515	3.594211274	3.613524576	3.634269996	3.656451641	3.680074206	3.705143027	3.731664144	3.759644356	3.78909129	3.820013462	3.85242035	3.886322459	3.921731396	3.958659944	3.997122141	4.037133358	4.078710386
3.531382258	3.545454462	3.560745039	3.577276439	3.595073212	3.6141621	3.634572162	3.656334889	3.679484348	3.704057334	3.730093533	3.75763571	3.78672991	3.817425677	3.849776306	3.883839104	3.919675695	3.957352337	3.996940292	4.03851622	4.082162619
3.536100143	3.549083055	3.563488215	3.579323886	3.596599313	3.615324779	3.635511657	3.657172471	3.680320964	3.704972167	3.731142475	3.758849729	3.788113305	3.818954207	3.851395166	3.885460751	3.92117748	3.958573945	3.997680943	4.038531616	4.081161605
3.530031044	3.544731273	3.560627202	3.577732937	3.596063849	3.615636627	3.63646934	3.658581508	3.681994166	3.706729952	3.732813183	3.760269957	3.789128244	3.819418	3.851171282	3.884422376	3.919207933	3.95556712	3.993541778	4.033176603	4.074519328
3.535171439	3.548382913	3.563014074	3.579066538	3.596542247	3.615443509	3.635773048	3.657534049	3.680730216	3.705365823	3.73144577	3.758975645	3.787961784	3.818411332	3.850332311	3.883733687	3.918625438	3.955018623	3.992925461	4.0323594	4.0733352
3.531149353	3.545170377	3.560432334	3.576955893	3.59476368	3.613880373	3.634332801	3.65615006	3.679363634	3.704007534	3.730118445	3.757735892	3.786902419	3.817663791	3.850069204	3.884171533	3.920027585	3.957698395	3.997249537	4.03875148	4.082279966
3.535046384	3.54850927	3.563331059	3.579522087	3.597093814	3.616058879	3.636431157	3.658225819	3.681459406	3.706149898	3.732316796	3.75998121	3.789165947	3.819895619	3.852196741	3.886097857	3.921629656	3.958825115	3.997719634	4.038351201	4.080760556
3.532798479	3.546762202	3.562032146	3.578620047	3.596538832	3.615802674	3.636427048	3.6584288	3.681826214	3.706639088	3.732888819	3.760598488	3.78979296	3.820498985	3.852745309	3.8865628	3.92198457	3.959046125	3.99778551	4.038243474	4.080463651
3.528395582	3.543189675	3.559167253	3.576344742	3.59474001	3.614372435	3.635262972	3.657434234	3.680910578	3.705718193	3.731885206	3.759441786	3.788420267	3.818855278	3.850783878	3.884245716	3.919283193	3.955941644	3.994269537	4.034318688	4.07614449
3.535986686	3.549355452	3.563998239	3.579935537	3.597189887	3.615785981	3.635750758	3.657113529	3.679906099	3.704162913	3.729921209	3.757221194	3.786106226	3.816623028	3.848821911	3.882757031	3.918486658	3.956073488	3.995584972	4.037093689	4.080677755
3.537418182	3.5503808	3.564745981	3.58052258	3.597720482	3.616350655	3.636425207	3.657957445	3.680961941	3.705454606	3.731452762	3.75897523	3.788042416	3.818676406	3.85090107	3.884742169	3.920227476	3.957386899	3.996252619	4.036859229	4.0792439
3.532070954	3.545969744	3.56119982	3.577771529	3.595696307	3.614986724	3.635656546	3.657720791	3.6811958	3.706099306	3.73245051	3.760270167	3.789580673	3.820406162	3.852772609	3.886707943	3.922242162	3.959407466	3.998238392	4.038771962	4.081047848
3.529220088	3.543534868	3.559078232	3.575868558	3.593925924	3.613272186	3.633931067	3.655928247	3.679291467	3.704050645	3.730237997	3.757888174	3.787038406	3.817728669	3.850001857	3.883903978	3.919484364	3.956795899	3.995895276	4.036843273	4.079705055
3.535937979	3.548275819	3.562160393	3.57759074	3.594565968	3.613085305	3.633148155	3.65475416	3.677903258	3.70259575	3.728832366	3.756614334	3.785943452	3.81682216	3.849253608	3.883241738	3.918791356	3.955908206	3.994599049	4.034871743	4.076735318
3.530541457	3.544531372	3.559910467	3.576679779	3.594840616	3.614394605	3.635343734	3.657690404	3.681437475	3.706588323	3.733146894	3.761117755	3.790506158	3.821318097	3.853560367	3.88724063	3.922367475	3.958950484	3.997000303	4.036528703	4.077548659
3.53456909	3.548435364	3.563511144	3.579816988	3.597375375	3.616210795	3.636349854	3.657821376	3.680656531	3.704888963	3.730554942	3.757693516	3.786346696	3.816559642	3.848380879	3.881862529	3.917060563	3.954035088	3.99285065	4.03357658	4.076287361
3.531232118	3.545695382	3.561419345	3.578412463	3.596684042	3.616244278	3.637104306	3.659276251	3.682773278	3.707609652	3.733800798	3.761363368	3.790315309	3.820675937	3.852466019	3.885707856	3.920425371	3.956644208	3.994391834	4.033697646	4.074593091
3.529752745	3.544429499	3.560301687	3.577383391	3.595689955	3.615238042	3.636045695	3.658132402	3.68151917	3.706228601	3.732284982	3.759714371	3.788544703	3.818805892	3.850529953	3.883751127	3.918506017	3.954833737	3.992776076	4.03237767	4.073686192
3.535884113	3.549527967	3.564480852	3.5807535	3.598357763	3.617306658	3.637614424	3.65929658	3.682369998	3.706852965	3.732765266	3.760128267	3.788965	3.819300264	3.851160728	3.88457504	3.919573951	3.956190442	3.994459865	4.034420095	4.07611169
3.531613342	3.544916547	3.559627359	3.575753333	3.593302919	3.612285509	3.632711491	3.654592308	3.677940519	3.702769867	3.72909535	3.756933299	3.78630146	3.817219077	3.849706994	3.883787747	3.91948567	3.956827013	3.995840055	4.03655524	4.079005304
3.537014443	3.549466439	3.563391669	3.578796989	3.595690142	3.614079811	3.633975668	3.655388442	3.678329978	3.702813312	3.728852741	3.756463909	3.78566389	3.816471279	3.848906292	3.882990868	3.91874878	3.956205753	3.995389591	4.036330308	4.07906027
3.529282536	3.543751022	3.55946969	3.576453266	3.594717857	3.614281019	3.635161816	3.657380901	3.680960593	3.705924967	3.732299947	3.760113412	3.789395305	3.820177762	3.852495237	3.88638465	3.921885542	3.959040246	3.99789407	4.038495496	4.080896404
3.535718629	3.549376265	3.564344466	3.580635599	3.5982633	3.617242531	3.637589643	3.659322441	3.68246026	3.707024041	3.733036426	3.760521841	3.789506608	3.820019049	3.85208961	3.885750985	3.92103826	3.957989064	3.996643729	4.037045472	4.079240585
3.535199432	3.549087959	3.56425692	3.580717697	3.59848281	3.617565974	3.637982148	3.659747602	3.682879981	3.707398375	3.733323399	3.760677276	3.789483928	3.81976907	3.851560323	3.884887316	3.919781821	3.956277872	3.994411917	4.034222967	4.075752762
3.531328478	3.545563909	3.561111144	3.57797404	3.596156948	3.615664749	3.636502896	3.658677464	3.68219519	3.707063528	3.733290695	3.76088573	3.789858549	3.820220001	3.851981935	3.885157259	3.919760008	3.955805415	3.993309983	4.03229156	4.072769419
3.528888524	3.543292591	3.558918918	3.575786677	3.593916796	3.613332044	3.634057121	3.656118758	3.679545822	3.704369438	3.730623117	3.7583429	3.787567513	3.818338536	3.850700592	3.884701549	3.920392741	3.95782922	3.997070014	4.038178429	4.08122237
3.53798433	3.550570438	3.564620818	3.580142343	3.59714276	3.615630747	3.635615968	3.657109129	3.680122044	3.704667705	3.730760358	3.758415579	3.787650361	3.818483204	3.850934213	3.885025197	3.920779782	3.958223525	3.997384035	4.038291113	4.080976882
3.528414623	3.54304422	3.558851291	3.575858078	3.59408879	3.613569696	3.634329235	3.65639813	3.679809524	3.704599113	3.730805312	3.758469419	3.787635804	3.818352118	3.850669516	3.884642909	3.920331232	3.957797749	3.997110381	4.038342071	4.081571187
3.537662748	3.550653224	3.565046429	3.580851239	3.59807756	3.616736382	3.636839836	3.658401253	3.681435233	3.705957712	3.731986042	3.759539076	3.78863725	3.819302685	3.851559289	3.885432861	3.920951216	3.958144306	3.997044356	4.037686013	4.080106497
3.538503548	3.551165491	3.565278047	3.580847321	3.597880206	3.616384431	3.636368615	3.657842321	3.680816122	3.70530166	3.73131172	3.758860307	3.787962718	3.818635633	3.850897203	3.88476714	3.920266826	3.957419408	3.996249923	4.036785406	4.079055028
3.529615547	3.543683347	3.559011545	3.575617998	3.593522272	3.612745722	3.633311573	3.655245018	3.678573319	3.703325916	3.729534559	3.757233435	3.78645932	3.817251736	3.849653134	3.883709081	3.919468471	3.95698376	3.996311212	4.037511184	4.080648419
3.529620348	3.544408837	3.560379599	3.577546612	3.595925088	3.61553153	3.636383791	3.658501138	3.681904324	3.706615666	3.732659122	3.760060388	3.788846992	3.819048398	3.85069612	3.883823849	3.918467582	3.954665765	3.992459457	4.031892489	4.073011657
3.535748089	3.548417936	3.562531834	3.578097278	3.595122698	3.613617501	3.633592131	3.655058128	3.678028192	3.702516254	3.728537551	3.756108702	3.785247802	3.815974507	3.848310136	3.882277771	3.917902375	3.955210907	3.994232447	4.034998339	4.077542329
3.528900415	3.543039929	3.558460178	3.57517604	3.593203843	3.612561424	3.633268205	3.655345269	3.678815441	3.703703386	3.730035706	3.757841052	3.787150237	3.817996374	3.850415008	3.884444277	3.920125071	3.957501214	3.996619665	4.037530726	4.080288281
3.528687648	3.543970961	3.560352716	3.577852712	3.596492343	3.616294683	3.637284561	3.659488658	3.682935609	3.707656104	3.733683018	3.761051531	3.789799277	3.819966494	3.851596195	3.884734351	3.919430092	3.955735925	3.993707978	4.033406257	4.074894942
3.527747281	3.542181712	3.557919658	3.574967703	3.593333151	3.613024065	3.634049314	3.656418621	3.680142616	3.705232891	3.731702055	3.759563804	3.788832977	3.819525633	3.85165912	3.885252157	3.920324909	3.956899083	3.99499801	4.034646752	4.075872199
3.538888578	3.551396168	3.565352158	3.580763998	3.597640073	3.615989752	3.635823443	3.657152654	3.679990055	3.704349548	3.730246346	3.75769705	3.786719736	3.817334045	3.849561286	3.883424535	3.918948754	3.956160903	3.995090073	4.035767617	4.078227304
3.53274082	3.545623948	3.559991015	3.575846297	3.593194681	3.612041713	3.632393653	3.65425753	3.677641202	3.702553424	3.729003913	3.757003422	3.786563816	3.817698151	3.850420757	3.88474733	3.92069502	3.958282527	3.997530207	4.038460175	4.081096417
3.538220064	3.55075892	3.56479903	3.580343369	3.597395391	3.615959083	3.636039014	3.657640391	3.68076912	3.705431868	3.73163613	3.759390298	3.788703737	3.819586857	3.852051194	3.886109496	3.921775804	3.959065546	3.997995628	4.038584532	4.08085242
];

FsArrD = [0.914456329	0.910906199	0.907077393	0.902978699	0.898619382	0.894009134	0.889158025	0.884076457	0.878775108	0.873264887	0.867556882	0.861662314	0.855592491	0.849358761	0.842972469	0.836444921	0.829787338	0.823010829	0.816126347	0.809144669
0.902575748	0.899025856	0.895190551	0.891078694	0.886699642	0.882063197	0.877179558	0.872059264	0.86671315	0.86115229	0.855387949	0.849431534	0.843294546	0.83698853	0.830525037	0.823915578	0.817171582	0.810304366	0.803325091	0.796244738
0.907069282	0.903514046	0.899675738	0.895563205	0.891185788	0.886553263	0.881675797	0.876563895	0.871228348	0.865680186	0.859930622	0.853991009	0.847872789	0.841587449	0.835146476	0.828561316	0.821843335	0.81500378	0.808053748	0.801004152
0.866111511	0.862850832	0.859335786	0.855573887	0.851573061	0.847341608	0.842888159	0.838221639	0.833351227	0.82828631	0.823036449	0.817611338	0.812020763	0.806274568	0.800382618	0.794354764	0.788200813	0.781930495	0.775553434	0.769079123
0.844979218	0.841870977	0.838521809	0.834938597	0.831128601	0.827099426	0.822858985	0.818415461	0.813777279	0.808953059	0.803951589	0.798781786	0.793452664	0.787973296	0.782352789	0.776600246	0.770724741	0.76473529	0.758640823	0.752450164
0.868538958	0.865274905	0.861757431	0.857994049	0.853992685	0.849761635	0.845309523	0.840645268	0.835778036	0.830717203	0.825472316	0.820053052	0.814469183	0.808730533	0.802846948	0.796828261	0.790684255	0.784424639	0.778059017	0.771596859
0.832887548	0.82993207	0.826750876	0.823350243	0.819736795	0.815917466	0.811899472	0.807690277	0.803297566	0.79872921	0.793993233	0.789097786	0.784051113	0.778861524	0.773537363	0.768086983	0.76251872	0.756840867	0.75106165	0.745189207
1.051175482	1.046670048	1.041786506	1.036538702	1.030941257	1.025009466	1.018759212	1.012206867	1.005369194	0.99826326	0.990906342	0.983315839	0.975509188	0.967503786	0.959316914	0.950965665	0.942466882	0.933837101	0.925092492	0.916248817
0.909755154	0.906216343	0.902395719	0.898302079	0.893944709	0.889333332	0.884478059	0.879389336	0.874077897	0.868554712	0.862830938	0.856917869	0.85082689	0.844569433	0.838156929	0.831600772	0.824912276	0.818102639	0.811182911	0.80416396
0.992651817	0.988555529	0.984120421	0.979358128	0.97428094	0.968901724	0.963233852	0.957291128	0.951087708	0.944638032	0.937956747	0.93105864	0.923958569	0.916671395	0.909211926	0.901594856	0.893834711	0.8859458	0.877942173	0.869837572
0.93150394	0.927822835	0.923846208	0.919583539	0.915044832	0.910240561	0.905181616	0.899879243	0.894344994	0.888590666	0.88262825	0.876469873	0.870127752	0.863614138	0.856941275	0.850121347	0.843166442	0.83608851	0.828899324	0.82161045
0.883799484	0.880543646	0.877041499	0.873300479	0.869328414	0.865133487	0.860724197	0.85610932	0.851297869	0.846299058	0.84112226	0.835776972	0.830272776	0.824619304	0.818826205	0.81290311	0.806859603	0.800705191	0.794449276	0.788101129
0.921273447	0.917720439	0.913893941	0.909802681	0.905455845	0.900863038	0.896034227	0.890979698	0.885710004	0.88023592	0.874568389	0.868718482	0.862697351	0.85651618	0.850186152	0.843718403	0.837123987	0.830413841	0.823598751	0.816689325
0.975403621	0.971438418	0.96715269	0.962557291	0.957663667	0.952483797	0.947030121	0.941315475	0.935353025	0.929156196	0.922738614	0.916114035	0.909296288	0.902299217	0.895136624	0.887822213	0.880369549	0.872792005	0.865102726	0.857314587
0.916497815	0.912937432	0.909099919	0.904994067	0.90062914	0.896014824	0.891161181	0.886078594	0.880777724	0.875269456	0.86956485	0.863675099	0.857611474	0.851385289	0.845007852	0.838490427	0.831844194	0.825080216	0.818209405	0.811242491
1.068616708	1.063989727	1.05896752	1.053564742	1.047796878	1.04168014	1.035231368	1.028467924	1.021407585	1.014068447	1.006468826	0.998627157	0.990561911	0.9822915	0.973834203	0.965208088	0.95643094	0.947520204	0.93849292	0.929365682
1.103162953	1.098295829	1.093003517	1.087302228	1.081209106	1.074742111	1.067919897	1.060761694	1.053287189	1.045516404	1.037469589	1.029167106	1.02062933	1.011876546	1.002928859	0.993806108	0.984527791	0.975112988	0.965580305	0.955947815
0.900601884	0.897204347	0.89354568	0.889633952	0.885477661	0.881085695	0.876467285	0.871631963	0.866589517	0.861349948	0.855923425	0.850320246	0.844550792	0.838625494	0.832554786	0.826349078	0.820018713	0.813573942	0.807024887	0.800381521
1.008842675	1.004624231	1.000045922	0.995120218	0.989860311	0.984280033	0.978393773	0.972216395	0.965763149	0.959049598	0.952091526	0.944904867	0.937505628	0.929909812	0.922133358	0.914192069	0.906101558	0.89787719	0.889534032	0.881086812
0.839475897	0.836691052	0.833701563	0.830513076	0.827131534	0.823563152	0.81981439	0.815891928	0.811802638	0.807553556	0.80315186	0.798604838	0.793919869	0.789104389	0.784165877	0.779111823	0.77394971	0.768686991	0.763331069	0.75788928
1.047742382	1.043263699	1.038394704	1.03314942	1.027542695	1.021590099	1.01530783	1.008712614	1.001821605	0.994652288	0.987222386	0.979549763	0.971652341	0.963548012	0.955254561	0.946789591	0.938170457	0.929414199	0.920537493	0.911556594
0.853168649	0.850053663	0.846700405	0.84311576	0.839306983	0.83528167	0.831047718	0.826613291	0.821986785	0.817176793	0.812192066	0.807041484	0.801734018	0.796278696	0.790684576	0.784960713	0.779116127	0.773159781	0.767100552	0.760947206
0.812340192	0.809711716	0.80688602	0.803868254	0.800663848	0.79727849	0.793718103	0.789988821	0.786096967	0.782049025	0.777851619	0.773511489	0.769035465	0.76443045	0.75970339	0.754861257	0.749911029	0.744859667	0.739714099	0.734481198
0.912376644	0.908818527	0.904980553	0.900871544	0.896500797	0.891878042	0.887013389	0.881917281	0.876600437	0.871073811	0.865348535	0.859435875	0.853347184	0.847093855	0.840687279	0.834138805	0.827459699	0.820661111	0.813754038	0.806749296
1.103049683	1.098179913	1.092882854	1.087174775	1.081072885	1.074595218	1.067760511	1.06058808	1.053097702	1.045309496	1.03724381	1.028921109	1.020361869	1.011586477	1.002615143	0.993467807	0.984164065	0.974723098	0.965163606	0.955503754
0.806233498	0.803771655	0.801127006	0.798304177	0.795308049	0.792143734	0.788816558	0.785332038	0.781695866	0.777913881	0.773992056	0.769936471	0.765753298	0.761448777	0.757029199	0.752500886	0.747870172	0.743143387	0.738326841	0.733426804
0.813576502	0.810942428	0.808111627	0.805089257	0.801880753	0.79849181	0.794928355	0.791196524	0.78730264	0.78325319	0.779054796	0.774714197	0.77023822	0.765633763	0.760907769	0.756067205	0.751119041	0.746070233	0.740927701	0.73569831
1.008865183	1.004648781	1.000072446	0.995148645	0.989890565	0.984312036	0.978427442	0.972251646	0.965799897	0.959087755	0.952131003	0.944945576	0.937547478	0.929952715	0.922177225	0.914236813	0.906147091	0.897923428	0.889580892	0.881134212
0.805633068	0.80317151	0.800526923	0.797703933	0.794707422	0.791542503	0.788214505	0.784728947	0.781091522	0.777308073	0.773384576	0.769327115	0.765141865	0.76083507	0.756413025	0.751882056	0.747248501	0.742518697	0.737698955	0.732795552
0.942407764	0.938711172	0.934729232	0.930471311	0.925947272	0.921167418	0.916142437	0.910883353	0.905401468	0.899708311	0.893815584	0.887735114	0.881478801	0.875058571	0.868486331	0.861773927	0.854933103	0.847975461	0.840912433	0.833755243

];

LsArrD = [3.532881065	3.546644136	3.561726633	3.578139485	3.595894755	3.615005695	3.635486806	3.657353893	3.680624138	3.705316174	3.731450158	3.759047865	3.78813277	3.818730158	3.850867221	3.884573181	3.919879408	3.956819558	3.995429715	4.035748545	4.077817466
3.533306105	3.546934179	3.561895125	3.578199919	3.595860697	3.614890803	3.635304855	3.657118801	3.680349993	3.705017258	3.731140981	3.758743192	3.787847662	3.81848	3.850667768	3.884440592	3.919830299	3.956871044	3.995599466	4.036054841	4.07827926
3.532916817	3.546610529	3.5616328	3.577994607	3.595708079	3.614786549	3.635244613	3.657098193	3.680364603	3.70506263	3.731212605	3.758836499	3.787958009	3.818602665	3.850797935	3.884573343	3.919960595	3.956993717	3.995709201	4.036146159	4.078346504
3.532815255	3.546580929	3.561655487	3.578051357	3.595782238	3.61486315	3.635310502	3.657142158	3.68037751	3.705037556	3.731144992	3.758724304	3.787801869	3.818406065	3.850567395	3.884318612	3.919694859	3.956733824	3.995475899	4.035964362	4.078245567
3.53303678	3.546812801	3.561889564	3.578280228	3.595999281	3.615062596	3.635487502	3.657292848	3.680499083	3.705128338	3.731204521	3.75875341	3.787802764	3.818382438	3.850524509	3.884263412	3.919636088	3.956682145	3.99544403	4.035967219	4.078300422
3.532651316	3.546453386	3.561562757	3.577991869	3.595754426	3.614865456	3.635341369	3.657200029	3.680460823	3.705144745	3.731274482	3.758874507	3.787971182	3.818592865	3.850770037	3.884535423	3.919924135	3.956973824	3.995724841	4.036220416	4.078506852
3.532717849	3.546566834	3.561705748	3.578148497	3.595910371	3.615008107	3.635459954	3.657285748	3.680506988	3.705146928	3.73123067	3.758785267	3.787839835	3.818425674	3.850576403	3.884328099	3.919719456	3.956791953	3.995590036	4.03616132	4.078556805
3.533217631	3.546725622	3.561615047	3.577891929	3.595563017	3.61463583	3.635118706	3.657020855	3.680352414	3.705124508	3.731349312	3.75904012	3.788211415	3.818878946	3.851059805	3.884772515	3.920037118	3.956875267	3.995310329	4.03536749	4.077073864
3.533505591	3.547167348	3.562154456	3.578477835	3.596149549	3.615182861	3.63559229	3.657393671	3.680604225	3.705242635	3.731329119	3.758885526	3.787935417	3.818504177	3.850619112	3.884309575	3.919607083	3.956545457	3.995160965	4.035492482	4.077581654
3.533529506	3.547051606	3.561935237	3.57818847	3.595820282	3.614840604	3.635260375	3.657091595	3.680347386	3.705042055	3.731191168	3.758811619	3.787921712	3.818541243	3.850691594	3.884395826	3.919678783	3.956567204	3.995089833	4.035277553	4.077163514
3.533452934	3.547086743	3.562054057	3.578365068	3.596031051	3.615064416	3.635478767	3.657288954	3.680511148	3.705162904	3.731263245	3.758832739	3.787893588	3.818469727	3.850586926	3.884272896	3.919557414	3.956472444	3.995052282	4.035333695	4.07735609
3.532288049	3.546224772	3.561453567	3.577986791	3.595838041	3.615022205	3.635555526	3.657455665	3.680741774	3.705434578	3.731556453	3.75913152	3.788185748	3.818747057	3.850845434	3.88451306	3.919784445	3.956696575	3.995289068	4.035604347	4.077687828
3.532519988	3.546394336	3.561577245	3.578079603	3.595913417	3.615091864	3.635629349	3.657541564	3.680845557	3.7055598	3.731704272	3.759300539	3.788371845	3.818943209	3.851041532	3.884695708	3.919936746	3.956797898	3.995314807	4.035525649	4.077471308
3.533494879	3.547130824	3.562111671	3.578446169	3.596144025	3.61521595	3.635673714	3.657530197	3.680799455	3.705496784	3.731638791	3.759243468	3.788330277	3.818920232	3.851035995	3.884701974	3.919944431	3.956791592	3.995273771	4.035423502	4.077275674
3.532223967	3.546052103	3.561196797	3.577668991	3.595480759	3.614645361	3.635177297	3.657092371	3.680407761	3.705142085	3.731315487	3.758949718	3.78806823	3.818696273	3.850861002	3.884591593	3.919919364	3.95687791	3.995503248	4.035833967	4.077911399
3.533292773	3.546717148	3.56153664	3.577756623	3.595383145	3.61442297	3.634883628	3.656773467	3.680101708	3.704878506	3.73111501	3.758823431	3.788017113	3.818710605	3.850919741	3.884661721	3.919955197	3.956820366	3.995279066	4.035354875	4.077073223
3.533774261	3.547104339	3.561849166	3.578012841	3.595600022	3.614615974	3.635066613	3.65695856	3.680299196	3.705096719	3.731360204	3.759099672	3.788326151	3.819051753	3.851289743	3.885054622	3.920362203	3.957229697	3.995675806	4.03572081	4.07738667
3.532681889	3.546530552	3.561681752	3.578147092	3.595939359	3.615072573	3.635562046	3.657424448	3.680677876	3.705341927	3.731437779	3.758988283	3.788018051	3.818553566	3.850623287	3.884257768	3.919489792	3.956354502	3.994889559	4.035135297	4.077134901
3.534106184	3.547469133	3.562215123	3.578351564	3.595886732	3.614829809	3.63519094	3.656981287	3.680213087	3.704899718	3.731055765	3.758697099	3.787840948	3.818505986	3.850712418	3.884482076	3.919838521	3.956807146	3.995415293	4.035692373	4.077669997
3.53231178	3.546350625	3.561654795	3.578238851	3.596118762	3.615311968	3.635837448	3.657715796	3.680969302	3.70562204	3.731699965	3.759231017	3.788245238	3.818774896	3.850854616	3.884521529	3.919815435	3.956778967	3.995457787	4.035900787	4.078160314
3.534137119	3.547422217	3.562110758	3.578208812	3.595723198	3.614661522	3.635032236	3.656844688	3.680109177	3.704837019	3.731040613	3.758733506	3.787930472	3.818647588	3.850902316	3.884713595	3.920101929	3.957089486	3.995700206	4.035959906	4.0778964
3.532427615	3.546302592	3.56147379	3.577954392	3.595758903	3.614903207	3.635404632	3.657282022	3.680555809	3.705248101	3.73138277	3.758985551	3.788084149	3.818708352	3.850890159	3.884663916	3.920066458	3.957137274	3.995918676	4.036455987	4.078797744
3.532670379	3.546568298	3.561738409	3.578196076	3.595958159	3.615043088	3.635470929	3.657263468	3.680444297	3.705038908	3.7310748	3.758581586	3.787591125	3.818137648	3.85025791	3.883991345	3.919380242	3.956469933	3.995308997	4.035949484	4.078447163
3.532855213	3.546624403	3.561715955	3.578140832	3.595911141	3.615040183	3.635542511	3.657433992	3.680731876	3.705454869	3.731623213	3.759258773	3.788385126	3.819027666	3.851213707	3.884972602	3.920335865	3.95733731	3.996013192	4.036402363	4.078546446
3.533513468	3.546826568	3.561558135	3.577712279	3.595293675	3.614307601	3.634759996	3.656657503	3.680007532	3.704818311	3.731098953	3.758859521	3.788111089	3.818865824	3.85113705	3.884939335	3.920288566	3.957202038	3.995698545	4.035798469	4.077523883
3.532633527	3.546567731	3.561765039	3.578241541	3.59601488	3.615104328	3.635530856	3.657317218	3.680488045	3.705069941	3.731091594	3.758583893	3.787580059	3.818115785	3.850229389	3.883961982	3.91935765	3.956463655	3.99533065	4.036012916	4.078568624
3.532422138	3.546363523	3.561576406	3.57807618	3.595879739	3.615005541	3.635473682	3.657305978	3.68052605	3.705159417	3.731233605	3.758778257	3.787825256	3.818408863	3.850565859	3.884335704	3.919760715	3.956886248	3.995760909	4.036436777	4.078969645
3.534035482	3.547390733	3.562128985	3.578257645	3.595784982	3.614720176	3.635073364	3.656855703	3.680079423	3.704757896	3.7309057	3.758538699	3.787674113	3.818330608	3.850528379	3.884289251	3.919636772	3.956596327	3.995195245	4.035462927	4.077430969
3.532681204	3.546606245	3.56179492	3.578263319	3.596029085	3.615111492	3.635531513	3.657311908	3.680477311	3.705054335	3.731071675	3.758560232	3.787553237	3.818086397	3.850198045	3.883929311	3.9193243	3.956430297	3.995297978	4.035981654	4.078539526
3.532409342	3.546278307	3.561460062	3.577964743	3.595803541	3.614988749	3.635533816	3.657453409	3.680763475	3.705481307	3.731625623	3.759216642	3.788276172	3.818827701	3.850896497	3.884509711	3.919696498	3.956488132	3.994918144	4.03502246	4.076839556
];

% 5 10 18 19 22 30 35 40

% 1 2 3 9 15 18 24
ii = 30;

Fsi = FsArrD(ii, :); 
Lsi = LsArrD(ii, :);

X = [];



for i=1:1:21
        
    T0 = hour(i) + minutes(i)/60;
    X(i) = T0;
end


plot(X, Ls, '*');
hold on;

plot(X, Lsi, 'red');
hold on;
plot(X(:, 1:20), Fs, '.');
hold on;
plot(X(:, 1:20), Fsi, 'blue');

xlabel('时间');
ylabel('影子长度');
title('306天 1.9M 36S 76E');