二阶常微分方程联立解有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/07/04 02:59:28
![二阶常微分方程联立解有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](/uploads/image/z/14343250-58-0.jpg?t=%E4%BA%8C%E9%98%B6%E5%B8%B8%E5%BE%AE%E5%88%86%E6%96%B9%E7%A8%8B%E8%81%94%E7%AB%8B%E8%A7%A3%E6%9C%89y1+y2+%E4%B8%A4%E5%8F%98%E6%95%B0%2C%E5%B9%B6%E4%B8%94%E6%9C%89%E4%B8%A4%E7%BB%84%E4%BA%8C%E9%98%B6%E5%B8%B8%E5%BE%AE%E5%88%86%E6%96%B9%E7%A8%8B%28%E5%AF%B9x%E5%BE%AE%E5%88%86%29%2C%E5%88%9D%E5%A7%8B%E5%80%BC%E5%B7%B2%E7%9F%A5%2Cy1%27%27+%3D+-a%2Aexp%28y1%29%2Ay1%27+-+b%2Ay1+%2B+c%2Ay2y2%27%27+%3D+-d%2Ay2+%2B+e%2Ay1+-+f%2Asin%28wx%29%E5%85%B6%E4%B8%ADa%2Cb%2Cc%2Cd%2Ce%2Cf%2Cw%E9%83%BD%E6%98%AF%E5%B7%B2%E7%9F%A5%E5%B8%B8%E6%95%B0%2C%E4%B8%94y1%280%29%3Dc1%2Cy1%27%280%29%3D0%2Cy2%280%29%3Dc2%2Cy2)
二阶常微分方程联立解有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对应的曲线
这个方程组不是线性的,很可能没有解析解吧
楼主如果把题中的数字具体写出来也许有可能做
一般是要考虑数值解