用matlab求微分方程 y’=-y2+x,y(0)=1的数值解,
来源:学生作业帮 编辑:作业帮 分类:综合作业 时间:2024/05/14 04:21:48
用matlab求微分方程 y’=-y2+x,y(0)=1的数值解,
function [x,y]=Classical_RK4(odefun,xspan,y0,h,varargin)
x=xspan(1):h:xspan(2);
y(1)=y0;
for k=1:length(x)-1
K1=feval(odefun,x(k),y(k),varargin{:});
K2=feval(odefun,x(k)+h/2,y(k)+h/2*K1,varargin{:});
K3=feval(odefun,x(k)+h/2,y(k)+h/2*K2,varargin{:});
K4=feval(odefun,x(k)+h,y(k)+h*K3,varargin{:});
y(k+1)=y(k)+h/6*(K1+K2+K3+K4);
end
x=x';y=y';
采用经典4阶龙格库塔方法计算,编写function函数,调用如下:
f=@(x,y)-y^2+x;
[x,y]=Classical_RK4(f,[0,5],1,0.1);
plot(x,y,'--*')
x=xspan(1):h:xspan(2);
y(1)=y0;
for k=1:length(x)-1
K1=feval(odefun,x(k),y(k),varargin{:});
K2=feval(odefun,x(k)+h/2,y(k)+h/2*K1,varargin{:});
K3=feval(odefun,x(k)+h/2,y(k)+h/2*K2,varargin{:});
K4=feval(odefun,x(k)+h,y(k)+h*K3,varargin{:});
y(k+1)=y(k)+h/6*(K1+K2+K3+K4);
end
x=x';y=y';
采用经典4阶龙格库塔方法计算,编写function函数,调用如下:
f=@(x,y)-y^2+x;
[x,y]=Classical_RK4(f,[0,5],1,0.1);
plot(x,y,'--*')
用matlab求微分方程 y’=-y2+x,y(0)=1的数值解,
MATLAB 求微分方程数值解 y''-y'+y=0 y(1)=2,y(2)=0,x∈[0,20] MATLAB代码
请问各位高手如何使用matlab解微分方程 y''(x)-y(x)sinx=0,y(0)=1,y'(0)=0;用数值解法
如何使用matlab用梯形法解常微分方程dy/dx=y/x-y^2/x,y(1)=0.5,求出其数值解
用MATLAB解方程(y+2)y'=x,当x=0:1:100时,求y的数值,并画出y的图形,怎么输入命令?
高手MATLAB 求微分方程的解 y''+4*y'+4*y=e^-2x
matlab 求线性方程的数值特解,并画出图形:y'''+8y'=0,y(0)=1,y'(0)=1y''())=2,x[
MATLAB求方程y’= y - 2x/y,y(0) = 1 (0≤x≤1,h = 0.1) 的数值解,要求编写程序
x*y'*sin(y/x)-y*sin(y/x)+x=0 求微分方程的解
如何处理matlab微分方程数值解,具体就是解微分方程组,得出一系列的(x,y,t)如何作sqrt(x^2+y^2)=z
MATLAB:求微分方程y*y'/5+x/4=0 的通解,并求任意常数为1时的解!
matlab解微分方程y“(x)+2y'(x)+2y(x)=xe^(-x),y(0)=0,y‘(0)=0,求过程及结果