查看: 79|回复: 7|关注: 0

[已答复] ode45求解微风方程组出现复数

[复制链接]

新手

7 麦片

财富积分


050


2

主题

6

帖子

0

最佳答案
发表于 2019-12-3 10:22:06 | 显示全部楼层 |阅读模式
新手自己编写的代码,求解结果出现复数,不知道是不是代码的问题,求解结果是否正确,求指教万分感谢!
function U11
[T,Y]=ode45(@rigid,[0 20],[1 0.8])
plot(T,Y(:,1),'-',T,Y(:,2),'--')
legend('G','q')
function dy=rigid(t,y)
dy=zeros(2,1);
dy(1)=0.1+0.2*y(2)-0.2*y(1);
dy(2)=(4.255*y(2)-3.38*y(1)+3.38*y(1)*y(2)-4.222)^(0.5)/2*(1-y(2))-0.01*y(2)-0.09;

论坛优秀回答者

中级

1331 麦片

财富积分


5001500


0

主题

1840

帖子

270

最佳答案
  • 关注者: 50
发表于 2019-12-3 10:41:48 | 显示全部楼层
你的dy(2)中存在开根号,出现复数是正常的吧

新手

7 麦片

财富积分


050


2

主题

6

帖子

0

最佳答案
 楼主| 发表于 2019-12-3 10:52:48 | 显示全部楼层
WarnerChang 发表于 2019-12-3 10:41
你的dy(2)中存在开根号,出现复数是正常的吧

所以代码没有问题?出现复数只是方程的问题吗?

论坛优秀回答者

中级

1331 麦片

财富积分


5001500


0

主题

1840

帖子

270

最佳答案
  • 关注者: 50
发表于 2019-12-3 11:44:46 | 显示全部楼层
牛奶贝拉糖53 发表于 2019-12-3 10:52
所以代码没有问题?出现复数只是方程的问题吗?

我是这么认为的呢

新手

7 麦片

财富积分


050


2

主题

6

帖子

0

最佳答案
 楼主| 发表于 2019-12-3 11:45:19 | 显示全部楼层
求解有复数,警告忽视了复数的虚部,请问做出的图是否受到影响,图像是不是会不正确?

出现警告,做出的图像是否正确

出现警告,做出的图像是否正确

新手

7 麦片

财富积分


050


2

主题

6

帖子

0

最佳答案
 楼主| 发表于 2019-12-3 11:47:19 | 显示全部楼层
WarnerChang 发表于 2019-12-3 11:44
我是这么认为的呢

好的谢谢您,以下的问题能不能再帮忙看看

论坛优秀回答者

中级

1331 麦片

财富积分


5001500


0

主题

1840

帖子

270

最佳答案
  • 关注者: 50
发表于 2019-12-3 12:13:16 | 显示全部楼层
牛奶贝拉糖53 发表于 2019-12-3 11:47
好的谢谢您,以下的问题能不能再帮忙看看

你这虚部这么小,而且虚部之间相差不多,我觉得影响不大,当然得看你的接受范围,哈哈

新手

7 麦片

财富积分


050


2

主题

6

帖子

0

最佳答案
 楼主| 发表于 2019-12-3 15:52:34 | 显示全部楼层
WarnerChang 发表于 2019-12-3 12:13
你这虚部这么小,而且虚部之间相差不多,我觉得影响不大,当然得看你的接受范围,哈哈 ...

非常感谢
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

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

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