검색결과 리스트
글
MATLAB을 이용해서 시뮬레이션을 하다보면 어떤 이유에서든 원신호에 약간의 오차를 주고 싶어질 때가 있습니다. 특히, 오차범위를 정해놓고 작업을 해야하는 경우가 간혹 발생하는데, 이부부을 구현하는 방법이야 많겠지만, 특별히 간단히 rand 함수를 이용해서 오차를 생성해보겠습니다.
먼저 대상 신호를 만들어보죠.
t=0:0.01:10;
y = sin(2*pi*t);
이렇게 하면
이런 신호가 하나 만들어 집니다.
이제 저의 목표는 저 신호를 10%이내 범위에서 파괴하는겁니다. 그럼 제목에서도 밝혔듯이 rand함수에 대해 간단히 알아보겠습니다.
rand함수는 랜덤변수를 생성합니다. 그런데 rand함수는 0과 1사이에서 임의의 변수를 생성합니다. 사용법은
rand(X,Y)
이렇게 하면 X*Y크기를 가진 0과 1사이의 임의의 실수로 만들어진 행렬을 생성해줍니다. 그래서 저는 방금전 코드에 이어서
e = rand(1,length(t))/10;
이렇게 에러를 생성합니다. 처음 데이터와 같은 크기를 가지도록하는 것이고, 뒤에 "/10"은 5%이내로 오차를 생성할 것이기 때문입니다. 왜 10%가 아니라 5%의 에러냐면, 센터를 기준으로 나눌것이기 때문입니다. 그래서 생성된 에러는
와 같습니다. 0부터 0.1까지의 값을 랜덤하게 가지게 된 것입니다. 이 값에서 0.95를 더하면
입니다.
이제 이 데이터를 원래 y값에다가 곱할 겁니다.
y1 = (e+0.95).*y;
이렇게 구현하면
이렇게 됩니다. 음... 확실히 하기 위해 피크치부분을 확대해 보겠습니다.
이렇게 해서 원하는 결과를 얻을 수 있습니다.
구현을 위한 전체코드는
t=0:0.01:10;
y = sin(2*pi*t);
e = rand(1,length(t))/10;
y1 = (e+0.95).*y;
plot(t, y)
axis([0 10 -1.5 1.5])
grid on
hold on
plot(t, y1, 'r')
입니다.
'The Robot > Prog.Lang.' 카테고리의 다른 글
| [MATLAB 연재] 9. 애니메이션 구현하기 (36) | 2010/04/07 |
|---|---|
| [MATLAB 연재] 8. 동역학 시뮬레이션 하기 (32) | 2010/04/07 |
| [MATLAB 연재] 7. Function 문의 사용 (32) | 2010/03/31 |
| RAND 함수를 이용하여 오차 생성하기 (22) | 2010/03/19 |
| [MATLAB 연재] 6. Simulink의 기초 2 (46) | 2010/03/17 |
| [MATLAB 연재] 5. Simulink의 기초 (60) | 2010/03/16 |
| [MATLAB 연재] 4. 조건문과 반복문, 에러대비 (36) | 2010/03/10 |
| [MATLAB 연재] 3. 기초연산법 및 데이터의 그래픽 출력 2 (26) | 2010/03/09 |
설정
트랙백
http://pinkwink.kr/trackback/273
전 공학은 다음 생애에서도 무리 ㅠ_ㅠ
저는 바람처럼님처럼 용감하게 배낭하나 들고 여행 다닐 용기가 없다는....^^
감사합니다 많은 도움이 되었습니다 이렇게 상세히 적어주시다니
이해가 잘 되네요 ㅎㅎ
네... 다행입니다...^^
근데.. 한가지 질문이 있는데요
파란색 선을 기준으로 빨간색(랜덤 변수) 선이 항상 크게 나타나나요? 오차라 하면은
피크치 부분이 파란색 선보다 0.05이내로 작게 나타날수는 없는 건가요?
저런... 음...
마지막 그림을 보시면 원래의 값보다 작은 값을 가지는 곳도 있고 높은 값을 가지는곳도 있는데요...ㅜㅜ
예를 들면 2.2초가 막 지난 부분이라던지...
3.2초 지점등등등등...^^
아 너무 친절한 설명!
머리에 쏙쏙 박히네용 감사합니당!!!
극찬에 몸둘바를 모르겠습니다... 감사합니다.^^
비밀댓글입니다
비밀댓글로 켄신님 이메일주소를 남겨주세요^^
비밀댓글입니다
발송했답니다..
아름다운 블로그 생활하세요^^
어렵지만 자꾸 보면 가깝게 느껴지는군요. ^^;
헉... 저..정말요??? ㅎㅎ
혹시 제가 설명을 잘해서는 아닐까요???? 헉...ㅜ.ㅜ(지성..^^)
좋은 한주 보내세요^^
허걱 영어보다 어렵네요,,, 원리는 알아듣게 설명은 잘해주셨는데...헐,
히히.. 전.. 영어가 가장 어렵답니다....ㅜ.ㅜ
신호에 error를 어떤식으로 주나 궁금해했던 적이 있었는데... 딱이거네요 ㅋㅋ
다른 여러가지 방법들이 있지만..
이 방법이 좀 심플하긴 하죠...
복잡한 신호에서의(대략위에서 세번째 그림같은) 양의 피크치는 어떻게 구하죠??
단순히 피크값을 구하고 싶으시다면, min, max함수를 사용해보시는 것도 괜찮을 겁니다. 만약 모든 피크치를 다 구하고 싶으시다면, 음... 약간 고민을 해봐야겠는데요.
하나 질문이 있습니다!!
위 5% 만큼의 에러를 생성하는 코드(e = rand(1, length(t))/10)에 추가적으로 -0.05를 해줘야 하지 않을까요?
그래야 -0.05~0.05 까지의 5% 에러를 낼 수 있는 uniform random code가 생성될 것 같습니다.
항상 좋은 자료 올려 주셔서 감사합니다!!
어떻게 구현하는가의 차이이지 않을까합니다.
그 후 저는 다시
y1 = (e+0.95).*y;
를 통해 조정을 했거든요
말씀하신것과 흡사한 방식으로 조절한것입니다.