简单的数值误差函数

2视图(30天)

显示旧的评论

Amit Kumar 2014年9月24日

0
链接

这个问题直接联系

//www.tianjin-qmedu.com/matlabcentral/answers/155969-error-in-simple-ode45-functions

评论道: 明星黾 2014年9月24日

答:接受明星黾

你好,我是试图解决non-stiff与数值函数二阶的颂歌。我的歌唱是d ^ ^ 2 = 2 y / dt - g + (4/15) * (1 / m) * (dy / dt) ^ 2

这是我ode函数文件

                         函数[ode_fun_vect] = ode_1_fun (t, z)
                        
                         全球mg
                        
                         m = 80;
                        
                         g = 9.81;
                        
                         ode_fun_vec = [z (1); - g + (4/15) * (z (1)) ^ 2 / m];
                        
                         结束

这是脚本解决歌唱:

                         clc;
                        
                         清晰的所有;
                        
                         关闭所有;
                        
                         initial_cond = (600,0);
                        
                         time_range = [0, 20];
                        
                         [t、y] =数值(@ode_1_fun, time_range initial_cond);
                        
                         图();
                        
                         次要情节(2,1,1);
                        
                         情节(t, y (: 1));
                        
                         包含(“时间”);
                        
                         ylabel (“位移”);
                        
                         次要情节(2,1,2);
                        
                         情节(t、y (:, 2));
                        
                         包含(“时间”);
                        
                         ylabel (“速度”);