看到一篇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的矩陣，即有多個特征值。應該取哪個為好。

求幫助，謝謝！

總結

以上是生活随笔為你收集整理的matlab 流程计算方法,吸波材料LLG公式计算复磁导率的过程及matlab程序的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。