Dodano użytkownika Krzysztof K. Aby usunąć, naciśnij Backspace.

Przejdź do treści
Korzystanie z usługi Gmail z czytnikami ekranu
Przenieś do Odebranych

Więcej
1 z 1 222

estymacja
Beata Mach <rasti.budyn@gmail.com>

Załączniki16 sty (2 dni temu)

do Paulina
have a fun

2 załączniki
Beata Mach <rasti.budyn@gmail.com>

Załączniki14:02 (2 minuty temu)

do Natalia, Krzysztof

---------- Forwarded message ----------
From: Beata Mach <rasti.budyn@gmail.com>
Date: 2018-01-16 20:23 GMT+01:00
Subject: estymacja
To: Paulina Krzyżak <paulinakrzyzak94@gmail.com>


have a fun


2 załączniki

Odpowiedz, odpowiedz wszystkim lub przekaż dalej
Używasz 12,97 GB (86%) z 15 GB
Zarządzaj
Warunki - Prywatność
Ostatnia aktywność konta: 20 godzin temu
Szczegóły



function [ metryki, mHat, yHat ] = data2poli( x, y ,st )
%UNTITLED3 Summary of this function goes here
%   Detailed explanation goes here


G = [ones(size(x),st)];
for i=2:st
    G(:,i)=x.^(i-1);
end
D = y;
%m hat
mHat = inv(G'*G) * G'*y;

% MDOEL REGRESJI
% yhat
yHat = G*mHat;
% y_est = mHat(1) * G(1) ;
% for i=2:st
%     y_est=+mHat(i) * G(i).*(i-1);
% end
% REsztkowa suma kwadratów odchyleń
% RSS = sum((D - y_est).^2);
RSS = sum((y - yHat).^2);

% Współczynnik determinacji modelu

% rss = sum( (y-y_est).^2);
% r2=1-(rss/tss);

ysr=mean(y);
TSS = sum((y-ysr).^2);
R2=1-RSS/TSS;
metryki = [RSS R2];
end

data2poli.m
Wyświetlanie data2poli.m.




=================



close all ; clear all; clc;

k2=load('E:\SemestrII\Modelowanie\pomiar_all.txt');
x=k2(:,1);
y=k2(:,2);
% wykres zależnosci
plot(x,y,'*')
st=15;
[ metryki1,mHat1,yHat1 ] = data2poli( x, y, st );

hold on;
plot(x,yHat1,'or')


%% zad3
close all; clear all; clc;
delete('wyniki.txt');
k2=load('E:\SemestrII\Modelowanie\pomiar_all.txt');
x=k2(:,1);
y=k2(:,2);
% wykres zależnosci



for k=1:15
    [ metryki,mHat,yHat ] = data2poli( x, y, k );
    save('wyniki.txt','metryki','-ASCII','-append')
end
subplot(2,2,1)
plot(y,'*')
hold on;
[ metryki,mHat,yHat ] = data2poli( x, y, 1 );
plot(yHat,'r')

subplot(2,2,2)
plot(y,'*')
hold on;
[ metryki,mHat,yHat ] = data2poli( x, y,2 );
plot(yHat,'.r')

subplot(2,2,3)
plot(y,'*')
hold on;
[ metryki,mHat,yHat ] = data2poli( x, y, 15 );
plot(yHat,'.r')
w=load('wyniki.txt');
subplot(2,2,4)

x1 = 0:1:12;
rss=w(:,1);
rss1=rss(3:15)
r2=w(:,2);
r21=r2(3:15)
for i=2:length(w)
    plotyy(x1,rss1,x1,r21)
    hold on;
end

estymacja_wielomianowa.m