ODE-matlab
이전 댓글 표시
Hello.
I want to solve this equation system:
x ̈ = q ( y ̇ B z − z ̇ B y ) m y ̈=−q (x ̇Bz−z ̇Bx) m z ̈= q (x ̇By −y ̇Bx) m x(0) = x0, y(0) = y0, z(0) = z0 x ̇(0) = v0,x, y ̇(0) = v0,y, z ̇(0) = v0,z.
and with B = [1 1 1]' and m=q=1.
When I look at the equation system it seems for me very complicated and I really don't have a clue how to solve this. But, here is my attempt so far:
function [dU]=ODE(t,U)
dU = zeros(size(U));
q=1;m=1;
x(0)=1; % Can I assume what x(0),y(0),z(0) , x1(0),y2(0) and z3(0) should be?
y(0)=1;
z(0)=1;
x1(0)=2;
y2(0)=2;
z3(0)=2;
v=[1;1;1];
x=[2;2;2];
B=[0;0;1];
Tt=100;
Tp=1;
B=Tt*file1(x)+Tp*file2(v);
x2=q./m*(diffy(y)*B(:,3)-diffy(z)*B(:,2))
y2=-q./m*(diffy(x)*B(:,3)-diffy(z)*B(:,1))
z2=q./m*(diffy(y)*B(:,3)-diffy(x)*B(:,1))
dy(1)= y(3)
dy(2)= y(4)
dy(3)= ?
dy(4)= ?
dy(5)= ?
dy(6)= ?
Can someone please, help me explain this?
Regards Cillian
댓글 수: 2
Cillian
2012년 6월 11일
Walter Roberson
2012년 6월 11일
Are you attempting to execute code that has the question marks in it ??
답변 (0개)
카테고리
도움말 센터 및 File Exchange에서 Mathematics에 대해 자세히 알아보기
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!