- 博客(5)
- 收藏
- 关注
RSS订阅原创 插值MATLAB实现（牛顿差商、插值误差、龙格现象、切比雪夫插值）
牛顿差商 function [c,y] = newtondd(a,b,x) n=length(a); for i=1:n v(i,1)=b(i); end for 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)); end end for j=1:n c(j)=v(1,j); end k=length(x); y=zeros(1,k); for w=1:k for i=n:-1:2
2021-11-25 21:30:35
6036
原创 非线性方程组求解Matlab实现 （多元牛顿方法、Broyden方法、Broyden方法2）
三元牛顿方法（非线性方程中三个未知数） function jie = multvarnewton3(g1,g2,g3,c) e=[inf,inf,inf]; syms x y z k1=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
4647
原创 方程组求解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); end end x(n)=b(n)/a(n,n); for i=n-1:-1:1
2021-11-15 20:50:17
701
原创 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
6240
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
706
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人