there are two coupled differential equation, I try to solve that differential for the i = 1 to 5, but i think my for loop is incorrect.
조회 수: 3 (최근 30일)
이전 댓글 표시
these are the differential equation I want to solve where the i = 1 to 5
I want to solve these differential equation but maybe the loop isn't correct or we can't apply loop here?
clc
ti = 0;
tf = 10E-6;
tspan = [ti tf];
a = 0.6;
n = 0.05;
tc = 10E-9;
r = 1.5;
F = 1;
k =1:1:5;
f = zeros(length(k),1) ;
for i = 1:length(k)
y(1) = A(i)
y(2) = O(i)
f = @(t,y) [
(-1/(2*tc)).*(1 - r/(1+ (A(i)^2)./F)).*A(i) + (n/(2*tc)).*((cos(O(i+1)- O(i))).*A(i+1) + (cos(O(i+1)- O(i))).*A(i-1)) ;
(a./(2*tc)).*(r/(1 + (A(i)^2)./F)) + (n/(2*tc)).*((A(i+1)./A(i)).*sin(O(i+1) - O(i)) - (A(i-1)./A(i)).*sin(O(i-1) - O(i)));
O(i+1) - O(i)
];
end
[time,Y] = ode45(f,tspan,[0;0;0]);
plot(time,Y(:,3))
댓글 수: 5
답변 (1개)
Walter Roberson
2022년 7월 31일
Note that because you use A(i+1) and O(i-1) you cannot produce plots for the first or last entries in A or O.
%number of data points
N = 5;
%put in your real A and O data here!!
A = randn(1,N);
O = rand(1,N) * 10;
%process data
ti = 0;
tf = 10E-6;
tspan = [ti tf];
a = 0.6;
n = 0.05;
tc = 10E-9;
r = 1.5;
F = 1;
k = 1:N;
f = zeros(N,1) ;
for i = 2:N-1
f = @(t,y) [
(-1/(2*tc)).*(1 - r/(1+ (A(i)^2)./F)).*A(i) + (n/(2*tc)).*((cos(O(i+1)- O(i))).*A(i+1) + (cos(O(i+1)- O(i))).*A(i-1)) ;
(a./(2*tc)).*(r/(1 + (A(i)^2)./F)) + (n/(2*tc)).*((A(i+1)./A(i)).*sin(O(i+1) - O(i)) - (A(i-1)./A(i)).*sin(O(i-1) - O(i)));
O(i+1) - O(i)
];
[time,Y] = ode45(f,tspan,[0;0;0]);
plot(time, Y(:,3), 'DisplayName', "i = " + i);
hold on
end
legend show
댓글 수: 5
Torsten
2022년 7월 31일
Usually, there are special "boundary" ODEs for the cases i=1 and i=N that differ from the ODEs for 1<i<N.
Walter Roberson
2022년 7월 31일
Note that once you have the definitions for the boundary ode, my recommendation would be to use the symbolic toolbox to set up the entire system of 15 coupled equations at the same time. Then follow the first example in the document for odeFunction() to convert to something that can be used with ode45
참고 항목
카테고리
Help Center 및 File Exchange에서 Ordinary Differential Equations에 대해 자세히 알아보기
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!