Departamento de produção e Sistemas – Métodos numéricos


Método dos mínimos quadrados para a função do primeiro grau em Matlab

 

>> x=[46.000000;

          20.000000;

          52.000000;

          30.000000;

          57.000000;

          25.000000;

          28.000000;

          36.000000;

          57.000000;

          44.000000;

          31.000000;

          31.000000;

          52.000000;

          23.000000;

          60.000000;

          48.000000;

          34.000000;

          51.000000;

          50.000000;

          34.000000;

          46.000000;

          23.000000;

          37.000000;

          40.000000;

          30.000000]

 

x =

 

    46

    20

    52

    30

    57

    25

    28

    36

    57

    44

    31

    31

    52

    23

    60

    48

    34

    51

    50

    34

    46

    23

    37

    40

    30

 

>> y=[354.000000;

         190.000000;

         405.000000;

         263.000000;

         451.000000;

         302.000000;

         288.000000;

         385.000000;

         402.000000;

         365.000000;

         209.000000;

         290.000000;

         356.000000;

         254.000000;

         395.000000;

         434.000000;

         220.000000;

         374.000000;

         308.000000;

         220.000000;

         311.000000;

         181.000000;

         274.000000;

         303.000000;

         244.000000]

 

y =

 

   354

   190

   405

   263

   451

   302

   288

   385

   402

   365

   209

   290

   356

   254

   395

   434

   220

   374

   308

   220

   311

   181

   274

   303

   244

 

>> p= polyfit (x,y,1)

 

p =

 

    5.4051   98.1605

 

>> f= polyval (p,x)

 

f =

 

  346.7934

  206.2618

  379.2238

  260.3124

  406.2491

  233.2871

  249.5023

  292.7428

  406.2491

  335.9833

  265.7175

  265.7175

  379.2238

  222.4769

  422.4643

  357.6036

  281.9327

  373.8187

  368.4137

  281.9327

  346.7934

  222.4769

  298.1478

  314.3630

  260.3124

 

>> r=sum((y-f).^2)

 

r =

 

  4.5797e+004

 

 

 

 

Anterior

 

Voltar ao inicio

Hosted by www.Geocities.ws

1