matlab怎么求解微分方程

作者&投稿:侨琼 (若有异议请与网页底部的电邮联系)

如何用matlab求解微分方程并画图
如何用matlab求解微分方程并画图,可以先用dsolve()或ode()求出其微分方程(组)的解析解或数值解,然后用plot()绘制其图形。例如:解微分方程 y'=y-2t\/y,y(0)=1,0<t<4 1、用dsolve()求解,代码及结果如下 2、用ode45()求解,代码及结果如下 3、当然喽,使用dsolve()或ode(...

matlab如何求解线性方程组?
拓展:1、MATLAB 是美国MathWorks公司出品的商业数学软件,用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境,主要包括MATLAB和Simulink两大部分。2、MATLAB是matrix&laboratory两个词的组合,意为矩阵工厂(矩阵实验室)。是由美国mathworks公司发布的主要面对科学计算、可视化以及交互...

如何用matlab求解微分方程?
分解的原则就是按照运算顺序一层一层的去掉,比如y=sin e^√x,他的最后一步运算时sin,所以先去sin,就设h(x)=sin(x),去掉sin后的最后一步运算是求指数e^√x,所以设g(x)=e^x,接着,去掉e后的最后一步运算是开方√x,所以设f(x)=√x,至此可以把h,g,f复合起来就是原函数了.第二题同理...

请问老师这个题用matlab怎么做?
题主给出的二阶微分方程组,用matlab求解其精确解和数值解,可以使用dsolve()函数得到其解析值,使用ode()函数得到其数值解。一、使用dsolve()函数求解 第一步,对x(t)、y(t)进行声明变量,即 syms x(t) y(t)第二步,使用diff()函数分别对x(t)、y(t)进行导数计算 Dy=diff(y,t,1...

matlab怎么计算方程组
1、对于比较简单的方程组,可以用solve()函数命令求解。如方程组 x + y = 1 ; x - 11y = 5 >>[x,y]=solve('x + y = 1','x - 11*y = 5')又如方程组 exp(x+1)-y²=10 ;ln(x)+3y=7 >>syms x y >>[x,y]=solve(exp(x+1)-y^2-10,log(x)+3*y-...

matlab中怎么求解一个式子中的某一个元素的值
2、如果我们想获取矩阵第2行第1列的值,输入a(2,1)。3、如果想获取矩阵第2行第3列的值,输入a(2,3)。4、如果想间隔的获取矩阵多个值,可以使用a(2,[13])方式,这是获取矩阵第2行第1列,第3列的值。5、如果我们想获取第2行的数据,可以使用a(2,:)。以上就是matlab中求解一个式子...

如何用Matlab求解如下一元四次方程,求指导
如何用Matlab求解如下一元四次方程,求指导 不带参数:solve('x^4-x^3+x^2=0')单引号内式子可以任意改变,但形式要与例子一致. 带参数:syms a b c x; solve('a*x^4-b*x^3+c*x^2=0',x),要解变量a就改为solve('a*x^4-b*x^3+c*x^2=0',a)一元四次方程求解 试根:...

如何用matlab求方程的所有整数解
如何用matlab求方程的所有整数解?其解决方案,可以通过下列途径解决:方法一、利用三循环语句和判断语句,求解其所有整数解 方法二、整数线性规划分支定界法,求解其所有整数解 下面按第一种方法求解,其过程为 m=50;n=50;l=50;disp('满足条件:2x+3y+4z=100的解')for i=1:m for j=1:n for...

使用matlab求解方程组
使用matlab求解线性方程组,可以这样解。首先,写出线性方程组的系数。即 A=[1 -9 -10;-9 1 -5;8 7 1];其二,写出线性方程组的常数项系数。即 B=[1;0;4];然后,用矩阵左除法,求出X=[x1,x2,x3] 的解。即 X=A\\B 运行结果 ...

如何利用matlab求微分方程通解
对于一般的微分方程,可以用dsolve()函数,直接求出其微分方程通解。例如:求下列微分方程的通解 求解代码:syms y(t) a %变量声明 eqn = diff(y,t,2) == a*y; %定义方程ySol(t) = dsolve(eqn) %方程求解 求解过程与结果 ...

隐晨14731514039问: matlab怎么对微分方程求解 -
海门市琥珀回答: 这个不用,如果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画图

隐晨14731514039问: matlab 求解微分方程 -
海门市琥珀回答: 编个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的图像

隐晨14731514039问: 如何使用MATLAB求解微分方程 -
海门市琥珀回答: 可以用desolve方程:r=dsolve('eqn1','eqn2',...,'cond1','cond2',...,'var').eqn1表示第1个微分方程,cond1表示第1个初始条件,var表示微分方程中的自变量,默认为t.

隐晨14731514039问: 怎样用matlab求解复杂的微分方程组 -
海门市琥珀回答: [x y]=dsolve('D2x=(x*x+3*x)*Dy','D2y=-(x*x+3*x)*Dx','x(0)=0','y(0)=0','Dx(0)=a','Dy(0)=a','D2x(0)=a','D2y(0)=a')%%%%%%%%%%%%%%%%%%%%%%%%%%% 设定的X(1)=dx/dt; x(2)=x; x(3)=dy/dt; x(4)=y;%%%%%%% function dx=shuzhijie(t,x) ...

隐晨14731514039问: 如何用matlab求解微分方程并画图 -
海门市琥珀回答: ^如何用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区间

隐晨14731514039问: 如何用MATLAB软件求微分方程? -
海门市琥珀回答: 微分方程的话可以用dsolve函数,格式可以总结为dsolve('equation','condition',自变量)的形式.自变量默认为t.这种方法解微分方程还好.如果是微分方程组的话建议用od...

隐晨14731514039问: 如何用MATLAB软件求微分方程?
海门市琥珀回答: 有符号法和数值法.一般使用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);

隐晨14731514039问: MATLAB中怎么解微分方程 -
海门市琥珀回答: 好像没什么特别好的办法: >> 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代入,再求解.

隐晨14731514039问: 如何利用matlab求解微分方程的表达式 -
海门市琥珀回答: dsolve()求解微分方程后,使用ezplot()作图,参考代码: >> y = dsolve('Dy=y*cos(x)/(1+2*y^2)','y(0)=1','x');>> ezplot(y,[-10,10]) % [-10,10]为x区间

隐晨14731514039问: 用MATLAB如何解微分方程,高手来救救
海门市琥珀回答: dsolve('D2y/Dx2=Dy/Dx-3*x^2+6*x') dsolve('(1+x)*D2y/Dx2+Dy/Dx=0') matlab里 一个一个复制进去,回车就行 第一个解 ans = [ x(t)=-1+(1-2*Int(diff(y(t),t)*diff(x2(t),t)/diff(y(t),$(t,2)),t)-2*C1)^(1/2), x(t)=-1-(1-2*Int(diff(y(t),t)*diff(x2(t),t)/diff(y(t),$(t,2)),t)-2*C1)^(1/2)] 第二个有些问题,自己看看吧


本站内容来自于网友发表,不代表本站立场,仅表示其个人看法,不对其真实性、正确性、有效性作任何的担保
相关事宜请发邮件给我们
© 星空见康网