数据下载地址在这里
Regularized Linear Regression
源码
clear;
clc;
%加载数据
x = load('ex5Linx.dat');
y = load('ex5Liny.dat');
%显示原始数据
plot(x,y,'o','MarkerEdgeColor','b','MarkerFaceColor','r','MarkerSize',10);
%将特征值变成训练样本矩阵
x = [ones(length(x),1),x,x.^2,x.^3,x.^4,x.^5];
[m,n]=size(x);
n=n-1;
%计算参数sidta,并且绘制出拟合曲线
rm=diag([0;ones(n,1)]);%lamda后面的矩阵
lamda=[0 1 10]';
colortype={
'g','b','r'};
theta=zeros(n+1,3);
xrange=linspace(min(x(:,2)),max(x(:,2)))';
hold on;
for i = 1:3
theta(:,i)=inv(x'*x+lamda(i).*rm)*x'*y;%计算参数sida
norm_sida=norm(theta)
yrange=[ones(size(xrange)) xrange xrange.^2 xrange.^3,...