用恒定系数求解一阶颂
2次观看(最近30天)
显示较旧的评论
我如何修复y = sym('y(t)),因为它总是显示错误并且无法正确运行程序。
%%以恒定系数求解一阶ODE
符号T S y
y = syms('y(t)');%写作y作为t的函数
y1 = diff(y,t);y(t)的符号形式的衍生物%
a =输入(“在方程式中输入DY/DT的系数”);
b =输入(“在方程式中输入Y的系数”);
f =输入(“输入方程的不孔子部分”);
y0 =输入(“初始条件为y(0) - ”);
eqnt = a*y1 + b*y -f;eqnt形式的%微分方程 - 0
eqns = laplace(eqnt,s);%采取DE的拉普拉斯变换
eqns = subs(eqns,{'拉普拉斯(y(t),t,s)',,,,'y(0)'},{y,y0});
y =简化(solve(eqns,y));%求解代数方程式eqns -0用于y和简化
y = ilaplace(y,s,t);%采取y的逆向转换y
disp([[“微分方程的解决方案是:y(t)”,char(y)])
0评论
答案(1)
Bjorn Gustavsson
2021年8月1日
编辑:Bjorn Gustavsson
2021年8月1日
应该像以下内容一样简单:
符号y(t)
另外:不是直接使用
DSOLVE
够了吗?
Hth