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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 材料区 » 功能材料 » 磁性材料 » 关于dM/dH退磁曲线
51/1返回列表
查看: 1553  |  回复: 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 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求材料调剂 +8 隔壁陈先生 2026-03-12 8/400 2026-03-18 22:19 by li123456789.
[考研] 354求调剂 +4 Tyoumou 2026-03-18 7/350 2026-03-18 21:45 by Tyoumou
[考研] 274求调剂 +5 S.H1 2026-03-18 5/250 2026-03-18 21:27 by guosr9609
[考研] 085600材料与化工调剂 324分 +8 llllkkkhh 2026-03-18 8/400 2026-03-18 21:01 by Catalysis25
[考研] 267一志愿南京工业大学0817化工求调剂 +8 SUICHILD 2026-03-12 8/400 2026-03-18 20:55 by winsuccess
[考研] 344求调剂 +6 knight344 2026-03-16 7/350 2026-03-18 20:13 by walc
[考研] 一志愿中国海洋大学，生物学，301分，求调剂 +4 1孙悟空 2026-03-17 4/200 2026-03-18 17:59 by fivewind
[教师之家] 焦虑 +8 水冰月月野兔 2026-03-13 12/600 2026-03-18 15:27 by 咪呜喵呜
[考研] 化学工程321分求调剂 +15 大米饭！ 2026-03-15 18/900 2026-03-18 14:52 by haxia
[考研] 0817 化学工程 299分求调剂 有科研经历 有二区文章 +7 rare12345 2026-03-18 7/350 2026-03-18 14:31 by laoshidan
[考研] 材料专硕306英一数二 +10 z1z2z3879 2026-03-16 13/650 2026-03-18 14:20 by 007_lilei
[考研] 0703化学调剂 +4 pupcoco 2026-03-17 7/350 2026-03-18 12:14 by djl2006
[考研] 11408 一志愿西电，277分求调剂 +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
[考研] [导师推荐]西南科技大学国防/材料导师推荐 +3 尖角小荷 2026-03-16 6/300 2026-03-16 23:21 by 尖角小荷
[基金申请] 今年的国基金是打分制吗？ 50+3 zhanghaozhu 2026-03-14 3/150 2026-03-16 17:07 by 北京莱茵润色
[考研] 326求调剂 +3 mlpqaz03 2026-03-15 3/150 2026-03-16 07:33 by Iveryant
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
[考研] 274求调剂 +3 S.H1 2026-03-12 3/150 2026-03-13 15:15 by JourneyLucky
[考研] 321求调剂（食品/专硕） +3 xc321 2026-03-12 6/300 2026-03-13 08:45 by xc321
信息提示
关闭
请填处理意见
关闭