y=1986:2015;
>> a=[512.73 521.23 534.69 544.45 550.05 553.56 555.38 555.42 559.44 562.82 567.53 571.91 576.88 582.47 583.19 587.09 595.46 601.76 610.38 620.92 631.00 637.36 641.74 646.84 650.12 656.62 660.62 662.81 671.41 680.36];
>> t1=[512.73 521.23 534.69 544.45 550.05 553.56 555.38 555.42 559.44 562.82 567.53 571.91 576.88 582.47 583.19 587.09 595.46 601.76 610.38 620.92 631.00 637.36 641.74 646.84 650.12 656.62 660.62 662.81 671.41 680.36 0];
>> a=a';a=a(:);
>> Rt=tiedrank(a);
>> n=length(a);t=1:n;
>> Qs=1-6/(n*(n^2-1))*sum((t-Rt).^2)
错误使用 -
矩阵维度必须一致。
我将Qs=1-6/(n*(n^2-1))*sum((t-Rt).^2)这一步拆开来写,发现是sum((t-Rt).^2)这的问题,但本人属于初次接触MATLAB,不知道怎么修改,我在工作列表中查到t是1x30 double
Rt是30x1 double。