matlab+微分方程组
@冯差1292：matlab中怎样写微分方程组 -
楚俗15281661366…… 以下列例子,来说明怎样写微分方程组.解 1、建立m-文件rigid.m如下: function dy=rigid(t,y) dy=zeros(3,1); dy(1)=y(2)*y(3); dy(2)=-y(1)*y(3); dy(3)=-0.51*y(1)*y(2);2、取t0=0,tf=12,输入命令: [T,Y]=ode45('rigid',[0 12],[0 1 1]); plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+')3、结果如图 图中,y1的图形为实线,y2的图形为“*”线,y3的图形为“+”线.
@冯差1292：如何用matlab求解微分方程组 -
楚俗15281661366…… 这说明此微分方程没有显示解,应用数值方法求解,例如ode45等函数 就一个例子(选自matlab帮助文档) 解入下微分方程组1,建立方程组函数 function dy = rigid(t,y)dy = zeros(3,1); % a column vectordy(1) = y(2) * y(3);dy(2) = -y(1) * y(3);dy(3) = -0.51 * y(1) * y(2);2,求解并绘图 [T,Y] = ode45(@rigid,[0 12],[0 1 1]);plot(T,Y(:,1),'-',T,Y(:,2),'-.',T,Y(:,3),'.') 结果
@冯差1292：如何用matlab求解一个二阶常系数微分方程组 -
楚俗15281661366…… 用命令:dsolve('s','s1','s2',…,'x') 其中s 为方程s1,s1,s3,…为初始条件x 为自变量方程s 中用d 表示求导 数d2,d3,…表示二阶三阶等高阶导数初始条件缺省时给出带任意常数 c1,c2,..的通解自变量缺省值为t 也可求解微分方程组 例 1、dsolve('dy=1+y^...
@冯差1292：怎样用matlab求解复杂的微分方程组
楚俗15281661366…… 用数值解法吧.matlab程序为: function hhh clear;clc; [TOUT,X] = ode23t(@odefun,[0,0.01],[1500 -500 0]) %t=TOUT,vx=X(:,1),vy=X(:,2),vz=X(:,3) plotyy(TOUT,X(:,1:2),TOUT,X(:,3)) function yhat=odefun(t,x) vx=x(1);vy=x(2);vz=x(3); dx=-(vx+1.4)*((vx+1...
@冯差1292：如何用matlab求解微分方程组的数值解
楚俗15281661366…… 在Matlab下输入:edit zhidao_feiying.m,然后将下面两行百分号之间的内容,复制进去,保存 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function y=zhidao_feiying(t,x) a=1; b=2; c=3; d=4; ...
@冯差1292：MATLAB 微分方程组求解 -
楚俗15281661366…… 除了有参数,还要有x,y,x',y'的初始值 看你的方程,类似有阻力的抛体运动 这里随便给了个初值,x(0)=0 y(0)=10 x'(0)=10 y(0)'=0 可以理解为物体从10米高的高度,以10米每秒的速度做平抛运动 m=1;k=0.2;g=10;%参数 ts=0:0.01:2;%时间区...
@冯差1292：matlab求解微分方程组的程序 -
楚俗15281661366…… >> syms x(t) y(t) r d a b r='1';d='0.5';a='0.1';b='0.02'; [x,y]=dsolve(diff(x)==(r-a*y)*x,diff(y)==-(d-b*x)*y,x(0)==25,y(0)==2) Warning: Explicit solution could not be found. > In dsolve at 194 x = [ empty sym ] y = [] 说明无解.
@冯差1292：matlab求解微分方程组,画出解函数图.x'= - x^3 - y,x(0)=1 y'=x - y^3,y(0)=0.5 0 - 作业帮
楚俗15281661366…… [答案] function dz=caonim(t,z) %保存为 caonim.m 文件 dz=[-z(1)^3-z(2);z(1)-z(2)^3]; ------------------------------------ 运行 [T,Z]=ode45('caonim',[0 30],[1,0.5]); plot(T,Z) legend('x-t','y-t') ------------------------
@冯差1292：用MATLAB编程求解微分方程组 -
楚俗15281661366…… dY=inline('[y(2);(y(1)*y(4)^2+9.8*(cos(y(3)-1)))/2;y(4);-(9.8*sin(y(3))+2*y(2)*y(4))/y(1)]','t','y');ode45(dY,1,[.5 0 0 .5])...
@冯差1292：MATLAB求解微分方程组dy(1)/ds=1.2*(0.16*(1 - 1*cosy(3)+0.5*1*siny(3)));dy(2)/ds=1.2*(0.16*(1 - 1*cosy(3)+0.5*1*siny(3)))*1*cosy(3)+(y(1)*y(4)*9.8*siny(3))/(300... - 作业帮
楚俗15281661366…… [答案] dy=inline('[1.2*(0.16*(1-1*cos(y(3))+0.5*1*sin(y(3))));1.2*(0.16*(1-1*cos(y(3))+0.5*1*sin(y(3))))*1*cos(y(3))+(y(1)*y(4)*9.8*sin(y(3)))/(300*y(2));-(1.2*(0.16*(1-1*cos(y(3))+0.5*1*sin(y(3))))*1*sin(y(...
楚俗15281661366…… 以下列例子,来说明怎样写微分方程组.解 1、建立m-文件rigid.m如下: function dy=rigid(t,y) dy=zeros(3,1); dy(1)=y(2)*y(3); dy(2)=-y(1)*y(3); dy(3)=-0.51*y(1)*y(2);2、取t0=0,tf=12,输入命令: [T,Y]=ode45('rigid',[0 12],[0 1 1]); plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+')3、结果如图 图中,y1的图形为实线,y2的图形为“*”线,y3的图形为“+”线.
@冯差1292：如何用matlab求解微分方程组 -
楚俗15281661366…… 这说明此微分方程没有显示解,应用数值方法求解,例如ode45等函数 就一个例子(选自matlab帮助文档) 解入下微分方程组1,建立方程组函数 function dy = rigid(t,y)dy = zeros(3,1); % a column vectordy(1) = y(2) * y(3);dy(2) = -y(1) * y(3);dy(3) = -0.51 * y(1) * y(2);2,求解并绘图 [T,Y] = ode45(@rigid,[0 12],[0 1 1]);plot(T,Y(:,1),'-',T,Y(:,2),'-.',T,Y(:,3),'.') 结果
@冯差1292：如何用matlab求解一个二阶常系数微分方程组 -
楚俗15281661366…… 用命令:dsolve('s','s1','s2',…,'x') 其中s 为方程s1,s1,s3,…为初始条件x 为自变量方程s 中用d 表示求导 数d2,d3,…表示二阶三阶等高阶导数初始条件缺省时给出带任意常数 c1,c2,..的通解自变量缺省值为t 也可求解微分方程组 例 1、dsolve('dy=1+y^...
@冯差1292：怎样用matlab求解复杂的微分方程组
楚俗15281661366…… 用数值解法吧.matlab程序为: function hhh clear;clc; [TOUT,X] = ode23t(@odefun,[0,0.01],[1500 -500 0]) %t=TOUT,vx=X(:,1),vy=X(:,2),vz=X(:,3) plotyy(TOUT,X(:,1:2),TOUT,X(:,3)) function yhat=odefun(t,x) vx=x(1);vy=x(2);vz=x(3); dx=-(vx+1.4)*((vx+1...
@冯差1292：如何用matlab求解微分方程组的数值解
楚俗15281661366…… 在Matlab下输入:edit zhidao_feiying.m,然后将下面两行百分号之间的内容,复制进去,保存 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function y=zhidao_feiying(t,x) a=1; b=2; c=3; d=4; ...
@冯差1292：MATLAB 微分方程组求解 -
楚俗15281661366…… 除了有参数,还要有x,y,x',y'的初始值 看你的方程,类似有阻力的抛体运动 这里随便给了个初值,x(0)=0 y(0)=10 x'(0)=10 y(0)'=0 可以理解为物体从10米高的高度,以10米每秒的速度做平抛运动 m=1;k=0.2;g=10;%参数 ts=0:0.01:2;%时间区...
@冯差1292：matlab求解微分方程组的程序 -
楚俗15281661366…… >> syms x(t) y(t) r d a b r='1';d='0.5';a='0.1';b='0.02'; [x,y]=dsolve(diff(x)==(r-a*y)*x,diff(y)==-(d-b*x)*y,x(0)==25,y(0)==2) Warning: Explicit solution could not be found. > In dsolve at 194 x = [ empty sym ] y = [] 说明无解.
@冯差1292：matlab求解微分方程组,画出解函数图.x'= - x^3 - y,x(0)=1 y'=x - y^3,y(0)=0.5 0 - 作业帮
楚俗15281661366…… [答案] function dz=caonim(t,z) %保存为 caonim.m 文件 dz=[-z(1)^3-z(2);z(1)-z(2)^3]; ------------------------------------ 运行 [T,Z]=ode45('caonim',[0 30],[1,0.5]); plot(T,Z) legend('x-t','y-t') ------------------------
@冯差1292：用MATLAB编程求解微分方程组 -
楚俗15281661366…… dY=inline('[y(2);(y(1)*y(4)^2+9.8*(cos(y(3)-1)))/2;y(4);-(9.8*sin(y(3))+2*y(2)*y(4))/y(1)]','t','y');ode45(dY,1,[.5 0 0 .5])...
@冯差1292：MATLAB求解微分方程组dy(1)/ds=1.2*(0.16*(1 - 1*cosy(3)+0.5*1*siny(3)));dy(2)/ds=1.2*(0.16*(1 - 1*cosy(3)+0.5*1*siny(3)))*1*cosy(3)+(y(1)*y(4)*9.8*siny(3))/(300... - 作业帮
楚俗15281661366…… [答案] dy=inline('[1.2*(0.16*(1-1*cos(y(3))+0.5*1*sin(y(3))));1.2*(0.16*(1-1*cos(y(3))+0.5*1*sin(y(3))))*1*cos(y(3))+(y(1)*y(4)*9.8*sin(y(3)))/(300*y(2));-(1.2*(0.16*(1-1*cos(y(3))+0.5*1*sin(y(3))))*1*sin(y(...
