lzth.net
当前位置:首页 >> mAtlAB中求微分方程组通解 >>

mAtlAB中求微分方程组通解

说明 求解常微分方程一般使用ode系列函数(其中最常用的是ode45),基本做法就是把微分方程用函数表示出来,然后调用ode函数求解即可。对于本题而言,使用匿名函数描述微分方程更为方便一些。 参考代码 k0=8; k1=1; K=1; k2=5; n=2; % 常数定义d...

可这样: >> 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...

这个方程要解出来还要初始条件吧 不妨设(x(0),y(0))=(A,B); 还有用ode45只能得到数值解,不能给出解析式表达 那么先建立一个函数文件 test_fun.m function dx=tsst_fun(x,t) dx(1)= c-b*x(1)-w*x(1)*x(2); dx(2)=w*x(1)*x(2)-(b+m)*x(2); 那么调...

没解析解你肯定得用数值解。以下是步骤。 第一步,在Matlab里写m-file,命名微分方程对应的函数: function f=myfun(t,x) f = [-x(1)+x(2)^2; -2*x(2)+x(1)^2]; 这里的x(1)就是你这儿的x, x(2)就是你这儿的y。写完后自动保存,不要对保存目录做...

在Matlab下输入:edit,然后将下面两行百分号之间的内容,复制进去,保存 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function dxdt=zhidao_rk4_5(t,x) global a b c d r x1=x(1); x2=x(2); x3=x(3); x4=x(4); dxdt=[ a*(x2-x1)+x4; d*x1-x...

>> syms x(t) y(t) r d a b r='1';d='0.5';a='0.1';b='0.02'; [x,y]=dsolve(diff(x)==(r-a*y)*x,diff(y)==-(d-b*x)*y,x(0)==25,y(0)==2) Warning: Explicit solution could not be found. > In dsolve at 194 x = [ empty sym ] y = [] 说明无解。

这说明此微分方程没有显示解,应用数值方法求解,例如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);d...

在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...

可以借助于嵌套函数或匿名函数实现附加参数的传递,例如 function main y0 = [1.4; 0.1; 0.1]; A = linspace(eps, 10, 20); Y = A * NaN; for ii = length(A) a = A(ii); y = ode45(@eq2, [0 a], y0); Y(ii) = y(end, 1); end plot(A, Y) functi...

对于不是很复杂的微分方程(最好不包含逻辑判断之类的),可以使用匿名函数(7.0之前使用内联函数inline)定义微分方程,然后用ode*系列函数求解。 例如,作为ode*函数文档中的例子: function dy = vdp1000(t,y) dy = zeros(2,1); % a column v...

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