알쓸유공

예제 1) 

%일단 3x3 행렬을 만들어 보도록 하겠습니다. 

>> A=[1 2 3 ; 4 5 6 ; 7 8 9]

A= 

     1     2     3 
     4     5     6 
     7     8     9 

%이때, 행렬 A의 특정 위치의 원소값을 행과 열을 위치로 확인이 가능하지만 행렬의 전체 원소의 개수에 따라서 확인이 가능합니다. 

%매트랩은 기본적으로 열(column)벡터를 기준으로 연산하기 때문에 A(1)=1, A(3)=7, A(4)=2, etc., 와 같은 순서대로 배치가 됩니다.

%그럼 이어서 A 행렬의 첫번째 열의 순서를 랜덤하게 재배치 해보도록 하겠습니다.

>>idx = randperm(3)

idx = 

      3 1 2 

% randperm(n) 함수는 1에서 n(구간의 끝점 포함)까지의 정수로 구성된 난수 순열을 행 벡터로 

% 반환합니다. 즉, 내가 재배치하고 싶은 벡터의 길이(n)에 맞춰 난수를 생성하여 사용하면 됩니다.

>> A( : , 1 ) = A( idx, 1)

A = 

     7     2     3

     1     5     6

     4     8     9 

% A 행렬의 첫번째 열벡터의 원소 순서가 바뀐걸 확인할 수 있습니다. 

예제 2)

% 예제1을 응용하면 아까 말씀드린 행렬 원소 전체를 랜덤하게 재배치할 수 있습니다.

>> A=[1 2 3 ; 4 5 6 ; 7 8 9]

A= 

     1     2     3 
     4     5     6 
     7     8     9 

>>idx2 = randperm( numel(A) )

% randperm함수를 이용해서 A 행렬의 원소 개수에 맞춰 난수를 생성해봅시다. 이때 numel(A) 함수는 A의 전체 원소개수를 출력해줍니다. 

%즉, 여기서는 A행렬의 원소가 9개이니 numel(A)의 값은 9가 되겠죠? 

idx2 =

     5     1     8     2     9     7     3     4     6

>>A = A(idx2)

A = 

      5     1     6     4     9     3     7     2     8

% A(idx2) = A(5), A(1), A(8),... A(6) 을 의미하게 됩니다. 그러면 각각의 원소값을 확인해보면 A(5)=5, A(1)=1, A(8)=6,... 이 되니까 이 순서대로 재배치가 되겠죠?

% 하지만 우리가 원하는건 행벡터가 아닌 처음 3x3 행렬의 형태이니 다시 원래대로 돌려 놓습니다.

>>A = reshape(A,3,3)

예제 3)

%예제 1과 2에서 배운 내용을 한꺼번에 정리해서 아래와 같이 사용하시면 됩니다. 

A = reshape( A( randperm( numel(A) ) ), 3, 3);

CVX is a modeling system for constructing and solving disciplined convex programs (DCPs). CVX supports a number of standard problem types, including linear and quadratic programs (LPs/QPs), second-order cone programs (SOCPs), and semidefinite programs (SDPs). CVX can also solve much more complex convex optimization problems, including many involving nondifferentiable functions, such as `1 norms. You can use CVX to conveniently formulate and solve constrained norm minimization, entropy maximization, determinant maximization, and many other convex programs. 

As of version 2.0, CVX also solves mixed integer disciplined convex programs (MIDCPs) as well, with an appropriate integer-capable solver. 

%% 원 그리기 

x1 = 0;  % 원 중심의 x좌표

y1 = 0;  % 원 중심의 y좌표

rc = 10; % 원의 반지름

[x,y] = cylinder(rc,200);  % z축의 파라미터를 입력할 경우 3차원의 원통형이 생성되지만 여기서는
                               % 2차원의 원만 표기하기 위해 z축 파라미터는 사용하지 않습니다. 

                               % 총 200개의 요소로 이루어진 원의 바운더리 좌표가 x, y에 각각 저장

figure(1)

plot(x(1,:)+x1, y(1,:)+y1, 'b', 'LineWidth', 2), hold on, grid on;

%앞에서 설정한 원의 중심값을 각 x, y좌표에 더해줌으로써 (x1,y1)을 중점으로 하는 반지름 10인 원을 그려주게 됩니다.'

title(['\fontsize{15}Radius=',num2str(rc)]);

xlabel('[m]'), ylabel('[m]');

1) 원의 중심을 기준으로 가장 멀리 떨어진 점의 거리(길이)를 반지름 이하가 되도록 한다. 

  -> rc값 이하의 임의의 숫자를 생성( r = rc*sqrt(rand) )

2) 원의 중심을 기준으로 임의의 라디안을 생성한다.

  -> gen_rad = 2*pi*rand;

3) 1)과 2)에서 생성한 값을 이용해 원의 중점(x1, y1)으로부터 임의의 위치에 point를 생성한다.

  -> x = r*cos(a)+x1, y = r*sin(a)+y1;

%% 임의의 점 생성

a = 2*pi*rand( 1, length(x) );

r = sqrt( rand( 1, length(x) ) );

x11 = rc.*r.*cos(a) + x1;

y11 = rc.*r.*sin(a) + y1;

plot(x11, y11, '.r')