qq_58208978 2021-05-21 13:22 采纳率: 33.3%
浏览 103
已采纳

为什么我这个程序出来y值全是NAN?

function dx= fun4(t,x) %建立函数 g=10; m=0.03; r=sqrt(x(1)^2+x(3)^2); r1=-x(1)/r; r2=-x(3)/r; dx=[x(2);g*m*r1/r^2;x(4);g*m*r2/r^2]; end x0=[0;0;0;0]; %调用函数 [t,y]=ode45(@fun4,[0,20],x0) plot(y(:,1),y(:,3))

  • 写回答

2条回答 默认 最新

  • CSDN专家-Matlab_Fans 2021-05-21 21:48
    关注

    初始值为全0的时候, r=sqrt(x(1)^2+x(3)^2)的取值为0了,r1=-x(1)/r; r2=-x(3)/r均为0/0,即为nan,nan与任何数运算返回值均为nan。解决办法是修改初始条件为非0值。

    x0=[1;-1;2;-2]; %调用函数 
    [t,y]=ode45(@fun4,[0,20],x0) 
    plot(y(:,1),y(:,3))

     

    本回答被题主选为最佳回答 , 对您是否有帮助呢?
    评论
查看更多回答(1条)

报告相同问题?

悬赏问题

  • ¥15 visionmaster启动失败,提示为“机器不满足授权而被禁用”
  • ¥50 用logisim设计16位单时钟周期cpu
  • ¥15 IDEA中圈复杂度如何具体设置
  • ¥50 labview采集不了数据
  • ¥15 请上面代码做什么处理或什么混淆
  • ¥15 英雄联盟自定义房间置顶
  • ¥15 W5500网线插上无反应
  • ¥15 如何用字典的Key,显示在WPF的xaml中
  • ¥15 weautomate读取Excel表格信息然后填写到网页一直报错,如何解决?
  • ¥15 C#如何在Webview2中获取网页验证码