how i can vectorize this code using bsxfun or meshgrid?

Question

0 个投票

Hi,

I am trying to vectorize this code, but I dont know how to start it, I read the help files for both bsxfun and meshgrid but no luck. Can someone please guide how to look at this rather than giving the code. My code is working for me but i want to improve it.

s1= 50;
s2=100;
E_=zeros(s2,s2);
grid_=linspace(-2.5,2.5,s2);
tau=linspace(-pi,pi,s1);
x=cos(tau);
y=sin(tau);
cur=myfun(x,y); %returns vector
for k=1:s2
    for m=1:s2
        for n=1:s1
            R= abs(sqrt((grid_(m)-x(n))^2+(grid_(k)-y(n))^2));
            u_(n)= besselj(0,R)*cur(n);
        end
        E_(k,m)=abs(trapz(tau,u_));
    end
end

Thanks

Turker

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

请先登录，再进行评论。

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

Follow Question

Answer 1

Jan 2017-3-14

编辑：Jan 2017-3-15

在 MATLAB Online 中打开

1 个投票

You can replace the inner loop by:

R  = abs(sqrt((grid_(m)-x) .^ 2 + (grid_(k)-y) .^ 2));  % [EDITED]: ^ to .^
u_ = besselj(0, R) .* cur;

I guess, that the abs() is not required.

4 个评论
显示 2更早的评论隐藏 2更早的评论

ttopal 2017-3-16

在 MATLAB Online 中打开

Hi Jan, it returns error as;

Assignment has more non-singleton rhs dimensions than 
non-singleton subscripts
Error in testfields (line 121)
        E_(k,m)=abs(trapz(tau,u_));

ttopal 2017-3-16

编辑：ttopal 2017-3-16

在 MATLAB Online 中打开

R  = abs(sqrt((grid_(m)-x) .^ 2 + (grid_(k)-y) .^ 2))'; %transpoze added.
u_ = besselj(0, R) .* cur;

Now it works! Thank you Jan!

请先登录，再进行评论。

how i can vectorize this code using bsxfun or meshgrid?

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

采纳的回答

4 个评论
显示 2更早的评论隐藏 2更早的评论

更多回答（0 个）

类别

标签

Community Treasure Hunt

how i can vectorize this code using bsxfun or meshgrid?

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

采纳的回答

4 个评论 显示 2更早的评论 隐藏 2更早的评论

更多回答（0 个）

类别

标签

另请参阅

Community Treasure Hunt

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

4 个评论
显示 2更早的评论隐藏 2更早的评论