数学建模学习总结..docxVIP

个人在参加数学建模时做到总结，用自己的语言写的，未必准确，但比较好懂吧。还包含一些方法的MATLAB实现。在浏览时推荐按Ctrl+F键，然后在导航窗口中浏览标题，快速转跳。数值和符号计算符号积分R=int(s,v):对符号表达式s中指定的符号变量v计算不定积分.表达式R只是表达式函数s的一个原函数,后面没有带任意常数C.?int(s):对符号表达式s中确定的符号变量计算计算不定积分.?int(s,a,b):符号表达式s的定积分，a,b分别为积分的上、下限?int(s,x,a,b):符号表达式s关于变量x的定积分，a,b分别为积分的上、下限?数值积分1变步长辛普生法数值积分的matlab实现[I,n]=quad(‘filename’,a,,b,tol,trace)Filename:被积函数名，a,b分别是积分上限和积分下限，tol是积分精度，缺省时取trace控制是否展现积分过程，非0时展现，缺省是为0。I为结果，n为函数调用次数。2被积函数由一个表格定义时求定积分用trapz(x,y):??梯形积分法，x时表示积分区间的离散化向量，y是与x同维数的向量，表示被积函数，z返回积分值。?3，二重定积分的数值求解，I=dblquad(‘f(x,y)’,a,b,c,d,tol,trace)表示f(x,y)在[a,b]*[c,d]上的二重积分，其他用法与quad相同/link?url=zcJ-sxrj0QDqzz8xCnHTnB7gxjoNRyOZzS4_4ZA22c8Bs9inYn6vVkqTVr_w-riLA_FOkWU47c9u03BlD56wA0BOLRAFPnDA98UAEr1uMI_求导数y=diff(fun),其中fun为被求导的函数例如 y1=diff(1/(exp(x)+x^2))-(2*x + exp(x))/(exp(x) + x^2)^2求n阶偏导y=diff(fun,x,n)使表达式变漂亮，可用pretty(y)变成易于数学习惯的形式微分方程dsolve(‘微分方程’，’自变量’)如dsolve(‘D2y+3*Dy+exp(x)=0’,’x’)求微分方程的特解dsolve(‘微分方程’,’初始条件’,’自变量’)如 在条件下dsolve(‘(x^2-1)*Dy+2*x*y-cos(x)=0’,’y(0)=1’,’x’)求微分方程组解通解dsolve(‘微分方程1，微分方程2’,’自变量’)特解dsolve(‘微分方程1，微分方程2’,’初始值’,’自变量’)求没有解析解的微分方程（方程组）的数值解这个解法有点难理解，主要是待解函数的写法有点奇怪1.在解n个未知函数的方程时，和x均为n维向量，m文件的待解方程可以x的分量形式写成，m文件中的函数式一个函数组，每个函数的左边都以该函数因变量数组中的一项来表示。2高阶微分方程组要等价的变成一阶微分方程组，/view/6367a003de80d4d8d15a4f3d.htmlHTTP:///view/2a3f2aea998fcc22bcd10d18.html符号计算技巧一般在符号计算过程syms x y;定义符号int diff?? 符号计算eval赋值的时候如果有表达式的话，eval可以将表达式当成语句运算，以达到赋值的目的subs是赋值函数，用数值替代符号变量替换函数，例如: subs(a+b,a,4) 意思就是把a用4替换掉，返回 4+b。也可以替换多个变量。例如：subs(cos(a)+sin(b),{a,b},[sym(alpha),2])分别用字符alpha替换a和2替换b，返回 cos(alpha)+sin(2)用法很灵活，仔细看帮助，会得到想要的形式的。插值与拟合插值就是要过已知点，拟合不一定要过已知点。一维插值：分段线性插值就是分段，在相邻的两个已知点之间用直线连接，虽然粗糙，但是很实用，工程现场用的很多。拉格朗日插值（Lagrange插值）就是在[a,b]已知n+1个已知点，过这n+1个点用至多n次的多项式去拟合原函数这就是Lagrange插值。又n次多项式有n+1个未知参数，刚好可以用这n+1个点唯一确定一个n次多项式。（不过个人感觉这是个不好的插值办法，第一计算量太大，第二n的次数一高，震荡太明显了，不符合原来的函数）Newton插值也是用n次多项式去拟合n+1个点，但是原理是用差商的累加来近似原函数的方法，Lagrange插值计算量太大，增加一个节点计算必须重新开始，牛顿法却只用在等式后面加一项就行了。埃尔米特插值（Hermite 插值）当知道已知点的导数，又要求插值函数和原函数在已知点的一阶甚至高阶导数相等时，可以采用埃米尔特插值，往往分段埃米尔特插值的效果会更好。样条线插值工程中以前常用有弹性的金属条（样条）将两点之间连起来，用金属条的位