cellfun with gradient [x,y]

1 次查看（过去 30 天）

Spyros Polychronopoulos 2017-12-4

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/371042-cellfun-with-gradient-x-y

评论： Spyros Polychronopoulos 2017-12-15

I have a structure array q(1:50,1).p(1:50,1:50). And I was taking the gradient:

for k1=1:100
[part_x(k1).p,part_y(k1).p]= gradient(q(k1).p(:,:)); %gradient
end

I was getting part_x(1:50,1).p(1:50,1:50) and part_y(1:50,1).p(1:50,1:50)

Now I have functions (z) so I created the same structure but with cells: q{1:50,1}.p{1:50,1:50}, and I am trying to do the same but it doesn't work

for k1=1:100
[part_x{k1}.p,part_y{k1}.p]=@(z) (cellfun(@gradient,q{k1}.p{:,:}(z)));
end

Any ideas? I tried arrayfun as well but not working either.

Here is a code example to clarify things:

for k1=1:10
q1(k1).p=rand(10,10);
end
for k1=1:10
[part_x1(k1).p,part_y1(k1).p]= gradient(q1(k1).p); %gradient
end
% Now with function handle the same as above
A=rand(10^3,1);
for k2=1:10^3
a{k2}=@(z) z*A(k2);
end
%preparing the array for the gradient
l=1;
for k3=1:10
for k4=1:10
    for k5=1:10
q2(k3).p{k4,k5}=@(z) a{l}(z);
l=l+1;
    end
end
end
clear k1 k2 k3 k4 k5 l
% now I want to do what I have done with the matrix above but I have
% functions
%this returns 1 value so doesn't work
for k1=1:10
[part_x2{k1}.p,part_y2{k1}.p]=@(z) (cellfun(@gradient,q2(k1).p(z)));
end
%this doesn't work either
for k1=1:10^3
[part_x2{k1}.p,part_y2{k1}.p]=arrayfun(@(z) gradient(q2(k1).p(z)));
end