Re: [問題] 微分方程是隱函數已回收
※ 引述《frank751012 (肉丸)》之銘言：
: dsolve('80*D2x+9110*x+(3763.45/(0.2-x))*x-cos(4*pi*t)=0','t')
: 請問要怎麼用matlab解隱函數
: 我用dsolve語法解
: 他只有顯示Warning: Explicit solution could not be found.
: 麻煩各位大大
: 謝謝~
先來個數值解吧
------------------------------程式內容----------------------------
function pttex121
[t x] = ode45(@noldeq,[0 1],[0;0])
plot(t,x)
legend('x','x''')
xlabel('t')
ylabel('x')
function f = noldeq(t,x)
f = [0;0];
f(1) = x(2);
f(2) = -9110/80.*x(1)-(3763.45./80./(0.2-x(1)))+cos(4.*pi.*t)./80;
----------------------------------------------------------------------
解為頻率函數.應該為三角函數.
不行就請出 MAPLE 或是 MATHEMATICA 大神來幫忙吧
--
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: 124.9.132.173
討論串 (同標題文章)
本文引述了以下文章的的內容：
完整討論串 (本文為第 2 之 2 篇)：
MATLAB 近期熱門文章
PTT數位生活區 即時熱門文章
