Equations to Matlab code
1 次查看(过去 30 天)
显示 更早的评论
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 ?
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 个评论
采纳的回答
Torsten
2014-11-20
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 个评论
更多回答(2 个)
Torsten
2014-11-18
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.
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Predictive Maintenance Toolbox 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!