PinkWink's Blog

안녕하세요 요즘 자동제어 공부를 하고 있는데 정말로 많은 도움 받고 있습니다 ㅎㅎ
그런데 한가지 궁금한게... 미방을 시뮬링크 말고 직접 풀순 없나요?
룬지 쿠타 알고리즘이나 매틀랩 내장함수를 이용해서요

위의 정리된 식들을 주욱 따라가다 보니, 역행렬 구하는 과정에서 계산 오류가 있는 것 같습니다. 1/(ad-bc) 에서 bc가 (m^2)(l^2) 인데, 따라서
1/{(M+m)(ml^2+I)-(m^2)(l^2)} 가 되어야 할것입니다.

200x200x200 크기의 매트릭스를 z축에대해서 100도 회전하고 x축에관해서 10회전시키고 싶은 경우에 이용할 수 있는 매트랩 내장함수나 다른 방법이있는지 알고 싶어서요... 부탁^^;;

앗..질문이 구체적이 못해서 죄송합니다. Matrix라는 말을 행렬이라고 표현했는데요. 제코드가 일단 텍스트화일로부터 데이터를 읽어들이는데, 그 데이타크기가 (x=200,y=200,z=200)이고, 데이타들이 x,y,z각각의 위치에 저장되어있습니다. 이 데이터를 meshgrid를 사용하면 3D그래픽화가능한데, 원점이 (0,0,0)인데요, 저는 회전축을을 (x=100,y=100,z=100)으로 이동시킨상태에서 z축에대해 100도회전하고, x축에관해 10도 회전시켜보고 싶어서요... 이렇게하려면 어떤 방법이있을지 여쭤보고싶습니다.

예. 그렇게 해보았는데, array index가 integer가 아니라는 오류때문에...

저는 잘 안되서 그러는데요. 계속 에러나오고요..
??? Subscript indices must either be real positive integers or logicals.

Error in ==> test at 31
val_voxel(updX, updY, updZ)=temp_val_dose(ix,jy,kz);
라구요.
제 코드좀 봐주시구요. 도움좀 주세요. 여기서 막혀서 진행이 안되서.. 도움을 바랍니다.

count = 1;
iteration = 200;
for kz = 1:iteration
for jy = 1:iteration
for ix = 1:iteration
val_test = count;
temp_val_dose(ix,jy,kz) = val_test;
count = count + 1;
end
end
end

the = 108*(2*pi)/360;
phi = 10*(2*pi)/360;

for kz = 1:iteration
for jy = 1:iteration
for ix = 1:iteration

cnt = iteration;
updX = ix*cos(the)*cos(phi)-jy*cos(the)*sin(phi)-kz*sin(the);
updY = jy*cos(phi);
updZ = ix*sin(the)*cos(phi)-jy*sin(the)*sin(phi)+kz*cos(the);

val_voxel(updX, updY, updZ)=temp_val_dose(ix,jy,kz);
cnt
end
end
end

예.. 코드전체를 봐달란 뜻은 아니였는데, 핑크윙크님께선 된다고 하시니까 제건 왜 안되는지 알고 싶었습니다. 제 생각은 array index를 좌표상의 공간으로 보고 회전변환을 시키면 되지않을까했었거든요...

정인 //

예전에 만들었던 코드인데, 지금 하고자 하시는 것과 비슷한 것 같아서 올려보입니다. 혹시나 도움 되실까 하여..

clear all; close all; clc;
d = randn(3,200);

theta_ = deg2rad(0);
phi_ = deg2rad(1);
psi_ = deg2rad(-2.5);

rotMat = [...
cos(theta_) * cos(psi_),...
-cos(phi_)*sin(psi_)+sin(phi_)*sin(theta_)*cos(psi_),...
sin(phi_)*sin(psi_)+cos(phi_)*sin(theta_)*cos(psi_);
cos(theta_)*sin(psi_),...
cos(phi_)*cos(psi_)+sin(phi_)*sin(theta_)*sin(psi_),...
-sin(phi_)*cos(psi_)+cos(phi_)*sin(theta_)*sin(psi_);
-sin(theta_),...
sin(phi_)*cos(theta_),...
cos(phi_)*cos(theta_);];

D = rotMat*d;

figure('color','w'),
plot3(d(1,:),d(2,:),d(3,:),'r*'); box on; hold on; axis image; grid on;
xlabel('X'); ylabel('Y'); zlabel('Z');
camproj('perspective'); view(gca,[-50 30]);
xlim([-max(d(:)),max(d(:))]);
ylim([-max(d(:)),max(d(:))]);
zlim([-max(d(:)),max(d(:))]);
plot3(D(1,:),D(2,:),D(3,:),'b^')

(핑크윙크님께서 기분 상해하시면 어쩌나 조마조마하며 ㅎㅎ;;)

보여주신 코드는 감사합니다. 그런데도 제가 아직 잘 이해를 못해서... 귀찮으실지모르지만, 제경우에 맞게 코드를 다시변형해 주시면 안될까요?
저는 단순히 3차원회전행렬변환을 z축변환 x축변환을 곱해서 다음과같은 변환행렬을 사용했는데요.
[ cos(the)*cos(phi)-cos(the)*sin(phi)-sin(the);
cos(phi);
sin(the)*cos(phi)-sin(the)*sin(phi)+cos(the); ];
그리고, 200x200x200인 matrix를 z축으로 100도 x축으로 10도 회전시키는 경우일때를 가정해서, 회전시켜보고싶은데요. 저 코드를 어떻게 하면좋을지 도움을 부탁드립니다.

따뜻하신 말씀에 감사합니다. 보여주신 코드를 실험해보았는데요. 핑크윙크님과 별밤님께서는 점(1,0,0)을 X,Z축 변환해서 점(0.707,-0.696,0.122)으로 이동하게되는데요. 거기까지는 이해되는데요. 제가 원하는 것은 점(1,1,1)에 어떤 데이타 값이 들어있습니다. 4차원데이타라고 부르던데... 해서 점(1,1,1)은 변환해서 위치가점(0.7,-0.6,0.1)로 변환된 위치에 점(1,1,1)에 들어있는 데이타값을 변화시키지않고 위치 시키려고하는데, array index가 정수만 가능한터라 데이터값을 점(0.7,-0.6,0.1)에 매핑하고 싶은데, 불가능한질문인가요?

네.. 감사합니다.

실행하고 나서 scope를 더블클릭 하면 마지막 캡쳐 부분의 그림이 뜨잖아요
근데 제가 그 그림의 axis를 조정하려고 그 그림에서 마우스 오른쪽클릭을
하니까 Axes properties 라구 조절하는 부분이 있더라구요. 근데 여기선 Y축만 조절이 되고 X축은 없던데... X축까지 axis를 조절 하려면 어떻게 해야되나요??

안녕하세요 정말 잘 배우고 있습니다.
질문이 있어요, 힘은 '0'이고, 일정 각도(초기치)에 자유진동을 하게끔 설정되었는데
결과 그래프를 보면, 왜 변위 'x'도 저렇게 진동하게 되는걸까요?
진자의 진동이 미세하게나마 카트에 진동을 일으킨다고 해석해야되나요?

안녕하세요 학교에서 시뮬링크를 이용한 프로젝트를 수행중에 꼭 필요한 부분이 생겨 질문드립니다.

혹시 input을 반복되는 step function으로 설정하는 방법 아시나요?

0에서 1초까지 3
1에서 2초까지 0
2에서 3초까지 3
3에서 4초까지 0
이런 식으로 무한 반복되는 input을 입력하는 방법을 아신다며 좀 알려주시겠어요?