24小时热门版块排行榜

>论坛更新日志 (1886)
>虫友互识 (178)
>论文投稿 (51)
>论文道贺祈福 (44)
>导师招生 (34)
>教师之家 (34)
>文献求助 (31)
>找工作 (30)
>考博 (23)
>基金申请 (22)
>硕博家园 (21)
>公派出国 (16)
>博后之家 (12)
>休闲灌水 (12)
>海外博后 (9)
>考研 (9)

返回列表

mom2Johnny

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 414.8
红花: 1
帖子: 9
在线: 28.2小时
虫号: 3313076
注册: 2014-07-08
性别: GG
专业: 金属功能材料

[交流] 吸波材料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的矩阵，即有多个特征值。应该取哪个为好。

求帮助，谢谢！

回复此楼