应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 材料区 » 功能材料 » 磁性材料 » 关于dM/dH退磁曲线
51/1返回列表
查看: 1654  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

zxlzqz

木虫 (正式写手)

[求助] 关于dM/dH退磁曲线 已有1人参与

大神们有谁知道如何得到dM/dH的图像，是M对H的整个求导，还是部分？M,和H都是磁滞回线就是上M和H 吗如何说明交换耦合
回复此楼

» 猜你喜欢
1楼 2016-03-17 22:02:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sam799824810

银虫 (小有名气)
牛人回复了

发自小木虫Android客户端
回复此楼
4楼2016-09-14 08:40:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

lxshisan

新虫 (初入文坛)

【答案】应助回帖

论文J-A磁化建模理论的质疑和修正方法里面有相关的公式推导。根据他的论文写的Matlab代码如下：
function [x,y]=lx17()
u0=4*pi*(1e-7);
alpha=7.092e-4;
x12=linspace(0,2e4,100);
% x12=linspace(3000,4000,100);
x13=linspace(2e4,-2e4,200);
x14=linspace(-2e4,2e4,100);
M0=[0.0001,0];
% 样本1
[x,y11]=ode45(@ode_fun11,x12,M0);
[x,y12]=ode45(@ode_fun12,x13,[y11(100,1),y11(100,2)]);
[x,y13]=ode45(@ode_fun11,x14,[y12(200,1),y12(200,2)]);

figure;
plot(x13',y12(:,2),x14',y13(:,2));
xlabel('H');
ylabel('M');
title('样本1');
figure;
plot(x13',u0*(x13'+y12(:,2)),x14',u0*(x14'+y13(:,2)));
xlabel('H');
ylabel('B');
title('样本1');

% 样本2
[x,y21]=ode45(@ode_fun21,x12,M0);
[x,y22]=ode45(@ode_fun22,x13,[y21(100,1),y21(100,2)]);
[x,y23]=ode45(@ode_fun21,x14,[y22(200,1),y22(200,2)]);

figure;
plot(x13',y22(:,2),x14',y23(:,2));
xlabel('H');
ylabel('M');
title('样本2');
figure;
plot(x13',u0*(x13'+y22(:,2)),x14',u0*(x14'+y23(:,2)));
xlabel('H');
ylabel('B');
title('样本2');

% 样本3；
[x,y31]=ode45(@ode_fun31,x12,M0);
[x,y32]=ode45(@ode_fun32,x13,[y31(100,1),y31(100,2)]);
[x,y33]=ode45(@ode_fun31,x14,[y32(200,1),y32(200,2)]);

figure;
plot(x13',y32(:,2),x14',y33(:,2));
xlabel('H');
ylabel('M');
title('样本3');
figure;
plot(x13',u0*(x13'+y32(:,2)),x14',u0*(x14'+y33(:,2)));
xlabel('H');
ylabel('B');
title('样本3');
figure
plot(x13',u0*(x13'+y12(:,2)),'g',x14',u0*(x14'+y13(:,2)),'g',x13',u0*(x13'+y22(:,2)),'r',x14',u0*(x14'+y23(:,2)),'r',x13',u0*(x13'+y32(:,2)),'b',x14',u0*(x14'+y33(:,2)),'b');
grid on
% [x,Man]=ode45(@ode_fun,x12,0.0001);
%
% figure;
% plot(x,Man(:,1))
end
%样本1磁场增大求解方程
function dy=ode_fun11(x,y)
ms=1.5743e6;
a=499;
alpha=7.092e-4;
k=1154.6;
c=0.0198;
deta=1;
km=k*(1-0.96*(y(2)/ms)^2);
dy(1,1)=ms*(-(csch((x+alpha*y(1))/a))^2/a+a/(x+alpha*y(1))^2)/(1-alpha*ms*(-(csch((x+alpha*y(1))/a))^2/a+a/(x+alpha*y(1))^2));
dy(2,1)=(-1/km/deta*(y(1)-y(2))-c/(1-c)*dy(1,1))/(alpha/km/deta*(y(1)-y(2))-1/(1-c));
end
%样本1磁场减小求解方程
function dy=ode_fun12(x,y)
ms=1.5743e6;
a=499;
alpha=7.092e-4;
k=1154.6;
c=0.0198;
deta=-1;
km=k*(1-0.96*(y(2)/ms)^2);
dy(1,1)=ms*(-(csch((x+alpha*y(1))/a))^2/a+a/(x+alpha*y(1))^2)/(1-alpha*ms*(-(csch((x+alpha*y(1))/a))^2/a+a/(x+alpha*y(1))^2));
dy(2,1)=(-1/km/deta*(y(1)-y(2))-c/(1-c)*dy(1,1))/(alpha/km/deta*(y(1)-y(2))-1/(1-c));
end
%样本2求磁场增大解方程
function dy=ode_fun21(x,y)
ms=1.5755e6;
a=1408.1;
alpha=2.4e-3;
k=2356.5;
c=0.0382;
deta=1;
km=k*(1-0.96*(y(2)/ms)^2);
dy(1,1)=ms*(-(csch((x+alpha*y(1))/a))^2/a+a/(x+alpha*y(1))^2)/(1-alpha*ms*(-(csch((x+alpha*y(1))/a))^2/a+a/(x+alpha*y(1))^2));
dy(2,1)=(-1/km/deta*(y(1)-y(2))-c/(1-c)*dy(1,1))/(alpha/km/deta*(y(1)-y(2))-1/(1-c));
end
%样本2求磁场减小求解方程
function dy=ode_fun22(x,y)
ms=1.5755e6;
a=1408.1;
alpha=2.4e-3;
k=2356.5;
c=0.0382;
deta=-1;
km=k*(1-0.96*(y(2)/ms)^2);
dy(1,1)=ms*(-(csch((x+alpha*y(1))/a))^2/a+a/(x+alpha*y(1))^2)/(1-alpha*ms*(-(csch((x+alpha*y(1))/a))^2/a+a/(x+alpha*y(1))^2));
dy(2,1)=(-1/km/deta*(y(1)-y(2))-c/(1-c)*dy(1,1))/(alpha/km/deta*(y(1)-y(2))-1/(1-c));
end
%样本3求增大解方程
function dy=ode_fun31(x,y)
ms=1.5827e6;
a=1017.9;
alpha=1.2e-3;
k=2735.8;
c=0.1051;
deta=1;
km=k*(1-0.96*(y(2)/ms)^2);
dy(1,1)=ms*(-(csch((x+alpha*y(1))/a))^2/a+a/(x+alpha*y(1))^2)/(1-alpha*ms*(-(csch((x+alpha*y(1))/a))^2/a+a/(x+alpha*y(1))^2));
dy(2,1)=(-1/km/deta*(y(1)-y(2))-c/(1-c)*dy(1,1))/(alpha/km/deta*(y(1)-y(2))-1/(1-c));
end
%样本3求减小解方程
function dy=ode_fun32(x,y)
ms=1.5827e6;
a=1017.9;
alpha=1.2e-3;
k=2735.8;
c=0.1051;
deta=-1;
km=k*(1-0.96*(y(2)/ms)^2);
dy(1,1)=ms*(-(csch((x+alpha*y(1))/a))^2/a+a/(x+alpha*y(1))^2)/(1-alpha*ms*(-(csch((x+alpha*y(1))/a))^2/a+a/(x+alpha*y(1))^2));
dy(2,1)=(-1/km/deta*(y(1)-y(2))-c/(1-c)*dy(1,1))/(alpha/km/deta*(y(1)-y(2))-1/(1-c));
end
赞一下 回复此楼
2楼2016-09-13 09:03:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1174895362

新虫 (小有名气)
Origin应该能拟合出来，你研究研究
赞一下 回复此楼
3楼2016-09-13 21:39:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] filecode +11 documentary 2026-08-10 12/600 2026-08-12 12:42 by kudofaye
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +20 Lanmanbaby 2026-08-09 33/1650 2026-08-12 11:26 by lantianda
[基金申请] FileCode能看出啥？ +7 要乐观耀哥 2026-08-10 19/950 2026-08-12 10:55 by 要乐观耀哥
[基金申请] 是这周出结果还是下周出结果？ +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[硕博家园] 读博的好处 +3 lnee 2026-08-11 3/150 2026-08-11 18:10 by 希望我好好的
[基金申请] 什么时候出结果，有咨询渠道？？？ +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 分享一下我之前已中青C的计划书的filecode +3 布布和一二 2026-08-11 4/200 2026-08-11 15:20 by aasahr
[基金申请] 我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人 +8 家与远方 2026-08-10 13/650 2026-08-11 09:49 by 家与远方
[基金申请] 关于代码变化问题，想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 静等基金结果 +5 gjjjzhong 2026-08-10 16/800 2026-08-10 17:19 by Tide man
[基金申请] 好奇怪的filecode +4 布布和一二 2026-08-08 5/250 2026-08-10 14:20 by 冰心玉壶晴
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] 关于filecode，很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 固定端突然变了，今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 好消息？这个有何含义？？？ +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
信息提示
关闭
请填处理意见
关闭