 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
小木蟲論壇-學術科研互動平臺 » 材料區(qū) » 功能材料 » 磁性材料 » 關于dM/dH退磁曲線
| 4 | 1/1 | 返回列表 |
| 查看: 1556 | 回復: 3 | ||
zxlzqz木蟲 (正式寫手)
|
[求助]
關于dM/dH退磁曲線 已有1人參與
|
| 大神們有誰知道如何得到dM/dH的圖像，是M對H的整個求導，還是部分？M,和H都是磁滯回線就是上M和H 嗎如何說明交換耦合 |
回復此樓» 猜你喜歡
296求調劑
已經有5人回復
-
0703化學調劑
已經有13人回復
-
085600材料與化工 求調劑
已經有15人回復
-
085700資源與環(huán)境308求調劑
已經有4人回復
-
生物學調劑招人！�。�
已經有4人回復
-
0817 化學工程 299分求調劑 有科研經歷 有二區(qū)文章
已經有17人回復
-
一志愿武漢理工材料工程專碩調劑
已經有5人回復
-
0856調劑，是學校就去
已經有7人回復
-
申博26年
已經有3人回復
-
一志愿吉林大學材料學碩321求調劑
已經有10人回復
【答案】應助回帖
|
論文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
1174895362
新蟲 (小有名氣)
- 應助: 5 (幼兒園)
- 金幣: 568.9
- 紅花: 7
- 帖子: 81
- 在線: 35.6小時
- 蟲號: 2343257
- 注冊: 2013-03-13
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
3樓2016-09-13 21:39:25
sam799824810
銀蟲 (小有名氣)
- 應助: 1 (幼兒園)
- 金幣: 837.9
- 帖子: 97
- 在線: 41小時
- 蟲號: 2607178
- 注冊: 2013-08-20
- 性別: GG
- 專業(yè): 有機高分子功能材料
4樓2016-09-14 08:40:28
相關版塊跳轉
我要訂閱樓主 zxlzqz 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 288求調劑 +14 | 于海海海海 2026-03-19 | 14/700 |
|
|---|---|---|---|---|
|
[考研] 招收調劑碩士 +4 | lidianxing 2026-03-19 | 10/500 |
|
|
|
[考研] 一志愿北京化工大學0703化學318分，有科研經歷，求調劑 +3 | 一瓶苯甲酸 2026-03-14 | 3/150 |
|
|
|
[考研] 一志愿天津大學化學工藝專業(yè)（081702）315分求調劑 +11 | yangfz 2026-03-17 | 11/550 |
|
|
|
[考研] 317求調劑 +3 | 申子申申 2026-03-19 | 6/300 |
|
|
|
[考研] 一志愿福大288有機化學，求調劑 +3 | 小木蟲200408204 2026-03-18 | 3/150 |
|
|
|
[考研] 一志愿天大材料與化工（085600）總分338 +5 | 蔡大美女 2026-03-13 | 5/250 |
|
|
|
[考研] 0817調劑 +3 | 沒有答案_ 2026-03-14 | 3/150 |
|
|
|
[考研] 085410人工智能專碩317求調劑（0854都可以） +3 | xbxudjdn 2026-03-18 | 3/150 |
|
|
|
[考研] 344求調劑 +6 | knight344 2026-03-16 | 7/350 |
|
|
|
[考研] 311求調劑 +11 | 冬十三 2026-03-15 | 12/600 |
|
|
|
[考研] 301求調劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|
|
[碩博家園] 湖北工業(yè)大學 生命科學與健康學院-課題組招收2026級食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
|
[考研] 材料專碩326求調劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
|
[考研] 274求調劑 +5 | 時間點 2026-03-13 | 5/250 |
|
|
|
[考研] 080500，材料學碩302分求調劑學校 +4 | 初識可樂 2026-03-14 | 5/250 |
|
|
|
[考研] 297一志愿上交085600求調劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|
|
|
[考研] 297求調劑 +4 | 學海漂泊 2026-03-13 | 4/200 |
|
|
|
[考研] 330求調劑 +3 | ?醬給調劑跪了 2026-03-13 | 3/150 |
|