N朵小红花-CSDN博客

原创牛顿迭代法的Matlab

function [k,X]=Newton(y,x0)format longsyms xm=zeros(1,100);for i=1:100 s=x0-subs(y,x,x0)/subs(diff(y),x,x0); if abs(s-x0)>5e-4 x0=s; m(i)=x0; else break; endendk=i-1;X=m(1,1:i-1);结果显示：> [k,X]=Newton(x

2020-12-27 15:18:31 566

原创改进欧拉的Matlab

改进欧拉function [x,y]=oula(fun,xspan,y0,h) x=xspan(1):h:xspan(2); y(1)=y0; for n=1:length(x)-1 k1=feval(fun,x(n),y(n)); y(n+1)=y(n)+h*k1; k2=feval(fun,x(n+1),y(n+1)); y(n+1)=y(n)+h*(k1+k2)/2; end x=x';y=y';

2020-12-27 15:02:30 5760

原创龙贝格的Matlab

function [RT,R,wugu,h]=romberg(fun,a,b, wucha,m)n=1;h=b-a; wugu=1; x=a;k=0;RT=zeros(4,4); RT(1,1)=h*(feval(fun,a)+feval(fun,b))/2;while((wugu>wucha)&(k<m)|(k<4)) k=k+1; h=h/2; s=0; for j=1:n x=a+h*(2*

2020-12-27 14:55:55 318

原创高斯-勒让德的Matlab

function [ql,Ak,xk]=guasslegendre(fun,a,b,n,tol)if nargin==1 a=-1;b=1;n=7;tol=1e-8;elseif nargin==3 n=7;tol=1e-8;elseif nargin==4 tol=1e-8;elseif nargin==2|nargin>5 error('The Number of Input Arguments Is Wrong!');endsyms xp=sym2poly(diff((x

2020-12-27 14:51:35 3361

原创杜利特尔（直接三角分解法）的Matlab

function [L,U,y,x] = Doolittle(A, b) [row_a, col_a] = size(A);for j = 1:col_a U(1,j) = A(1,j);end L(1,1) = 1;for i = 2:row_a L(i,1) = A(i,1)/A(1,1); end for i = 2:row_a for j = i:col_a temp_sum = 0; for k = 1:i-1

2020-12-27 14:43:04 7680 1

原创最佳平方逼近的Matlab

最佳平方逼近是先求积得到希尔伯特矩阵，再根据矩阵求解方程。function [A,B,p]=zjpfbj_zzy(weightfun,primaryfun,fun,a,b)% ---wfun：权函数% ---phifun：基函数% ---fun：逼近函数% ---a,b：逼近区间的端点% ---n：最佳逼近的次数，默认值为1% 输出参数：% ---A：法方程组的系数矩阵% ---b：法方程组的右端向量% --

2020-12-27 14:25:08 5863 3

原创梯形公式和辛普森的Matlab

梯形公式：function y = TXGS(a,b,n,f)h = (b - a) / n;x = linspace(a,b,n+1);yy = h * feval(f,x);yy(1) = yy(1) / 2;yy(n+1) = yy(n+1) / 2;y = sum(yy);运行结果：>> f(x)=inline('sqrt(4-sin(x).^2)');>> TXGS(0,pi/6,6, f(x))ans = 1.0356219003136

2020-12-27 14:14:10 6718 2

原创三次样条插值的Matlab

function [m,H,lambda,mu,D,A,dY,s]=ClampedSp_1 (X,Y,dyo,dyn)m=length(X);A=zeros(m,m); n=m-1;H=zeros(1,n);lambda=zeros(1,n);mu=zeros(1,n);lambda(1)=1;A(1,1)=2;A(1,2)=lambda(1);D=zeros(1,n);mu(n)=1;for k=1:n hk=X(k+1)-X(k);H(k)=hk;endfor k=1:n-1

2020-12-27 14:10:36 2738 6

原创牛顿插值的Matlab

牛顿插值先求出差商表，再去利用差商表的对角线的值去求多项式。function [A,C,L,wcgs,Cw]= newploy_final(X,Y)n=length(X);%向量X中元素xi的个数，此处n=5A=zeros(n,n); %把变量A定义为一个n行n列的零矩阵用来存储各阶差商A(:,1)=Y'; %将Y转置赋值给零矩阵A的第一列 s=0.0; p=1.0; q=1.0; c1=1.0; for j=2:n for i=j:n

2020-12-27 14:08:00 825

qq_37233416的博客