필터 지우기
필터 지우기

better method for evaluating matrix

조회 수: 2 (최근 30일)
JaeSung Choi
JaeSung Choi 2017년 11월 21일
편집: Stephen23 2017년 11월 21일
For given vector such as following, I want to make square matrix s such that
%constant, original form
N = 500;
x = 2*pi*linspace(0,1,N);
for i = 1:N
for j = 1:N
s(i,j) = sin(x(i)-x(j))
end
end
But it was too slow, so recently I edited it to following, but it's still too slow!!
Can anybody please help me??
%constant
N = 500;
x = 2*pi*linspace(0,1,N);
for i = 1: N;
s(:,i) = x-x(i);
end
s=sin(s);

채택된 답변

Stephen23
Stephen23 2017년 11월 21일
편집: Stephen23 2017년 11월 21일
You really need to learn how to write vectorized code. Solving every task using lots of ugly loops is not an efficient way to write MATLAB code. Try this:
N = 500;
vec = 2*pi*linspace(0,1,N);
mat = sin(bsxfun(@minus,vec(:),vec))

추가 답변 (1개)

Andrei Bobrov
Andrei Bobrov 2017년 11월 21일
편집: Andrei Bobrov 2017년 11월 21일
N = 500;
x = 2*pi*linspace(0,1,N);
s = sin(x(:)' - x(:));
for old versions of MATLAB:
s = sin( bsxfun(@minus, x(:)',x(:)) );
  댓글 수: 1
Stephen23
Stephen23 2017년 11월 21일
편집: Stephen23 2017년 11월 21일
Note that the output is transposed compared to the code given in the question.

댓글을 달려면 로그인하십시오.

카테고리

Help CenterFile Exchange에서 Creating and Concatenating Matrices에 대해 자세히 알아보기

Community Treasure Hunt

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

Start Hunting!