作者zupo (外派越南工作者)
标题Re: [讨论] 非线性微分方程
时间Thu Nov 28 23:54:11 2013
※ 引述《lachese (不想抓老鼠的猫)》之铭言：
: 在下为MATLAB新手
: 想请问各位高手
: 关於解非线性的问题
: 如果我的方程式包含exp的t次方
: 那该怎麽解呢?
: 我有试过ODE，可是对於exp的t次方项显示描述错误??
: 我的非线性如下：
: x1'=-10*x1+x2*exp^.3t
: x2'=-x2
: 是说这样的方程式不能用ODE的解法吗?
: 谢谢各位的指教!
本题用 dsolve 好像还是解不出解析解
毕竟MATLAB在解数学函数解析解领域不强
所以我用数值解来解
------------------------------------------------
function pttex139
[t x] = ode45(@nonlfun,[0 2],[1;1]);
plot(t,x(:,1),t,x(:,2))
legend('x1','x2')
title('联立方程式解')
xlabel('t')
ylabel('Value')
grid
function f = nonlfun(t,x)
f = zeros(2,1);
f(1) = -10.*x(1)+x(2).*exp(3.*t);
f(2) = -x(2);
------------------------------------------------
本题有零解 x1 = 0 , x2 = 0;
故初始值不选用以上条件
另外本题受两个初始值正负号不同
趋势线会有较大的变化
可以自己试看看
--
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.163.21.205