实验十.欧拉法MATLAB实现-创新互联

实 验 步 骤、过 程

1.

求出其特解为

(1)欧拉法的向前、向后、改进欧拉法

编制函数文件:

function [y1,y2,y3,y4]=su(a,b,k)
% a起始点，b迭代终止点，h=(b-a)/k
h=(b-a)/k;x=a:h:b;[m,n]=size(x);
y1=zeros(m,n);y2=zeros(m,n);y3=zeros(m,n);
y1(1)=1;y2(1)=1;y3(1)=1;
for i=1:n-1 %向前欧拉法
    x1=x(i);x2=x(i+1);
    y1(i+1)=y1(i)+h*((3*(x1).^2-2)*y1(i)/(x1).^3+1);
end
for i=1:n-1 %向后欧拉法
    x1=x(i);x2=x(i+1);
    yy=y2(i)+h*((3*(x1).^2-2)*y2(i)/(x1).^3+1);
    y2(i+1)=y2(i)+h*((3*(x1)^2-2)*yy/(x1).^3+1);
end
for i=1:n-1 %改进欧拉法
    x1=x(i);x2=x(i+1);
    z1=(3*(x1).^2-2)*y1(i)/(x1).^3+1;
    yy=y3(i)+h*z1;
    z2=(3*(x1).^2-2)*yy/(x1).^3+1;
    y3(i+1)=y3(i)+h*0.5*(z1+z2);
end
for i=1:n %精确值
    y4(i)=((x(i))^3*(1-exp((x(i))^(-2)-1)))/2;
end

命令窗口下执行：

>>[y1,y2,y3,y4]=su(pi,2*pi,5)

(2)二阶、四阶龙格-库塔方法

题目可转化为：

编制函数文件:

function f=su(x,y)
f=(3*x^2-2)*y/(x^3)+1

命令窗口下执行：

>>[x,y]=ode23('su',[pi,2*pi],1);[x',y']
>>[x,y]=ode45('su',[pi,2*pi],1);[x',y']

答案供参考,如有错误请自行改正，本人只是随意分享之前的作业，代码均来源于网络+个人稍微修改，认为有用可参考哈～

你是否还在寻找稳定的海外服务器提供商？创新互联www.cdcxhl.cn海外机房具备T级流量清洗系统配攻击溯源，准确流量调度确保服务器高可用性，企业级服务器适合批量采购，新人活动首月15元起，快前往官网查看详情吧