Solution to second order differential equation

Vin
2013 10 18
1 个回答
回答已采纳
更新时间：2014 2 5
4 次查看（30 天）

0 个投票

I want to solve a second order differential equation which passes parameters to the ode solver. The parameters are column arrays and hence I used for loop in the function. I don't know if what I wrote is right or wrong. I am getting the solution. But I have no idea if it is right or wrong..
Is the code within the function right? I mean the syntax...
Thankyou...
The function is given below
function r = f1(tspan,y2,I,J,K,L,W)
r = zeros(2,1);
c = zeros(length(tspan));
r(1) = y2(2);
for n =2:length(tspan)
c(n) = -(I(n).*y2(2) + J(n).*y2(1) + K(n) + L(n).*cos(W*tspan(n)));
end
r(2)=c;

请先登录，再回答此问题。

 采纳的回答

Azzi Abdelmalek
Azzi Abdelmalek 2013-10-18

1 个投票

function r = f1(t,y2,tspan1,I,J,K,L,W)
r = zeros(2,1);
c = zeros(length(tspan1));
r(1) = y2(2);
for n =2:length(tspan1)
c(n) = -(I(n).*y2(2) + J(n).*y2(1) + K(n) + L(n).*cos(W*tspan1(n)));
end
r(2)=c;
% How to call your function
% Define Tspan1, Ic and your parameters
[T Y] = ode45(@(t,y2) f1(t,y2,Tspan1,I,J,K,L,W),Tspan1,IC)

3 个评论
显示 1更早的评论 隐藏 1更早的评论

Vin
Vin 2013-10-18
编辑：Vin 2013-10-18
This is my function call ( from part of my code)
for n = store2:numel(t)
store1=n
m=n;
if (theta(n)<thetaD)
tspan = store2:0.1:n-1;
init2 = [y1(store2) y1dot(store2)];
[tspan,y2] = ode45(@(tspan,y2) f1(tspan,y2,I,J,K,L,W),tspan,init2);
c2 = vertcat(c2,y2(:,1));
break
end
end
Azzi Abdelmalek
Azzi Abdelmalek 2013-10-18
Look at the above answer for how to write the function f1, and call it like below:
[t,y2] = ode45(@(t,y2) f1(t,y2,tspan,I,J,K,L,W),tspan,init2);
Vin
Vin 2014-2-5
I think I got it Thanks

请先登录，再进行评论。

更多回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 Programming 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by