二阶常微分方程联立解有y1 y2 两变数,并且有两组二阶常微分方程(对x微分),初始值已知,y1'' = -a*exp(y1)*y1' - b*y1 + c*y2y2'' = -d*y2 + e*y1 - f*sin(wx)其中a,b,c,d,e,f,w都是已知常数,且y1(0)=c1,y1'(0)=0,y2(0)=c2,y2
来源:学生作业帮助网 编辑:六六作业网 时间:2024/11/28 01:22:34
二阶常微分方程联立解有y1 y2 两变数,并且有两组二阶常微分方程(对x微分),初始值已知,y1'' = -a*exp(y1)*y1' - b*y1 + c*y2y2'' = -d*y2 + e*y1 - f*sin(wx)其中a,b,c,d,e,f,w都是已知常数,且y1(0)=c1,y1'(0)=0,y2(0)=c2,y2
二阶常微分方程联立解
有y1 y2 两变数,并且有两组二阶常微分方程(对x微分),初始值已知,
y1'' = -a*exp(y1)*y1' - b*y1 + c*y2
y2'' = -d*y2 + e*y1 - f*sin(wx)
其中a,b,c,d,e,f,w都是已知常数,且y1(0)=c1,y1'(0)=0,y2(0)=c2,y2'(0)=0
c1,c2也都是已知常数,求y1,y2的解
图片中最后一项sin(wt)应改正为sin(wx)
解完再增加100的奖励!
很感谢诸位用MATLAB帮在下解题,然而先前已有使用MATLAB跑波形了,
主要是希望有常微分的高手,能解出CROSS-FORM,也就是Y1与Y2的通解!
另外y2'(0)的初值不是0了,改为c3,真抱歉!
二阶常微分方程联立解有y1 y2 两变数,并且有两组二阶常微分方程(对x微分),初始值已知,y1'' = -a*exp(y1)*y1' - b*y1 + c*y2y2'' = -d*y2 + e*y1 - f*sin(wx)其中a,b,c,d,e,f,w都是已知常数,且y1(0)=c1,y1'(0)=0,y2(0)=c2,y2
用ode45解决
题目中的所有已知常数我都用1代替.即a=b=c=d=e=f=w=c1=c2=c3=1.
在x从1到10之间求解
编写m文件odex2.m储存
function yy=odex2(x,z)
yy=zeros(4,1)
yy=[z(2);-exp(z(1))*z(2)-z(1)+z(3);z(4);-z(3)+z(1)-sin(x)]
在matlab命令窗口中输入
[x,z]=ode45('odex2',[0,10],[1;0;1;1])
plot(x,z(:,1),'r',x,z(:,3))
运行之后得到2条曲线:
红色的曲线就是y1对应的曲线,蓝色的是y2对应的曲线
这个方程组不是线性的,很可能没有解析解吧
楼主如果把题中的数字具体写出来也许有可能做
一般是要考虑数值解