interpolarea_functiilor.m

Function interpolation

clear all
close all
clc

P=0:0.1:1;
T=0.2+0.8*(P+0.7*sin(2*pi*P));
Ptest=0:0.01:1;
Ttest=0.2+0.8*(Ptest+0.7*sin(2*pi*Ptest));
net1=newff(minmax(P),[15 1], {'tansig' 'logsig'}, 'traingdm');
net1.trainParam.lr=2;
net1.trainParam.mc=0.5;
net1.trainParam.min_grad=1e-10;
net1.trainParam.show=1;
net1.trainParam.epochs=300;
net1.trainParam.goal=0.001;

[net1,tr1]=train(net1,P,T);
y1test=sim(net1,Ptest);
y1antr=sim(net1,P);

er1antr=(sum(sum(abs(T-y1antr))))/length(T);
er1test=(sum(sum(abs(Ttest-y1test))))/length(Ttest);
er=[er1antr er1test];

clf;
plot(Ptest,Ttest,'+r');
title('Aproximarea functiei f(x)');
hold

plot(P,T,'.')
plot(Ptest,y1test,':');
xlabel('x');
ylabel('Antrenament . Test + Iesire MLP --')
hold off
figure
bar(er*100);
ylabel('Eroare aproximare [%]');
xlabel('Antrnament                Test')