matlab求解微分方程
答:第一种方法:利用dsolve函数求微分方程的符号解(通解):对于一些不是很难,要求出通解的微分方程,用dsolve函数求解。1、 打开Matlab软件-->点击新建脚本菜单,新建一个脚本文件用于编写微分方程求解程序。2、 输入微分方程求解程序-->点击保存-->点击运行。3、在matlab的命令窗口即可看到求解结果,是...
答:第一步、分析微分方程的特点确定模块 我们将y"+4y'+4y=cos2t方程改写成y"=cos2t-4y'-4y,从右边可知,式子是由cos2t外加载信号,一阶导数项和一次变量组成。所以我们,可以用下列功能模块来搭建其仿真模型 cos2t:使用Fcn模块和Clock模块;y':使用lntegrator模块;y':使用lntegrator模块;-4:...
答:用matlab求解微分方程可以按下列思路来求解。1、用dsolve函数,求解该微分方程的通解,T表达式 T=dsolve(diff(T)==-k*(T-21.1));2、求当t=0时的T表达式,eq1=C1/10 + 211/10 3、求当t=1时的T表达式,eq2=(C1*exp(-k))/10 + 211/10 4、用solve函数,求解由2和3组成的方程组的解...
答:1. 准备工作 在开始求解微分方程之前,需要先安装MATLAB软件,并掌握一些基本的MATLAB语法和数学知识。2. 定义微分方程 在MATLAB中,可以使用syms命令定义符号变量,并使用diff命令对变量进行求导。例如,假设要求解微分方程dy/dx=x^2,可以使用以下代码:syms y(x)eqn = diff(y,x) == x^2;3. 求...
答:1、首先得介绍一下,在matlab中解常微分方程有两种方法,一种是符号解法,另一种是数值解法。在本科阶段的微分数学题,基本上可以通过符号解法解决。2、用matlab解决常微分问题的符号解法的关键命令是dslove命令。该命令中可以用D表示微分符号,其中D2表示二阶微分,D3表示三阶微分,以此类推。值得注意的...
答:如何用matlab求解微分方程并画图,可以先用dsolve()或ode()求出其微分方程(组)的解析解或数值解,然后用plot()绘制其图形。例如:解微分方程 y'=y-2t/y,y(0)=1,0<t<4 1、用dsolve()求解,代码及结果如下 2、用ode45()求解,代码及结果如下 3、当然喽,使用dsolve()或ode(...
答:如何用matlab求解微分方程的数值解和解析解?一、微分方程的数值解可以ode函数来求解。其求解方法:1、自定义微分方程组函数,odefun(x,y)2、定义x【0,10】间若干等份的数值,如n=50 3、定义y的初值,即 y0=[2,7];4、使用ode45函数求y(x),z(x)数值解,即 [x,y]=ode45(@odefun,x,y0)...
答:1、找到关于求解常微分方程的习题。2、这里用matlab求解,主要用到的就是“dsolve”,具体的格式如下,但是不仅仅局限下面两种,其中eq就是代表方程式,而con为初始条件。3、这里的“x”关于x的函数。知道语法就可以进行计算了。4、在matlab中输入如下,dsolve('D2y+2*Dy+exp(x)=0','x'),并直接...
答:1、首先看一下matlab求解方程的方法,首先指明所解方程的变量,然后指明方程,未知数和限制条件,最后求解方程。2、在matlab命令行窗口中输入symsx[x,params,conds]=solve(sin(x)==1,'ReturnConditions', true),按回车键可以得到方程解。3、转换一下,可以看到sin(x)=1方程的解是如下图所示 。4、...
答:MATLAB提供了7个常微分方程求解器(solver),分别是ode45,ode23,ode113,ode15s,ode23s,ode23t,ode23tb,其中前3个适用于求解非刚性(Nonstiff)问题,后4个适用于刚性问题。所谓刚性问题,简单点说,就是系统包含多个相互作用但变化速度相差十分悬殊的子过程。这些求解器并不是可以互相取代的,它们分别...
网友评论:
钮帘17279817643:
matlab 求解微分方程 -
19926蔺姣
: 编个M文件,命名成rigid: %x(1)指原微分方程中的x,x(2)指x' function dx = rigid(t,x) dx=zeros(2,1);%两行 dx(1)=x(2); dx(2)=-0.06*x(1)+0.35*sin(2*t+0.07); %令w=2,必须是个给定常数在命令窗口里输入命令: [T,Y]=ode45(@rigid,[0,500],[0.03,0]); plot(T,Y(:,1)) 就输出x的图像
钮帘17279817643:
matlab怎么对微分方程求解 -
19926蔺姣
: 这个不用,如果f是数,直接 dsolve('m*D2y+c*Dy+k*y=f') 就可以了 如果f是f(t) 需要编写程序用ode45 首先编写m-file%%%%%%%%%%%%%function dotx=myfun(t,y) dotx=[y(2);1/m*(-c*y(2)-k*y(1)+f(t))];%f(t)要写出具体形式%%%%%%%%%%%%%%%M-file 以上要单存一个文件在工作路径下%%%%%%%%%%%%% 在command 窗口 [t,y]=ode45('myfun',[t_0,t_final],[y_0,y'_0])%[t_0,t_final]和[y_0,y'_0]分别问求解区间和初值 再用plot画图
钮帘17279817643:
MATLAB 微分方程组求解 -
19926蔺姣
: 除了有参数,还要有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;%时间区...
钮帘17279817643:
如何用matlab求解微分方程并画图 -
19926蔺姣
: ^如何用matlab求解程并画图 dsolve()求解微分方程后,使用ezplot()作图,参考代码:1 2>> y = dsolve('Dy=y*cos(x)/(1+2*y^2)','y(0)=1','x'); >> ezplot(y,[-10,10]) % [-10,10]为x区间
钮帘17279817643:
如何使用MATLAB求解微分方程 -
19926蔺姣
: 可以用desolve方程:r=dsolve('eqn1','eqn2',...,'cond1','cond2',...,'var').eqn1表示第1个微分方程,cond1表示第1个初始条件,var表示微分方程中的自变量,默认为t.
钮帘17279817643:
matlab分段微分方程组求解 -
19926蔺姣
: 用matlab求解分段微分方程组,可以这样来考虑. 1、建立分段函数,ode_fun1(t,y);ode_fun2(t,y);ode_fun3(t,y) 2、分三次用ode45求解,其【t,x,y】的数值 [t1,y1]=ode45(@ode_fun1,tspan,X0) [t2,y2]=ode45(@ode_fun2,tspan,X0) [t3,y3]=ode45(@ode_fun3,tspan,X0) 3、最后用plot函数绘出微分方程组解的图形.plot(x1,y1,x2,y2,x3,y3) 运行结果如下图.
钮帘17279817643:
如何利用matlab求解微分方程的表达式 -
19926蔺姣
: dsolve()求解微分方程后,使用ezplot()作图,参考代码: >> y = dsolve('Dy=y*cos(x)/(1+2*y^2)','y(0)=1','x');>> ezplot(y,[-10,10]) % [-10,10]为x区间
钮帘17279817643:
如何用MATLAB软件求微分方程?
19926蔺姣
: 有符号法和数值法.一般使用dsolve和ode命令,比如: 符号法 dsolve('Df = f + sin(t)') ans = C4*exp(t) - sin(t)/2 - cos(t)/2 数值法 sol = ode45(@vdp1,[0 20],[2 0]);x = linspace(0,20,100);y = deval(sol,x,1); plot(x,y);
钮帘17279817643:
MATLAB中怎么解微分方程 -
19926蔺姣
: 好像没什么特别好的办法: >> syms t q y; q=4*t^3;>> syms Dy>> eq=Dy-(q-2*t);>> y=dsolve(char(subs(eq)),'y(0)=5','t') y = t^4-t^2+5 也就是要先移项把微分方程写成等号右端为零的形式(需额外定义Dy),然后subs代入,再求解.
钮帘17279817643:
matlab求解微分方程的通解问题解某微分方程:x^2*D2y+x*Dy+(x^2 - 1/2)*y=0,初值:y(pi/2)=2,Dy(pi/2)= - 2/pi,书本标准答案为:ans=2^(1/2)*pi*^(1/2)/x^(1/2)*sin... -
19926蔺姣
:[答案] syms x y y=sqrt(2*pi/x)*sin(x);y1=diff(y,'x');y2=diff(y,'x',2); h=x^2*y2+x*y1+(x^2-1/2)*y;h=simplify(h) 结果: h = -1/4/x*2^(1/2)*sin(x)*pi^(1/2)/(1/x)^(1/2) 可见不等于0,=y,说明原答案有问题