怎样用matlab写出解微分方程x x=1 ,其中 x(0)=1的程序语句.
来源:学生作业帮助网 编辑:作业帮 时间:2024/07/03 08:01:46
[x,y]=dsolve('Dx=x-y-x*(x^2+y^2)','Dy=x+y-y*(x^2+y^2)','x(0)=2','y(0)=1')得到的结果是解析解没有找到.用数值解.在Matlab下
在Matlab下输入:edit,然后将下面两行百分号之间的内容,复制进去,保存%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
告诉我邮箱,发你代码
昨天不是回了一个吗?在Matlab下输入:edit,然后将下面两行百分号之间的内容,复制进去,保存%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
注意,12.19t这种写法不对,应该写成12.19*t.结果如下>>symsS1S2t>>[S2,S1]=dsolve('DS2=12.19*t-24209-0.3*S2-0.01*S2','DS1=
我改了你的一些内容,不知道有没有被我改错.在Matlab下输入:edit,然后将下面两行百分号之间的内容,复制进去,保存%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
functiondx=appollo(t,x)mu=1/82.45;mustar=1-mu;r1=sqrt((x(1)+mu)^2+x(3)^2);r2=sqrt((x(1)-mustar)^2+x(
把方程和初值条件都列出来这一大堆东西几个人看得懂?再问:公式(1)(2)就是方程啊,我写的那两个就是对应的方程。倒数第二张图片有实验数据。还需要什么参数吗?非常抱歉,我学的专业不太扎实,问题可能有点弱
常微分方程的解有标准格式,需要根据议程形式假设其通解,比如假设y=c1x+c2x'+c3,代入方程解通解,再加边界条件x(0)=y(0)=1解特解.
怎么说呢...可能是因为所求系数(整型)太大了,也可能是你初始量设置的有问题.具体这会儿我没参考书,要不然就告诉你怎么改进.那句"Error,(incombine/radical)integertoo
用solve命令就可以了,具体用法你查看一下MATLAB的帮助文档吧比如:solve('2*x+y-z=1','2*y+3*z=-2','x-z=-1');上面的方程是我胡乱写的,意思是有几个方程就往
ode23,ode45,ode113,ode15s,ode23s,ode23t,ode23tb这些函数都是的,用法都是一样的,但是针对不同的问题,最常用的是ode45例子:见ode45的帮助文档的后半
人家问的就是解析解,不知道上边几个在说什么你这个方程高版本的matlab可以解,只是算出来之后带积分项目:symsabQSTzZmcAPMeESS=dsolve('Dx=a*x','Dy=(b*Q*x
symsaxyy=int(x^2+a,x,0,a)然后运行solve('1/3*a^3+a^2=3','a');vpa(a)这结果出来了,三个.自己根据需要选择合适解.一般是第一个,后面俩个是复数解
1、你贴出来的报错信息和代码对不上号:前面显示错误的那行代码和你贴出来的完全不一样;而后面的错误(Inputargument'u1'isundefined)也不可能是目前的代码所导致的
f1=sym('x1+x01=sqrt(x2^2+y2^2)*cos(a+atan(y2/x2))+qrt(x01^2+y01^2)*cos(a+atan(y01/x01))');f2=sym('y1
这个好像是常微分方程组把,先建立一个调用式子的函数:functionf=hanshu(y,x)f=[a(x(3)-x(2));b(x(1)-0.5*sqrt(x(2)-c));d(2.5-x(3))]
[xy]=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',
角度用sind、cosd,例如>>sind(30)ans=0.5000%真有意思,用不同ode,得到的结果不同~!图形见:http://hi.baidu.com/chemical%5Fengineer
这下子就可以了:functionhh[T,Y]=ode45(@rigid,[012],[011])plot(T,Y)functiondy=rigid(t,y)dy=zeros(3,1);dy(1)=y