24小时热门版块排行榜

返回列表

当前主题已经存档。

【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 positron 的 10 个金币

positron

木虫 (职业作家)

应助: 0 (幼儿园)
贵宾: 0.025
金币: 705.6
散金: 91
红花: 1
帖子: 4541
在线: 310.1小时
虫号: 274184
注册: 2006-08-26
专业: 核物理

[交流] 【求助】matlab解3D 偏微分方程能这样处理成二维的不？

matlab解三维偏微分方程（HF方程）时，由于在晶体中，势场是三维的，3D偏微分方程解着有点难（用ode15s求解，老提示内存不足，过程中引入的参数及时的clear掉都不成），想换一直思路，这时可否这样考虑：

由于晶体的势场是在一系列Z值（晶胞的Z方向）对应的二维势场分布（有限元分成的网格上的）这样我就可以轻松的用matlab的eig函数求二维矩阵的特征向量和特征值
每一组特征向量归一化---》波函数平方，得到空间几率分布
即A=n*n*n
for k=1:n
AA=A(:,:,k);
[V,D]=eig(AA);
给每组V特征向量归一化
有n组特征向量，再除以n,相当于对整个空间归一化
VV(:,:,zz)=V
end

还是一定要解3D 偏微分方程，得到空间的特征向量和特征值，得到空间分布？

回复此楼