※ 引述《newlymoon (NEWLYMOON)》之銘言:
: dx1=x2
: dx2=(-3+cost)x1
: 請問怎求X1 X2
: 解析解找不出來 求數值解
回文版上討論即可.我回你信箱是提醒你我回你了
這個問題 當初始值 X1 X2 皆為 0 時 為奇異解 (零解)
系統的解會隨時間上下震盪.
================================= 程式 ====================================
function pttex25
[t x]=ode15s(@pttex,[0 10],[1;1])
plot(t,x)
function f=pttex(t,x)
f=zeros(2,1);
f(1)=x(2);
f(2)=(-3+cos(t)).*x(1);
===========================================================================
你可以試看看不同的初始值會有什麼結果 改 [1;1] 那兩個值
--
Basic -> Advance -> Master -> Expert -> Great-Master
Skill: AspenPlus AspenDynamic MATLAB Numerical-Method
Level: Master Advance Expert Master
-----------......什麼時候升級呢? 囧> --------------
Fortran Tecplot Sigmaplot Corldraw
Basic Basic Basic Basic
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.118.160.163