| 查看: 1653 | 回复: 3 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
关于dM/dH退磁曲线 已有1人参与
|
||
| 大神们有谁知道如何得到dM/dH的图像，是M对H的整个求导，还是部分？M,和H都是磁滞回线就是上M和H 吗如何说明交换耦合 |
回复此楼» 猜你喜欢
综述论文作为代表作会不会影响评审专家的印象分？
已经有13人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有4人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有29人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有9人回复
-
是这周出结果还是下周出结果？
已经有3人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有9人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有4人回复
-
帮忙看看fileCode
已经有13人回复
-
为什么网上很多人说本周 12号出结果
已经有7人回复
-
读博的好处
已经有3人回复
1174895362
新虫 (小有名气)
- 应助: 5 (幼儿园)
- 金币: 568.9
- 红花: 7
- 帖子: 81
- 在线: 35.6小时
- 虫号: 2343257
- 注册: 2013-03-13
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
3楼2016-09-13 21:39:25
【答案】应助回帖
|
论文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
sam799824810
银虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 837.9
- 帖子: 97
- 在线: 41小时
- 虫号: 2607178
- 注册: 2013-08-20
- 性别: GG
- 专业: 有机高分子功能材料
4楼2016-09-14 08:40:28
