matlab求微分的命令
车炉力748如何使用MATLAB求解微分方程 -
潘官琰19352535348 ______ 可以用desolve方程:r=dsolve('eqn1','eqn2',...,'cond1','cond2',...,'var').eqn1表示第1个微分方程,cond1表示第1个初始条件,var表示微分方程中的自变量,默认为t.
车炉力748matlab 矩阵间的微分如何求? -
潘官琰19352535348 ______ 这实际上是求Jacobian矩阵 设有向量 v(t)=[v1(t),v2(t),……,vm(t)]' 则有向量函数 f(v)=[f1(v),f2(v),……,fn(v)]' 则有f(v)的obian矩阵: Jacobian=[df1/dv1…df1/dvm … … dfn/dv1…dfn/dvm] 在MATLAB中用函数命令jacobian()来实现.
车炉力748如何用matlab求解微分方程并画图? 求解的话 是输入deslove这我知道,那么作图呢? -
潘官琰19352535348 ______ 求解微分方程的函数是dsolve().如要画微分方程的解函数,可以用plot()函数.例如:求解dy/dx=1+y^2,并且当y(0)=1的解,画出y=f(x)的图形. y=dsolve('Dy=1+y^2','y(0)=1','x') x=-2*pi:pi/10:2*pi; y1=eval(y); plot(x,y1,'k-') xlabel('x'),ylabel('y(x)');grid on
车炉力748matlab分段微分方程组求解 -
潘官琰19352535348 ______ 用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) 运行结果如下图.
车炉力748如何用matlab求解微分方程组的数值解 -
潘官琰19352535348 ______ 在Matlab下输入:edit zhidao_feiying.m,然后将下面两行百分号之间的内容,复制进去,保存 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function y=zhidao_feiying(t,x) a=1; b=2; c=3; d=4; ...
车炉力748如何运用matlab求解非线性偏微分方程组 -
潘官琰19352535348 ______ 比如:y''-y'+y=x+1 y是关于x的函数 命令如下: syms x t;dsolve('D2x-Dx+x=t+1') 结果:ans = t + C5*exp(t/2)*cos((3^(1/2)*t)/2) + C6*exp(t/2)*sin((3^(1/2)*t)/2) + 2 因为系统默认的变量为t 所以那样编程 你可以百度的
车炉力748如何用matlab求解微分方程组 -
潘官琰19352535348 ______ 这说明此微分方程没有显示解,应用数值方法求解,例如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),'.') 结果
车炉力748用matlab命令求微分方程d^2y/dx^2+2*dy/dx+2y=o,满足初始条件y(0)=1,dy/dx(0)=0的解,并绘制出方程解y(t)的时间曲线图 -
潘官琰19352535348 ______[答案] y=dsolve('D2y+2*Dy+2*y=0','y(0)=1,Dy(0)=0'); t=0:0.01:20; y=eval(y); plot(t,y) 祝你学习愉快!
车炉力748下列是MATLAB中求函数零点的命令的是 - 上学吧普法考试
潘官琰19352535348 ______ 用数值解法吧.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...