%?程序P7-1：?logistic方法Matlab實(shí)現(xiàn)程序
%?《MATLAB數(shù)學(xué)建模方法與實(shí)踐》（《MATLAB在數(shù)學(xué)建模中的應(yīng)用》升級(jí)版），北航出版社，卓金武、王鴻鈞編著.?
%%?數(shù)據(jù)準(zhǔn)備
clear?all
clc
X0=xlsread（‘logistic_ex1.xlsx‘?‘A2:C21‘）;?%?回歸數(shù)據(jù)X值
XE=xlsread（‘logistic_ex1.xlsx‘?‘A2:C26‘）;?%?驗(yàn)證與預(yù)測(cè)數(shù)據(jù)
Y0=xlsread（‘logistic_ex1.xlsx‘?‘D2:D21‘）;?%?回歸數(shù)據(jù)P值
%--------------------------------------------------------------------------
%%?數(shù)據(jù)轉(zhuǎn)化和參數(shù)回歸
n=size（Y01）;
for?i=1:n
????if?Y0（i）==0
????????Y1（i1）=0.25;
????else
????????Y1（i1）=0.75;
????end
end
X1=ones（size（X01）1）;?%?構(gòu)建常數(shù)項(xiàng)系數(shù)
X=[X1?X0];
Y=log（Y1./（1-Y1））;
b=regress（YX）;
%--------------------------------------------------------------------------
%%?模型驗(yàn)證和應(yīng)用
for?i=1:size（XE1）
????Pai0=exp（b（1）+b（2）*XE（i1）+b（3）*XE（i2）+b（4）*XE（i3））/（1+exp（b（1）+b（2）*XE（i1）+b（3）*XE（i2）+b（4）*XE（i3）））;
????if?Pai0<=0.5
???????P（i）=0;
????else?
???????P（i）=1;
????end
end
%--------------------------------------------------------------------------
%%?顯示結(jié)果
disp（[‘回歸系數(shù):‘?num2str（b‘）]）;
disp（[‘評(píng)價(jià)結(jié)果:‘?num2str（P）]）;

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件????????1118??2018-06-18?14:58??logical\P7_1_logisctic_ex1.m
?????文件???????13150??2015-02-26?12:16??logical\logistic_ex1.xlsx
?????目錄???????????0??2018-09-18?22:01??logical\