Matlab 常微分方程数值解法

矩阵与数值分析实验 用Matlab编写的 常微分方程数值解法 程序

function [ output_args ] = Untitled1（ input_args ）
%UNtitleD1 Summary of this function goes here
%  Detailed explanation goes here
clc;
clear;
h=0.1;
a=0;
b=1;
ya=1;
N=（b-a）/h;
T=linspace（abN+1）;
Y=zeros（1N+1）;
Y（1）=ya;

%Euler法
for j=1:N
    Y（j+1）=Y（j）+h*（Y（j）-T（j）*Y（j）^2）;
end
plot（TY‘g‘）;grid on;
hold on;

%Euler改进法
for j=1:N
   Y（j+1）=Y（j）+（h/2）*（Y（j）-T（j）*Y（j）^2+  （Y（j）+h*（Y（j）-T（j）*Y（j）^2））-T（j+1）*（Y（j）+h*（Y（j）-T（j）*Y（j）^2））^2  ）;
end
plot（TY‘r‘）;
hold on;

%精确解
Y=1./（T-1+2*exp（-T））;
plot（TY）;
xlabel（‘x‘）;
ylabel（‘f（x）‘）;
title（‘Euler（绿色）、Euler改进法（红色）及精确解（蓝色）的比较‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        611  2010-01-03 22:41  常微分方程数值解法\常微分方程数值解法.m

     目录          0  2010-01-09 23:06  常微分方程数值解法

----------- ---------  ---------- -----  ----

                  611                    2