matlab中方程组数小于未知数,求助,方程个数多于未知数个数的微分方程组该怎样求解...

function dy=fangcheng(t,y)

dy=zeros(11,1);

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);

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);

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);

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;

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;

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;

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;

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;

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;

dy(8)=(1/(3.6e-4))/(-y(1)-y(3));

dy(9)=(1/(1.8e-4))/(y(3)-y(4));

dy(9)=(1/(1.8e-4))/(-y(5)+y(7));

dy(10)=(1/(1.8e-4))/(-y(2)+y(4));

dy(10)=(1/(1.8e-4))/(-y(6)-y(7));

dy(11)=(1/(3.6e-4))/(y(5)+y(6));

end

%x=zeros(11,1);

%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;

%[t,y]=ode45(@fangcheng,[0:1e-4:0.01],x);

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值