![](https://img-blog.csdnimg.cn/20201014180756925.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
数值分析
NCEPU-除尘
环境工程学生一枚!
展开
-
非线性方程组求解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 · 4143 阅读 · 0 评论 -
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 · 591 阅读 · 0 评论 -
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 · 5955 阅读 · 1 评论 -
方程组求解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 · 567 阅读 · 0 评论