用Matlab实现apriori算法关联规则的挖掘程序，完整有详细注解

下面这段是apriori算法中由2频繁项集找k频繁项集的程序，程序中有两个问题：

1、似乎while循环的K永远都是固定的，也就是都是频繁2项集的个数。得到频繁3项集后K的个数不是要变吗？如何体现呢？

2、程序中有两个for的大循环，但是发兆族现结果是只要找到一个频繁3项集第二个for循环就会结束，但是其实还应该有其它的频繁3项集。for循环不是应该无条件执行到参数k结束吗？当时k值是15，可是程序结束的时候i=2，j=3,然后j就不执行4以及一直到k的部分了。是什么原因呢？麻烦高手指点一下。急啊……

while( k>0)

le=length(candidate{1})

num=2

nl=0

for i=1:k-1

for j=i+1:k

x1=candidate{i} %candidate初始值为频繁2项集，这个表示频繁项集的第i项

x2=candidate{j}

c = intersect(x1, x2)

M=0

r=1

nn=0

l1=0

if (length(c)==le-1) & (sum(c==x1(1:le-1))==le-1)

houxuan=union(x1(1:le),x2(le))

%树剪枝,若一个候选项的某个K-1项子集为非频繁，则剪枝掉

sub_set=subset(houxuan)

%生成该候选项的所有K-1项子集

NN=length(sub_set)

%判断这些御猜段K-1项自己是否都为频繁的

while(r &M<NN)

M=M+1

r=in(sub_set{M},candidate)

end

if M==NN

nl=nl+1

%候选k项集

cand{nl}=houxuan

%记录每个候选k项集出现的次数

le=length(cand{1})

for i=1:m

s=cand{nl}

x=X(i,:)

if sum(x(s))==le

nn=nn+1

end

end

end

end

%从候选集中找频繁项集

if nn>=th

ll=ll+1

candmid{nl}=cand{nl}

pfxj(nl).element=cand{nl}

pfxj(nl).time=nn

disp('得到的频繁项集为:'镇誉)

result=(candmid{nl})

disp(result)

end

end

end

end

主要问题包括：

1、变量x未定义。

2、函数MUSIC里面：

S=[S(257:512)S(1:256)]

应为

S=[S(257:512) S(1:256)]

另外，clearR未定义，不知道干什么用的，可闭告以直枣庆接删掉。

3、函数ARMA里面，调用的Burg未定凳态握义。

function R=main_Dj()

clcclear

G=[1 2 51 4 12 3 12 4 62 5 5.82 6 5.72 7 5.63 7 2 3 11 1.5 3 12 4...

    4 5 0.54 8 3 5 6 15 9 36 7 0.66 10 2.5 7 11 2.78 9 18 12 6...

    9 10 1.59 12 510 11 0.510 12 411 12 3]

opt=0

route=sroute(G,opt)

R=[]

r=route(3,end)

R=[r,R]

while r~=1

  r=route(3,r)

  R=[r,R]

end

R=char(R+64)

R=[R,'O']

end

function route=sroute(G,opt)

% 求图的最短路的Dijkstra算法，规定1是起液漏点

% G是给定图的邻接矩阵或弧表矩阵，程序能够自动识别

% 当opt=0（或缺省）时求无向图的最短路，opt=1时求有向图的最短路

% d——标记最短距离 

% route是一个矩阵，第一行标记顶点，第2行标记1到该点的最短距离，

% 第3行标记最短路上该点的先驱顶点

if (nargin==1) opt=0 end

while 1      % 此循环自动识别或由弧表矩阵生成邻接矩阵

   if G(1,1)==0

      A=G

      n=size(A,1)

      M=sum(sum(A))break

   else

      e=G

      n=max([e(:,1)e(:,2)])          % 顶点数

      m=size(e,1)                     % 边数

      M=sum(e(:,3))                   % 代表无穷大

      A=M*ones(n,n)

      for k=1:m

          A(e(k,1),e(k,2))=e(k,3)  

          if opt==0 

               A(e(k,2),e(k,1))=e(k,3) % 形成无向图的邻接矩阵

          end         

      end

      A=A-M*eye(n)                      % 形成图的邻接矩阵      

   end

   break

end

pb(1:length(A))=0pb(1)=1               % 永久标号点记为1

index1=1                                % 依次记录永久标号顶点

index2=ones(1,length(A))                % 标记最短路上各点的先驱顶点

d(1:length(A))=Md(1)=0                 % 标记距离

temp=1                                  % 标记最近一个永久标号点

while sum(pb)<length(A)

   tb=find(pb==0)                       % 找出临时标号点

   d(tb)=min(d(tb),d(temp)+A(temp,tb))  % 更新距离

   tmpb=find(d(tb)==min(d(tb)))         % 确定新最小距离点

   temp=tb(tmpb(1))                     % 其中之一记为新永久标号点

   pb(temp)=1                           % 增加新永久标号点

   index1=[index1,temp]                 % 记录新永久标号点

   index=index1(find(d(index1)==d(temp)-A(index1,temp)')) % 确定前驱顶点

   if length(index)>=2                   % 前驱顶点多于1个时取第一个

      index=index(1)

   end

   index2(temp)=index         碧枝          % 记录前驱顶点

end

route=[1:n d index2]

end

运行结果

R =

ADEFJKO

代码自己看，不闹慧烂解释，也别叫我解释了，很麻烦的。

---