"Arrays have incompatible sizes for this operation" with for loop in integral2

1 次查看（过去 30 天）
ho man
ho man 2024-3-13
评论： ho man 2024-3-13
I am going to find a fixed point for some intergal represented function.
(this is just an example, I do not know if the fixed point exists or not) by the iteration:
My code is like that:
input = @(y,s) y;
N_it=10;
output=cell(1, N_it);
output{1}=input;
density = @(z,t,s) exp( -(z.^2)./(2.*(s-t)) ) ./ (sqrt(2*pi).*sqrt(s-t));
for i=1:N_it
input=output{i};
integral = @(y,t) integral2( @(z,s) input(z,s).*density(z-y,t,s),0,100,@(t)t,T,'RelTol',0,'AbsTol',1e-5);
n_output=@(y,t)integral(y,t);
output{i+1}=n_output;
end
For N_it =1, I can get the answer, there is no error. For N_it>1, when I find the value of n_output(1,1), there is an error "Arrays have incompatible sizes for this operation." and then a very long list indicating my problematic lines.

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

回答（1 个）

VBBV
VBBV 2024-3-13
Ran in:
input = @(y,s) y;
N_it=10;
output=cell(1, N_it);
output{1}=input;
density = @(z,t,s) exp( -(z.^2)./(2.*(s-t)) ./ (sqrt(2*pi).*sqrt(s-t)));
for i=1:N_it
input=output{i};
integral = @(y,t) integral2( @(z,s) input(z,s).*density(z-y,t,s),0,100,@(t)t,T,'RelTol',0,'AbsTol',1e-5);
n_output=@(y,t)integral(y,t);
output{i+1}=n_output;
end
output(1,1)
ans = 1×1 cell array
{@(y,s)y}
  8 个评论
显示 6更早的评论
ho man
ho man 2024-3-13
Yes, your code is correct. Thank you so much. My code is the same as yours but the last line "output(1,1)" is replaced by "n_output(1,1)". What is the mistake here?

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Agriculture 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by