机器学习(周志华) 参考答案 第三章 线性模型 3.5

机器学习(周志华西瓜书) 参考答案 总目录

• http://blog.csdn.net/icefire_tyh/article/details/52064910

机器学习(周志华) 参考答案 第三章 线性模型

• http://blog.csdn.net/icefire_tyh/article/details/52069025

---

5.编程实现线性判别分析，并给出西瓜数据集3.0α上的结果。

    线性判别法没啥好说的，由于西瓜数据集数据分布不太友好，所以LDA结果也不是那么

西瓜数据集3.0

LDA在线性不可分的数据集下确实确实结果不理想

%线性判别分析（LDA） 西瓜数据集3.0ɑx = xlsread('C:\Users\icefire\Desktop\ml\西瓜3.0.xlsx', 'sheet1', 'A1:Q2');y = xlsread('C:\Users\icefire\Desktop\ml\西瓜3.0.xlsx', 'sheet1', 'A4:Q4');%更改y的值用来适合程序 1好瓜 2坏瓜y=2-y;u= zeros(2);%计算均值for i=1:17    u(:,y(i)) = u(:,y(i)) + x(:,i);endu(:,1) = u(:,1) / 8;u(:,2) = u(:,2) / 9;%计算两类协方差矩阵和Sw=zeros(2);for i=1:17    Sw=Sw+(x(:,i)-u(:,y(i)))*(x(:,i)-u(:,y(i))).';end%使用奇异值分解计算w [U,S,V] = svd (Sw);w=V/S*U.'*(u(:,1) - u(:,2));%画出散点图以及计算出的直线%逐点画  分别表示是否好瓜for i=1:17         if y(i)==1       plot(x(1,i),x(2,i),'+r');       hold on;    else if y(i)==2          plot(x(1,i),x(2,i),'og');              hold on;        end    endend%计算出直线边界点 并绘制直线ply=-(0.1*w(1)-0.01)/w(2);pry=-(0.9*w(1)-0.01)/w(2);line([0.1 0.9],[ply pry]);xlabel('密度');ylabel('含糖率');title('线性判别分析（LDA）');