[討論] 利用數值求解簡諧振盪已回收
初學Matlab,對於寫程式語言仍然生疏
以下是我自己解的,請問錯誤在哪? 該如何修正?謝謝
<Ques.>
slove x"+wx=0
where w^2=k/m
for given k=1 and m=1
x(0)=0;
v(0)=0;
<sol>:
w=1;
x(1)=0;
v(1)=0;
t=10*pi;
n=256;
h=t/n;
x(2)=x(1)+v*h;
t=0:h:10*pi;
for i=1:255;
x(i+1)=x(i)+h;
x(i+2)=2*(x(i)+v*h)-(1+((w^2)*(h^2))x(i);
end
plot(t, x(i+2));
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 123.195.250.167
→
12/03 07:58, , 1F
12/03 07:58, 1F
討論串 (同標題文章)
完整討論串 (本文為第 1 之 2 篇):