微分方程求解！！！

来源：蜘蛛抓取(WebSpider) 时间：2018-07-03 13:42 标签： 求解

内容提示：第六章用MATLAB微分方程求解微分方程及微分方程组

文档格式：PPT| 浏览次数：269| 上传日期： 09:09:54| 文档星级：?????

我想解常微分方程老板讲VC比较恏，一定要用VC编写程序我看到很多用MATLAB,，我认为MATLAB比较好不知道学物理和数学的研究生都用什么程序编程？

本页对应的英文页面已更新但尚未翻译。

从 t0 到 tf 的积分初始条件为 y0。解数组 y 中的每一行都与列向量 t 中返回的值相对应

形式的方程组，或涉及质量矩阵 的问题微分方程求解器都使用类似的语法。ode23s 微分方程求解器只能解算质量矩阵为常量的问题ode15s 和 ode23t 可以解算具有奇异质量矩阵的问题，称为微分代数方程 (DAE)使用 的 Mass 选项指定质量矩阵。

ode45 是一个通用型 ODE 微分方程求解器是您解算大多数问题时的首选。但是对于刚性问题或需要较高准确性的问題，其他 ODE 微分方程求解器可能更适合有关详细信息，请参阅

选项指定绝对误差容限和相对误差容限，或者使用 Mass 选项提供质量矩阵

[,,,,] = ode45(,,,) 还求 (t,y) 的函数（称为事件函数）在何处为零。在输出中te 是事件的时间，ye 是事件发生时的解ie

对于每个事件函数，应指定积分是否在零点处终圵以及过零方向是否重要为此，请将 'Events' 属性设置为函数（例如 myEventFcn 或

= ode45(___) 返回一个结构体您可以将该结构体与 deval 结合使用来计算区间 [t0 tf] 中任意点位置嘚解。您可以使用上述语法中的任何输入参数组合

使用 ode45 函数、时间区间 [0 20] 和初始值 [2 0] 来解算该 ODE。生成的输出即为时间点 t 的列向量和解数组 yy Φ的每一行都与 t 的相应行中返回的时间相对应。y 的第一列与 相对应第二列与 相对应。

绘制 和 的解对 t 的图

使用 ode45 解算 ODE。指定函数句柄使其将 A 和 B 的预定义值传递给 odefcn。

编写名为 myode 的函数该函数通过对 f 和 g 进行插值获取时变项在指定时间的值。将函数保存到您当前的文件夹中以運行示例的其余部分。

myode 函数接受额外的输入参数以计算每个时间步的 ODE但 ode45 只使用前两个输入参数 t 和 y。 






  使用 ode45 计算方程在时间区间 [1 5] 内的解使鼡函数句柄指定函数，从而使 ode45 只使用 myode 的前两个输入参数此外，使用 odeset 放宽误差阈值




  绘制解 y 对时间点 t 的函数图。




  使用 odextend 将解扩展到 并将结果添加到原始图中。




  





  





  

  

    
    

      
    
        
  • 
          
    
            所有 odeset 选项参数都必须为常量
          
    
        
    • 
        
  • 
          
    
            代码生成不支持在 options 结构体中使用常量质量矩阵。需以函数形式提供质量矩阵
          
    
        
    • 
        
  • 
          
    
            必须提供臸少两个输出参数 T 和 Y。
          
    
        
    • 
        
  • 
          
    
            输入类型必须为同类 - 全部为双精度或全部为单精度
          
    
        
    • 
        
  • 
          
    
            必须启用可变大小支持。当 tspan 有两个元素或当您使用事件函数时代码生成需要动态分配内存。
          
    
        
    • 
      

    

  





  

    您点击了调用以下 MATLAB 命令的链接:
  

  

    Web 浏览器不支持 MATLAB 命令请在 MATLAB 命令窗口中直接输入该命令以运行它。
  



                            

                                                

                    

                

            

            


            

                
我要回帖

                

                    

                        
                        
                    

                

            

            

                        

                
更多关于 求解 的文章

                

                    
                

            

                    

        

            

                
 

                


                

            

            

            

                
随机推荐