数值分析
花树堆雪
这个作者很懒,什么都没留下…
展开
-
【复合五点高斯-勒让德公式】
首先先看五点高斯-勒让德公式,其求积系数和求积节点可以在数值积分的教材上查到,matlab代码如下function I=Gauss_Legendre_5(fun,a,b)w=[0.2369269,0.2369269,0.4786287,0.4786287,0.5688889];x1=[-0.9061798,0.9061798,-0.5384693,0.5384693,0];x=(b+a)...原创 2018-12-03 23:48:39 · 9813 阅读 · 0 评论 -
matlab实现牛顿下山法
说起牛顿下山法,首先要提牛顿法,牛顿法是求解非线性方程的一个重要方法,具体可以点击牛顿法;虽然牛顿法作为一个二阶的迭代收敛方法,但是其对于函数和初始点的要求都比较高,而牛顿下山法则是有效降低这些要求的一种技巧。牛顿下山法的迭代公式如下xk+1(λ)=xk−λf(xk)f(xk)′x_{k+1}(\lambda)=x_k-\lambda\frac{f(x_k)}{f(x_k)^\prime}xk...原创 2018-12-23 00:22:09 · 14249 阅读 · 1 评论 -
MATLAB实现多元非线性回归
简单多元非线性回归算例现有以下数据ix1x2x3y11.123.210.12123.210.231.21.831041.11.92.910.150.92.12.910假如有以下模型:y=ax1+bx2+cx32y=ax_1+bx_2+cx_3^2y=ax1+bx2+cx32可以基于最小二乘...原创 2019-01-10 16:25:36 · 61094 阅读 · 15 评论 -
MATLAB实现列主元高斯消去法
列主元高斯消去法function x=gauss_column(A,b) %输入矩阵A和列向量b,返回解向量x[ni,nj]=size(b);if rank(A)~=rank([A,b]) %若系数矩阵秩和增广矩阵秩不相等,则无解 fprintf('无解\n') returnelse if rank(A)<ni %若系数矩阵秩和增广矩阵秩相等,但...原创 2019-01-10 17:13:11 · 23100 阅读 · 6 评论