24小时热门版块排行榜    

查看: 1983  |  回复: 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朵)

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zrzneu

新虫 (初入文坛)


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

发自小木虫Android客户端
3楼2018-02-08 12:16:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答

板凳儿

新虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
2楼2016-03-07 09:42:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 314求调剂 +8 无懈可击的巨人 2026-03-12 8/400 2026-03-18 14:50 by haxia
[考研] 材料与化工一志愿南昌大学327求调剂推荐 +8 Ncdx123456 2026-03-13 9/450 2026-03-18 14:40 by haxia
[考研] 一志愿西南交大,求调剂 +4 材化逐梦人 2026-03-18 4/200 2026-03-18 14:22 by 007_lilei
[考研] 288求调剂,一志愿华南理工大学071005 +4 ioodiiij 2026-03-17 4/200 2026-03-18 12:36 by Linda Hu
[考博] 环境领域全国重点实验室招收博士1-2名 +3 QGZDSYS 2026-03-13 5/250 2026-03-18 11:13 by QGZDSYS
[考研] 293求调剂 +11 zjl的号 2026-03-16 16/800 2026-03-18 08:10 by zhukairuo
[考研] 296求调剂 +5 大口吃饭 身体健 2026-03-13 5/250 2026-03-17 21:05 by 不惑可乐
[考研] 268求调剂 +8 一定有学上- 2026-03-14 9/450 2026-03-17 17:47 by laoshidan
[考研] 材料工程专硕调剂 +5 204818@lcx 2026-03-17 5/250 2026-03-17 17:27 by Little-xue
[考研] 274求调剂 +5 时间点 2026-03-13 5/250 2026-03-17 07:34 by 热情沙漠
[考研] 304求调剂 +5 素年祭语 2026-03-15 5/250 2026-03-16 17:00 by 我的船我的海
[考研] 283求调剂 +10 小楼。 2026-03-12 14/700 2026-03-16 16:08 by 13811244083
[考研] 0703化学调剂,求各位老师收留 +8 秋有木北 2026-03-14 8/400 2026-03-16 15:21 by 哦哦123
[基金申请] 有必要更换申报口吗 20+3 fannyamoy 2026-03-11 3/150 2026-03-14 00:52 by zhanghaozhu
[考研] 材料工程调剂 +4 咪咪空空 2026-03-11 4/200 2026-03-13 19:57 by JourneyLucky
[考研] 求调剂 +5 一定有学上- 2026-03-12 5/250 2026-03-13 18:31 by ms629
[考研] 考研调剂 +4 芬达46 2026-03-12 4/200 2026-03-13 16:04 by ruiyingmiao
[考研] 工科材料085601 279求调剂 +8 困于星晨 2026-03-12 10/500 2026-03-13 15:42 by ms629
[考研] 308求调剂 +3 是Lupa啊 2026-03-12 3/150 2026-03-13 14:30 by 求调剂zz
[论文投稿] 投稿问题 5+4 星光灿烂xt 2026-03-12 6/300 2026-03-13 14:17 by god_tian
信息提示
请填处理意见