方程求解程序清单 a=-1,b=2,c=-1; w=1; m=2; n=1; h = 0.02; t=0:h:30; s1=dsolve('a*D2y+b*Dy+c*y=sin(w*t)','y(0)=m,Dy(0)=n','t'); s1_n = eval(s1); hold on plot(t,s1_n,'ko'); EulerOED(a,b,c,w,m,n,h); hold off function EulerOED(a,b,c,w,x0,x1,h) A = [x0;x1]; t=0:h:30; for i = 1:1:length(t)-1 A(:,i+1) = [1,h;(-(c/a)*h),(1-(b/a)*h)]*A(:,i) + [0;(h/a)]*sin(w*t(i)); end plot(t,A(1,:),'r*'); 对于二阶全微分方程a*y''(t)+b*y'(t)+c=sin(wt) ,不同的a,b,c,w取值和初始条件会求出不同的解,通解又是由齐次解和特解组成。其中,齐次解由特征方程决定,而特解的决定因素则比较复杂。 讨论思路 (1)通解随初始条件变化情况 (2)通解随a,b,c变化情况 b^2-4ac>0(两个不同的实根) b^2-4ac=0(两个相同的重根) b^2-4ac<0(两个不同的复数根) 1).b>0 2).b=0 3).b<0 (3)通解随w变化情况 b^2-4ac=0情况 b^2-4ac<0情况 (3)通解随w变化的规律 W属于(0,1)时,随w的增大在齐次解的旁边波动 w属于(1,+),随w的增大逐渐趋近于齐次解。 Matlab解二阶常微分方程 方程:a*y''(t)+b*y'(t)+c=sin(wt) 求解:1.解析解 2.数值解(欧拉方法) 目的:1.比较两种求解方式的拟合情况 2.通解随w变化的规律 (1)通解随初始条件变化情况 Ex1: a=2,b=3,c=1,y(0)=0;y'(0)=0,w=1 Ex2: a=2,b=3,c=1,y(0)=2;y'(0)=0,w=1 Ex3: a=2,b=3,c=1,y(0)=2;y'(0)=4,w=1 (2)通解随a,b,c变化情况 Ex1: a=2,b=3,c=1,y(0)=0;y'(0)=0,w=1 Ex4: a=-2,b=3,c=1,y(0)=0;y'(0)=0,w=1 Ex5: a=2,b=-3,c=1,y(0)=0;y'(0)=0,w=1 Ex6: a=2,b=3,c=-1,y(0)=2y'(0)=1,w=1 ? EX: a=2 ,b=2*sqrt(2) ,c=1,y(0)=0;y'(0)=0,w=1 (3).b^2-4ac<0 EX:a=4,b=-1,c=2,y(0)=0;y'(0)=0,w=1 EX:a=4,b=1,c=2,y(0)=3,y'(0)=0,w=1 EX:a=4,b=0,c=1,y(0)=2;y'(0)=0,w=1

matlab求二阶微分方程的通解,二阶常微分方程matlab的数值解和解析解分析总报告.ppt...相关推荐

  1. matlab求分段函数的值.,如何用MATLAB求分段函数的最小值和最大值?

    7.1.1 分段线性插值 所谓分段线性插值就是通过插值点用折线段连接起来逼近原曲线,这也是计算机绘制图形的基本原理.实现分段线性插值不需编制函数程序,MATLAB自身提供了内部函数interp1其主要 ...

  2. 【Matlab】求解微分方程{上}(通解和特解)

    求解微分方程 desolve函数 实例1 实例2 实例3 实例4 求解有条件的微分方程 微分方程显示隐式解 未找到显式解决方案时查找隐式解决方案 求微分方程级数解 为具有不同单边限制的函数指定初始条件 ...

  3. matlab求方程在X附近的根,matlab 实验03 求代数方程的近似根(解)

    实验三 求代数方程的近似根(解) 求代数方程的根是最常见的数学问题之一(这里称为代数方程,主要是想和后面的微分方程区别开.为简明起见,在本实验的以下叙述中,把代数方程简称为方程),当是一次多项式时,称 ...

  4. MATLAB中输入微分方程dy表示,怎么用MATLAB求解如Dy = y+1/y 的微分方程

    怎么用MATLAB求解如Dy = y+1/y 的微分方程 关注:239  答案:2  mip版 解决时间 2021-01-28 19:40 提问者妳熄滅叻菸,説啓従偂 2021-01-27 19:41 ...

  5. matlab 求向量的交集_从零开始的matlab学习笔记——(16)函数绘图

    matlab应用--求极限,求导,求积分,解方程,函数绘图....更多内容尽在个人专栏:matlab学习 这一节开始我们正式进入函数绘图,内容比较繁杂,我尽量把它整合的简洁一些. 单一图像: 首先我们 ...

  6. 用Matlab求拐点和凹凸区间需要使用MATLAB

    1.启动MATLAB,输入如下代码,求解函数曲线的一阶导数.二阶导数.拐点和拐点处的切线方程. close all; clear all; clc;format compact % MATLAB之考研 ...

  7. matlab求系统稳定时k的范围,Matlab大作业

    一.通过举例说明运用MATLAB 判别控制系统稳定的所有方法 稳定是控制系统是否能进行工作的首要条件.一般来说,稳定性成为区分系统是否有用的标志.从实际应用的角度来看,可以认为只有稳定的系统才有用.而 ...

  8. matlab求不等式的方法,求解变分不等式的matlab程序我需要

    2008-02-21 不等式的解法关于一元三次不等式的解法, 全日制普通高级中学教科书(试验修订本.必修)第一章第四节是"含绝对值的不等式解法",第五节是"一元二次不等式 ...

  9. matlab 求系统码的校验矩阵,数字通信系统matlab仿真

    课程设计报告 题目:基于MATLAB的通信系统仿真 ---信道编码对通信系统性能的影响 专业: 通信工程 姓名: XXX 学号: 0730xxxx 基于MATLAB的通信系统仿真 ---信道编码对通信 ...

最新文章

  1. 设置VSCode快速切换多个项目窗口的快捷键Alt+E
  2. 我的Debian 8.0 (jessie)配置文档
  3. 浅析企业网站应该如何选择适合自己的服务器?
  4. 关于路径搜索的算法, 可能用到
  5. PHP实现图片压缩的两则实例(转)
  6. Vue中通过Axios向SpringBoot发送get和post请求
  7. 数据结构与算法:已知二叉树两种遍历序列,求第三种遍历序列
  8. Java verts_Java-JOGL绘制VBO黑屏
  9. 面试官、女朋友都满意系列 - 决策树
  10. Oracle date日期类型使用
  11. 微信小程序如何快速生成wxml,wxss,json,js目录文件
  12. (转)马云的江湖 PK 史玉柱的兵法
  13. java 数组溢出异常,Java数组超出范围时如何处理多个异常?
  14. eclipse混淆打包出错
  15. 面试问到处理过什么棘手问题_为什么调节人工智能如此棘手?
  16. 带时间窗的车辆路径规划问题(VRPTW)
  17. c语言编木马程序,5分钟教会!C语言远程控制木马:“控制端”制作,附送源码!...
  18. mac版android sdk安装手机模拟器
  19. Python-pandas-Excel排序、筛选
  20. 基于 Petri 网的软件过程支撑环境设计

热门文章

  1. iPad /iPhone无法开机、白苹果怎么办?教你快速解决iPad /iPhone开不了机的难题
  2. 2020北大信科计算机考研公示,2020年北京大学信息科学技术学院硕士研究生拟录取名单.pdf...
  3. 禁止复制服务器文件夹,远程桌面服务器 禁止复制文件夹
  4. python实现AI换脸
  5. 【转】贾佳亚港中文团队冠军技术分享:最有效的COCO物体分割算法
  6. 基于金字塔模型的自动化测试设计
  7. 互联网名称与数字地址分配机构ICANN简介
  8. 【flutter】type ‘Null‘ is not a subtype of type ‘bool‘
  9. 带中文字库的12864驱动笔记
  10. 参考文献格式自动生成器