matlab改进平方根算法,改进平方根请教

该楼层疑似违规已被系统折叠隐藏此楼查看此楼

function [x]=gaijin(A,b,n) %用改进平方根法求解Ax=b

L=zeros(n,n); %L为n*n矩阵

D=diag(n,0); %D为n*n的主对角矩阵

S=L*D;

for i=1:n %L的主对角元素均为1

L(i,i)=1;

end

for i=1:n

for j=1:n %验证A是否为对称正定矩阵

if (eig(A)<=0)||(A(i,j)~=A(j,i)) %A的特征值小于0或A非对称时，输出wrong

disp('wrong');break;

end

D(1,1)=A(1,1); %将A分解使得A=LDLT

for i=2:n

for j=1:i-1

S(i,j)=A(i,j)-sum(S(i,1:j-1)*L(j,1:j-1)');

L(i,1:i-1)=S(i,1:i-1)/D(1:i-1,1:i-1);

end

D(i,i)=A(i,i)-sum(S(i,1:i-1)*L(i,1:i-1)');

end

y=zeros(n,1); % x，y为n*1阶矩阵

x=zeros(n,1);

for i=1:n

y(i)=(b(i)-sum(L(i,1:i-1)*D(1:i-1,1:i-1)*y(1:i-1)))/D(i,i); %通过 LDy=b解得y的值

end

for i=n:-1:1

x(i)=y(i)-sum(L(i+1:n,i)'*x(i+1:n)); %通过LTx=y解得x的值

end

>> A=[2 -1 1,-1 -2 3,1 3 1]

b=[4 5 6]

n=3

gaijin(A,b,n)

错误使用 eig

对于标准特征值问题 EIG(A)，A 必须为方阵。

出错 gaijin (line 15)

if (eig(A)<=0)||(A(i,j)~=A(j,i))

%A的特征值小于0或A非对称时，输出wrong

matlab改进平方根算法,改进平方根请教相关推荐

【预测模型】基于Logistic混沌映射改进麻雀算法改进BP神经网络实现数据预测matlab源码
1 模型简介见这里. 2 部分代码 function [FoodFitness,FoodPosition,Convergence_curve]=SSA(N,Max_iter,lb,ub,dim,f ...
matlab算法改进,pso算法改进含MATLAB代码
%% No-dealing LOO RBF % 程序运行完以后,观察一下数据 %QPSO输出 % Validation_err(curCount,Swarmsize) ...
《混沌映射与比特重组的图像加密》（平萍等）一文的算法改进和展望
索引传送门算法改进算法展望 (初步) 传送门文献原文可前往其收录处<混沌映射与比特重组的图像加密>. 文献的算法流程 (包括算法勘误, 算法整理) 可以参见博文<<混沌 ...
实战28：数字图像可视化水印系统的设计与实现(LSB算法、DCT算法、随机间隔算法、区域校验位算法、图像降级算法、图像降级算法改进等6种数字水印算法的实现)
基于数字图像的可视化水印系统按照水印算法的不同主要分为空间域水印和变换域水印两大类.空间域水印以 LSB 算法--最低有效位算法为代表,变换域水印以 DCT 算法--离散余弦变换算法为代表[10]. ...
高光谱图像分类python语言编写改进lstm算法
高光谱图像分类python语言编写改进lstm算法改进的思想是高光谱图像存在同物易谱和异物同谱现象,导致原始的光谱信息在反映地物类别上不够准确,因此通过结合像元的类别信息,定义一个误差损失函数,求 ...
光谱分类算法 matlab,Matlab K-means聚类算法对多光谱遥感图像进行分类（一）
Matlab K-means聚类算法对多光谱遥感图像进行分类作者: 白艺亭测试了下matlab自带kmeans函数,作者编写函数,以及ENVI下的Kmeans方法,对比其效果,代码及结果图展示见下 ...
kmeans改进 matlab,基于距离函数的改进k―means 算法
摘要:聚类算法在自然科学和和社会科学中都有很普遍的应用,而K-means算法是聚类算法中经典的划分方法之一.但如果数据集内相邻的簇之间离散度相差较大,或者是属性分布区间相差较大,则算法的聚类效果十分有 ...
【自适应盲均衡8】基于分数间隔的复数改进常模算法（FSE-MMA）的瑞利衰落信道盲均衡MATLAB仿真
关注公号[逆向通信猿]更精彩!!! 知识回顾关于最基本的盲均衡算法,即CMA盲均衡的原理.推导与MATLAB仿真,可以参考 [自适应盲均衡2]多径衰落信道的复数常模算法(CMA)的理论推导与MATL ...
MATLAB基于小波变换的语音信号去噪算法改进
MATLAB基于小波变换的语音信号去噪算法改进概述 0. 需要调用的子函数 0.1 Gnoisegen函数 0.2 snrr函数 1. 语音信号输入和加噪 1.1 语音信号输入 1.2 语音信号加噪 ...
基于模拟退火算法改进粒子群SA-PSO优化shubert函数，测试函数shubert(十)，MATLAB编程实现
测试函数shubert(十) shubert函数属于周期性多峰函数,如图1所示拥有多个全局最优值,如图二所示在一个周期内只有一个全局最优值,局部最优解较多,适合测试算法的收敛性能,粒子群算法是一种收敛 ...

matlab改进平方根算法,改进平方根请教

matlab改进平方根算法,改进平方根请教相关推荐

最新文章

热门文章