| 查看: 2032 | 回复: 2 | ||||
mom2Johnny新虫 (初入文坛)
|
[交流]
吸波材料LLG公式计算复磁导率的过程及matlab程序 已有2人参与
|
|
看到一篇paper，利用Landau-Lifshitz-Gilbert 公式计算片状颗粒的复磁导率。(JAP 107,033913, 2010) http://scitation.aip.org/content/aip/journal/jap/107/3/10.1063/1.3291129 找到原始文献 (JAP 99, 083905, 2006) http://scitation.aip.org/content/aip/journal/jap/99/8/10.1063/1.2190719 想应用在自己的实验上，取阻尼系数a=0.1, 旋磁比r=2.8, 4*pi*Ms=20000, Ha=600, w=2*pi*f, 写了如下matlab程序求磁导率张量u clc; clear; a=0.1; % alpha damping coeeficient r=2.8; % gamma gyromagnetic ratio (MHz/Oe) Ha=600; % anisotropy field (Oe) Ms=20000/(4*pi); % Ms c=1; % count the num of input array for f=1:0.1:18 % frequency array w=2*pi*f*10^3; % omega angular frequency wm=r*Ms; % omega m w0=r*Ha; % omega 0 w1=w0+j*w*a; % omega 0' u=1+((w1*wm)/((w1)^2-w^2)); % mu k=(-1)*w*wm/((w1)^2-w^2); % kappa u_tansor=[u, j*(-1)*k, 0; j*k, u, 0; 0, 0, 1;]; % permeability tensor ui_eig=eig(u_tansor); % eigenvalue ui ui_real=real(ui_eig); % real permeability ui_imag=-1*imag(ui_eig); % imagnary permeability D=[f; ui_real; ui_imag;]; U_data(:,c)=D; u_data=U_data'; c=c+1; end 求助各位两个问题： 磁导率张量是可以求出来，但是看起来数量级不太对。不知是哪里出了问题。 文献说计算的复磁导率是张量的特征值。然而其特征值是一个3*1的矩阵，即有多个特征值。应该取哪个为好。 求帮助，谢谢！ |
回复此楼» 收录本帖的淘帖专辑推荐
 吸波材料专辑 |
» 本帖已获得的红花（最新10朵）
zrzneu» 猜你喜欢
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有0人回复
-
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有110人回复
-
【方案分享】单环磁场+轴心控制+偏转导出电子束约束系统（可行性实验）
已经有6人回复
-
【修正版】单环用磁约束低速电子实验方案（简化版）
已经有0人回复
-
桂林理工大学物理学专业招收调剂，还有三个名额！！！
已经有22人回复
-
考博自荐
已经有4人回复
-
山东大学第二批博士研究生招生
已经有0人回复
-
中国科学院东莞材料科学与技术研究所-2026年博士招生-吴昊研究员-磁学与自旋电子学
已经有0人回复
-
《电磁学》教材推荐
已经有1人回复
-
【急招】合肥工大核聚变材料计算方向2026级工程博士生
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
2楼2016-03-07 09:42:05
3楼2018-02-08 12:16:18