matlab中高精度计算函数vpa与非厄米矩阵本征值的求解

bash 复制代码
clear;clc;close all
tic
%并行设置%
delete(gcp('nocreate'));%关闭之前的并行
c=parcluster('local');
c.NumWorkers=50;%手动设置线程数(否则默认最大线程为12)
parpool(c, c.NumWorkers);
%并行设置%

w=1;
u=2.5;

N=30;


valsto=zeros(2*N+2,100);
v10=linspace(-3,3,100).';


parfor jj=1:length(v10)


    v1=v10(jj);

    [~,val]=eig(vpa(H(u,v1,w,N),30));%设置30位高精度计算
    %[~,val]=eig(H(u,v1,w,N));
    tval=diag(val);
    [~,index]=sort(real(tval));
    val=tval(index);
    valsto(:,jj)=val;

end

figure
hold on
plot(linspace(-3,3,100),abs(valsto(:,:)),'k.')

%%%%%%%%%%=========%%%%%%%%%%%%%%
%%%%%%%%%%=========%%%%%%%%%%%%%%



clearvars -except N
syms u v w real


val2=eig(H(u,v,w,N));
val2=subs(val2,{u,w},{2.5,1});

valsto2=zeros(2*N+2,100);
v10=linspace(-3,3,100).';


parfor jj=1:length(v10)


    v1=v10(jj);
    tval=val2;
    tval=double(subs(tval,v,v1));
    [~,index]=sort(real(tval));
    tval=tval(index);
    valsto2(:,jj)=tval;

end

plot(linspace(-3,3,100),abs(valsto2(:,:)),'bo')
hold off
saveas(gcf,"VAPtest0.fig")
toc

function [H0]=H(u,v,w,N)

T1=[];
T2=[];
for hh=1:N
    T1=[T1;v+u;w+u];
    T2=[T2;v-u;w-u];
end
T1=[T1;v+u];
T2=[T2;v-u];
H0=diag(T1,1)+diag(T2,-1);


end

上面代码没用高精度计算函数vpa之前数值计算的特征值能谱有部分和符号计算(篮圈)是对不上的,如下:

使用了vpa高精度函数后,数值计算和符号计算的结果完全一致:

相关推荐
豆沙沙包?16 小时前
2025年--Lc186--64. 最小路径和(多维动态规划,矩阵)--Java版
java·矩阵·动态规划
lingchen190619 小时前
矩阵的除法
人工智能·算法·矩阵
机器学习之心HML19 小时前
MATLAB基于GWO-BP神经网络对某拨叉件锻造金属流动性的参数分析
开发语言·神经网络·matlab
珞瑜·21 小时前
MATLAB2025B版本新特点
matlab
豆沙沙包?1 天前
2025年--Lc188--931. 下降路径最小和(多维动态规划,矩阵)--Java版
java·矩阵·动态规划
熬了夜的程序员1 天前
【LeetCode】74. 搜索二维矩阵
线性代数·算法·leetcode·职场和发展·矩阵·深度优先·动态规划
点云SLAM1 天前
矩阵奇异值分解算法(SVD)的导数 / 灵敏度分析
人工智能·线性代数·算法·机器学习·矩阵·数据压缩·svd算法
坚持编程的菜鸟1 天前
LeetCode每日一题——矩阵置0
c语言·算法·leetcode·矩阵
信息快讯2 天前
“COMSOL+MATLAB光子学仿真:从入门到精通,掌握多物理场建模
开发语言·matlab·comsol·光学工程
北山太湖2 天前
Matlab安装硬件支持包
开发语言·matlab