错误:矩阵尺寸必须在matlab中一致，才能建立混淆矩阵。

Question

我知道以前有人问过这个问题，但我没有从他们那里得到答案。我也没有从mathworks.com那里得到多少，我正试图为ANN构建混淆矩阵。testSize是我拥有的testSet文件的数量。

output = sim(net,Yt);%testing network
    confusionMatrix(10,10);%building confusion matrix
for i=1:(testSize-2)*10
    for j=1:10
        m(j) = abs(output(i,i)-1);
end
minimum = find(m == min(m));
confusionMatrix(floor((i-1)/(testSize-2)) + 1,minimum) = confusionMatrix(floor((i-1)/(testSize-2))+1);
end

下面是我的全部密码。它用于基于testSets和trainSets的语音识别。

TrainSet=dir('TrainSet');

maxLength=0;


for i=3:length(TrainSet)
File=strcat('TrainSet\',TrainSet(i).name);
sub_direction=dir(File);
    for j=3:length(sub_direction)

        path=strcat(File,'\',sub_direction(j).name);
        sample=wavread(path); % Reading Files
        if(length(sample) > maxLength)
            maxLength = length(sample);
        end     
    end
end

TestSet=dir('TestSet');

for i=3:length(TestSet)
File=strcat('TestSet\',TestSet(i).name);
sub_direction=dir(File);
    for j=3:length(sub_direction)

        path=strcat(File,'\',sub_direction(j).name);
        sample=wavread(path);
        if(length(sample) > maxLength)
            maxLength = length(sample);
        end
    end

end


%Equalizing Trainset data

for i=3:length(TrainSet)
File=strcat('TrainSet\',TrainSet(i).name);
sub_direction=dir(File);
trainSize = length(sub_direction);
    for j=3:length(sub_direction)

        path=strcat(File,'\',sub_direction(j).name);
        trainSetMatrix(i-2,j-2,:) = zeros(1,maxLength);
        sample=wavread(path); % Reading Files
        for k=1:length(sample)
            trainSetMatrix(i-2,j-2,k) = sample(k);

        end 





    end
end


%Equalizing Testset data

for i=3:length(TestSet)
File=strcat('TestSet\',TestSet(i).name);
sub_direction=dir(File);
testSize = length(sub_direction);
    for j=3:length(sub_direction)

        path=strcat(File,'\',sub_direction(j).name);

        testSetMatrix(i-2,j-2,:) = zeros(1,maxLength);
        sample=wavread(path); % Reading Files
        for k=1:length(sample)
            testSetMatrix(i-2,j-2,k) = sample(k);

        end     
    end
end

for i=3:length(TrainSet)
File=strcat('TrainSet\',TrainSet(i).name);
sub_direction=dir(File);
    for j=3:length(sub_direction)
    frameSize = 160;%framing: each frame equals to 8khz*20msec = 160samples
    frameOverlap = 80;%overlapping is 50% of each frame
    framedTrainSetMatrix(1,:) = trainSetMatrix(i-2,j-2,1:160);
    for k=1:floor(maxLength/80)-2
        framedTrainSetMatrix(k+1,:)= trainSetMatrix(i-2,j-2,k*80:(k*80+frameSize-1));
    end
    end
end


for i=3:length(TrainSet)
File=strcat('TrainSet\',TrainSet(i).name);
sub_direction=dir(File);
    for j=3:length(sub_direction)
    window = hamming(frameSize);
    for k=1:floor(maxLength/80)-1
        windowedFrame(k,:) = framedTrainSetMatrix(k,:).*window';
        LinearPredictiveCoding(k,:) = lpc(windowedFrame(k,:),12);
        lpcResult(k,:) = LinearPredictiveCoding(k,2:13);
    end
    coefficient=12*(floor(maxLength/80)-1);
    X((i-3)*(trainSize-2)+(j-2),:)=reshape(lpcResult,1,coefficient);

    end
end




for i=3:length(TestSet)
File=strcat('TestSet\',TestSet(i).name);
sub_direction=dir(File);
    for j=3:length(sub_direction)
    frameSize = 160;%framing: each frame equals to 8khz*20msec = 160samples
    frameOverlap = 80;%overlapping is 50% of each frame
    framedTestSetMatrix(1,:) = testSetMatrix(i-2,j-2,1:160);
    for k=1:floor(maxLength/80)-2
        framedTestSetMatrix(k+1,:)= testSetMatrix(i-2,j-2,k*80:(k*80+frameSize-1));
    end
    end
end


for i=3:length(TestSet)
File=strcat('TestSet\',TestSet(i).name);
sub_direction=dir(File);
    for j=3:length(sub_direction)
    window = hamming(frameSize);
    for k=1:floor(maxLength/80)-1
        windowedFrame(k,:) = framedTestSetMatrix(k,:).*window';
        LinearPredictiveCoding(k,:) = lpc(windowedFrame(k,:),12);
        lpcResult(k,:) = LinearPredictiveCoding(k,2:13);
    end
    coefficient=12*(floor(maxLength/80)-1);
    Y((i-3)*(testSize-2)+(j-2),:)=reshape(lpcResult,1,coefficient);

    end
end


Xt = transpose(X);
T(10,1900);
hiddelLayer=1158;

net=newff(Xt,T,hiddelLayer);%building network

net.divideParam.trainRatio=0.2;

net.efficiency.memoryReduction=60;
net=train(net,Xt,T);%training network

Yt = transpose(Y);
output = sim(net,Yt);%testing network


%building confusion matrix
confusionMatrix(10,10);



for i=1:(testSize-2)*10
    for j=1:10
        m(j) = abs(output(i,i)-1);
    end
minimum = find(m == min(m));
confusionMatrix(floor((i-1)/(testSize-2)) + 1,minimum) = confusionMatrix(floor((i-1)/(testSize-2))+1);
end

for i=1:10
    for j=1:10
        confusionMatrix(i,j) = confusionMatrix(i,j)/(testSize-2);
    end
end

完全错误是：

?？?使用==> eq出错 矩阵尺寸必须一致。 ==>测试在167处的误差 最小=查找(m == min(m))；

请帮帮忙。

Answer 1

请参阅min的文档：

如果A是矩阵，则min(A)是包含每个列的最小值的行向量。

您的min(m)调用正在返回一个向量，而MATLAB无法在==中使用该向量( ==的输入A和B必须相同大小，除非其中一个是标量)，如此测试用例所示，这会产生相同的错误：

a = rand(4);
test = find(a == min(a));

如果要查找整个矩阵中最小值的索引，请使用find(a == min(a(:)));

或者，更简单地说：[minval idx] = min(a(:))

Answer 2

如果您想找到最小参数的索引，最好使用

[mn mIdx] = min( m(:) ); 

而不是

mIdx = find( m == min(m(:)) );

有关Matlab的argmax和argmin的更多信息，请参见this thread。