clzm.net
当前位置:首页 >> mAtlAB解微分方程 >>

mAtlAB解微分方程

dsolve('D3y+D2y/x-Dy/(x^2)=x','y(500)=0,Dy(500)=0','x'),用这个能够得到y的一个含参数的解,然后根据x=0时,y为有限值来确定参数的值,之后就能解出相应的y’’ 和y’ ,后面那个貌似用MATLAB不太好实现,最好是自己来分析。

你没给初值,我随便射了 clear allclc f=@(t,y)([y(2);(3*y(1)^2+cos(y(1)))/(3*y(2)^2)]);[t,Y]=ode45(f,[0 10],[1 1]);plot(t,Y(:,1),t,Y(:,2)),legend('X','dX/dt')

dsolve('Dg=k1*r/V-k2*g','g(0)=g0') 结果: ans = -(exp(-k2*t)*(k1*r - V*g0*k2) - k1*r)/(V*k2)

可以解。但这种非线性微分方程能求出解析解的可能性微乎其微,一般用数值法求解。 有以下几个问题需要明确: 1、需要初值,包括x(0)和x'(0)。 2、需要知道常数ξ的值。 3、按照惯例,一般以t或x作为自由变量,y作为因变量,像这个微分方程这样dx/...

ezplot(y)

syms x(t) y(t) z(t) [x,y,z]=dsolve(diff(x)==4*x+3*y+2*z,diff(y)==6*x-3*y+2*z,diff(z)==6*x+5*y+4*z) 按上面试试

如果是常微分方程,可以用dsolve函数。该函数可以解单变量常微分方程或者多变量常微分方程组,所以5个变量也不在话下。 调用格式如下: [y1,...,yN] = dsolve(eqns) solves the system of ordinary differential equations eqns and assigns the...

可这样: >> syms Vx Vy Vz b m y1(t) y2(t) y3(t) y4(t) y5(t) y6(t) y1=dsolve(diff(y1)==Vx) y2=dsolve(diff(y2)==-b/m*(Vx^2+Vy^2+Vz^2)*Vx) y3=dsolve(diff(y3)==Vy) y4=dsolve(diff(y4)==-b/m*(Vx^2+Vy^2+Vz^2)*Vy) y5=dsolve(diff(y5)==V...

参考代码: % 参数定义m1 = 1; m2 = 2; m3 = 3;k1 = 4; k2 = 5; k3 = 6;c1 = 7; c2 = 8; c3 = 9;Qy = 10; % 微分方程K = [k1 -k1 0; -k1 k1+k2 -k2; 0 -k2 k2+k3];C = [c1 -c1 0; -c1 c1+c2 -c2; 0 -c2 c2+c3]dy = @(t,y)[y(4:6); (-K*y(1:3)-C*...

在Matlab下输入:edit zhidao_feiying.m,然后将下面两行百分号之间的内容,复制进去,保存 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function y=zhidao_feiying(t,x) a=1; b=2; c=3; d=4; %%比方说 %f=a*x+y-z; %g=b*sin(x*y)-c*cos...

网站首页 | 网站地图
All rights reserved Powered by www.clzm.net
copyright ©right 2010-2021。
内容来自网络,如有侵犯请联系客服。zhit325@qq.com