4--1.m

本程序是利用神经网络来计算五个城市的最短距离

clear all
% 初始化
x=-1:0.2:1;
y=0.5+0.4*sin(2*pi*x);
for i=1:length(x)
    noise=sqrt(0.05)*randn(1);
    yn(i)=y(i)+noise;
end

% M阶多项试拟合
M=10;
n=length(x);
B=zeros(1:M+1);
F=zeros(n,M+1);
for k=1:M+1
    F(:,k)=x'.^(k-1);
end
A=F'*F;
B=F'*yn';
C=A\B;


% 代入测试集进行测试
x=-1:0.2:1;
yc=0;
for k=1:M+1
    yc=yc+C(k).*(x'.^(k-1));
   
end

%显示结果
plot(x,y,'b-o',x,yn,'g-o',x,yc,'r--*');
legend({'待拟合函数','加了噪声的函数','逼近结果'});