首页 >>  正文

matlab求微分方程

来源:baiyundou.net   日期:2024-09-20

祁鸿到1046matlab求解微分方程的通解问题解某微分方程: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... -
时爸枯15176475325 ______[答案] 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,说明原答案有问题

祁鸿到1046怎么用matlab解微分方程 dy/dx=A - By - Cxy -
时爸枯15176475325 ______[答案] syms A B C x y; >> dsolve('Dy=A-B*y-C*x*y','x') ans = -1/2*A*erf(1/2*(-2*C)^(1/2)*x-B/(-2*C)^(1/2))*exp(-1/2*x*(2*B+C*x)-1/2*B^2/C)/C*2^(1/2)*(-pi*C)^(1/2)+exp(-1/2*x*(2*B+C*x))*C1

祁鸿到1046matlab微分方程编程求解 -
时爸枯15176475325 ______ 直接解就可以了>> S=dsolve('Dx+5*x+y=exp(t)','Dy-x-3*y=0','x(0)=1','y(0)=0') S = x: [1x1 sym] y: [1x1 sym]>> S.x,S.y ans = -4*exp((-1+15^(1/2))*t)*(1/22+13/330*15^(1/2))+exp((-1+15^(1/2))*t)*(1/22+13/330*15^(1/2))*15^(1/2)-4*exp(-(1+15^(1/2))*t)*(-...

祁鸿到1046用MATLAB求微分方程,并画出图形解y=dsolve('D3y=D2y - Dy - y+t^2','y(0)=0','Dy(0)=1','D2y(0)= - 1')这样解了之后怎么画出图形啊! -
时爸枯15176475325 ______[答案] y=dsolve('D3y=D2y-Dy-y+t^2','y(0)=0','Dy(0)=1','D2y(0)=-1'); tt=linspace(0,10,1001); %表示t从0到10 for i=1:1001 t=tt(i); yy(i)=eval(y); end plot(tt,yy)

祁鸿到1046如何运用matlab求解非线性偏微分方程组 -
时爸枯15176475325 ______ 比如: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 所以那样编程 你可以百度的

祁鸿到1046matlab求解微分方程组,画出解函数图.x'= - x^3 - y,x(0)=1 y'=x - y^3,y(0)=0.5 0 -
时爸枯15176475325 ______[答案] 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') ------------------------

祁鸿到1046matlab怎么求微分方程的变量表达式 -
时爸枯15176475325 ______ 对于一般的微分方程用matlab的dsolve()函数命令求得其解析值(即函数表达式).例如1、求 的解析值.syms y(t) a Dy = diff(y); D2y = diff(y,2); dsolve(D2y == -a^2*y, y(0) == 1, Dy(pi/a) == 0) 运行结果 y=exp(-a*t*1i)/2 + exp(a*t*1i)/2 式中li(x)=logint(x) 对数积分2、求 的解析值.syms x(s) a x = dsolve(diff(x) == -a*x, x(0) == 1) 运行结果 x = exp(-a*s)

祁鸿到1046怎样用matlab求解复杂的微分方程组
时爸枯15176475325 ______ 用数值解法吧.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...

祁鸿到1046matlab里面求微分方程dsolve('Df=exp( - i*s*t)*f^2+sin(f)','t')为什么解不出来? -
时爸枯15176475325 ______[答案] 因为它其实不存在解析解(除了平凡解0以外),而Matlab一直在试图找解.

祁鸿到1046matlab中符号法求常微分方程
时爸枯15176475325 ______ 把自变量x变为t >> y=dsolve('(t^2+1)*Dy-(y^2+1)=0') y =tan(atan(t)+C1) 即 y =tan(atan(x)+C1) 给10分吧,好吗?看看我的回答 http://passport.baidu.com/?center&tpl=ik&aid=7&default_tab=4#4,0 总分:16464 ;回答总数: 682;题均: 16464/682=24.1408

(编辑:自媒体)
关于我们 | 客户服务 | 服务条款 | 联系我们 | 免责声明 | 网站地图 @ 白云都 2024