tic; % this method is transform from Ritz method %is used for solving two point BVP %this code was writen by HU.D.dong in February 11th 2017 %MATLAB 7.0 clear clc N=50; h=1/N; X=0:h:1; f=inline('pi^2/2*sin(pi/2*x)'); %以下是右端向量: for i=2:N fun1=@(x) f(X

考拉有n个字符串字符串,任意两个字符串长度都是不同的.  考拉最近学习到有两种字符串的排序方法:   1.根据字符串的字典序排序.例如: "car" < "carriage" < "cats" < "doggies < "koala" 2.根据字符串的长度排序.例如: "car" < "cats" < "koala" &l

一般而言,方程没有能够普遍求解的silver bullet,但是有几类方程的求解方法已经非常清晰确凿了,比如线性方程.二次方程或一次分式.一次方程可以直接通过四则运算反解出答案,二次方程的求根公式也给出了只需要四则运算和开根号的符号表达式.而一次分式的分子即为一次函数.更多的方程并没有普适的符号表达式,但通过用便于求零点的函数模仿.代替之也可以估计零点的位置.插值方法可以实现这一思路. 插值迭代方法包括割线法.二次插值法等多项式插值方法,反插法以及线性分式插值法等等,其核心是用几个点及其函数值信

#先上代码后补笔记# #可以直接复制粘贴调用的MATLAB函数代码!# 1. 亚当斯(Adams)预测-修正算法 由亚当斯-巴什福特(Adams-Bashforth)显式预测公式和亚当斯-莫顿(Adams-Moulton)隐式修正公式组成的预测-修正(PECE)对. function [ YMat ] = Adams( func, tvec, y_init, order ) % Adams预测-修正算法,用于求解常微分初值问题 % 输入四个参数:函数句柄func(接收列向量.返回列向量),积分时

#先上代码后补笔记# #可以直接复制粘贴调用的MATLAB函数代码!# 1. 朗格-库塔(Runge-Kutta)方法族 目前只实现了四阶Runge-Kutta方法. function [ YMat ] = Runge_Kutta( func, tvec, y_init, order ) % 朗格-库塔(Runge-Kutta)算法解常微分方程(组): % 输入四个参数:函数句柄func(接收列向量.返回列向量),积分时间列向量tvec,初值行向量y_init,阶数order: % 输出一个参数

#先上代码后补笔记# #可以直接复制粘贴使用的MATLAB函数!# 1. 定步长牛顿-柯茨积分公式 function [ integration ] = CompoInt( func, left, right, step, mode ) % 分段积分牛顿-柯茨公式: % 输入5个参数:被积函数(FUNCTIONHANDLE)'func',积分上下界'left'/'right',步长'step', % 模式mode共三种('midpoint','trapezoid','simpson'): % 返

计算科学(也称科学计算 scientific computation 或 SC)是一个快速增长的多学科领域,使用先进的计算能力来理解和解决复杂的问题. 计算科学包括三个不同的方面: 1. 开发用于解决科学(例如生物,物理和社会).工程.人文问题的算法(数值和非数值),或建模软件,或仿真软件:    2. 计算机和信息科学,开发和优化高级系统硬件,软件,网络和数据管理组件,解决计算要求严苛的问题:    3. 计算基础设施,支持科学和工程问题解决.发展计算机和信息科学: 在实际使用中,通常是从数值

FESTUNG - 3. 采用 HDG 方法求解对流问题[1] 1. 控制方程 线性对流问题控制方程为 \[\begin{array}{ll} \partial_t c + \nabla \cdot f = h, & \mathrm{in} \; J \times \Omega \\ c(x, 0) = c_0(x), & \mathrm{on} \; {0} \times \Omega \\ c(x, t) = c_D, & \mathrm{on} \; J \times \pa

选修了cs的AI课,开始有点不适应,只能用matlab硬着头皮上了,不过matlab代码全网仅此一份,倒有点小自豪. 一.练习题目 分别用宽度优先.深度优先.贪婪算法和 A*算法求解"罗马利亚度假问题".具体地图我这里不给出了,有兴趣的可以去搜索.即找到从初始地点 Arad到 目的地点 Bucharest 的一条路径. 要求:分别用文件存储地图和启发函数表,用生成节点数比较以上四种算法在同一问题求解时的效率,列表给出结果. 附:罗马尼亚度假问题图(图1.1) 图1.1 罗马尼亚度假问题

最近建立了一个网络流模型,是一个混合整数线性规划问题(模型中既有连续变量,又有整型变量).当要求解此模型的时候,发现matlab优化工具箱竟没有自带的可以求解这类问题的算法(只有bintprog求解器,但是只能求解不含连续变量的二值线性规划问题).于是在网上找了一些解决问题的途径,下面说说我试过的几种可能的解决方案,包括cplex.GLPK.lpsolve 和 yalmip. cplex 首先想到的是IBM公司大名鼎鼎的cplex.cplex是IBM公司一款高性能的数学规划问题求解器,可以快速.

只是用来求解函数的部分一个根,具体算法没查询,只是调用自带的函数 代码如下 % 求函数零点和极小值 xiszero=func(0) x=-1:0.1:1.5; y=func(x); plot(x,y); hold on h1=refline(0,0); set(h1,'color','r')%先画图帮助选择合理的初值 xsolve=fzero('func',-0.5);%第一次参数是函数名称,第二个参数是初值 mpoint = fminbnd('func',0.5,1);%求极小值,注意只能用小

matlab中有专门的solve函数来解决方程组的(a-x)^2+(b-y)^2=e^2(C-x)^2+(D-y)^2=v^2已知a,b,c,d,e,v 值求解 X,Y 请问用 matlab 如何写,就是求2个园的交点问题.仿真程序为:global a b c d e v;>> a=1;b=0;c=-1;d=0;e=1.5;v=1.5;%设定你这几个未知数的值>> syms x y;%%%%%%x,y是变量>> [x,y]=solve('x^2+y^2-2*a*x-2*

曾攀老师的<有限元分析基础教程>第三章有二维杆单元的推导,并结合一个例题进行了解析解和基于Matlab的程序求解.但是我感觉书中的MATLAB代码有点罗嗦,而且一些实现方法也比较麻烦,比如已经知道了杆单元的起点和终点坐标,仍然需要另外给出单元局部坐标与整体坐标的夹角,这完全没必要.于是我就用Python重构了这段程序,当然并不是把书中的MATLAB代码翻译成python(事实上完全可以这么干,而且很快!).比如我使用了面向对象的思想,把杆单元写成了一个类,这样思路比较清晰. #! /usr/b

小白往往听到微分方程就觉得害怕,其实数学建模中的微分方程模型不仅没那么复杂,而且很容易写出高水平的数模论文. 本文介绍微分方程模型边值问题的建模与求解,不涉及算法推导和编程,只探讨如何使用 Python 的工具包,零基础求解微分方程模型边值问题. 通过 3个 BVP 案例层层深入,手把手教你搞定微分方程边值问题. 欢迎关注『Python小白的数学建模课 @ Youcans』系列,每周持续更新 1. 常微分方程的边值问题(BVP) 1.1 基本概念 微分方程是指含有未知函数及其导数的关系式. 微分

之前转载了一篇博客http://blog.sina.com.cn/s/blog_6163bdeb0102dvay.html,讲Matlab网格划分程序Distmesh,看了看程序,感觉程序写得有很多值得学的地方,所以又自己重新看了一看,加了一些注释,最后再总结一下学到的东西吧. 源代码的地址已经改变,如下http://people.sc.fsu.edu/~jburkardt/m_src/distmesh/distmesh.html.程序给出了20个划分的例子,文件名为p1_demo.m~p20_

原文:[原创]开源Math.NET基础数学类库使用(06)数值分析之线性方程组直接求解 开源Math.NET基础数学类库使用系列文章总目录:   1.开源.NET基础数学计算组件Math.NET(一)综合介绍    2.开源.NET基础数学计算组件Math.NET(二)矩阵向量计算    3.开源.NET基础数学计算组件Math.NET(三)C#解析Matlab的mat格式   4.开源.NET基础数学类库使用Math.NET(四)C#解析Matrix Marke数据格式   5.开源.NET基

遗传算法 1.前言 遗传算法是一种基于生物界自然群体遗传进化机制的自适应全局优化概率搜索算法.它与传统算法不同,不依赖梯度信息,而是通过模拟自然进化过程来搜索最优解. 例子:兔子的遗传进化 有人说,现代医学阻碍了人类的进化?你怎么看? 2．发展历程 遗传算法由密歇根大学的约翰·霍兰德和他的同事于二十世纪六十年代在对细胞自动机(英文:cellular automata)进行研究时率先提出.在二十世纪八十年代中期之前,对于遗传算法的研究还仅仅限于理论方面,直到在匹兹堡召开了第一届世界遗传算法大会.随

灰色预测的主要特点是只需要4个数据,就能解决历史数据少,序列的完整性以及可靠性低的问题,能将无规律的原始数据进行生成得到规律性较强的生成序列,易于检验 但缺点是只适合中短期的预测,且只适合指数级增长的预测. 在建立灰色预测模型之前,需先对原始时间序列进行数据处理,经过数据预处理后的数据序列称为生成列.对原始数据进行预处理,不是寻找它的统计规律和概率分布,而是将杂乱无章的原始数据列通过一定的方法处理,变成有规律的时间序列数据,即以数找数的规律,再建立动态模型. 灰色预测通过鉴别系统因素之间发展趋势

前言 最近小编学了运筹学中的单纯形法.于是,很快便按奈不住跳动的心.这不得不让我拿起纸和笔思考着,一个至关重要的问题:如何用单纯形法装一个完备的13? 恰巧,在我坐在图书馆陷入沉思的时候,一位漂亮的小姐姐靠过来,说:“同学,你是在看线性规划吗?你能帮我看看这道题该怎么解好吗?” 纳尼?还真是瞌睡来了送枕头.但是,尽管心里万马奔腾,还是要装作若无其事的样子,蛋蛋一笑.“这个啊,简单!让我来算算.” 但是一拿到题目之后,扫了一眼.惊得差点没把笔吞下去.这……城里人都这么会play的吗?我*,25个变

Eigen库矩阵运算使用方法 Eigen这个类库,存的东西好多的,来看一下主要的几个头文件吧: ——Core 有关矩阵和数组的类,有基本的线性代数(包含 三角形 和 自伴乘积 相关),还有相应对数组的操作. ——Geometry 几何学的类,有关转换.平移.进位制.2D旋转.3D旋转(四元组和角轴相关) ——LU 逻辑单元的类,有关求逆,求行列式,LU分解解算器(FullPivLU,PartialPivLU) ——Cholesky 包含LLT和LDLT的乔里斯基因式分解法. (小科普:Chole

一起来学matlab-matlab学习笔记13函数 13_2 匿名函数 觉得有用的话,欢迎一起讨论相互学习~Follow Me 参考文献 https://ww2.mathworks.cn/help/matlab/matlab_prog/anonymous-functions.html?searchHighlight=%E5%8C%BF%E5%90%8D%E5%87%BD%E6%95%B0&s_tid=doc_srchtitle 什么是匿名函数 匿名函数是不存储在程序文件中.但与数据类型是 fun