Errores (problema 3 hoja 4)

3_4.m

>> A = [1 0 0; 0 1 0; 0 0 1; 1 -1 0; 1 0 -1; 0 1 -1;]

A =

     1     0     0
     0     1     0
     0     0     1
     1    -1     0
     1     0    -1
     0     1    -1

>> longs = [2.49 1.74 3.25 4.65 2.63 4.40]

longs =

    2.4900    1.7400    3.2500    4.6500    2.6300    4.4000

>> des = [-809.4 -177.4 -5763.3 634.4 -4953.5 -5584.8]

des =

  1.0e+003 *

   -0.8094   -0.1774   -5.7633    0.6344   -4.9535   -5.5848

>> sum(longs)

ans =

   19.1600

>> longs[1]/sum(longs)
??? longs[1]/sum(longs)
         |
Error: Unbalanced or unexpected parenthesis or bracket.
 
>> longs(1)/sum(longs)

ans =

    0.1300

>> longs*des
??? Error using ==> mtimes
Inner matrix dimensions must agree.
 
>> longs'*des

ans =

  1.0e+004 *

   -0.2015   -0.0442   -1.4351    0.1580   -1.2334   -1.3906
   -0.1408   -0.0309   -1.0028    0.1104   -0.8619   -0.9718
   -0.2631   -0.0577   -1.8731    0.2062   -1.6099   -1.8151
   -0.3764   -0.0825   -2.6799    0.2950   -2.3034   -2.5969
   -0.2129   -0.0467   -1.5157    0.1668   -1.3028   -1.4688
   -0.3561   -0.0781   -2.5359    0.2791   -2.1795   -2.4573

>> 
>> longs*des'

ans =

 -5.5706e+004

>> (longs*des')/sum(longs)

ans =

 -2.9074e+003

>> lp = (longs*des')/sum(longs)

lp =

 -2.9074e+003

>> sigma_c = 1/le
??? Error using ==> le
Not enough input arguments.
 
>> length(10)

ans =

     1

>> length([10 2])

ans =

     2

>> m = length(longs)

m =

     6

>> %Observaciones independientes -> Matriz diagonal
>> %Varianzas de los desniveles observados -> proporcionales a las longitudes de los itinerarios
%var(tria h_j) = sigma_c^2*s_ij
>> % s_ij no da unidades (var(triag h_ij)) -> puesto que se utilizan como un factor de ponderacion
>> %lambda_iv*W = 6*6 diag
>> %sigma_y = sigma_0^2*p^-1 (peso unidad)
>> %Asigno al desnivel triag h_ij un peso: 10/s_ij
>> % ej (s_12 = 2.49 -> w_11 ~ 4.0, s_13... -> w_22 = ..., ...)
>> %W = diag(4.0, ...)
>> %x^p (A^t*W*A)^-1*A^t*w*y
>> 10/longs
??? Error using ==> mldivide
Matrix dimensions must agree.
 
>> longs*10

ans =

   24.9000   17.4000   32.5000   46.5000   26.3000   44.0000

>> longs*1/10

ans =

    0.2490    0.1740    0.3250    0.4650    0.2630    0.4400

>> for i=1:m,
 s(i) = 10/longs(i)
end

s =

    4.0161


s =

    4.0161    5.7471


s =

    4.0161    5.7471    3.0769


s =

    4.0161    5.7471    3.0769    2.1505


s =

    4.0161    5.7471    3.0769    2.1505    3.8023


s =

    4.0161    5.7471    3.0769    2.1505    3.8023    2.2727

>> W = diag(s)

W =

    4.0161         0         0         0         0         0
         0    5.7471         0         0         0         0
         0         0    3.0769         0         0         0
         0         0         0    2.1505         0         0
         0         0         0         0    3.8023         0
         0         0         0         0         0    2.2727

>> var_des = var(des)

var_des =

  8.7612e+006

>> sigma_c = var_des/10

sigma_c =

  8.7612e+005

>> sigma_c_0 = 10*var_des

sigma_c_0 =

  8.7612e+007

>> sigma_c_0 = 10*sigma_c

sigma_c_0 =

  8.7612e+006

>> Q_y = sigma_c_0*W^-1

Q_y =

  1.0e+006 *

    2.1815         0         0         0         0         0
         0    1.5244         0         0         0         0
         0         0    2.8474         0         0         0
         0         0         0    4.0739         0         0
         0         0         0         0    2.3042         0
         0         0         0         0         0    3.8549

>> x = (A'*W*A)^-1*A'*W

x =

    0.5433    0.2505    0.2062    0.1972    0.2595   -0.0533
    0.1750    0.6790    0.1460   -0.1603   -0.0147    0.1607
    0.2692    0.2727    0.4581    0.0421   -0.3113   -0.2306

>> x = (A'*W*A)^-1*A'*W*des
??? Error using ==> mtimes
Inner matrix dimensions must agree.
 
>> x = (A'*W*A)^-1*A'*W