数值分析
花了心的大萝卜
这个作者很懒,什么都没留下…
展开
-
Matlab实现——to find the fixed point
欢迎前往个人博客驽马点滴和视频空间 哔哩哔哩-《挨踢日志》fixpt.mfunction [k,p,err,P]=fixpt(g,p0,tol,max1)% Input - g is the iteration function input as a string ??g??% - p0 is the initial guess for the fi...原创 2011-09-13 18:44:28 · 1391 阅读 · 0 评论 -
Matlab实现——捕食者-被捕食者模型
欢迎前往个人博客驽马点滴和视频空间 哔哩哔哩-《挨踢日志》问题重述:捕食者-被捕食者模型非线性微分方程的一个例子是捕食者-被捕食者模型。设x(t)和y(t)分别表示兔子和狐狸在时刻t的数量,捕食者-被捕食者模型表明,x(t)和y(t)满足:一个典型的计算机模拟可作用系数:A=2, B=0.02, C=0.0002, D=0.8且满足:x(0)=30...原创 2011-09-13 20:06:53 · 22093 阅读 · 7 评论 -
Matlab实现——Euler's Method
欢迎前往个人博客驽马点滴和视频空间 哔哩哔哩-《挨踢日志》euler.m%Program9.1 (Euler's Method)function E=euler(f,a,b,ya,M)%Input - f is the function entered as a string ??f??% - a and b are the left and right ...原创 2011-09-13 20:05:55 · 3957 阅读 · 0 评论 -
Matlab实现——Recursive Trapezoidal Rule
欢迎前往个人博客驽马点滴和视频空间 哔哩哔哩-《挨踢日志》rctrap.m%Program 7.3 (Recursive Trapezoidal Rule).function T=rctrap(f,a,b,n,tol)%Input - f is the integrand input as a string ’f’% - a and b are upper an...原创 2011-09-13 19:40:43 · 2140 阅读 · 0 评论 -
Matlab实现——Lagrange approximation
欢迎前往个人博客驽马点滴和视频空间 哔哩哔哩-《挨踢日志》lagran.m%Program 4.1 (Lagrange approximation)function [C,L]=lagran(X,Y)%Input - X is a vector that contains a list of abscissas% - Y is a vector that...原创 2011-09-13 19:39:23 · 1463 阅读 · 0 评论 -
Matlab实现——Differentiation Using Limits
欢迎前往个人博客驽马点滴和视频空间 哔哩哔哩-《挨踢日志》difflim.m%Program 6.1 (Differentiation Using Limits). To approximate %f'(x) numerically by generating the sequencef'(x)≈Dk=%[f(x+(10^-k)*h)-f(x-(10^-k)*h]/[2*(...原创 2011-09-13 19:38:12 · 842 阅读 · 0 评论 -
Matlab实现——Clamped Cubic Spline
欢迎前往个人博客驽马点滴和视频空间 哔哩哔哩-《挨踢日志》%csfit.m%Program 5.3 (Clamped Cubic Spline ).function S=csfit(X,Y,dx0,dxn)%Input - X is the 1xn abscissa vector% - Y is the 1xn ordinate vector% ...原创 2011-09-13 19:30:07 · 2819 阅读 · 0 评论 -
Matlab实现——Adaptive Quadrature Using Simpson’s Rule
欢迎前往个人博客驽马点滴和视频空间 哔哩哔哩-《挨踢日志》%Input - f is the integrand input as a string ¡¯f¡¯% - a and b are upper and lower limits of integration% - tol is the tolerance% Output - SRmat is the table...原创 2011-09-13 19:29:30 · 1683 阅读 · 0 评论 -
Matlab实现——Adams-Bashforth-Method
欢迎前往个人博客驽马点滴和视频空间 哔哩哔哩-《挨踢日志》abm.mfunction A = abm(f,T,Y)%Input - f is the function entered as a string 'f'% - T is the vector of abscissas% - Y is the vector of ordina...原创 2011-09-13 19:18:36 · 5489 阅读 · 1 评论 -
Matlab实现——Jacobi Iteration and Gauss-Seidel Iteration
欢迎前往个人博客驽马点滴和视频空间 哔哩哔哩-《挨踢日志》Use Jacobi Iteration and Gauss-Seidel Iteration to solve the following linear system with several different initial guesses for p.x+z=2-x+y=0x+2y-3z=0下面是这两个程...原创 2011-09-13 19:15:34 · 2330 阅读 · 0 评论 -
Matlab实现——Upper Triangularization Followed by Back Subtitution
欢迎前往个人博客驽马点滴和视频空间 哔哩哔哩-《挨踢日志》uptrbk.m%Program 3.2 (Uppr Triangularization Followed by Back Subtitution)function X = uptrbk(A,B)%Input - A is an N x N nonsingular matrix% ...原创 2011-09-13 19:13:24 · 7302 阅读 · 0 评论 -
Matlab实现——Forward-substitution
欢迎前往个人博客驽马点滴和视频空间 哔哩哔哩-《挨踢日志》%Program 3.1.1 (Forward-substitution)function X=forwardsub(A,B)%Input - A is an n x n lower-triangular nonsingular matrix% - B is an n x 1 matrix%Out...原创 2011-09-13 19:09:35 · 3007 阅读 · 0 评论 -
Matlab实现——严格对角占优三对角方程组求解(高斯赛尔德Gauss-Seidel迭代、超松弛)
欢迎前往个人博客驽马点滴和视频空间 哔哩哔哩-《挨踢日志》严格对角占优三对角方程组求解对中等规模的n阶的(n<100)线性方程组,直接法的准确性和可靠性,所以常采用直接法对于较高阶的方程组,特别是地于某些偏微分方程离散化后得到的大型稀疏方程组(系统矩阵绝大多数为零元素),由于直接解法的计算代价较高,使得迭代法更具有竞争力。于是设计以下的2种算法:...原创 2011-09-13 18:54:49 · 9892 阅读 · 0 评论 -
Matlab实现——Newton-Raphson Iteration
欢迎前往个人博客驽马点滴和视频空间 哔哩哔哩-《挨踢日志》Newton.m%Program 2.5 (Newton-Raphson Iteration) function [p0,err,k,y]=newton(f,df,p0,delta,epsilon,max1)% Input - f is the object function input as a string '...原创 2011-09-13 18:45:30 · 6745 阅读 · 2 评论 -
Matlab实现——求矩阵的逆(LU分解)
欢迎前往个人博客驽马点滴和视频空间 哔哩哔哩-《挨踢日志》Program ( :Solve By Factorization with Pivoting )思路及原理:就得到:程序:function X=Ni(A)%Input - A is an N x N matrix%Output - I is an N x N inverse matrix of ...原创 2011-09-13 20:07:46 · 25723 阅读 · 8 评论