MathWorks - Mobile View
  • MathWorks 계정에 로그인합니다.MathWorks 계정에 로그인합니다.
  • Access your MathWorks Account
    • 내 계정
    • 나의 커뮤니티 프로필
    • 라이선스를 계정에 연결
    • 로그아웃
  • 제품
  • 솔루션
  • 아카데미아
  • 지원
  • 커뮤니티
  • 이벤트
  • MATLAB 다운로드
MathWorks
  • 제품
  • 솔루션
  • 아카데미아
  • 지원
  • 커뮤니티
  • 이벤트
  • MATLAB 다운로드
  • MathWorks 계정에 로그인합니다.MathWorks 계정에 로그인합니다.
  • Access your MathWorks Account
    • 내 계정
    • 나의 커뮤니티 프로필
    • 라이선스를 계정에 연결
    • 로그아웃

비디오 및 웨비나

  • MathWorks
  • 비디오
  • 비디오 홈
  • 검색
  • 비디오 홈
  • 검색
  • 영업 상담
  • 평가판 신청
  Register to watch video
  • Description
  • Full Transcript
  • Code and Resources

Solving ODEs in MATLAB, 10: Tumbling Box

From the series: Solving ODEs in MATLAB

Cleve Moler, MathWorks

Throw a rectangular object with sides of three different lengths (such as a cereal box), into the air. You can get the box to tumble stably about its longest axis, or about its shortest axis. But if you try to make it tumble about it middle axis, you will find the motion is unstable. The model of the angular momentum is a nonlinear system of three differential equations. There are six critical points: the four corresponding to the long and short axes are stable; the two corresponding to the middle axis are unstable.

Related MATLAB code files can be downloaded from MATLAB Central

 

Here are the differential equations for the angular momentum of a tumbling box. Try throwing a book, or a box, or any rectilinear object whose three dimensions are all different, into the air with a twist, to make a tumble.

You could go to rotate about its longest axis, or about its shortest axis. But you can't get to rotate about its middle axis. Let's examine that phenomena numerically.

Here's the anonymous function defining those system of three first order differential equations. Now I'm going to start with an initial condition that's near the first critical point. 1, 0, 0 is a critical point. And I'm going to take 0.2 times a random number, to sort of be near the critical point, and then normalize it, so that it has length 1.

So the largest component is the first component. And the other two are small, but not too small. This is an easy problem numerically. There's no stiffness involved here. So I'm going to use ODE 23, integrate from 0 to 10, and here's the solution.

The blue component is the first one, and it stays near 1. And the other two are periodic, rotating around the 0. Let's go back and take another starting condition. Here it is again.

Now the other two components are really small. And when we integrate that, the blue component stays flat near 1. And the other two guys hardly move at all.

Now I'm going to go to the third critical point, 0, 0, 1. Do the same thing. Take a random number near there. Use ODE 23. Now the yellow component stays near one. And the other two move periodically around 0.

Run that again. The third component is near 1. The other two are not too big. And run the ODE 23. The other component stays near 1. And the other two rotate periodically around 0.

Now we're going to go to the middle critical point. We're going to try and get the box to rotate around its middle axis. The second component is the one near 1. And now we see completely different behavior.

This sienna component doesn't stay near 1. It goes down near -1, and comes back up. Let's integrate over a longer period, so we can see that behavior.

So it's periodic. But it goes down to -1 and comes back to 1. And the other two move in large amplitudes around 0. So this is the instability of that middle critical value.

Let's doing again. Same thing. 1 down to -1, and back up. It's periodic. These solutions are all periodic. But that middle critical point is unstable. Now I want to view these in a different way, graphically.

The differential equations have these three critical points. Any solutions to start exactly in these initial conditions stay there. But what happens if you start near these initial conditions?

Well it turns out, that x and z are stable critical points. But y is an unstable critical point. If the angular momentum is near x or near z, it stays near there. But if it starts near y, it moves away quickly.

You can think of x as the short axis, and z is the long axis. Rotation near the short axis is stable. And rotation near the long axis a stable. But rotation near the middle axis is unstable.

We can see that in the following graphic. It turns out that if a solution starts with an initial condition that has norm 1, it stays with norm 1. So the solution lives on the unit sphere.

Here's our unit share with our three critical points, x, y, and z. If this were the Earth, z will be the North Pole. Axis where the 0-th meridian crosses the equator. That's in the east Atlantic, a little off of West Africa. y would be where the 90th meridian crosses the equator. That's in the Indian Ocean, west of Sumatra.

If we start with an initial condition near x, the solution orbits around x. That's stable rotation around the short axis. If we start with an initial condition near z, the solution orbits around z. That's stable rotation around the long axis.

But if we start near y, the solution takes off, goes over to near -y, turns around, and comes back to y. Periodic, but goes clear around the globe. That turns out, that's a circle actually, an orbit around x.

If we come up a little above y, we get an orbit around z. Go down a little below y, we get an orbit around -z. Go to the right of y, we get an orbit around -x.

Let's zoom in a little bit. And we can see that y is a classic unstable critical point. Let's conclude by drawing a few orbits.

Download Code and Files

Download the code in this video

Related Products

  • MATLAB

Learn More

Explore more resources from Cleve Moler
Related Information
Learn differential equations

Feedback

Featured Product

MATLAB

  • Request Trial
  • Get Pricing

Up Next:

predprey
14:16
11: Predator-Prey Equations
View full series (12 Videos)

Related Videos:

5:34
Solving ODEs in MATLAB, 9: The MATLAB ODE Suite
15:21
Solving ODEs in MATLAB, 1: Euler, ODE1
14:16
Solving ODEs in MATLAB, 11: Predator-Prey Equations
10:24
Solving ODEs in MATLAB, 12: Lorenz Attractor and Chaos
9:37
Solving ODEs in MATLAB, 3: Classical Runge-Kutta, ODE4

View more related videos

MathWorks - Domain Selector

Select a Web Site

Choose a web site to get translated content where available and see local events and offers. Based on your location, we recommend that you select: .

Select web site

You can also select a web site from the following list:

How to Get Best Site Performance

Select the China site (in Chinese or English) for best site performance. Other MathWorks country sites are not optimized for visits from your location.

Americas

  • América Latina (Español)
  • Canada (English)
  • United States (English)

Europe

  • Belgium (English)
  • Denmark (English)
  • Deutschland (Deutsch)
  • España (Español)
  • Finland (English)
  • France (Français)
  • Ireland (English)
  • Italia (Italiano)
  • Luxembourg (English)
  • Netherlands (English)
  • Norway (English)
  • Österreich (Deutsch)
  • Portugal (English)
  • Sweden (English)
  • Switzerland
    • Deutsch
    • English
    • Français
  • United Kingdom (English)

Asia Pacific

  • Australia (English)
  • India (English)
  • New Zealand (English)
  • 中国
    • 简体中文Chinese
    • English
  • 日本Japanese (日本語)
  • 한국Korean (한국어)

Contact your local office

  • 영업 상담
  • 평가판 신청

제품 소개

  • MATLAB
  • Simulink
  • 학생용 소프트웨어
  • 하드웨어 지원
  • File Exchange

다운로드 및 구매

  • 다운로드
  • 평가판 신청
  • 영업 상담
  • 가격 및 라이선스
  • MathWorks 스토어

사용 방법

  • 문서
  • 튜토리얼
  • 예제
  • 비디오 및 웨비나
  • 교육

지원

  • 설치 도움말
  • 사용자 커뮤니티
  • 컨설팅
  • 라이선스 센터
  • 지원 문의

회사 정보

  • 채용
  • 뉴스 룸
  • 사회적 미션
  • 영업 상담
  • 회사 정보

MathWorks

Accelerating the pace of engineering and science

MathWorks는 엔지니어와 과학자들을 위한 테크니컬 컴퓨팅 소프트웨어 분야의 선도적인 개발업체입니다.

활용 분야 …

  • Select a Web Site United States
  • 특허
  • 등록 상표
  • 정보 취급 방침
  • 불법 복제 방지
  • 매스웍스코리아 유한회사
  • 주소: 서울시 강남구 삼성동 테헤란로 521 파르나스타워 14층
  • 전화번호: 02-6006-5100
  • 대표자 : 이종민
  • 사업자 등록번호 : 120-86-60062

© 1994-2021 The MathWorks, Inc.

  • Naver
  • Facebook
  • Twitter
  • YouTube
  • LinkedIn
  • RSS

대화에 참여하기