How to create a function as the convolution of two functions? Thanks

Question

Tmat 2020-3-18

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/511627-how-to-create-a-function-as-the-convolution-of-two-functions-thanks

评论： Tmat 2020-3-19

Hi guys, I want to compute the convolution of two probability densities, and as an experiment I tried two uniform [0,1] densities:

fune=@(e)(e>=0).*(e<=1);
funz=@(z)(z>=0).*(z<=1);
funezt=@(ez,t)fune(ez-t).*funz(t);
funez=@(ez)integral(@(t)funezt(ez,t),0,2); % calculation using the definition of convolution
vecez=0:0.01:2;
vecfunez=funez(vecez);
plot(vecez,vecfunez)

However, this didn't work and the error message I got is:

Matrix dimensions must agree.
Error in convolutionexperiment>@(ez,t)fune(ez-t).*funz(t) (line 4)
funezt=@(ez,t)fune(ez-t).*funz(t);

Need some insignts and guidance on fixing this. Thanks!

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

David Goodmanson 2020-3-19

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/511627-how-to-create-a-function-as-the-convolution-of-two-functions-thanks#answer_420813

编辑：David Goodmanson 2020-3-19

在 MATLAB Online 中打开

Hello James,

fune=@(e)(e>=0).*(e<=1);
funz=@(z)(z>=0).*(z<=1);
funezt=@(ez,t)fune(ez-t).*funz(t);
funez=@(ez)integral(@(t)funezt(ez,t),0,2,'ArrayValued',true);   % <--
vecez=0:0.01:2;
vecfunez=funez(vecez);
plot(vecez,vecfunez)