查看: 120|回复: 2|关注: 0

[已答复] 求助,方程个数多于未知数个数的微分方程组该怎样求解

[复制链接]

新手

5 麦片

财富积分


050


1

主题

2

帖子

0

最佳答案
发表于 2019-3-15 15:48:56 | 显示全部楼层 |阅读模式
  1. function dy=fangcheng(t,y)
  2. dy=zeros(11,1);
  3. dy(1)=(203.227*y(1)+200.227*y(3)+200*(y(5)+y(6))+0.6837*dy(3)+0.667*(dy(5)+dy(6))-y(8))/(-0.7437);
  4. dy(2)=(200*y(1)+3.303*y(2)+200*y(3)-0.303*y(4)+200.227*y(5)+206.227*y(6)+0.667*dy(1)+0.667*dy(3)-0.06*dy(4)+0.6837*dy(5)+0.8037*dy(6)+y(11)-y(10))/(-0.0793);
  5. dy(3)=(y(11)+200*y(1)+3*y(2)+200*y(3)+200.227*y(5)+206.53*y(6)+0.303*y(7)+0.667*dy(1)+0.06*dy(2)+0.6837*dy(5)+0.823*dy(6)+0.0193*dy(7)-y(10))/(-0.667);
  6. dy(4)=(y(9)-y(11)+0.227*y(1)+4.53*y(3)-0.303*y(4)-4.227*y(5)-0.227*y(6)+0.0167*dy(1)+0.256*dy(3)-0.2367*dy(5)-0.0167*dy(6)-y(8))/0.0193;
  7. dy(5)=(y(9)-y(11)+0.227*y(1)+4.227*y(3)-4.53*y(5)-0.227*y(6)+0.303*y(7)+0.0167*dy(1)+0.2367*dy(3)-0.0167*dy(6)+0.0193*dy(7)-y(8))/0.256;
  8. dy(7)=(y(10)-0.303*y(2)-0.303*y(3)+3.606*y(4)-3*y(7)-0.0193*dy(2)-0.0193*dy(3)+0.1286*dy(4)-y(9))/0.09;
  9. dy(6)=(y(10)+3*y(4)+0.303*y(5)-0.303*y(6)-3.606*y(7)+0.09*dy(4)+0.0193*dy(5)-0.1286*dy(7)-y(9))/0.0193;
  10. dy(6)=(y(10)-0.303*y(3)+3.303*y(4)-0.303*y(6)-3.303*y(7)-0.0193*dy(3)+0.1093*dy(4)-0.1093*dy(7)-y(9))/0.0193;
  11. dy(7)=(y(10)-0.303*y(2)+3.303*y(4)+0.303*y(5)-3.303*y(7)-0.0193*dy(2)+0.1093*dy(4)+0.0193*dy(5)-y(9))/0.1093;
  12. dy(8)=(1/(3.6e-4))/(-y(1)-y(3));
  13. dy(9)=(1/(1.8e-4))/(y(3)-y(4));
  14. dy(9)=(1/(1.8e-4))/(-y(5)+y(7));
  15. dy(10)=(1/(1.8e-4))/(-y(2)+y(4));
  16. dy(10)=(1/(1.8e-4))/(-y(6)-y(7));
  17. dy(11)=(1/(3.6e-4))/(y(5)+y(6));
  18. end
  19. %x=zeros(11,1);
  20. %x(1)=0.748;x(2)=-0.747;x(3)=0.219;x(4)=1.195;x(5)=-0.219;x(6)=-0.748;x(7)=-1.195;x(8)=204.654;x(9)=406.997;x(10)=399.355;x(11)=204.654;
  21. %[t,y]=ode45(@fangcheng,[0:1e-4:0.01],x);

复制代码
我编写的程序无法得出正确的结果,特此求助各位大神


论坛优秀回答者

权威

3223 麦片

财富积分



2

主题

3509

帖子

724

最佳答案
  • 关注者: 148
发表于 2019-3-16 10:28:09 | 显示全部楼层
方程个数多于未知数个数的微分方程组??

不要用一般的方程来说。微分方程有这样的说法吗?
怎么才是正确的?你知道答案了?
检查自己的方程,检查初值。

新手

5 麦片

财富积分


050


1

主题

2

帖子

0

最佳答案
 楼主| 发表于 2019-3-16 10:42:30 来自手机 | 显示全部楼层
我是根据电路关系推出的微分方程组,我也有这些未知量的仿真值,理论上来讲的话,计算值和仿真值应该不会相差太大的
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

站长推荐上一条 /3 下一条

快速回复 返回顶部 返回列表