matlab解二阶微分方程,解出的结果如下,请问怎样画出y和时间t的函数图像?

y=dsolve('101*D2y+100*Dy+1000*y=102*sin(20*pi*t)','y(0)=0,Dy(0)=0','t')
解出的结果是
y=exp(-(50*t)/101)*sin((10*985^(1/2)*t)/101)*((51*985^(1/2)*exp((50*t)/101)*((50*sin(t*(20*pi - (10*985^(1/2))/101)))/101 - cos(t*(20*pi - (10*985^(1/2))/101))*(20*pi - (10*985^(1/2))/101)))/(9850*((20*pi - (10*985^(1/2))/101)^2 + 2500/10201)) + (51*985^(1/2)*exp((50*t)/101)*((50*sin(t*(20*pi + (10*985^(1/2))/101)))/101 - cos(t*(20*pi + (10*985^(1/2))/101))*(20*pi + (10*985^(1/2))/101)))/(9850*((20*pi + (10*985^(1/2))/101)^2 + 2500/10201))) - (51*cos((10*985^(1/2)*t)/101)*(99485*pi^2*sin(20*pi*t
- (10*985^(1/2)*t)/101) - (4925*sin(20*pi*t + (10*985^(1/2)*t)/101))/2 - (4925*sin(20*pi*t - (10*985^(1/2)*t)/101))/2 + 99485*pi^2*sin(20*pi*t + (10*985^(1/2)*t)/101) + 4925*pi*cos(20*pi*t - (10*985^(1/2)*t)/101) + 4925*pi*cos(20*pi*t + (10*985^(1/2)*t)/101) + (25*985^(1/2)*cos(20*pi*t - (10*985^(1/2)*t)/101))/2 - (25*985^(1/2)*cos(20*pi*t +

继续输入:
ezplot(y,[0,2]) % t∈[0,2]时的图像.
就出来了. 想要得到某一点的y的值, 输入
subs(y,'t',1) 就得到t=1时的y的值.追问

敢问ezplot是什么意思啊?这个公式没见过啊。。

追答

ezplot是画符号表达式的函数,输入付号表达式会自动画图。subs是带入求值,令变量取某值,得出符号表达式的结果。

温馨提示:内容为网友见解,仅供参考
无其他回答

matlab 求解二阶微分方程并画出y与t函数曲线
[t,x] = ode45(@(t,x)dy(t, x, m, alpha, beta, U, k, p, A),[0,1e-3],[0,0]);subplot(1,2,1)plot(t,x(:,1))xlabel('$$t$$','interpreter','latex','fontsize',16);ylabel('$$y$$','interpreter','latex','fontsize',16)subplot(1,2,2)plot(t,x(:,2)...

如何用matlab求二阶微分方程的图像解
tspan=[0,10];[t,y]=ode45(@ode_fun,tspan,y0)3、利用plot函数绘出二阶微分方程的图像解 polt(t,y)xlabel('t'),ylabel('y(t)','dy\/dt');4、执行结果

matlab:利用Simulink求解二阶微分方程?
1、打开Simulink的模块库 (1)选择Commonly User Blocks打开常用模块插入Integrator模块,插入Mux模块,插入Out1模块;(2)选择User-defined Functions打开自定义函数模块插入Fcn模块 2、根据运算流程,画连接线 3、在Integrator模块中,设置y'(0),y(0)的初始值 4、在Fcn模块中,输入微分方程的关系式,...

如何用matlab求解二阶微分方程,以及程序实例
2、在matlab命令行窗口中输入symsx[x,params,conds]=solve(sin(x)==1,'ReturnConditions', true),按回车键可以得到方程解。3、转换一下,可以看到sin(x)=1方程的解是如下图所示 。4、也可以求解下面的一个方程。5、输入symsabcyx;[x,y]=solve([a*x^2+b*y+c==0,a*x+2*y==4],[...

matlab如何求解二阶微分方程数值解的图象
题主的问题可以使用ode45函数求解,得到【x,y,y'】的数值解。其求解方法:1、确定初值,y0=[pi,0];2、确定x的范围,xspan=[0,10];3、自定义函数,odefun=@(x,y)[y(2);sin(y(1))*cos(y(1))*((2\/9*pi)\/(2+0.4*(sin(y(1)))^2))^2-5*sin(y(1))];4、求解,[x,...

二阶微分方程组的MATLAB解法
过程:先降阶为一阶微分方程组 y'=z z'=-z-y+sin(t)编制如下函数m文件 function dy=weifen(t,x)dy=zeros(2,1);y=x(1)z=x(2)dy(1)=x(2);dy(2)=sin(t)-x(2)-x(1);然后用ode45解方程 [t,y]=ode45(@weifen,[0 20],[5 6])plot(t,y)就可以画出y和y'的图像 ...

matlab二阶微分方程求解及画图
参考代码:m=10600; k=31986497; c=8.9E4;dy = @(t,y)[y(2); (k*y(1)-c*y(2)\/abs(y(2))*y(2)^2)\/m];ode45(dy,[0 5],[0 3.66])其中,蓝线为x,绿线为x'。

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

如何用matlab求解下述的二阶变系数常微分方程,并且画出图像
>> y=(C1*exp(gamma^(1\/2)*x)+C2*exp(-gamma^(1\/2)*x))\/x;>> p=diff(y);>> simpliyf(diff(p)+p*2\/x-gamma*y)ans = 0 表明确实是方程的解.然后求解初始条件:>> p_x_equ_0=subs(p,x,0)错误使用 symengine Division by zero.提示错误,因为x=0时解无意义,因此需要在...

matlab 二阶微分方程 求解
1、建立微分方程的自定义函数,odefun(容x,y)2、当a=10、θ=π\/6时,执行下列代码 theta=pi\/6;[x,y]= ode45(@odefun,[0,0.18],[0.5,theta])plot(x,y),grid on legend('y (x)','y’(x)')xlabel('x'),ylabel('y (x),y’(x)')figure(2)plot(y(:,1),y(:,2))...

相似回答