(2) 设方程式为y’=-40y,y(0)=2用欧拉数值方法与ode45求方程式的数值解 请大神用MATLAB编写啊
来源:学生作业帮 编辑:作业帮 分类:综合作业 时间:2024/10/06 09:41:30
(2) 设方程式为y’=-40y,y(0)=2用欧拉数值方法与ode45求方程式的数值解 请大神用MATLAB编写啊
欧拉方法
第一步:建立并保存名为 funfcn1.m 的 m 文件函数;
function f=funfcn1(x,y)
f=-40*y;
第二步:建立并保存名为Qeuler1.m的m文件函数;
function[X,Y,n]=Qeuler1(funfcn,x0,y0,b,h)
x=x0;n=(b-x0)/h;X=zeros(n,1);y=y0;Y=zeros(n,1);k=1;X(k)=x;Y(k)=y;
for k=2:n+1
fxy=feval(funfcn,x,y);
x=x+h;y=y+h*fxy;
X(k)=x;Y(k)=y';
end
plot(X,Y,'rp');grid,xlabel('自变量X'),ylabel('因变量Y')
title('用向前欧拉公式计算dy/dx=f(x,y),y(x0)=y0在[x0,b]上的数值解')
第三步:建立如下程序并保存名为jisuan1_1.m的m文件
clc,clear
x0=0;y0=3;b=2;h1=0.1;%赋初始值
[x1,Y1 ,n1]=QEuler1(@funfcn1,x0,y0,b,h1);%调用函数
第一步:建立并保存名为 funfcn1.m 的 m 文件函数;
function f=funfcn1(x,y)
f=-40*y;
第二步:建立并保存名为Qeuler1.m的m文件函数;
function[X,Y,n]=Qeuler1(funfcn,x0,y0,b,h)
x=x0;n=(b-x0)/h;X=zeros(n,1);y=y0;Y=zeros(n,1);k=1;X(k)=x;Y(k)=y;
for k=2:n+1
fxy=feval(funfcn,x,y);
x=x+h;y=y+h*fxy;
X(k)=x;Y(k)=y';
end
plot(X,Y,'rp');grid,xlabel('自变量X'),ylabel('因变量Y')
title('用向前欧拉公式计算dy/dx=f(x,y),y(x0)=y0在[x0,b]上的数值解')
第三步:建立如下程序并保存名为jisuan1_1.m的m文件
clc,clear
x0=0;y0=3;b=2;h1=0.1;%赋初始值
[x1,Y1 ,n1]=QEuler1(@funfcn1,x0,y0,b,h1);%调用函数
(2) 设方程式为y’=-40y,y(0)=2用欧拉数值方法与ode45求方程式的数值解 请大神用MATLAB编写啊
matlab ode45 求方程(1+x^2)y''=2xy',y(0)=1,y'(0)=3 的数值解,我的运行不行,
MATLAB求方程y’= y - 2x/y,y(0) = 1 (0≤x≤1,h = 0.1) 的数值解,要求编写程序
用matlab求微分方程 y’=-y2+x,y(0)=1的数值解,
用MATLAB解方程(y+2)y'=x,当x=0:1:100时,求y的数值,并画出y的图形,怎么输入命令?
matlab 求线性方程的数值特解,并画出图形:y'''+8y'=0,y(0)=1,y'(0)=1y''())=2,x[
MATLAB 求微分方程数值解 y''-y'+y=0 y(1)=2,y(2)=0,x∈[0,20] MATLAB代码
matlab 数值解比如有方程 y^3+x^2+3x-2y=1 其中x的区间为0到10 如何让matlab求解数值解y呢
用matlab求这个方程的数值解:y''-(1-y^2)y'+y=0,y(1)=2,y(2)=0,[0,20],麻烦把程
求大神支招,我想用ode45来求解一个微分方程的数值解,但是MATLAB老是说有错,怎么办?
y=1/(x^4 + x^3 + x^2 + x + 1)用matlab怎么求(1,2)区间的定积分?要数值解
MATLAB中求Y=F(X)数值积分