NCEPU-除尘-CSDN博客

原创插值MATLAB实现（牛顿差商、插值误差、龙格现象、切比雪夫插值）

牛顿差商function [c,y] = newtondd(a,b,x)n=length(a);for i=1:n v(i,1)=b(i);endfor j=2:n for i=1:n-j+1 v(i,j)=(v(i+1,j-1)-v(i,j-1))/(a(i+j-1)-a(i)); endendfor j=1:nc(j)=v(1,j);endk=length(x);y=zeros(1,k);for w=1:kfor i=n:-1:2

2021-11-25 21:30:35 5953

原创非线性方程组求解Matlab实现（多元牛顿方法、Broyden方法、Broyden方法2）

三元牛顿方法（非线性方程中三个未知数）function jie = multvarnewton3(g1,g2,g3,c)e=[inf,inf,inf];syms x y zk1=g1(x,y,z);k2=g2(x,y,z);k3=g3(x,y,z);while(norm((e-c),inf)>0.5*10^-6)f1=diff(k1,x);f2=diff(k1,y);f3=diff(k1,z);f4=diff(k2,x);f5=diff(k2,y);f6=diff(k2,z)

2021-11-18 20:01:00 4551

原创方程组求解matlab实现（朴素高斯求解、LU分解、雅可比迭代方法、高斯-塞德尔方法、连续过松弛（SOR））

朴素高斯求解function x= oridinarygauss(a,b)n=length(b);for j=1:n-1 for i=j+1:n mult=a(i,j)/a(j,j); for k=j+1:n a(i,k)=a(i,k)-mult*a(j,k); end b(i)=b(i)-mult*b(j); endendx(n)=b(n)/a(n,n);for i=n-1:-1:1

2021-11-15 20:50:17 643

原创 Matlab 牛顿方法求解方程

function [jie ,n]= newton(f,g,x0,tol)%输出：jie为牛顿法函数解；n为迭代次数；%输入：f为求解函数，g为其导数，x0为初始估计，tol为控制精度。if g(x0)==0 error('初始值不合适');end%避免出现导致牛顿法导数分母为零情况。i=1;x(1)=x0;x(i+1)=x(i)-f(x(i))/g(x(i));while (abs(x(i+1)-x(i))>tol) i=i+1; x(i+1)=x(i)-f.

2021-11-01 21:20:19 6168 1

原创 Numerical analysis 数值分析编程问题 Matlab 实现--二分法计算希尔伯特矩阵问题

希尔伯特Hilbert矩阵是一个n×n矩阵，其对应的第ij个元素的值是1/（i+j-1）.令A表示5×5希尔伯特矩阵.它最大的特征值是1.567.使用二分法确定如何改变左上元素A11，使得A的最大特征值等于π.确定A11，精确到小数点后6位.可以使用MATLAB命令hilb、pi、eig和max简化计算任务.function a11 = Hilbert(n,a,b,tol)k=hilb(n);while(b-a)/2>tol c=(a+b)/2; k(1,1)=c;fc=max...

2021-10-28 19:45:57 657