二、关于求解变参数微分方程 看到有不少人问过如何处理"变参数微分方程", 所以抽了一点时间,写了一个例子, 希望能帮助到那些需要求解此类问题的人。 %%%================================%%% clear all fun=inline('[y(2);sin(w*t)-2*y(1)-3*y(2)]','t','y','flag','w'); tsp=[0 10]; y0=[1 1]; xlim(tsp) for w=1:10 [t,y]=ode45(fun,tsp,y0,[],w); plot(t,y) title(['w = ',num2str(w)]); pause(0.5); end %%%================================%%% 讨论、回复专贴:“关于求解变参数微分方程”回复专贴 [ 本帖最后由 eight 于 2007-6-1 00:15 编辑 ] |
GMT+8, 2025-4-12 02:45 , Processed in 0.034947 second(s), 16 queries , Gzip On.
Powered by Discuz! X3.4
Copyright © 2001-2021, Tencent Cloud.