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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 材料区 » 功能材料 » 磁性材料 » 关于dM/dH退磁曲线
41/1返回列表
查看: 1551  |  回复: 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的回帖

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的回帖

sam799824810

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

发自小木虫Android客户端
回复此楼
4楼2016-09-14 08:40:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zxlzqz 的主题更新
41/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085700资源与环境308求调剂 +3 墨墨漠 2026-03-18 3/150 2026-03-18 22:35 by bingxueer79
[考研] 求材料调剂 +8 隔壁陈先生 2026-03-12 8/400 2026-03-18 22:19 by li123456789.
[考研] 一志愿985，本科211，0817化学工程与技术319求调剂 +7 Liwangman 2026-03-15 7/350 2026-03-18 20:08 by walc
[考研] 0703化学调剂 +7 妮妮ninicgb 2026-03-15 11/550 2026-03-18 19:48 by macy2011
[考研] 一志愿武理材料305分求调剂 +5 想上岸的鲤鱼 2026-03-18 6/300 2026-03-18 17:53 by 无际的草原
[考研] 085600材料与化工 +5 安全上岸！ 2026-03-16 5/250 2026-03-18 15:33 by cmz0325
[考研] 材料与化工一志愿南昌大学327求调剂推荐 +8 Ncdx123456 2026-03-13 9/450 2026-03-18 14:40 by haxia
[考研] 311求调剂 +11 冬十三 2026-03-15 12/600 2026-03-18 14:36 by 星空星月
[考研] 材料专硕274一志愿陕西师范大学求调剂 +6 薛云鹏 2026-03-13 6/300 2026-03-18 14:14 by 脱颖而出
[考研] 303求调剂 +4 睿08 2026-03-17 6/300 2026-03-18 11:01 by Iveryant
[考研] 0703化学调剂 +3 妮妮ninicgb 2026-03-17 3/150 2026-03-18 10:29 by macy2011
[考研] 268求调剂 +6 简单点0 2026-03-17 6/300 2026-03-18 09:04 by 无际的草原
[考研] 301求调剂 +9 yy要上岸呀 2026-03-17 9/450 2026-03-18 08:58 by 无际的草原
[考研] 考研求调剂 +3 橘颂. 2026-03-17 4/200 2026-03-17 21:43 by 有只狸奴
[考研] 一志愿华中师范071000，325求调剂 +6 RuitingC 2026-03-12 6/300 2026-03-16 14:50 by 可淡不可忘
[考研] 材料与化工求调剂一志愿 985 总分 295 +8 dream…… 2026-03-12 8/400 2026-03-13 22:17 by 星空星月
[考研] 329求调剂 +3 miaodesi 2026-03-12 4/200 2026-03-13 20:53 by 18595523086
[考研] 一志愿211化学学硕310分求调剂 +8 努力奋斗112 2026-03-12 9/450 2026-03-13 15:41 by JourneyLucky
[考研] 283求调剂，材料、化工皆可 +8 苏打水7777 2026-03-11 10/500 2026-03-13 09:06 by Linda Hu
[考研] 333求调剂 +3 152697 2026-03-12 4/200 2026-03-13 07:08 by Iveryant
信息提示
关闭
请填处理意见
关闭