声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1188|回复: 9

[编程技巧] 请教一个几组数据的方差排序的问题?

[复制链接]
发表于 2007-11-8 22:24 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
请教一个几组数据的方差排序的问题?
比如,有4组数据,A=[36,40,45,46];B=[40,48,52,36]
C=[37,41,47,45];D=[43,44,45,38]
需要在每个数组中选一个数,然后对选取的这四个数计算其方差,目的是找到方差最小、次小和第三小的三组数据,并且还要显示找到的这三组数据所对应原数组的数据!
我的初步想法是用循环来计算,即取A(i),B(j),C(k),D(m),然后循环求var(A(i),B(j),C(k),D(m)),最后找到方差最小的三组数据,并显示原A(i),B(j),C(k),D(m);
但是找到方差最小的三组数据,并且显示原数据是难题,恳请熟手指点,感激不尽!
回复
分享到:

使用道具 举报

发表于 2007-11-9 06:53 | 显示全部楼层
1.大概也只能如此
2. 这个命题的运算量太大, 具体意义何在?
 楼主| 发表于 2007-11-9 12:41 | 显示全部楼层
自然有这个问题的具体意义,急需高人指点
发表于 2007-11-9 22:10 | 显示全部楼层
第一:不要用循环,用meshgrid产生穷举数据
第二:使用方差函数后,sort函数排序;最后三个就是最小的三个

[ 本帖最后由 donkeyxu 于 2007-11-9 22:44 编辑 ]
发表于 2007-11-11 17:05 | 显示全部楼层
%程序如下,复制保存为haha.m,然后在主窗口输入haha即可,脱离了循环,速度很快;
function haha
A=[36,40,45,46];B=[40,48,52,36];C=[37,41,47,45];D=[43,44,45,38];
[x,y,z,w]=nmeshgrid(A,B,C,D);
x=x(:)';y=y(:)';z=z(:)';w=w(:)';
U=[x;y;z;w];
K=var(U);
[T,R]=sort(K);
disp('the fist 3 smaller answer is:');
disp(num2str(T(1:3)'));
disp('and the elements is:');
disp('A   B   C   D')
disp(num2str(U(:,R(1:3))'));
function [xx1,xx2,xx3,xx4]=nmeshgrid(x1,x2,x3,x4)
nx1 = numel(x1); nx2 = numel(x2); nx3 = numel(x3);nx4 = numel(x4);
xx1 = reshape(full(x1(:)),[1 nx1 1 1]);
xx2 = reshape(full(x2(:)),[nx2 1 1 1]);
xx3 = reshape(full(x3(:)),[1 1 1 nx3]);
xx4 = reshape(full(x4(:)),[1 1 nx4 1]);
xx1 = xx1(ones(nx2,1),:,ones(nx4,1),ones(nx3,1));
xx2 = xx2(:,ones(1,nx1),ones(nx4,1),ones(nx3,1));
xx3 = xx3(ones(nx2,1),ones(nx1,1),ones(nx4,1),:);
xx4 = xx4(ones(nx2,1),ones(nx1,1),:,ones(nx3,1));

评分

1

查看全部评分

 楼主| 发表于 2007-11-12 12:45 | 显示全部楼层
程序运行很好,谢谢donkeyxu 倾囊相助,在下佩服的很:lol ,现在正在学习中!
 楼主| 发表于 2007-11-20 21:25 | 显示全部楼层
希望donkeyxu 能给点程序的注解,不胜感激!
发表于 2007-11-20 22:36 | 显示全部楼层
 楼主| 发表于 2007-11-30 10:40 | 显示全部楼层
请问donkeyxu,如果我的数组变为5个数组,比如A=[36,40,45,46];B=[40,48,52,36]
C=[37,41,47,45];D=[43,44,45,38],E=[43,46,45]
我修改了一下程序,运行的结果不对,很着急!请你帮忙看一下,谢谢了
function haha
A=[36,40,45,46];B=[40,48,52,36];C=[37,41,47,45];D=[43,44,45,38];E=[43,46,45]
[x,y,z,w,v]=nmeshgrid(A,B,C,D,E);
x=x(:)';y=y(:)';z=z(:)';w=w(:)';v=v(:)';
U=[x;y;z;w;v];
K=var(U);
[T,R]=sort(K);
disp('the fist 3 smaller answer is:');
disp(num2str(T(1:3)'));
disp('and the elements is:');
disp('A   B   C   D   E')
disp(num2str(U(:,R(1:3))'));
function [xx1,xx2,xx3,xx4,xx5]=nmeshgrid(x1,x2,x3,x4,X5)
nx1 = numel(x1); nx2 = numel(x2); nx3 = numel(x3);nx4 = numel(x4);nx5 = numel(x5);
xx1 = reshape(full(x1(:)),[1 nx1 1 1 1]);
xx2 = reshape(full(x2(:)),[nx2 1 1 1 1]);
xx3 = reshape(full(x3(:)),[1 1 1 nx3 1]);
xx4 = reshape(full(x4(:)),[1 1 nx4 1 1]);
xx5 = reshape(full(x5(:)),[1 nx5 1 1 1]);
xx1 = xx1(ones(nx2,1),:,ones(nx5,1),ones(nx4,1),ones(nx3,1));
xx2 = xx2(:,ones(1,nx1),ones(nx5,1),ones(nx4,1),ones(nx3,1));
xx3 = xx3(ones(nx2,1),ones(nx1,1),ones(nx4,1),:,ones(nx5,1));
xx4 = xx4(ones(nx2,1),ones(nx1,1),:,ones(nx3,1),ones(nx5,1));
xx5 = xx5(ones(nx1,1),:,ones(nx4,1),:,ones(nx3,1),ones(nx2,1),);
发表于 2007-12-3 14:23 | 显示全部楼层

回复 #10 tigersyu 的帖子

相信看懂donkeyxu程序,应该就可以解决

只是调整了个参数而已,要多尝试

相信自己可以解决
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-29 10:00 , Processed in 0.061431 second(s), 21 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表