How to adjust the "isnan" to make it work for cell?

1 次查看（过去 30 天）

Ismail Qeshta 2018-2-20

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/383681-how-to-adjust-the-isnan-to-make-it-work-for-cell

评论： Ismail Qeshta 2018-2-21

Hi I have the below for loop. I would like to replace the NaN data output with 0.07. I keep getting the below error message. Could anyone please advise me how to stop the error message? Thank you.

for i=1:size(A,2)%number of columns
x{i} = interp1(y3, x3, A(:,i), 'linear');
k=1:1;
temp=x(k,:);
temp(isnan(temp))=0.07;
x(k,:)=temp;
 fid=fopen(['result_' num2str(1) '.txt'],'w');
 fprintf(fid,'%f\n',x); 
 fclose(fid);
end

The error message:

Undefined function 'isnan' for input arguments of type 'cell'.
Error in InterpolaeMFJavadCompleted (line 24)
temp(isnan(temp))=0.07;

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Birdman 2018-2-20

Check Stephen's answer.

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

Stephen23 2018-2-20

2
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/383681-how-to-adjust-the-isnan-to-make-it-work-for-cell#answer_306063

编辑：Stephen23 2018-2-20

在 MATLAB Online 中打开

Don't waste your time with cell arrays and cell2mat and the like, you don't need them. Just use a numeric variable:

N = size(A,2) %number of columns
%C = cell(1,N);
for k = 1:N
    tmp = interp1(y3, x3, A(:,k), 'linear');
    tmp(isnan(tmp))=0.07;
    %C{k} = tmp;
    fnm = sprintf('result_%d.txt',k);
    [fid,msg] = fopen(fnm,'wt');
    assert(fid>=3,msg)
    fprintf(fid,'%f\n',tmp); 
    fclose(fid);
end

I also made some other small changes to make your code more robust.

25 个评论
显示 23更早的评论隐藏 23更早的评论

Ismail Qeshta 2018-2-20

编辑：Ismail Qeshta 2018-2-20

在 MATLAB Online 中打开

I cannot thank you enough Stephen for your great help.

Actually you are right, the code produces 400 file each has 10 column, but it has only one row.

I am using the following code that combines my plots and your provided code above.

clear; clc;
Folder = cd;
N=100;
for k = 1:N;
  Driftt  = sprintf('Drift%d.out', k);
  Reactt  = sprintf('React%d.out', k);
  matDrift     = importdata(fullfile(Folder, Driftt));
  matReact = importdata(fullfile(Folder, Reactt));
  x1= matDrift(:,2);
  y1= -sum(matReact(:,2:11),2);
[x3, ix] = unique(x1);
y3       = y1(ix);
A=dlmread('AllPiers2.txt');
Nr = size(A,1) % number of rows
Nc = size(A,2) % number of columns
Np = numel(100) % number of plots
fmt = repmat(',%f',1,Nr);
fmt = [fmt(2:end),'\n'];
for kc = 1:Nc
    fnm = sprintf('result_%d.txt',kc);
    [fid,msg] = fopen(fnm,'wt');
    assert(fid>=3,msg)
    for kp = 1:Np
        ...
        tmp = interp1(y3, x3, A(:,kc), 'linear');
        tmp(isnan(tmp))=0.07;
        fprintf(fid,fmt,tmp);
    end
    fclose(fid);
end
end