24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1998  |  回复: 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的矩阵,即有多个特征值。应该取哪个为好。

求帮助,谢谢!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

板凳儿

新虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
2楼2016-03-07 09:42:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zrzneu

新虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
送红花一朵
你好,我也遇到这个问题了,计算出的特征值是三个,该怎么办

发自小木虫Android客户端
3楼2018-02-08 12:16:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 mom2Johnny 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 349学科化学045106求调剂,化学类都可以 +5 保好懂懂 2026-04-08 5/250 2026-04-08 23:32 by cheerful9622
[考研] 293调剂 +20 yj1221 2026-04-08 21/1050 2026-04-08 22:37 by 猪会飞
[考研] 318求调剂 +13 ykyhsa 2026-04-05 15/750 2026-04-08 21:37 by wj165256
[考研] 一志愿厦大生物学332求调剂 +5 池池池池池池 2026-04-08 5/250 2026-04-08 21:23 by zhyzzh
[考研] 本科211,293分请求调剂 +12 莲菜就是藕吧 2026-04-03 13/650 2026-04-08 20:30 by 背对大海出发
[考研] 288求调剂,一志愿华南理工大学071005 +11 ioodiiij 2026-04-08 11/550 2026-04-08 16:48 by tjzhao
[考研] 求调剂 +6 chenxrlkx 2026-04-05 8/400 2026-04-08 16:03 by screening
[考研] 材料专硕322分 +12 哈哈哈吼吼吼哈 2026-04-02 12/600 2026-04-08 11:43 by 1753564080
[考研] 338求调剂 +8 wxygxsaaaaa 2026-04-06 8/400 2026-04-08 06:58 by 无际的草原
[考研] 287求调剂 +3 通信学硕081000 2026-04-03 4/200 2026-04-06 21:03 by going home
[考研] 生物学求调剂 一志愿沪9,326分 +6 刘墨墨 2026-04-06 6/300 2026-04-06 19:36 by lijunpoly
[考研] 一志愿江南大学085501机械工程专硕326分,本科佳木斯大学 +5 顾若浮生 2026-04-03 9/450 2026-04-05 09:57 by 1753564080
[考研] 341求调剂 +3 学无止境,冲 2026-04-05 3/150 2026-04-05 09:40 by lbsjt
[考研] 26考研调剂0710 0860 +9 补补不补 2026-04-03 14/700 2026-04-04 23:32 by 果冻大王
[考研] 一志愿0817化学工程与技术,求调剂 +24 我不是只因 2026-04-02 28/1400 2026-04-04 15:15 by dongzh2009
[考研] 293求调剂 +5 末未mm 2026-04-02 6/300 2026-04-03 15:20 by 王保杰33
[考研] 求调剂机会 +5 意染ivy 2026-04-03 5/250 2026-04-03 15:13 by qoooooo614
[考研] 11408,284分,二战真诚求调剂 +4 12.27 2026-04-02 4/200 2026-04-03 14:14 by dxiaoxin
[考研] 一志愿北京交通大学材料工程总分358 +4 cs0106 2026-04-03 4/200 2026-04-03 13:41 by 百灵童888
[考研] 一志愿武汉理工0856,初试334 +3 26考研材料 2026-04-02 3/150 2026-04-02 21:22 by dongzh2009
信息提示
请填处理意见