谁有MATLAB求解偏微分方程的模板程序吗,急求
发布网友
发布时间：2024-09-08 21:37
共1个回答
热心网友
时间：2024-11-15 05:35
ode45 是解微分方程的函数，你在帮助里的说明书点ode45就能查到。
我这给你写一段，假设有微分方程
你要新建一个函数 写出这个微分方程，然后在主程序里调用这个函数，然后加入初始值，就是0 1 1。 然后用ode解。
这是定义的函数的样子：
function dy = rigid(t,y)
dy = zeros(3,1); % a column vector
dy(1) = y(2) * y(3);
dy(2) = -y(1) * y(3);
dy(3) = -0.51 * y(1) * y(2);
然后在主函数里，也就是另一个M文件里，解微分方程，
options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-5]);
[T,Y] = ode45(@rigid,[0 12],[0 1 1],options);
plot(T,Y(:,1),'-',T,Y(:,2),'-.',T,Y(:,3),'.')；%这个是绘制图形
最后那句就是画图，出来的图形就是y的值，y1 y2 y3 .如下图：
希望解决了你的疑惑。谢谢。
