24小时热门版块排行榜    

Znn3bq.jpeg
查看: 362  |  回复: 0

542950171

金虫 (小有名气)

[求助] matlab求助

下面的一段代码为啥只实现了排序,期待解答
%基于熵的属性离散化方法,参见《一种基于熵的连续属性离散化方法》
function [result,dds]=shang(data)
[m,n]=size(data);
result=zeros(10,5);
for i=1:n
    %%%%%%%%%%%%%%%%%计算划分点
    setpoints=[];
    seqs=[];
    locs=[];
    %%%%%%%%%%%%%%%当出现有相通属性值的样本点时,去掉其余的只保留一个,为后面
    %为后面计算断点做准备
    [seqs,locs]=sort(data(:,i));
    se=seqs;
    for bbs=1:m-1
    if se(bbs)==se(bbs+1)
        qs=bbs+1;
        while se(bbs)==se(qs)
            se(qs)=inf;
        end
    end
    end
    se(find(se==inf))=[];
    newlen=length(se);
    %%%%%%%%%%%%%%%%%
    for j=1:newlen-1       %计算断点的值
        setpoints(j)=(se(j)+se(j+1))/2;
    end
    %%%%%%%%%%%%%%%%%
    nn=length(setpoints);%计算现有内部断点数
    for t1=1:nn+1   %计算各个断点处的卡方统计值
   %%%%%%%%%%计算卡方分布值
    if t1==1
    s1=length(find(-inf     elseif t1==nn+1
    s1=length(find(setpoints(t1-1)<=seqs&seqs     else
    s1=length(find(setpoints(t1-1)<=seqs&seqs     end
    num(t1)=s1;
    P(t1)=s1/m;
    end%接 for t1=1:nn+1
    Hmin(1)=-sum(P.*log2(P));%初始熵值
    k0=length(setpoints);
    H0=Hmin(1);
    Ck=zeros(1,2);%初始化Ck,另Ck(2)=0是为了它能够进入循环while Ck(2)-Ck(1)>=0
    %Ck(2)是合并区间后的值,在while循环中会设置
    %%%%%%%%%%%%%%%%%%循环
    symbol=0;
    while Ck(2)-Ck(1)>=0
    Ck(1)=Ck(2);
    if symbol==1
        setpt(ad)=[];
        setpoints=[];
        setpoints=setpt;
    end
    Q=[];
    H=[];
    num=[];
    setpt=[];
    nn=length(setpoints);%计算现有内部断点数
    setpt=setpoints;
    for c1=1:nn%计算去掉每个断点所得的熵值,找出其中熵差最小的一个作为合并对象
        setpoints=[];
        setpoints=setpt;
        setpoints(c1)=[];
        tt=length(setpoints);
        for t2=1:tt+1   %计算各个断点处的卡方统计值
        %%%%%%%%%%计算卡方分布值
        if t2==1
        s1=length(find(-inf         elseif t2==tt+1
        s1=length(find(setpoints(t2-1)<=seqs&seqs         else
        s1=length(find(setpoints(t2-1)<=seqs&seqs         end
        num(t2)=s1;
        Q(t1)=s1/m;
        end%接 for t1=1:nn+1
        H(c1)=-sum(Q.*log2(Q));%后续计算所得熵值
    end
    symbol=1;
    [Hmin(2),ad]=min(H);
%     Ck(2)=(nn)*Hmin(2)-Hmin(1)*(nn-1);
    Ck(2)=(k0)*Hmin(2)-H0*(nn-1);
   
    Hmin(1)=Hmin(2);
    end
    nt=length(setpt);
    result(1:nt,i)=setpt;
end
pointsnum=[];
    %将各个区间的值离散化
    dds=zeros(m,n);
    for i=1:n
        locd=[];
        locd=result(:,i);
        locd(find(locd==0))=[];
        pointsnum(i)=length(locd)+1;%区间个数
        if length(locd)~=0
        for j=1:pointsnum(i)
            if j==1
                dds(find(-inf             elseif j==pointsnum(i)
                dds(find(locd(j-1)<=data(:,i)&data(:,i)             else
                dds(find(locd(j-1)<=data(:,i)&data(:,i)<=locd(j)),i)=j;
            end
        end
        end
    end
附件中是那篇论文[ 来自科研家族 英乐盒 ]
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : 一种基于熵的连续属性离散化算法.pdf
  • 2012-04-19 09:06:28, 159.71 K

» 猜你喜欢

低头要有勇气,抬头要有底气,做人要有志气。吃得苦中苦,方为人上人。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 542950171 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 又一批高校组建人工智能学院 师资行吗 不是骗人吗 +6 yexuqing 2026-04-19 7/350 2026-04-23 12:32 by yexuqing
[基金申请] 国自然面上和省基金B类撒花 +18 花田半亩~白 2026-04-21 18/900 2026-04-23 11:31 by 12021227
[考研] 有没有学校收留 +3 蒋昌鹏qtj 2026-04-20 3/150 2026-04-22 20:25 by 学员JpLReM
[考研] 312求调剂 +3 山河似你温柔 2026-04-22 3/150 2026-04-22 20:17 by 学员JpLReM
[论文投稿] 急需审稿人!!! +3 陆小果画大饼 2026-04-21 3/150 2026-04-21 23:54 by jzy_123456
[考博] 申博/考博 +4 啃面包的小书虫 2026-04-17 8/400 2026-04-21 16:26 by 啃面包的小书虫
[考研] 295分求调剂 +6 ?要上岸? 2026-04-17 6/300 2026-04-21 08:18 by Equinoxhua
[考研] 085600材料与化工调剂 5+3 孜孜不倦2002 2026-04-19 6/300 2026-04-20 21:25 by babero
[论文投稿] 有没有接收比较快的sci期刊呀,最好在一个月之内的,研三孩子求毕业 20+4 之护着 2026-04-16 7/350 2026-04-20 15:45 by 豆豆7758
[考研] 337求调剂 +3 jyz04 2026-04-18 3/150 2026-04-20 12:24 by 研可安
[考博] 申博 +3 Xyyx. 2026-04-18 3/150 2026-04-20 10:44 by YuY66
[考博] 湖南大学刘巧玲课题组2026年第二批次博士研究生招生信息 +3 南风观火 2026-04-18 5/250 2026-04-20 10:13 by 南风观火
[考研] 求计算机方向调剂 +3 Toffee2 2026-04-16 6/300 2026-04-19 22:37 by ll叶
[考研] 294求调剂 +8 淡然654321 2026-04-17 9/450 2026-04-19 19:51 by Equinoxhua
[考研] 304求调剂 +8 castLight 2026-04-16 8/400 2026-04-19 17:14 by 中豫男
[考研] 求调剂 +6 苦命人。。。 2026-04-18 7/350 2026-04-19 16:27 by 中豫男
[考研] 接受任何调剂 +6 也就是栗子 2026-04-17 7/350 2026-04-18 17:20 by 涵竹刘
[考研] 260求调剂 +4 Zyt1314520.. 2026-04-17 5/250 2026-04-18 08:28 by babysonlkd
[有机交流] 二苯甲酮酸类衍生物 50+3 小白爱主人 2026-04-17 6/300 2026-04-17 18:47 by kf2781974
[考研] 322求调剂 +6 tekuzu 2026-04-17 6/300 2026-04-17 13:48 by Espannnnnol
信息提示
请填处理意见