?? 新建 文本文檔2.txt
字號:
load d:\exercise\iris\train.txt;
load d:\exercise\iris\test.txt;
a=train(:,1:4);
a0=(a-ones(75,1)*mean(a))./(ones(75,1)*std(a));
b=test(:,1:4);
a1=diag((a0-ones(75,1)*mean(a0(1:25,:)))*inv(cov(a0(1:25,:)))*(a0-ones(75,1)*mean(a0(1:25,:)))');
a1(:,2:2)=diag((a0-ones(75,1)*mean(a0(26:50,:)))*inv(cov(a0(26:50,:)))*(a0-ones(75,1)*mean(a0(26:50,:)))');
a1(:,3:3)=diag((a0-ones(75,1)*mean(a0(51:75,:)))*inv(cov(a0(51:75,:)))*(a0-ones(75,1)*mean(a0(51:75,:)))');
a2=exp(-a1/2);
a2(:,4:4)=train(:,8:8);
a2(:,1:1)=a2(:,1:1)/((2*pi)*(2*pi)*sqrt(det(cov(a(1:25,:)))));
a2(:,2:2)=a2(:,2:2)/((2*pi)*(2*pi)*sqrt(det(cov(a(26:50,:)))));
a2(:,3:3)=a2(:,3:3)/((2*pi)*(2*pi)*sqrt(det(cov(a(51:75,:)))));
fid=fopen('d:\exercise\iris\resulttr.txt','w');
fprintf(fid,'%g %g %g %g %g\n',a2');
b0=(b-ones(75,1)*mean(a))./(ones(75,1)*std(a));
b1=diag((b0-ones(75,1)*mean(a0(1:25,:)))*inv(cov(a0(1:25,:)))*(b0-ones(75,1)*mean(a0(1:25,:)))');
b1(:,2:2)=diag((b0-ones(75,1)*mean(a0(26:50,:)))*inv(cov(a0(26:50,:)))*(b0-ones(75,1)*mean(a0(26:50,:)))');
b1(:,3:3)=diag((b0-ones(75,1)*mean(a0(51:75,:)))*inv(cov(a0(51:75,:)))*(b0-ones(75,1)*mean(a0(51:75,:)))');
b2=exp(-b1/2);
b2(:,1:1)=b2(:,1:1)/((2*pi)*(2*pi)*sqrt(det(cov(a(1:25,:)))));
b2(:,2:2)=b2(:,2:2)/((2*pi)*(2*pi)*sqrt(det(cov(a(26:50,:)))));
b2(:,3:3)=b2(:,3:3)/((2*pi)*(2*pi)*sqrt(det(cov(a(51:75,:)))));
b2(:,4:4)=test(:,8:8);
fid=fopen('d:\exercise\iris\resultte.txt','w');
fprintf(fid,'%g %g %g %g %g\n',b2');?? 快捷鍵說明
復(fù)制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -