模糊推理系统

时间:2011-04-20 05:58:13

标签: matlab fuzzy-logic

我试图通过使用evalfis()函数提供iput来使用.fis文件。我正在以下列形式提供意见:

        jo= myImage(2);   % jo contains the graysacle value of pixel at location 2         

        a=readfis('terrain_classification_fuzzy');
        Result = evalfis(jo,a);

但我收到以下错误:

???使用==>时出错evalfismex 第一个输入必须是定义的DOUBLE矩阵

==>中的错误evalfis在84 [输出,IRR,ORR,ARR] = evalfismex(输入,fis,numofpoints);

==>中的错误see_movie_segmentationTry4_incorrect at 113             结果= evalfis([jo],a);

请帮助我...

2 个答案:

答案 0 :(得分:2)

检查这个例子和参数,看看你是否有类似的参数:

function trn_2in(mf_n, epoch_n)    
% This script requires the Fuzzy Logic Toolbox from the MathWorks.
if nargin < 2, epoch_n = 100; end;    
if nargin < 1, mf_n = 4; end;    
% ====== collect training data    
point_n = 11;    
x = linspace(-10, 10, point_n);    
y = linspace(-10, 10, point_n);    
[xx, yy] = meshgrid(x, y);    
tmp1 = sin(xx)./(xx);    
index = find(isnan(tmp1)==1);    
tmp1(index) = ones(size(index));  
tmp2 = sin(yy)./(yy);    
index = find(isnan(tmp2)==1);    
tmp2(index) = ones(size(index));  
zz = tmp1.*tmp2;    
trn_data = [xx(:) yy(:) zz(:)];       
% ====== training options    
ss = 0.1;    
ss_dec_rate = 0.9;    
ss_inc_rate = 1.1;
mf_type = 'gbellmf';
% ====== generate the initial FIS     
in_fismat = genfis1(trn_data, mf_n, mf_type);
% ====== start training    
[trn_out_fismat trn_error step_size] = ...    
    anfis(trn_data, in_fismat, [epoch_n nan ss ss_dec_rate ss_inc_rate]);    
% ====== compute ANFIS output     
z_hat = evalfis([xx(:) yy(:)], trn_out_fismat);  
% ====== plot of training data
genfig('training data');    
subplot(221);    
mesh(xx, yy, zz);    
limit = [min(xx(:)) max(xx(:)) min(yy(:)) max(yy(:)) ...    
    min(zz(:)) max(zz(:))];
    axis(limit); set(gca, 'box', 'on');    
xlabel('X'); ylabel('Y'); title('Training data');    
    zz_hat = evalfis([xx(:) yy(:)], trn_out_fismat);    
subplot(222);    
mesh(xx, yy, reshape(zz_hat, point_n, point_n));    
axis(limit); set(gca, 'box', 'on');    
xlabel('X'); ylabel('Y'); title('ANFIS Output');    
subplot(223)    
plot(1:epoch_n, trn_error);
xlabel('epoch number'); ylabel('root mean squared error');    
title('error curve');
subplot(224)    
plot(1:epoch_n, step_size);    
xlabel('epoch number'); ylabel('step size');    
title('step size curve');    
% ====== plot MFs    
figH = genfig('MFs for SINC function training');    
% plot initial MFs on x and y    
subplot(221); plotmf(in_fismat, 'input', 1);    
subplot(222); plotmf(in_fismat, 'input', 2);    
subplot(223); plotmf(trn_out_fismat, 'input', 1);
subplot(224); plotmf(trn_out_fismat, 'input', 2);    
delete(findobj(figH, 'type', 'text'));    
subplot(221); title('Initial MFs on X');    
subplot(222); title('Initial MFs on Y');    
subplot(223); title('Final MFs on X');    
subplot(224); title('Final MFs on Y');

答案 1 :(得分:2)

我发现了错误。只是“jo”最初设置为int8而不是double类型。