clear all;
close all;
clc;

randn('seed',0);
mu1=[0 0];
S1=[0.3 0;0 0.35];
cls1_data=mvnrnd(mu1,S1,1000);
plot(cls1_data(:,1),cls1_data(:,2),'+');
hold on;

mu2=[4 0];
S2=[1.2 0;0 1.85];
cls2_data=mvnrnd(mu2,S2,1000);
plot(cls2_data(:,1),cls2_data(:,2),'r+');
axis([-8 8 -8 8]);

for i=-4:0.1:8
for j=-6:0.1:6
d1=([i,j]-mu1)*inv(S1)*([i,j]-mu1)';
d2=([i,j]-mu2)*inv(S2)*([i,j]-mu2)';
d=d1-d2;
if d<0.1
plot(i,j);
end

end
end
grid on;

figure;
mu1=[0 0];
S1=[0.1 0;0 0.75];
cls1_data=mvnrnd(mu1,S1,1000);
plot(cls1_data(:,1),cls1_data(:,2),'+');
hold on;

mu2=[3.2 0];
S2=[0.75 0;0 0.1];
cls2_data=mvnrnd(mu2,S2,1000);
plot(cls2_data(:,1),cls2_data(:,2),'r+');
axis([-8 8 -8 8]);

for i=-8:0.1:8
for j=-8:0.1:8
d1=([i,j]-mu1)*inv(S1)*([i,j]-mu1)';
d2=([i,j]-mu2)*inv(S2)*([i,j]-mu2)';
d=d1-d2;
if d<0.1
plot(i,j);
end

end
end
grid on;

figure;
mu1=[0 3];
S1=[0.3 0;0 0.35];
cls1_data=mvnrnd(mu1,S1,1000);
plot(cls1_data(:,1),cls1_data(:,2),'+');
hold on;

mu2=[4 0];
S2=[0.3 0;0 0.35];
cls2_data=mvnrnd(mu2,S2,1000);
plot(cls2_data(:,1),cls2_data(:,2),'r+');

axis([-8 8 -8 8]);

for i=-8:0.1:8
for j=-8:0.1:8
d1=([i,j]-mu1)*inv(S1)*([i,j]-mu1)';
d2=([i,j]-mu2)*inv(S2)*([i,j]-mu2)';
d=d1-d2;
if d<0.1
plot(i,j);
end

end
end
grid on;

figure;
mu1=[0 -3];
S1=[0.5 1;1 2.5];
cls1_data=mvnrnd(mu1,S1,1000);
plot(cls1_data(:,1),cls1_data(:,2),'+');
hold on;

mu2=[4 0];
S2=[0.5 1;1 2.5];
cls2_data=mvnrnd(mu2,S2,1000);
plot(cls2_data(:,1),cls2_data(:,2),'r+');

axis([-8 8 -8 8]);

for i=-8:0.1:8
for j=-8:0.1:8
d1=([i,j]-mu1)*inv(S1)*([i,j]-mu1)';
d2=([i,j]-mu2)*inv(S2)*([i,j]-mu2)';
d=d1-d2;
if d<0.1
plot(i,j);
end

end
end
grid on;

http://www.cnblogs.com/tiandsp/category/348031.html

最新文章

1. node.js基础 1之 HTTP流程实例
2. invalid byte 1 of 1-byte UTF-8 sequence
3. VBS学习：流程控制语句判断结构
4. 重构14-Break&#160;Responsibilities
5. C内存分配函数
6. bzoj1005
7. VB 中ListView 某一列的颜色添加不上去的解决方法
8. 编译SASS
9. js自定义的简易滚动条
10. getScrollX()理解
11. Interface Comparable&lt;T&gt;
12. canvas动态修改宽高问题
13. ES 7 async／await Promise
14. CISCO 过载NAT配置（小型网络)
15. JSP页面用&lt;a&gt;标签访问 Action 出错
16. Android-仿“抖音”的评论列表的UI和效果
17. Visual Studio中设置Nuget程序包源
18. Asp.Net导出文件名中文乱码
19. XAF-如何实现自定义权限系统用户对象
20. Python中协程Event()函数

热门文章

1. [CF1454] Codeforces Round #686 (Div. 3) solution
2. Linux 学习笔记04丨Linux的用户和用户组管理
3. VUE：组件总结
4. angular中datetime-local属性使用ng-model报错
5. Spring Boot + RabbitMQ 使用示例
6. Mac用brew更新完python2.7后无法找到虚拟环境
7. moviepy音视频开发：使用credits1给视频加片头片尾字幕
8. python-字典dict、去除重复set
9. postman学习网址
10. 团队作业part4--项目冲刺