多项式逻辑回归代码美赛思路即将发布科研交流

1、logisticy=1/(1+exp(-在matlab 中可以画出其x =x = -y= 1./(exp(-title(logisticfunction); 以上是一维的情况。对于多维变量，可以定义一个超平面 代入原来的变量x中，得到：对于任意变量x，可以代入上式计算出 y 值并与 0.5 比较进行分类，其中sgn(x) 为符号函数。为了演示logistic 函数是版怎样用于分类的，定义平方和（或L2norm)代价函数为：% generaterandom shift=% generaterandom shift=n= 2;%2 dimN = x=randn(n,N/2)-shift,randn

2、(n,N/2)*2+shift; y = zeros(N/2,1);ones(N/2,1);%showthedata hold on; title(2dtraining codesegmentfunction w,b,cost=function w,b,cost=fprintf(trainingstarted.n); n = size(x,1);N= w=ones(n,1)/n; b = 1;cost= count= 1;while %find partial_w=zeros(n,1); partial_b = 0;for a= partial_w=partial_w%findstep old

3、_cost=logistic_cost(x,y,w,b); step = 1;while step 1e-w1= w - partial_b=partial_bb1= b - new_costnew_cost=if new_coststep=step* if step= 1e-fprintf(finishedserachingthestepsizen); w=w1; b=cost= if new_costfprintf(trainingstopedafter%diterates,notconvergedto desired precision!,count);count=count+ 调用上面

4、logistic_train 函数就能得到w,b。 如前述，训练方法为梯度下降法。%training.%training.w,b,cost = logistic_train(x,y,1e-6,100);%show convergence curevedisp(training ended);title(convergencecurve);cost训练结束后，得到的w30.199430.4356b2.7303. （由于样本是随机人工生成的，每次得到的w,b 都会有些不同。）最后，有了模型（即logistic function)和模型参数 w,b，% plot% plot thetraininghold on;%visualizetheclassificationarea hold on;for x = -shift*5:0.5:shift*5for y