matlab ode45 求方程(1+x^2)y''=2xy',y(0)=1,y'(0)=3 的数值解,我的运行不行,我编的函数文件:function dy=sdre(x,y) dy=[y(2);2*x/(1+x^2)*y(2)]; 程序 [x,y]=ode45('sdre',[0 10],[1 3]);plot(x,y)

来源:学生作业帮助网 编辑:作业帮 时间:2024/05/06 14:06:07
matlab ode45 求方程(1+x^2)y''=2xy',y(0)=1,y'(0)=3 的数值解,我的运行不行,我编的函数文件:function dy=sdre(x,y) dy=[y(2);2*x/(1+x^2)*y(2)]; 程序 [x,y]=ode45('sdre',[0 10],[1 3]);plot(x,y)

matlab ode45 求方程(1+x^2)y''=2xy',y(0)=1,y'(0)=3 的数值解,我的运行不行,我编的函数文件:function dy=sdre(x,y) dy=[y(2);2*x/(1+x^2)*y(2)]; 程序 [x,y]=ode45('sdre',[0 10],[1 3]);plot(x,y)
matlab ode45 求方程(1+x^2)y''=2xy',y(0)=1,y'(0)=3 的数值解,我的运行不行,
我编的函数文件:function dy=sdre(x,y)
dy=[y(2);2*x/(1+x^2)*y(2)];
程序 [x,y]=ode45('sdre',[0 10],[1 3]);
plot(x,y)

matlab ode45 求方程(1+x^2)y''=2xy',y(0)=1,y'(0)=3 的数值解,我的运行不行,我编的函数文件:function dy=sdre(x,y) dy=[y(2);2*x/(1+x^2)*y(2)]; 程序 [x,y]=ode45('sdre',[0 10],[1 3]);plot(x,y)
1.建立M文件vdp1001a.m如下:function dy=vdp1001a(x,y)dy=zeros(2,1);dy(1)=y(2);dy(2)=2*x*y(2)/(1+x^2); 2.取x0=0,xf=10,输入命令:[X,Y]=ode45('vdp1001a',[0 10],[1 3]);plot(X,Y(:,1),'-')3.结果如图



你写的程序只有一处错误
ode45中要取函数句柄,
[x,y]=ode(@sdre,[0 10],[1 3]);
我已仿真成功。有问题再问我

.建立M文件vdp1001a.m如下:
function dy=vdp1001a(x,y)
dy=zeros(2,1);
dy(1)=y(2);
dy(2)=2*x*y(2)/(1+x^2);
2.取x0=0,xf=10,输入命令:

[X,Y]=ode45('vdp1001a',[0 10],[1 3]);
plot(X,Y(:,1),'-'

在matlab中用 ode45解方程 y'=x-y,y(0)=1,0 matlab ode45 求方程(1+x^2)y''=2xy',y(0)=1,y'(0)=3 的数值解,我的运行不行,我编的函数文件:function dy=sdre(x,y) dy=[y(2);2*x/(1+x^2)*y(2)]; 程序 [x,y]=ode45('sdre',[0 10],[1 3]);plot(x,y) 用matlab求y=f(x),y''=y^2的数值解,待解方程为:y=f(x),y''=y^2,即y的二阶导等于y的平方边界条件:y(0)=1,y'(100)=0用matlab函数ode45求解function dy=test(x,y)dy=zeros(2,1);dy(1)=y(2);dy(2)=[y(1)]^1;end[X,Y]=ode45('test',[0 100 matlab中ode45 及odeset的使用.在蝴蝶效应的解方程组中,方程是%x'=a*(X(2)-X(1))%y'=X(1)*(c-X(3))-X(2)%z'=X(1)*X(2)-b*X(3)使用ode45解方程.语句为lorenz=@(t,X) [a*(X(2)-X(1));X(1)*(c-X(3))-X(2);X(1)*X(2)-b*X(3)];[t X]=ode45(lorenz 求用Matlab的ode45编一个双线性方程X''+20X'+40πX=sin(t);X>=0X''+20X'+20πX=sin(t);X matlab用ode45 解微分方程如何用ODE45解微分方程 其中初值给出方程如下,请帮忙写出方程在程序中如何表达,我用function dy=func(x,y,r)dy = zeros(size(y));dy(1) = y(1)dy(2) = sqrt((1+r^2)/r^2*log((1+r^2)/(1+y(1)^2))+y(1 用matlab ode45求微分方程组 dx/dt+x+y=0 dy/dt+x-y=0 x(0)=0 y(0)=1 t=[0,0.5] matlab中,ode45能解决含有(dy/dt)^2和d2y/d2t的方程么 matlab ode45 解方程方程就是上面那样.m,c,k,f和φ可以单独赋值.求ω给一个值,x和t的关系,x一阶导数和t的关系.(它们的关系要画两幅图)然后,ω再给值,那两者之间的关系.这个有人会吗?谢谢! MATLAB 求微分方程数值解 y''-y'+y=0 y(1)=2,y(2)=0,x∈[0,20] MATLAB代码我的代码:funtion dy=fabc(x,y)dy=zeros(2,1);dy(1)=y(2);dy(2)=y(2)-y(1);[x,y]=ode45('fabc',[0 20],[2 0]);plot(x,y(:,1),'-') matlab解二阶微分方程组大致方程为方程1:y1+y2+y1'+y1=c方程2:y1+y2+y1'+y2'+y2=0有2阶导数和1阶导数混合的,y1和y2都是关于x的函数,初值确定,c是一个常数,是用ode45加龙格库塔法解吗?网上提到的 matlab解微分方程组,这个程序报错,求指导如何改~function x_dot=cha721lotka_volterra(t,x)x_dot=zeros(2,1);x=zeros(2,1);x_dot=[4*x(1)-2*x(1)*x(2);x(1)*x(2)-3*x(2)];t_final=100;x0=[2,3];[t,x]=ode45('cha721lotka_volterra',[0,t_final],x0);p matlab ode45 可以解二阶微分方程吗 如何用ode45解这个微分方程 这个方程的matlab程序又该如何写?ay-b(y'^2)+c*siny=0 自变量为x a,b,c均为正实数 matlab 求一元方程数值解请问如何用matlab求方程exp(-x)=1-x/2的数值解? Matlab 函数M文件 转化 匿名函数以及函数句柄例 ode45解常微分方程组一般解法:[t,x]=ode45(fun,[t0 td],x0);function dx=fun(t,x)dx(1)=2x(1)+x(2);dx(2)=-3x(1)-2x(2);end若写到一个M文件里,怎么表达? 跪求matlab大神指导!错在哪里 如何修改啊?>> figure;for s0=[0.9 0.7 0.5 0.3 0.1] x0=[1-s0,s0]; [t,x]=ode45(@sirfun,ts,x0); s=x(:,2);i=x(:,1); plot(s,i);hold on;endxlabel('S');ylabel('R')axis([0 1 0 1])plot([0 1],[1,0],'r')tit matlab怎么对微分方程求解知道是用ode45,方程是my ̈+cy ̇+ky=f