x = load(‘ex3x.dat‘);y = load(‘ex3y.dat‘);x = [ones(size(x,1),1) x];%每一行是一个样本，在这里每个样本增加一维1，原因在前面课说了（讲wx+b变成w‘x齐次的）meanx = mean(x);%求均值           接下来四行是让样本的每一维度（除第一维1外）的值标准化。sigmax = std(x);%求标准偏差        但是前面不是说线性的不用进行feature scale吗（第一课讲的）x(:,2) = (x(:,2)-meanx(2))./sigmax(2);x(:,3) = (x(:,3)-meanx(3))./sigmax(3); figureitera_num = 100; %尝试的迭代次数sample_num = size(x,1); %训练样本的个数alpha = [0.01, 0.03, 0.1, 0.3, 1, 1.3];%因为差不多是选取每个3倍的学习率来测试，所以直接枚举出来plotstyle = {‘b‘, ‘r‘, ‘g‘, ‘k‘, ‘b--‘, ‘r--‘};%建了一个包，每一个值代表画出的曲线样式不同，b是blue蓝色，%r是red ，g是green..b--是blue颜色--代表的是虚线，而前面那些不加的是实现。theta_grad_descent = zeros(size(x(1,:)));for alpha_i = 1:length(alpha) %alpha_i是1,2，...6，表示的是学习速率向量和曲线格式向量的坐标：alpha(alpha_i)，plotstyle(alpha_i)    theta = zeros(size(x,2),1); %theta是cost function的参数，初始值赋值为0向量（3*1的向量，x有几维theta就是几维的参数向量）    Jtheta = zeros(itera_num, 1);%Jthete是个100*1的向量，第n个元素代表第n次迭代cost function的值（预测与真实y的总均方误差）    for i = 1:itera_num %计算出某个学习速率alpha下迭代itera_num次数后的参数               Jtheta(i) = (1/(2*sample_num)).*(x*theta-y)‘*(x*theta-y);%Jtheta是个100*1的列向量。(x*theta-y)‘*(x*theta-y)代表的就是        %cost function 公式的那个平方，因为在向量水平上平方没有直接平方，所以就是这种转置后内积的形式。并且得到的是        %一个标量，所以再与前面的系数相乘可以直接用*，而不用.*  还有一点是前面的系数 我还是不明白为什么        %是(1/(2*sample_num))        grad = (1/sample_num).*x‘*(x*theta-y);        theta = theta - alpha(alpha_i).*grad;    end    plot(0:49, Jtheta(1:50),char(plotstyle(alpha_i)),‘LineWidth‘, 2)%此处一定要通过char函数来转换因为包用（）索引后得到的还是包cell，    %所以才要用char函数转换，也可以用{}索引，这样就不用转换了。    %一个学习速率对应的图像画出来以后再画出下一个学习速率对应的图像。        hold on    if(1 == alpha(alpha_i)) %通过实验发现alpha为1时效果最好，则此时的迭代后的theta值为所求的值        theta_grad_descent = theta    endendlegend(‘0.01‘,‘0.03‘,‘0.1‘,‘0.3‘,‘1‘,‘1.3‘);xlabel(‘Number of iterations‘)ylabel(‘Cost function‘)%下面是预测公式price_grad_descend = theta_grad_descent‘*[1 (1650-meanx(2))/sigmax(2) (3-meanx(3)/sigmax(3))]‘