Euler method RLC series circuit,will loop forever

2015 1 18
3 个回答
回答已采纳
更新时间：2015 1 18
1 次查看（30 天）

0 个投票

I am trying to integrate for the electrical current but when i run the script it keeps going forever an no plot is shown... Here is the script:
t(1)=0;
y=[0,0];
h=0;001;
tf=10;
n=(tf-t(1))/h;
for k=1:n
y(k+1,:)=y(k,:)+h*ord5(t(k),y(k,:));
t(k+1)=t(k)+h;
end
here is the ord5 function that needs to be integrated
function yd = ord4(t,y)
R=1,4;
L=0;
C=0;
Vs=5;
yd(1)=1/C*y(2);
yd(2)=(Vs-R*y(2)-y(1))/L;
end
Vs is the total circuit voltage and y(1)=Uc(the capacitor voltage) and y(2)=I(the circuit electric current).yd(1) is the derivative of Uc on time and yd(2) is the derivative of I on time.

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

 采纳的回答

Star Strider
Star Strider 2015-1-18

0 个投票

I haven’t run your code yet, but I see one problem straight away:
h=0;001;
I’m not certain what you intend, but here ‘h=0’, so ‘n=Inf’.
Might this be typo?
You may instead want:
h = 0.001;

更多回答（2 个）

John Etkin
John Etkin 2015-1-18

0 个投票

Thats a typing error but thats not the problem. Odd is when i run another version of Euler's method it works fine
clear all;close all;
x=[0 0];
h=0.01;
t0=0;
tf=10;
t(1)=0;
i=1;
for m=0.1:h:tf,
x(i+1,:)=x(i,:)+h* func(t,x(i,:));
t(i+1)=m;
i=i+1;
end
This one works fine,but i like the first one better cause its easier to memorize. Also,any reason why my plot is so small ?I have to give 10x the values that my teacher uses to get to see the plot else its just a small dot..
John Etkin
John Etkin 2015-1-18

0 个投票

Thanks ,found the error i used , instead of .

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

Star Strider
Star Strider 2015-1-18
My pleasure!
I didn’t see it in the code you posted, but so long as you found the error, great!
I did see:
R=1,4;
but since the comma is a legitimate decimal separator in other locations (not mine), I ignored it. (That assigns ‘R=1’ where I am.)

请先登录，再进行评论。

类别

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

标签

Community Treasure Hunt

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

Start Hunting!

Translated by