逻辑回归拟合的MATLAB代码（数据完整，）

clc clear close all load maydata1.mat n=size(y,1); X1=ones(size(x,1),1); % 构建常数项系数 X=[X1, x]; Y=log(y./(1-y)); b=regress(Y,X); XE =x; %-------------------------------------------------------------------------- %% 模型验证和应用 for i=1:size(XE,1) Pai0=exp(b(1)+b(2)*XE(i,1)+b(3)*XE(i,2)+b(4)*XE(i,3)+... b(5)*XE(i,4)+b(6)*XE(i,5)+b(7)*XE(i,6)+... b(8)*XE(i,7)+b(9)*XE(i,8)+b(10)*XE(i,9))/(1+... exp(b(1)+b(2)*XE(i,1)+b(3)*XE(i,2)+b(4)*XE(i,3)+... b(5)*XE(i,4)+b(6)*XE(i,5)+b(7)*XE(i,6)+... b(8)*XE(i,7)+b(9)*XE(i,8)+b(10)*XE(i,9))); % if Pai0<=0.5 % P(i)=0; % else % P(i)=1; % end P(i)=Pai0; end %-------------------------------------------------------------------------- %% 显示结果 disp(['回归系数:' num2str(b')]); figure plot(7*y,'bo') hold on plot(7*P,'r^') yz1 = round(y*7); yz2 = round(P*7)'; zhunquelv = sum(yz1==yz2)/length(yz1) [MSE, RMSE ,MBE, MAE ] =MSE_RMSE_MBE_MAE(7*y,7*P'); R_2 = R_2(7*y,7*P'); figure boxplot([ones(1,5);[MSE, RMSE ,MBE, MAE R_2]]) set(gca,'XTick',1:5) set(gca,'XTickLabel',{'均方误差','均方根误差','平均绝对误差','平均偏差','相关系数'}) set(gca,'fontsize',12) legend('进化策略优化小波软阈值去噪')