Equations to Matlab code
2 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
Maria
le 18 Nov 2014
Réponse apportée : Warren Dennis
le 17 Fév 2020
hi
I want to convert the two equations below into Matlab code , I wrote the function below for that and I don't know if that is correct or not , specially the integral form ? any help ?
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/146144/image.png)
function [VI,Vw]= Modelfun(TL,TR)
dt=0.01; %sampling time
r=0.052;
L=0.306;
M=9;
I=5;
VI=(TL+TR)/ M * r;
VI=VI*dt/2; %first integral
Vw=(L*(TL-TR)) / 2 * r * I;
Vw=Vw*dt/2; %first integral
Vw=Vw*dt/2; %second integral
end
0 commentaires
Réponse acceptée
Torsten
le 20 Nov 2014
It is not possible simply to output (VI,Vw) given (TL,TR).
As input you need (VIold,Vwold,Vwdotold,dt,TL,TR).
The output VI,Vw and Vwdot is then calculated according to
VI=VIold+dt*(TL+TR)/(M*r)
Vw=Vwold+dt*Vwdotold
Vwdot=Vwdotold+dt*L*(TL-TR)/(2*r*I).
This would mean that your function "Modelfun" advances VI and Vw from their values at time t to time t+dt.
Best wishes
Torsten.
0 commentaires
Plus de réponses (2)
Torsten
le 18 Nov 2014
If the right-hand sides of your equations does not depend on t, the solution to your differential equations is explicitely given by
VI=VI0+(TL+TR)/(M*r)*t
Vw=Vw0+Vw0dot*t+0.5*L*(TL-TR)/(2*r*I)*t^2
where
VI0=VI at time t=0
Vw0=Vw at time t=0
Vw0dot = d/dt(Vw) at time t=0.
Best wishes
Torsten.
5 commentaires
Voir également
Catégories
En savoir plus sur Characters and Strings dans Help Center et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!