Re: [問題] 解多次微分方程式已回收
※ 引述《tarepanda275 (趴趴熊)》之銘言：
: 我想問一下 我有一個方程式為 m*D2u+c*Du^a+k*u=0
: 我在matlab寫
: eq='m*D2u+c*(Du)^a+k*u=0'
: int='u(0)=0,Du(0)=0'
: dsolve(eq,int,'t')
: 請問我還缺少啥，還是說要改成怎樣寫才求的出u??
: 感謝
: 另外 我若想求得對u微分兩次的數值
: 我的方程式假如為 20*D2u+2*Du+100*u=60
^^^^
這裡不是多次方?
: 在M檔案裡我寫成
: function du=mck(t,u)
: du(1)=u(2);
: du(2)=(-100*u(1)-2*u(2)+60)/20;
: 在主程式裡輸入
: 還要再輸入啥 才能得到微分兩次 也就是du(2)的值
: 感謝
---------------------------------------------
function pttex141
clc
[t1 u1] = ode45(@mck,[0 10],[0;0]);
[t2 u2] = ode15s(@mck,[0 10],[0;0]);
u11 = (-100.*u1(:,1)-2*u1(:,2).^2+60)/20;
u22 = (-100.*u2(:,1)-2*u2(:,2).^2+60)/20;
data1 = [t1 u1 u11]
data2 = [t2 u2 u22]
figure(1)
plot(t1,data1(:,2:4))
figure(2)
plot(t2,data2(:,2:4))
function du=mck(t,u)
du = zeros(2,1);
du(1)=u(2);
du(2)=(-100*u(1)-2*u(2).^2+60)/20;
--------------------------------------------
我假設一次微分項為二次方
--
1.MATLAB programming 2.ASPEN process simulation package
3.FORTRN programming 4.Advance Engineering Mathematics
5.Process Control Theory
6.Chemical Engineering Basic Theory(Kinetic.thermodynamics.transport)
7.Numerical Method and Analysis
8.MATLAB Toolbox.Simulink system basic design
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 113.160.148.219
討論串 (同標題文章)
MATLAB 近期熱門文章
PTT數位生活區 即時熱門文章
