I don't understand what is going wrong with my code!!, I get error of a matrix singular
2 visualizaciones (últimos 30 días)
Mostrar comentarios más antiguos
function dy = rocketequat( t,y )
%UNTITLED6 Summary of this function goes here
g=9.81; % constant of gravity m/s^2
Ve=2800;%celerity of inject gas m/s
m0=500000;%initial mass =500000kg
dmdt=-2300;%kg/s
dy(1)= -g*sin(y(2)) - Ve * dmdt/( m0 - dmdt*t );%dy(1)=dvdt acceleration
dy(2)= ( g/y(1) )*cos( y(2) ); %dy(2)=dbetadt variation in angle of inclination
end
%in a script i wrote
[t,y]=ode23s(@rocketequat,[0 160],[0 pi/2]);
% i gut this message:
In rocket (line 1)
Warning: Matrix is singular, close to singular or badly scaled. Results may be inaccurate. RCOND = NaN.
3 comentarios
Respuestas (1)
Stephan
el 30 de Dic. de 2018
Editada: Stephan
el 30 de Dic. de 2018
Hi,
you do a divide by zero at the initial point - use another x0 for y(1) - for example try:
[t,y]=ode45(@rocketequat,[0 160],[0.0001 pi/2]);
plot(t,y(:,1),t,y(:,2))
function dy = rocketequat(t,y)
%UNTITLED6 Summary of this function goes here
g=9.81; % constant of gravity m/s^2
Ve=2800;%celerity of inject gas m/s
m0=500000;%initial mass =500000kg
dmdt=-2300;%kg/s
dy(1)= -g.*sin(y(2)) - Ve .* dmdt./( m0 - dmdt*t );%dy(1)=dvdt acceleration
dy(2)= ( g./y(1) ).*cos( y(2) ); %dy(2)=dbetadt variation in angle of inclination
dy=dy';
end
gives:
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/199941/image.png)
There is also no need for the useage of a stiff solver - ode45 will do the job.
Best regards
Stephan
7 comentarios
Ver también
Categorías
Más información sobre Resizing and Reshaping Matrices en Help Center y File Exchange.
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!